Accumulation capacity of soluble rocks as the main stability index in covered karst areas

Full Text

ВВЕДЕНИЕ

Какими бы слабыми ни были грунты, перекрывающие закарстованный массив, каким бы сильным ни было техногенное воздействие, проявления карста в основании сооружений и на земной поверхности исключаются, если аккумуляционная емкость растворимых пород, отражающая их способность принимать и накапливать обломочный материал, недостаточно велика. Другими словами, трещиноватость и закарстованность подстилающих отложений определяют саму возможность образования провалов. Но именно эти показатели не только состояния массива, но и провальной опасности до сих пор остаются наименее изученными, что не в последнюю очередь связано с высокой стоимостью глубокого бурения.

В отличие от горного дела, которое оперирует заданными геометрическими характеристиками полого подземного пространства, инженерная геология сталкивается с высокой степенью их неопределенности. Математически строгое описание не просто трещинно-порового, но и полого пространства закарстованного массива – задача чрезвычайно сложная. Если же кроме его геометрии нужно учитывать растворимость горных пород, режим подземных вод и техногенные изменения последнего – сложная вдвойне.

Решение этой проблемы требует долгой и кропотливой работы, в основе которой лежит совместный анализ данных инженерно-геологических, гидрогеологических, гидрохимических, геофизических, спелеологических и других исследований. Несмотря на серьезные трудности и высокие трудозатраты, ставить и решать задачу оценки аккумуляционной емкости закарстованных отложений необходимо.

ОБЪЕМ ВЫНЕСЕННЫХ ГРУНТОВ И РАЗРЫХЛЕНИЕ ПОКРОВНОЙ ТОЛЩИ

Самым общим результатом выноса несвязных и раздробленных связных грунтов объемом DV из покровной толщи с не очень высокой исходной пористостью n₀ является увеличение последней до среднего значения n в области деформирования-разрушения объемом V₀ (рис. 1). И уравнение баланса массы грунта до и после разрыхления записывается в виде [5, 7]:

r_s(1 - n₀)V₀ = r_s(1 - n)V = r_s(1 - n)(V₀ + DV),(1)

где r_s- плотность минеральной части грунта; V = V₀ + DV.

Из уравнения (1) следует, что

DV = V₀(n - n₀)/(1 - n) = kV₀,(2)

а из опытов – что провал земной поверхности образуется только после того, как высота области V₀ достигнет мощности толщи h [5, 17, 24]. При этом для канонических схем области деформирования-разрушения, показанных на рис. 1, минимальный объем грунтов DV = V_кр, удаление которых обеспечивает появление заметной глазу воронки, равен:

V_кр = kpR²h/3,(3)

V_кр = kpR²h, (4)

V_кр = 4kpa²h/3, (5)

где a – малая полуось “эллипсоида”, крайние значения которой x₁, x₂ (x₁ ≤ a ≤ x₂) определяются моделью зонального строения области влияния ослабленного участка [8].

Рис. 1. Схемы к расчету критического объема грунта V_кр = ∆V, вынос которого из области деформирования V₀, аппроксимированной круговым конусом, цилиндром и эллипсоидом, приведет к провалам земной поверхности: 1 – грунтовая толща, 2 – закарстованный массив, z_0,1,2,3 – характерные отметки мощности в модели области влияния отверстия [8], позволяющей определить объем V₀.

Уравнение (2) и вытекающие из него зависимости (3-5) могут служить для оценки устойчивости территории и прогноза возможности появления на ней воронок, если сопоставить значение V_кр с величиной аккумуляционной емкости V_а закарстованного массива. Из них также следует, что V_кр→ 0 при h → 0 и особенно при R → 0, так как объем V_кр прямо пропорционален h в первой степени и R = D/2 – во второй [7]. Этим хорошо объясняется тот факт, что деформирование мощной грунтовой толщи и оседание большой площади земной поверхности происходят, как правило, в форме медленного прогиба с небольшой амплитудой, тогда как провалы появляются внезапно, могут быть очень глубокими и иногда возобновляются. Кроме того, V_кр→ 0 при n₀→ n, что указывает на возможность просадок грунтов при отсутствии закарстованных отложений, например при увлажнении (лессы) или вибрации (рыхлые пески).

Большой опыт горных инженеров свидетельствует, что три основные формы разрыхления (см. на рис. 1) наблюдаются в самых разных породах – скальных, полускальных и дисперсных. А переход пород из монолитного в раздробленное состояние характеризуется коэффициентом разрыхления [1, 9, 11, 15, 18, 20-22, 25, 28]:

К_р = r_м/r_р, (6)

где r_м, r_р – плотности пород в массиве исходная и в раздробленном состоянии соответственно. При первоначальном обрушении К_р может изменяться в пределах 1.08-2.10. Однако чаще всего его значения, при которых породы, разрушаясь, заполняют полое пространство и воспринимают на себя давление вышележащих пластов, составляют К_р = 1.15-1.80. Используется также коэффициент остаточного разрыхления К_ор, которым характеризуются раздробленные породы по истечении времени или в результате нескольких циклов обрушения. Он заметно меньше коэффициента первичного разрыхления: К_ор= 1.10-1.30 [15, 25, 28].

Нетрудно видеть, что коэффициент К_р = r_м/r_р = (1 – n₀)/(1 – n) в формуле (6) связан с коэффициентом k = (n - n₀)/(1 - n₀) в выражениях (2-5) соотношением:

К_р = k + 1. (7)

Поэтому уравнение (2) записывается также в виде:

DV = (К_р - 1)V₀. (2а)

Таким образом, рассмотренный выше подход к прогнозу провалов справедлив для районов не только покрытого, но и голого карста, которые занимают соответственно 70% и 30% площади распространения на земной суше растворимых отложений (46.9×10⁶ км²). Однако первым, кто предложил оценивать устойчивость карбонатного массива, сравнивая объем карстовой полости с объемом заполняющих ее при обрушении скальных пород, был, по-видимому, Г.М. Шахунянц [27].

АККУМУЛЯЦИОННАЯ ЕМКОСТЬ ЗАКАРСТОВАННОГО МАССИВА

Использование обсуждаемого в статье, по сути дела, горно-геологического подхода к оценке устойчивости закарстованных территорий требует нахождения аккумуляционной емкости V_а растворимых отложений − объема их полого пространства, которое способно принимать и накапливать обломочный материал. Понятие “аккумуляционная емкость”, введенное достаточно давно [5-7], предполагает зависимость этой характеристики от многих условий и факторов, в том числе от минерального состава, трещиноватости и растворимости закарстованных пород, режима и химического состава подземных вод. В первом приближении аккумуляционная емкость карстующихся отложений равна:

V_а= V_м(К_к +v_к•T/Dh),(8)

где V_м – объем массива пород или области аккумуляции, в пределах которой может откладываться обломочный материал; Dh – высота этой области (V_м/Dh – ее средняя площадь); К_к – коэффициент открытой трещинно-карстовой пустотности; v_к – скорость карстовой денудации; Т – рассматриваемый промежуток времени (в инженерной геологии обычно срок службы зданий и сооружений). Если v_к→ 0 и v_к•T/Dh << К_к, то:

V_а = V_м•К_к. (8а)

Так, для инженерно-геологических условий, показанных на рис. 2, область аккумуляции V_м представляет собой усеченный конус, и выражение (8а) принимает вид:

V_а= pmК_к[3D² + 2m•ctgj_o´(3D + 2m•ctgj_o´)]/12,(9)

где m = Dh – мощность карбонатных пород, D = 2R – диаметр или пролет ослабленного участка (входного отверстия в трещинно-поровое пространство), j_o´– угол откоса песков под водой. При малом пролете входного отверстия (D << m•ctgj_o´) уравнение (9) упрощается:

V_а=pm³К_к•ctg²j_o´/3. (9а)

Если пьезометрическая поверхность трещинно-карстовых вод имеет заметный уклон, то их горизонтальная скорость велика, угол при нижнем основании конуса V_м (см. рис. 2) в направлении движения воды будет значительно меньше, а в противоположном направлении – несколько больше j_o´, и значение V_м увеличится.

Рис. 2. Схема к определению аккумуляционной емкости растворимых пород: 1 – пески, 2 – закарстованные известняки с маломощным глинистым элювием в кровле, 3 – глины тугопластичные, H ≥ H₀ – уровни подземных вод, D = 2R – диаметр (пролет) участка ослабления массива – входного (в область аккумуляции V_м) отверстия; j_о´ – угол откоса песков под водой.

При интенсивном же отборе трещинно-карстовых вод объем V_м и, следовательно, емкость V_а возрастают очень сильно (рис. 3).

Рис. 3. Вид физической модели сверху после зачистки верхней песчаной толщи и экранирующего глинистого слоя, разрушенных в ходе опытов (а), и увеличение области аккумуляции V_м в результате интенсивной откачки трещинно-карстовых вод (б).

а: 1, 2 – лицевая и задняя стенки лотка; 3 – выполненный песком “палеоврез” (ослабленная зона); 4 – суффозионные полости; 5 – трещинно-поровое пространство “закарстованных пород”, запечатанное песчано-глинистым материалом; 6 – правая дрена для снижения напора воды в нижнем слое модели (левая дрена на фотографии не видна); стрелками показано направление движения напорной воды.

б: 4 – водозаборные скважины; остальные обозначения приведены на рис. 2.

На рис. 3а видно, что после принудительного разрушения экранирующего слоя над системой суффозионных полостей и выноса вышележащих песков в поровое пространство нижнего слоя, имитирующего закарстованные породы, это пространство оказывается полностью закольматированным. При этом размеры области аккумуляции в массиве зависят от положения эксплуатационных скважин в плане, интервала и интенсивности откачек (рис. 3б). В пользу этого вывода свидетельствуют и результаты, полученные Ф.В. Котловым в ходе комплексного изучения суффозии на участке Замоскворецкой дренажной системы [16].

Обращает на себя внимание, что в бортах ослабленной зоны 3, несмотря на одинаково задаваемые там условия, возникают не сплошные линейные, а отдельные (локальные) суффозионные полости (см. рис. 3а). Причем наибольших размеров они достигают у прозрачных стенок лотка 1, 2, где и образовались провалы в модели. Очевидно, это связано с тем, что на контактной поверхности работает лишь трение скольжения грунтов по стенке, которое заметно меньше угла их внутреннего трения. Именно снижением трения во многом объясняется образование воронок над подземными коллекторами или у стенок колодцев и буровых скважин. Поэтому модельную ситуацию, показанную на рис. 3а, с позиций оценки геологического строения можно трактовать как пересечение погребенного эрозионного вреза (3) разрывными тектоническими нарушениями (1, 2).

ЧАСТНЫЙ СЛУЧАЙ АККУМУЛЯЦИИ ГРУНТОВ В МАССИВЕ

Из ф. (8а) следует, что при К_к= 1 оценка аккумуляционной емкости сводится к определению объема одиночной полости-приемника V_а = V_м = DV, которой в массиве может и не быть. Но даже если она есть, то чаще всего не настолько большая, чтобы обеспечить образование воронки на земной поверхности: V_а = V_м < DV = V_кр (см. рис. 1).

Объемы изолированной полости DV и области деформирования-разрушения (разрыхления) V₀ прямо пропорциональны их высотам – Dh и h₀ соответственно. Причем область V₀ может воспринимать на себя давление перекрывающих грунтов и “подбучивать”, т.е. подпирать их, не достигая земной поверхности, поэтому в общем случае h₀ ≠ h. Тогда из выражений (2) и (2а) имеем:

Dh ~ h₀k = h₀(К_р - 1), (10)

h₀ ~ Dhk = Dh/(К_р - 1). (10а)

В работе В.Н. Андрейчука предложено следующее прогностическое уравнение [3]:

h_см = DhК_р/(К_р- 1), (11)

где, по образному выражению автора, h_см – высота смыкания обвального свода и догоняющей его осыпи, при которой провалообразование прекращается. Вывод ф. (11) не приводится, но поскольку h_см = Dh + h₀ (рис. 4), то легко видеть, что это уравнение является еще одной формой записи условий (10) и (10а), если высоту области разрыхления V₀ считать не от кровли, а от подошвы полости.

Рис. 4. Схемы обрушения горных пород над полостью, для которых в выражениях (10) справедлив знак равенства (а, б), и схема наблюдаемого в массиве обрушения кровли пещер (в).

Подход, основанный на сравнении h₀ (или h_см) и Dh, применяет и Г. Адерхолд [2], реконструируя механизм карстовых провалов в Земле Гессен (Германия). Того же подхода придерживается В.П. Хоменко [26], выполняя ретроспективный прогноз известного в Нижегородской области Пивоваровского провала. А в статье И.А. Южанина и А.М. Терлецкого [28] рассмотрена методика послойного (снизу вверх) суммирования мощности отдельных слоев после разрушения h_i∙К_р и сопоставления значения ∑h_i∙К_р с общей высотой возможной зоны обрушения m_в + ∑h_i, где m_в – вынимаемая мощность пласта очистной выработки (m_в≡ Dh).

Все вышесказанное свидетельствует о том, что одним из самых важных, если не главным, критерием устойчивости закарстованных территорий является глубина залегания ослабленных участков растворимых пород и, следовательно, мощность перекрывающих их грунтов. Это прямо следует из уравнений (2-5). И огромная заслуга В.Н. Андрейчука [3] заключается в том, что он наглядно показал справедливость этого утверждения на примере районов развития Кунгурской и Ледяной пещер на Урале и ряда пещер в Западной Украине.

К сожалению, в настоящее время величина мощности покровной толщи закарстованных массивов является критерием устойчивости только для территории Республики Башкортостан¹. В нормативных документах, регулирующих изыскания и строительство в других карстовых районах ^2,3,4,5,этот критерий не прописан.

Возвращаясь к уравнениям (10, 10а), отметим, что заменить в них знак пропорциональности “~” на знак равенства мы имеем право тогда, когда форма карстовой полости и области разрыхления над ней примерно одинакова (рис. 4 а, б). На это обстоятельство обратил внимание еще Г.М. Шахунянц [27].

На рис. 4в показана схема развития процесса над карстовой пещерой в зрелой стадии ее развития. Она отражает результаты спелеологических наблюдений и обследования гротов Кунгурской и Ледяной пещер на Урале [3, 4, 13], а также Новомурадымовской пещеры в Башкирии. Очевидно, что зависимость h_см (или h₀) от Dh при этом будет более сложной, чем та, что предусмотрена ф. (11). Этот вопрос подробно рассматривается в последнем разделе статьи на конкретном примере. Ниже остановимся на некоторых особенностях оценки устойчивости закарстованных и подработанных территорий.

НЕКОТОРЫЕ ОСОБЕННОСТИ ОЦЕНКИ ПРОВАЛЬНОЙ ОПАСНОСТИ В КАРСТОВЫХ РАЙОНАХ

Из трех канонических схем (см. на рис. 1) условию равенства форм областей выноса (деформирования, разрыхления покровной толщи) и аккумуляции удовлетворяет только схема на рис. 1б. Но если трубообразная форма канала течения или обрушения присуща многим воронкам, то вертикальная цилиндрическая полость в растворимых породах – явление довольно редкое. Крупные горизонтальные полости, протяженность которых много больше их высоты, встречаются гораздо чаще.

Еще чаще приходится сталкиваться с ситуацией, когда трещинно-карстовые коллекторы представляют не вполне ясную пространственную систему связанных между собой небольших полостей. Размеры же полого пространства подработанных массивов, как отмечалось во введении, заранее известны. И в этом состоит, пожалуй, главное отличие инженерной геологии карстовых районов от горного дела.

Другое серьезное отличие заключается в выборе и оценке величины коэффициента разрыхления. Известно [15], что с ростом отношения глубины выработки h к ее высоте m_в он возрастает и максимальных значений достигает при h/m_в ≈ 3-4, а затем снижается. Уменьшается он, как уже говорилось, и с течением времени, а также при увеличении количества циклов обрушения, что объясняется соответственно статическим и динамическим уплотнением грунтов. Эти зависимости учитываются при размещении выработок, определении порядка отработки пластов и управлении горным давлением [25, 28]. И горные инженеры вправе ориентироваться на промежуточные, между К_ор и К_р, ближе к первому, значения коэффициента разрыхления: К_ор ≤ К_р′ < К_р.

Тех же горно-геологических позиций, по-видимому, придерживается и В.Н. Андрейчук [3], не затрагивающий этот вопрос, но рекомендующий достаточно низкие значения К_р некоторых типов пород, встречающихся в районе Кунгурской и Ледяной пещер. Так, для песчаных и гравийно-галечниковых отложений назначаемый им коэффициент разрыхления равен К_р = 1.0-1.1.

Действительно, в плотных моделях, приготовленных из сухого кварцевого мелкого полуокатанного песка, значения коэффициента разрыхления получаются равными К_р = 1.05-1.12 [5-8]. Однако заметим, что даже при самой тщательной укладке и трамбовке плотность этих моделей не превышала r₀ = 1.60-1.67 г/см³, т.е. начальная пористость была больше либо равной n₀ = 0.37-0.40. В массиве же значения n₀ ≥ 0.37-0.40 характерны для песчаной толщи средней плотности сложения [12]; у плотных мелких, а тем более средних и крупных песков начальная пористость заметно меньше.

Средняя конечная пористость в области деформирования моделей составляла n = 0.42-0.44, тогда как для максимально рыхлого сложения идеальных сыпучих тел n = 0.48. В несвязных грунтах, частицы которых угловатые, и в пылеватых песках она может быть больше n = 0.48. Кроме того, в грунтовой толще пески, как правило, водонасыщенные или влажные. А в них, хотя и по разным причинам, могут существовать, не разрушаясь, достаточно крупные пустоты – от 2-10 мм до первых дециметров [7]. И их общая пористость может быть существенно выше, чем у воздушно-сухих песков. Именно поэтому для оценки устойчивости несвязных грунтов предложено пользоваться величиной конечной пористости n = 0.5 [5, 6].

С учетом всего вышесказанного значения коэффициента разуплотнения песчаной толщи, равные К_р = 1.15-1.25, не кажутся чрезмерно большими. Для гравийно-галечниковых отложений они будут еще больше хотя бы по причине высоких значений угла дилатансии [22].

Если не рассматривать обрушение пород над крупной карстовой полостью, то в отличие от подработанных массивов маловероятно, что вынос обломочного материала в трещинно-поровое пространство растворимых отложений будет сопровождаться его последующим уплотнением в этом пространстве – статическим или динамическим. Скорее, следует ожидать неограниченного увеличения пористости вынесенных грунтов на фоне роста области аккумуляции V_м (см. рис. 2, 3). В области выноса V₀ (см. рис. 1) уплотнение возможно, однако появление воронки на земной поверхности или под сооружением, на наш взгляд, намного опасней ее последующей эволюции, так как мы знаем об этом и можем принять соответствующие тому или иному случаю меры защиты. О редких крупных провалах с катастрофическими последствиями речь сейчас не идет, а примеры успешного выполнения противокарстовых мероприятий и предотвращения большого экономического и социального ущерба после того, как воронка образовалась в основании здания, есть [19].

Таким образом, используя обсуждаемый подход для оценки возможности возникновения карстово-суффозионных воронок в покровной толще закарстованного массива и на ее поверхности, необходимо учитывать форму не только области выноса (деформирования-разрушения, разрыхления) грунтов (см. рис. 1), но и области аккумуляции (см. рис. 2-4). И, проводя расчеты, следует ориентироваться на значения коэффициентов k или К_р, но не К_ор. При этом, конечно, надо быть уверенным, что собственно карстовый процесс (коррозия, выщелачивание и отчасти эрозия) не сильно отразится на величине аккумуляционной емкости растворимых пород в течение прогнозируемого отрезка времени Т. Другими словами, нужно показать, что эти породы закарстованные, но не карстующиеся.

УСТОЙЧИВОСТЬ ОСНОВАНИЯ НИЖЕГОРОДСКОЙ АТОМНОЙ ЭЛЕКТРОСТАНЦИИ (НИАЭС)

Участок строительства НИАЭС площадью 20 км² выбран на границе Нижегородской и Владимирской областей по ряду причин энергетического и социально-экономического характера с учетом [10] и требований НП-032-01⁶. Учет инженерно-геологических условий при этом отошел на второй план, и Нижегородская АЭС оказалась первой в мире атомной станцией, размещенной на территории развития активного сульфатно-карбонатного карста.

В ходе скрупулезных комплексных изысканий, выполненных проектным институтом ОАО “НИАЭП” с привлечением ряда производственных и научно-исследовательских организаций [10], с большим трудом выделена площадка размещения основных сооружений (3.25 км²), в границах которой проявления карста на земной поверхности отсутствуют. И, естественно, встал вопрос, не образуется ли “завтра” провал на этой площадке, пока еще свободной от поверхностных карстовых форм? Чтобы ответить на этот вопрос, нужно хотя бы вкратце остановиться на инженерно-геологических условиях строительства.

На рис. 5 видно, что площадка размещения энергоблоков удалена и от р. Оки, и от палеовреза р. Бол. Кутра. Мощность слабых четвертичных грунтов здесь минимальная, а твердых глинистых и полускальных отложений уржумского яруса – максимальная. Гидрогеологические условия участка характеризуются наличием двух водоносных комплексов. Верхний, напорно-безнапорный слабоводоносный комплекс с абсолютными отметками установившегося уровня от 92-118 до 124-126 м приурочен к четвертичным пескам, прослоям и линзам уржумских песчаных пород, известняков, доломитов и мергелей. Нижний комплекс состоит из казанско-сакмарского (верхнего) и сакмарского (нижнего) водоносных пластов. Они разделены водоупорными ангидритами и монолитными гипсами общей мощностью до 20-28 м. Но на площадке АЭС верхний водоносный (казанско-сакмарский) пласт практически отсутствует, а приуроченные к сакмарским доломитам воды нижнего пласта залегают на глубине 97-98 м. Их уровень устанавливается на 46-50 м выше кровли доломитов, на абсолютных отметках 80-97 м.

Рис. 5. Геологический разрез участка строительства Нижегородской АЭС (по материалам ОАО “НИАЭП”). Пермская система, нижний отдел: P₁s – сульфатно-карбонатная толща сакмарского яруса; средний отдел: P₂kz – доломиты, мергели и глины казанского яруса; P₂ur₁, P₂ur₂ – мергелисто-глинистая толща нижнего (P₂ur₁¹, P₂ur₂² – нижняя и верхняя пачки) и верхнего подъярусовуржумского яруса. Плиоцен – нижний плейстоцен: acN₂–Q₁, dpN₂–Q₁ – песчано-глинистые аллювиально-коллювиальные и оползневые отложения. Четвертичная система: gIkš² – моренные суглинки донского оледенения; f,lgIkš¹, f,lgIkš³ – водноледниковые пески времени наступания и отступания ледника; aIIck-ms –песчаные отложения III н/п террасы (чекалинский и московский горизонты); dlII-III – делювиальные суглинки и супеси; aIIImk-kl, aIIIln-os, aIV – аллювиальные песчано-глинистые грунты II (микулинский и калининский горизонты) и I (ленинградский и осташковский горизонты) н/п террас и поймы. НИАЭС – площадка строительства.

Из результатов инженерно-геологических изысканий следует также, что водопроницаемость карбонатных пород массива-основания АЭС компенсируется их слабой растворимостью, а средняя растворимость сульфатных пород – их непроницаемостью. Начиная с глубины 18-19 м, подземные воды становятся слабоагрессивными и к сульфатным, и к карбонатным породам. Их минерализация увеличивается с глубиной и на отметках подошвы уржумской толщи настолько высока (3-5 г/л), что следует ожидать не растворения пород, а выпадения вещества из раствора.

Конечно, гипс может растворяться, даже если минерализация воды достигает 7 г/л, и его растворимость зависит не только от минерализации подземных вод [14]. Однако в данных условиях наблюдается именно осаждение гипса. Об этом свидетельствуют результаты комплексных геологических, гидрогеологических и геохимических исследований – натурных и лабораторных.

Изменение напора трещинно-карстовых вод в окрестности промплощадки крайне мало: I = DH/DL = 0.003–0.004, где DH = 8 м – падение напора на пути движения подземных вод DL = 2.0-2.7 км, мала и скорость их горизонтального движения. Таким образом, из четырех условий развития карста, сформулированных Д.С. Соколовым [23], не выполняются три: 1 – для пород и 2 – для подземных вод, что позволяет оценить устойчивость массива пород, исходя из определенных предпосылок следующим образом.

Максимальная мощность сильнотрещиноватых и раздробленных сакмарско-казанских отложений, вскрытых скважинами в их кровле, равна (h_о.з)_max = 4.1 м за пределами стройплощадки. Тем не менее допускаем существование на площадке АЭС (3.25 км²) не только этой ослабленной (раздробленной) зоны, но и достаточно большой в плане полости высотой Dh ≈ (h_о.з)_max = 4 м.

Средняя мощность покровной толщи на площади 3.25 км² составляет h = 69.6 м, и почти вся эта толща представлена уржумскими отложениями, которые относятся в основном к твердым связным и полускальным грунтам (К_р = 1.3-1.5). Поэтому принимаем: h = 70 м, (К_р)_ср = 1.4 или, в соответствии с ф. (7), (k)_ср = 0.4, а также цилиндрическую форму канала обрушения уржумских пород.

Учитывая низкую скорость движения трещинно-карстовых вод, считаем на основании уравнения Бернулли, что скорость движения воды в крупной изолированной полости будет пренебрежимо малой. Это позволяет выбрать коническую форму области аккумуляции с углом откоса j_о = 45° (рис. 6).

Рис. 6. Схематизация массива пород на площадке строительства энергоблоков и стадии обрушения покровной толщи: ∆h – высота гипотетической карстовой полости на контакте сакмарско-казанских (P₁s–P₂kz) растворимых и уржумских (P₂ur) нерастворимых отложений; L_с – пролет свода обрушения и области разрыхления грунтов; h₀ – максимальная высота канала обрушения, h_н – мощность ненарушенной толщи, h₀ + h_н = h – мощность покровной толщи, с поверхности сложенной маломощными четвертичными отложениями (Q).

Полагая в ф. (4) h ≤ h₀ и R = L_c/2, подставляем в нее k = k_ср = 0.4 и заменяем V_кр на объем усеченного круглого конуса высотой Dh, диаметр верхнего и нижнего оснований которого равен соответственно L_c и L_о = L_c + 2Dh (см. рис. 6). Тогда после несложных преобразований ф. (4) получим условие выхода провала на земную поверхность в виде:

(Dh- 0.4h₀)L_c² + 2Dh²L_c + 4Dh³/3 ≥ 0. (12)

Неравенство (12) по-разному выполняется в 3-х случаях: при: 1) Dh - 0.4h₀ > 0; 2) Dh - 0.4h₀ = 0, и 3) Dh - 0.4h₀ < 0. Поскольку Dh = 4 м, h₀ = 70 м и Dh - 0.4h₀ = -24 м, то мы имеем дело с последним, третьим случаем, значения L_c для которого лежат в интервале 0-2.67 м.

Столь малая ширина области деформирования (разуплотнения, обрушения) характерна для каналов течения несвязных грунтов. Иными словами, проявления карста на земной поверхности при условии 0 ≤ L_c ≤ 2.67 м возможны лишь в случае реализации процесса массовой суффозии [8]. Но доля суффозионно-неустойчивых грунтов в общей мощности пород уржумского яруса средней перми (P₂ur) очень мала. Выдержанный по простиранию прослой мелких песков мощностью 2.9-4.6 м наблюдается лишь в подошве отложений верхнего подъяруса (P₂ur₂). В нижней пачке нижнего подъяруса (P₂ur₁¹) линзы и прослои алевритов и песков вообще отсутствуют (см. рис. 5). Отсюда следует, что образование провалов вследствие развития суффозии массы на стройплощадке исключается.

Поэтому, исходя из свойств и состояния уржумских пород, полагаем, что пролет свода их обрушения (диаметр верхнего основания усеченного конуса, в пределах которого накапливается обломочный материал) будет в 3-5 раз больше значения L_c = 2.67 м. Решая уравнение (12) относительно h₀ подстановкой Dh = 4 м, L_c= 8.0-13.4 м, находим, что высота области деформирования-разрушения равна h₀ = 23.3–17.2 м, т.е. примерно в 3-4 раза меньше h.

Таким образом, в природных условиях площадка размещения основных сооружений НИАЭС характеризуется 3-4-кратным запасом устойчивости: К_у = h/h₀ = 3.0-4.1. При этом, учитывая, что мощность ненарушенной толщи равна h_н = h - h₀ ≈ 47-53 м (см. рис. 6), трудно представить себе приложенное с поверхности техногенное воздействие (статическое, динамическое или фильтрационное), которое могло бы вывести нижнюю разрушенную часть покровной толщи из равновесия. Другое дело, если оно будет приложено изнутри массива, например в виде откачек трещинно-карстовых вод, но согласно проектным решениям такой вид воздействия исключается.

ЗАКЛЮЧЕНИЕ

Для инженера-геолога районы покрытого карста представляют наибольший интерес не только потому, что имеют на Земле самое широкое распространение. Важнее то, что в отличие от районов голого карста они наиболее благоприятны, если не сказать пригодны, для хозяйственного освоения и жизни человека в целом. И чем больше мощность покровной толщи закарстованных массивов, тем более справедливо это утверждение.

В инженерной геологии карстовых районов или, что то же самое, в инженерном карстоведении наиболее актуальной, да и сложной, пожалуй, проблемой является локальный прогноз провалообразования. Особенно остро она стоит там, где воронки провала и оседания отсутствуют, не зарегистрированы или плохо выражены в рельефе, т.е. в районах нереализованной (потенциальной, ожидаемой) провальной опасности. Один из возможных подходов к решению этой проблемы, изложенный в статье, базируется на сравнении объема удаленных из покровной толщи несвязных и раздробленных связных грунтов с аккумуляционной емкостью закарстованного массива – объемом полого карстового пространства, способного принимать и накапливать обломочный материал.

Понятие “аккумуляционная емкость” предполагает зависимость этой характеристики от целого комплекса условий и факторов, в том числе минерального состава, трещиноватости и растворимости пород, режима трещинно-карстовых вод и интенсивности техногенного воздействия. Ее отождествление с объемом отдельной карстовой полости справедливо далеко не всегда, но главное – оно не способствует постановке и решению проблемы количественной оценки этого очень важного, если не основного показателя устойчивости закарстованных территорий.

Определение аккумуляционной емкости растворимых пород – сложная и далеко не решенная в настоящее время задача, которую, тем не менее, нужно и ставить, и решать. При этом необходимо учитывать кардинальные отличия инженерной геологии карстовых районов от горного дела, которое давно использует эффект увеличения объема горных пород при их обрушении в своих целях.

Разрыхление грунтов – второй краеугольный камень обсуждаемого в статье подхода к оценке устойчивости закарстованных территорий. Этот процесс служит одним из значимых механизмов самоорганизации геологической среды, которая “стремится” противостоять внешним воздействиям и сохранить равновесие. Здесь в явном виде проявляется справедливость принципа Ле Шателье-Брауна. Но в отличие от подработанных территорий, используя эффект разуплотнения грунтов при оценке провальной опасности в карстовых районах, следует ориентироваться на максимальные значения коэффициента первичного разрыхления пород. В противном случае прогнозируемая устойчивость массива окажется неоправданно заниженной.

Если инженерно-геологические условия поддаются достаточно простой схематизации, то рассмотренный в статье подход к прогнозу провалообразования уже сейчас можно использовать на практике. Об этом свидетельствуют результаты оценки устойчивости грунтового основания энергоблоков Нижегородской атомной электростанции

Особенность строительства НИАЭС заключается в том, что все ее объекты расположены на территории развития активного сульфатно-карбонатного карста. В результате инженерно-геологических изысканий выделена площадка размещения основных сооружений, где поверхностные формы карста отсутствуют, а мощность твердых глинистых грунтов, перекрывающих закарстованные породы, максимальна. Установлено также, что скорость растворения пород в зоне затрудненного водообмена крайне мала. В этих природных условиях, как следует из выполненных расчетов, основанных на сопоставлении аккумуляционной емкости закарстованного массива с объемом области разрыхления покровной толщи, площадка главных сооружений НИАЭС характеризуется 3-4-кратным запасом устойчивости.

Источник финансирования. Статья подготовлена в рамках выполнения государственного задания и плана НИР по теме №г.р. АААА-А19-119021190077-6.

Примечания:

¹ ТСН 302-50-95. РБ. Инструкция по изысканиям, проектированию, строительству и эксплуатации зданий и сооружений на закарстованных территориях / Госстрой Республики Башкортостан. Уфа, 1996. 41 с. https://files.stroyinf.ru/Index2/1/4294814/4294814948.htm

² Инструкция по проектированию зданий и сооружений в районах г. Москвы с проявлением карстово-суффозионных процессов. М.: Картолитография, 1984. 15 с. https://docplan.ru/Index1/45/45188.htm

³ Рекомендации по проведению инженерных изысканий, проектированию, строительству и эксплуатации зданий и сооружений на закарстованных территориях Нижегородской области / Департамент градостроительного развития Нижегородской обл. Нижний Новгород, 2012. 140 с. https://files.stroyinf.ru/Index2/1/4293758/4293758833.htm

⁴ ТСН 22–308–98 НН. Инженерные изыскания, проектирование, строительство и эксплуатация зданий и сооружений на закарстованных территориях Нижегородской области / Администрация Нижегородской обл. Комитет архитектуры и градостроительства. Нижний Новгород, 1999. 72 с. https://mooml.com/d/spravochnye-posobiya-k-snip/12969/

⁵ ТСН 11-301-2004 По. Инженерно-геологические изыскания для строительства на закарстованных территориях Пермской области / Администрация Пермской обл. Пермь, 2004. 122 с. https://studopedia.org/4-177843.html

⁶ НП-032-01. Размещения атомных станций. Основные критерии и требования по обеспечению безопасности / Госатомнадзор РФ. М., 2002. 11 с. https://internet-law.ru/stroyka/doc/41356

About the authors

A. V. Anikeev

Author for correspondence.
Email: anikeev_alex@mail.ru
Russian Federation, Ulanskii per., 13, str. 2, Moscow, 101000

References

Supplementary files

Supplementary Files

Action

1. JATS XML

Download

2. Fig. 1. Schemes for calculating the critical volume of soil Vcr = ∆V, the removal of which from the deformation region V0, approximated by a circular cone, cylinder and ellipsoid, will lead to failures of the earth's surface: 1 - soil thickness, 2 - karst massif, z0,1,2, 3 - characteristic power marks in the model of the region of influence of the hole [8], which makes it possible to determine the volume V0.

Download (150KB)

Indexing metadata

3. Fig. 2. Scheme for determining the accumulative capacity of soluble rocks: 1 - sands, 2 - karst limestones with low-power clay eluvium in the roof, 3 - refractory clay, H ≥ H0 - groundwater levels, D = 2R - diameter (span) of the area of mass attenuation - input (in the accumulation region Vm) hole; jо´ is the angle of repose of the sands under water.

Download (314KB)

Indexing metadata

4. Fig. 3. Top view of the physical model after cleaning the upper sandy layer and the shielding clay layer destroyed during the experiments (a), and the increase in the accumulation region Vm as a result of intensive pumping of fractured-karst waters (b).

Download (482KB)

Indexing metadata

5. Fig. 4. Schemes of rock collapse over the cavity, for which the equal sign (a, b) is valid in expressions (10), and the scheme of cave roof collapse observed in the cave massif (c).

Download (621KB)

Indexing metadata

6. Fig. 5. Geological section of the construction site of the Nizhny Novgorod nuclear power plant (based on the materials of NIAEP OJSC). Perm system, lower section: P1s - sulfate-carbonate stratum of the Sakmara stage; middle department: P2kz - dolomites, marls and clays of the Kazan tier; P2ur1, P2ur2 - marl-clay stratum of the lower (P2ur11, P2ur22 - lower and upper packs) and upper podjarusovurzhumskogo layer. Pliocene - Lower Pleistocene: acN2 – Q1, dpN2 – Q1 - sandy-clay alluvial-colluvial and landslide deposits. Quaternary system: gIkš2 - moraine loams of Don glaciation; f, lgIkš1, f, lgIkš3 - water-ice sands of the time of glacier advance and retreat; aIIck-ms - sandy deposits of the III n / a terrace (Chekalinsky and Moscow horizons); dlII-III - deluvial loam and sandy loam; aIIImk-kl, aIIIln-os, aIV - alluvial sandy-clay soils II (Mikulin and Kalinin horizons) and I (Leningrad and Ostashkov horizons) n / a terraces and floodplains. NIAES - construction site.

Download (959KB)

Indexing metadata

7. Fig. 6. Schematization of the rock mass at the construction site of power units and the stage of collapse of the cover stratum: ∆h is the height of the hypothetical karst cavity at the contact of Sakmara-Kazan (P1s – P2kz) soluble and Urzhum (P2ur) insoluble deposits; Lс - span of the collapse arch and soil loosening area; h0 is the maximum height of the collapse channel, hн is the thickness of the undisturbed thickness, h0 + hн = h is the thickness of the integument, from the surface composed of low-power Quaternary deposits (Q).

Download (280KB)

Indexing metadata