ГИДРОГЕОЛОГИЧЕСКИЕ ИССЛЕДОВАНИЯ НА ОРОШАЕМЫХ ТЕРРИТОРИЯХ
ГЕОФИЛЬТРАЦИОННАЯ СХЕМАТИЗАЦИЯ
Как уже отмечалось, в процессе геофильтрационной схематизации гидрогеологические закономерности представляются в гидродинамической (фильтрационной) постановке с полной математической формулировкой рассматриваемой задачи; при этом устанавливаются характер протекания процесса во времени, пространственная структура потока, граничные условия и внутреннее строение потока, определяемые закономерностями изменения геофильтрационных параметров.
Обычно рациональным оказывается решение рассматриваемой задачи в два этапа. На первом этапе задача носит региональный характер и направлена на обоснование изменения основных элементов питания подземных вод (площадного или напорного питания, расходов потока по контурам орошаемого массива); в ней рассматриваются основные региональные закономерности формирования подземных вод с учетом взаимодействия орошаемых массивов и очередности орошения и дренажа; при этом инженерные мероприятия (орошение и дренаж) задаются ориентировочно (их действие осредняется в пределах орошаемой территории).
На втором этапе рассматриваются локальные задачи, направленные на обоснование конкретных мелиоративных мероприятий и особенно основных элементов дренажа; при этом обычно могут изолированно рассматриваться отдельные участки орошаемой территории, в частности, участки между дренами, где, исходя из решения региональных задач, выбираются данные для задания площадного питания или бокового притока.
В соответствии с постановкой прогнозных расчетов следует выделить, по крайней мере, два типа фильтрационных схем, различающихся по масштабности рассматриваемых задач: региональные и локальные.
Региональной фильтрационной схемой охватывается вся область формирования подземных вод или отдельные ее части в пределах зоны влияния орошаемого массива.
При обосновании локальных фильтрационных схем рассматривается структура фильтрационного потока в зоне влияния конкретных мелиоративных мероприятий (в частности, дренажных сооружений) с учетом региональных закономерностей формирования, развития и разгрузки подземных вод. В этом случае наиболее важными представляются анализ вертикальной неоднородности потока и оценка значимости неравномерности инфильтрационного питания (во времени и пространстве).
Формирование геофильтрационного потока на орошаемых территориях происходит главным образом в пределах осадочных отложений, представленных чередованием песчаных и песчано-граве - листых водоносных пластов и глинистых разделяющих слоев, включая возможное наличие довольно пестрых покровных отложений, причем мощность этих отложений может быть весьма значительной (до 1000 м и более). В этих условиях важно найти оптимальные пути для наиболее четкого графического представления (картирования) этих элементов строения массива (потока), которые могли бы дать основу для гидрогеологического обоснования необходимости дренажа и принципиального обоснования его типа.
При постановке такой задачи следует иметь в виду, что основным фактором, определяющим выбор типа систематического дренажа и эффективность его действия, является наличие в разрезе хорошо проницаемых водоносных горизонтов [28, 35]. В связи с этим при картировании геофильтрационного строения для обоснования локальных расчетных схем систематического дренажа должны быть прежде всего установлены участки орошаемого массива, однотипные по развитию водоносных пластов, так что наиболее крупной таксонометрической единицей рассматриваемого районирования может быть выбрана фильтрационная система водоносных пластов.
Очевидно, детальность геофильтрационной схематизации так же, как и гидрогеологического районирования в целом, определяется решаемыми практическими задачами. Остановимся прежде всего на методике мелкомасштабной геофильтрационной схематизации (1 :500ООО-—1 : 1 000 000), предназначенной для ориентировочного обоснования гидрогеологических предпосылок развития мелио - раций.
Применительно к рассматриваемой детальности районирования могут быть даны только самые общие представления о геофильтрационных схемах. В частности, целесообразно выделить три случая: а — наличие в разрезе над региональным водоупором хорошо проницаемых водоносных пластов с проводимостью более 20 м2/сут; b — наличие в разрезе над региональным водоупором только относительно слабопроницаемых пород с суммарной проводимостью от 2—5 до 20 м2/сут; с—наличие в разрезе мощной толщи весьма слабопроницаемых пород, коэффициент фильтрации которых соизмерим с инфильтрационным питанием при поливах (менее 0,01—• 0,001 м/сут).
В случаях «а» и «Ь» возхможны многообразные геофильтрационные схемы, поэтому следует дополнить эти обозначения следующей индексацией: а' — суммарная проводимость 20—100 м2/сут; а"— „ „ 100—500 „ ;
А!" — „ „ более 500 „ ;
Аг — водоносный пласт представлен трещиноватыми породами;
А0 — нет слабопроницаемых покровных отложений (с поверхности
Хорошо проницаемый водоносный пласт); аг — в верхней части разреза могут быть встречены относительно водоупорные прослои глин или тяжелых суглинков, осложняющие действие горизонтального дренажа; Ь' — мощность слабопроницаемых отложений менее 50 м; Ь" — то же, более 50 м; Ьч и 6Г — то же, что и ат и аг.
Например, водоносный пласт в трещиноватых породах, проводимость меняется в широких пределах, в покровном слое суглинков могут быть встречены прослои глин с k<\0~2 м/сут —
Или водоносный пласт песчано-гравийных пород, залегающих с поверхности, проводимость 200—400 м2/сут — а"; или толща слабопроницаемых трещиноватых песчаников мощностью 80 м (подстилаются глинами), в верхней части разреза нередко встречаются прослои водоупорных глин — Ь" .
Используя данные геофильтрационной схематизации и другие сведения мелиоративно-гидрогеологического районирования (глубины до воды, сведения о подземном стоке), можно в первом приближении судить о прогнозных изменениях в положении уровней подземных вод и в мелиоративной обстановке. Наряду с этим можно следующим образом выделить рациональные типы дренажа в различных геофильтрационных условиях:
А' — горизонтальный и комбинированный, а"-'" — вертикальный и комбинированный, Ь — горизонтальный.
При более детальной разработке геофильтрационной схематизации для Схемы генерального плана развития мелиорации часто возникает необходимость более определенной количественной информации для создания доказательных представлений о наиболее рациональных технических решениях. В этих случаях в зависимости от конкретных задач предпроектных проработок требуется как более надежное построение региональных прогнозов и выбор наиболее благоприятных с гидрогеологической точки зрения объектов первоочередного развития мелиорации, так и обоснование начальных геофильтрационных схем для приближенных расчетов сооружений мелиоративных систем. Вместе с тем методика геофильтрационной схематизации наиболее полно может быть представлена для конкретного проектирования, когда ставится определенная задача количественного обоснования мелиоративных расчетов.
В этом случае результаты геофильтрационной схематизации удобнее всего выносить на специальные карты, отражающие основные закономерности строения массива применительно к обоснованию конкретных инженерных решений.
Как уже отмечалось, для построения локальных расчетных схем систематического дренажа целесообразно прежде всего выделить участки массива, однотипные по развитию водоносных пластов. В связи с этим наиболее крупной таксономической единицей при рассматриваемом картировании является фильтрационная система
водоносных пластов. Исходя из реальных особенностей строения водоносного комплекса, можно выделить следующие четыре основных типа фильтрационных систем: I — однопластовая система (один водоносный пласт однородного или двухслойного строения); II — двухпластовая система (два водоносных пласта разделены слабопроницаемым слоем, верхний пласт однородного или двухслойного строения); III — «беспластовая» система (неоднородная или сравнительно однородная толща слабопроницаемых пород); IV — водоупорная система (развитие в разрезе относительно водоупорных пород).
В реальных условиях встречаются и многопластовые фильтрационные системы, представленные чередующимися в разрезе водоносными пластами и слабопроницаемыми слоями, однако при построении расчетных моделей мелиоративного дренажа такая система может быть схематизирована и реализована как двухпластовая или однопластовая [29]. Анализ результатов такого районирования позволяет судить о целесообразности применения различных типов дренажа. Так, в случае I и II систем рационально применять систематический вертикальный дренаж. При двух и более водоносных пластах следует сопоставлять несколько вариантов размещения в разрезе фильтров скважин вертикального дренажа. На участках развития систем III наиболее эффективен горизонтальный дренаж (реже комбинированный). Система IV характеризует случаи, когда искусственный дренаж неприменим.
Рис. 2. Пример картирования водоносных пластов. |
А — карта, б — литологический разрез; 1 — отроги хребта; 2 — граница конуса выноса: 3 — границы фильтрационных схем; 4 — границы типовых схем |
Рассмотрим в качестве примера один из конусов выноса предгорного шлейфа Туркестанского хребта в Джизакской степи (рис. 2). Здесь отмечается характерная для флювиальных отложе-
Водоносных пластов
Ний предгорных районов дифференциация обломочного материала в направлении от горных сооружений к базису эрозии. В верхней части конуса выноса развиты мощные толщи галечников, которые по направлению сноса постепенно расслаиваются глинистыми породами. На периферии предгорного шлейфа последние преобладают в разрезе (рис. 3).
На начальных стадиях исследования территории наряду со схематизацией по мор- фогенетическим признакам целесообразно прибегать к количественным обобщениям, позволяющим выявить общие закономерности в изменениях основных показателей строения и гидрогеологических условиях территории. В частности, для Джизакского массива характерно закономерное уменьшение суммарной мощности хорошопроницаемых слоев в направлении сноса материала. На рис. 4 показан график изменения отношения мощности песчаных слоев к общей вскрытой мощности, построенный на основании обобщения материала по 40 скважинам, находящимся на различных профилях. Количество хорошопроницаемых прослоев в разрезе определяется законом распределения Пуассона. В частности, вероятность Рп встречи п прослоев в среднем может быть определена по формуле
Р„ = 0,09-М!, (Ы)
Расчетные вероятности встречи того или иного числа горизонтов
В разрезе следующие:
П. " 0 1 2 3 4 5 6 Ра. 0,09 0,21 0,26 0,2 0,13 0,05 0
Как видно из приведенных данных, для большей части площади территории (более 60%) характерно распространение от 1 до
Т
Рис. 4. Изменение относительной мощности песчаиых прослоев вдоль склона (точки, обозначенные кружками, треугольниками и квадратами, соответствуют различным профилям) |
3 слоев хорошопроницаемых прослоев. На 10 % площади хорошо - проницаемые слои до глубины 200 м вообще не встречены. Причем большая их часть (более 60 %) сосредоточена в верхней части разреза до глубины 50 м.
Характер распределения слоев в разрезе виден из графика кумулятивной кривой вероятности встреч прослоев на различных глубинах (рис. 5). Этот график позволяет правильно поставить дальнейшие изыскания, ограничив их глубиной 100 м и сосредоточив главным образом в верхней 50-метровой толще.
Характерной особенностью строения предгорий Туркестанского хребта является широтное простирание генетических комплексов, формировавшихся в процессе сноса материала. Эти процессы определяют масштаб геологической неоднородности, выразившейся
Кт(х) |
Рис. 5. Кумулятивная кривая вероятности встречи скважиной хорошопроницае - мого прослоя на различной глубине
Рис. 6. Корреляционная функция изменения суммарной мощности песчаных слов в зависимости от расстояния.
1 — при направлении вдоль склона; 2 — поперек склона
в чередовании зон, для которых характерны большие мощности хорошопроницаемых пород, с зонами, где они практически не встречаются. Зоны с достаточно развитыми линзами и прослоями песка имеют характер русел, вытянутых в направлении сноса материала. Эта общая тенденция хорошо видна из графиков корреляционной функции, построенных для суммарной мощности песчаных прослоев (рис. 6).
Корреляционная функция, показывающая характер изменения мощности песчаных слоев вдоль склона, имеет плавный характер и указывает на хорошую коррелированность разрезов по скважинам, расстояние между которыми около 5 км. Корреляционная функция, характеризующая изменение мощности поперек склона, имеет ярко выраженную периодичность и указывает на слабую коррелированность разрезов по скважинам, удаленным друг от друга уже на 2 км. Полученные характеристики позволяют обоснованно подойти к построению карты суммарной водопроводимости, необходимой при решении региональных задач. При этом целесообразно приведение сложной толщи к одному водоносному горизонту, эквивалентному ей по проводимости Гэ
Гэ = г1 + Г2+... +Та, (1.2)
Где Т\, Т2, Тп — проводимости различных водоносных горизонтов. Этот водоносный горизонт (назовем его «эквивалентным») будет иметь напор #э
//э_ Н1ті + Н2т2+ ... +нптп ^ {х
Т э
Где Ні, #2, Нп — напоры в соответствующих горизонтах.
При решении региональных задач линейные размеры рассматриваемой области, как правило, значительно превышают величину фактора перетекания, определяемого сопротивлением слабопроницаемых слоев, что позволяет считать Нс = Подобное допущение оказывается справедливым и при рассмотрении естественного потока. В области транзита напоры в различных горизонтах отличаются незначительно и в большинстве случаев соответствуют отметкам свободной поверхности. Последнее обусловлено тем, что верхняя толща отложений конусов выноса (до глубины 40—50 м) является основной водопроводимой частью разреза. Таким образом, величина напора Нэ может быть получена либо непосредственно по скважине, вскрывающей всю толщу отложений, либо с некоторыми допущениями по скважинам, вскрывающим только свободную поверхность. В последнем случае необходимо иметь достаточно данных, обосновывающих возможность принятия в качестве эквивалентного напора напор на свободной поверхности. Значения эквивалентной проводимости Тв можно получить путем численного анализа напорной поверхности эквивалентного пласта исходя из соотношения
(І.4)
Ч
Где Ті и іі — водопроводимость и градиент на г-м участке линии тока. Для получения абсолютных значений Т необходимо иметь данные опытных откачек, позволяющие непосредственно рассчитать локальные значения Т и распространить их на всю рассматриваемую область. Подобного рода построения профиля проводимости можно проводить только по линии тока естественного потока при отсутствии на рассматриваемом участке испарения или питания. Карта проводимости эквивалентного пласта строится интерполяцией между отдельными профилями. Такая карта была составлена для той части конуса выноса, где отсутствует разгрузка подземных вод за счет испарения. Для построения была использована карта гидроизогипс, отражающая характер свободной поверхности потока. Предварительно был проведен корреляционный анализ, позволяющий судить о соотношении градиентов свободной и пьезометрической поверхностей эквивалентного пласта на тех участках, где она известна. Результаты анализа показали полное соответствие (коэффициент корреляции 0,95) градиентов. Изучение распределения напоров в различных водоносных горизонтах на электрической модели выявило, что существенное различие (более 1 — 2 м) в них наблюдается только в том случае, когда коэффициент фильтрации разделяющих слоев меньше 10~4 м/сут.
Изучение региональных схем потока, помимо рассмотренных выше карт параметров, требует обоснования граничных условий рассматриваемой территории. Верхние области конуса выноса являются областями с обеспеченным питанием, что дает возможность принять здесь границу с постоянным расходом. Границы, идущие вдоль склона, совпадают с линиями тока естественного потока, изменение направления которого при орошении незначительно. В связи с этим их можно рассматривать как непроницаемые.
В нижней части конуса выноса происходит разгрузка подземных вод за счет испарения со свободной поверхности потока там, где она находится на глубине меньше критической. Здесь в качестве граничного условия нужно задавать условие третьего рода, определяемое характером испарения и глубиной залегания свободной поверхности подземных вод. Это условие должно задаваться на всей площади, где в естественных условиях происходит испарение. В тех областях, где в результате работы дренажа уровень понизился на глубину, большую критической, граничное условие меняется на условие второго рода, а при решении задачи сложением решений, питание, задаваемое в этом случае, должно соответствовать интенсивности испарения в естественных условиях.
Сложность гидрогеологических условий и их разнообразие обусловливают необходимость по-разному подходить к обоснованию дренажа в различных частях конуса выноса. Следуя закономерностям геологического строения и гидрогеологических условий в пределах рассматриваемого конуса выноса, можно выделить три района. I район характеризуется глубоким (до 50 м) залеганием свободной поверхности, которая проходит в отложениях галечника с высокими фильтрационными свойствами (коэффициент фильтрации до 30 м/сут). В этих условиях дренаж может происходить либо естественным путем, имея в виду, что при подъеме уровня проводимость водоносного горизонта увеличится почти в три раза, либо рядом скважин (заградительный дренаж). II район охватывает область средних частей конуса выноса. Свободная поверхность здесь находится на незначительной, в некоторых случаях близкой к критической, глубине и проходит в отложениях слабопроницаемых супесей и суглинков. Проводимость эквивалентного пласта в этой области от 100 до 500 м2/еут. Участки орошения в этом районе можно рассматривать изолированно с учетом дополнительного притока из первого района. Наиболее сложной представляется работа дренажа в нижних частях конуса выноса (III район). Здесь уровни подземных вод уже в естественном состоянии находятся на глубине выше критической, что обусловливает интенсивное испарение и, таким образом, разгрузку всего комплекса водоносных горизонтов. Эта область характеризуется преобладанием в разрезе суглинков и супесей. Песчаные горизонты встречаются редко и обладают незначительной мощностью (до 5 м). Суммарная проводимость в этой области уменьшается до 50 м2/сут.
Анализ результатов геофильтрационной схематизации показывает, что для системы 1-а (см. рис. 2, 3) целесообразно прежде всего рассмотреть вариант заградительного вертикального дренажа; для систем 1-6 и II-6, в—различные формы систематического дренажа; для системы ИІ-а — горизонтальный и комбинированный дренаж, а для системы Ш-б — преимущественно горизонтальный дренаж.
Следует отметить, что в природных условиях нередко встречаются водоносные комплексы с тонкослоистым строением. В таких случаях формальный подход к фильтрационной схематизации может привести к существенным просчетам в оценке параметров дренажа. Это положение прежде всего относится к обоснованию гидродинамической роли прослоев хорошо - и слабопроницаемых пород при определении фильтрационной системы. Необходимость дифференцированного рассмотрения таких прослоев в расчетной схеме для наиболее простых условий может быть установлена по предварительным расчетам, а в более сложных случаях — решением тестовых задач на аналоговых моделях.
Особенности строения покровных отложений влияют главным образом на гидродинамическое несовершенство горизонтального дренажа, которое существенно зависит от проницаемости этих отложений и сравнительно слабо от мощности. Поэтому основным параметром покровных отложений, который должен отражаться при картировании, является проницаемость слоя под уровнем заложения горизонтального дренажа. С точки зрения эффективности дренажа можно выделить три типа проницаемости покровных отложений: легкие, средние и тяжелые, причем в тяжелых практически исключается эффективная работа горизонтального дренажа, в средних требуется расчетное обоснование выбора типа дренажа, а в легких горизонтальный дренаж обычно оказывается предпочтительным.
Для оценки перечисленных условий целесообразно составлять карту-срез поверхности, соответствующей уровню заложения горизонтального дренажа с выделением зон тяжелых, средних и легких отложений. В зоне развития легких отложений можно, кроме того, показать участки, где покровный слой отсутствует или дрена полностью его прорезает, В остальных случаях целесообразно выделять слои более проницаемых пород над уровнем заложения дренажа и слои менее проницаемых пород*—ниже этого уровня. На рис. 7 приводится пример картирования покровных отложений одного из участков совхоза № 13 в Каршинской степи (дельтовая равнина), для которого характерно разнообразное строение покровных отложений (рис. 8).
При оценке условий просачивания инфильтрационных вод необходимо иметь представление о наличии в покровных отложениях слоев с низкой проницаемостью (<0,001 м/сут), поскольку наличие таких горизонтов в почвенном слое ставит под сомнение возможность орошения этих земель, а присутствие их в нижней части разреза затрудняет применение вертикального дренажа. В ряде случаев целесообразно составлять отдельную карту распространения таких слоев с указанием глубины их залегания.
При обобщении гидрогеологических предпосылок применения различных типов систематического дренажа синтезируются результаты картирования водоносных
Рис. 7. Картирование покровных отложений в совхозе № 13 Каршинской степи. / — глины и тяжелые суглинки; 2- легкие и средние суглинки; 3 — супеси; 4 — пески основйого водоносного пласта; 5—-скважины; 6 — изолинии мощности покровных отложений |
Рис. 8. Схемы строения покровной толщи в совхозе Jsfe 13 Каршинской степи в зависимости от условий дренирования.
А — тяжелые; Б — средине; В — легкие отложения
Систем I и II, содержащих водоносные пласты, следует рассмотреть ряд случаев. Так, при небольшой мощности покровных отложений горизонтальный дренаж, прорезающий покровный слой, может оказаться наиболее эффективным. В случае низкой проницаемости и существенной мощности покровных отложений нужно рассмотреть варианты применения комбинированного или вертикального дренажа, а при средней и высокой проницаемости — все отмеченные типы дренажа. Участки с различными рекомендуемыми формами дренажа, установленные в результате картирования, могут занимать сравнительно небольшие площади. Часто существенные различия в условиях осуществления дренажа определяются плановой неоднородностью строения покровных отложений, которая может проявляться значительно резче, чем изменения проводимости подстилающего водоносного пласта. Например, для условий совхоза № 13 в Каршинской степи (см. рис. 8) расчетные расстояния между горизонтальными дренами меняются от 20 до 800 м и зависят от величины фильтрационного сопротивления покровных отложений. Последние могут существенно изменяться на сравнительно небольших расстояниях.
Следует отметить, что эффективность горизонтального дренажа в описываемых условиях практически не зависит от строения покровных отложений между дренами, так как при значении kB, превышающем интенсивность инфильтрации, потери напора на вертикальное просачивание несущественны. Поэтому учет особенностей строения покровных отложений целесообразно проводить только по линиям дренажа. Использование карты проницаемости покровных отложений наряду с топоосновой при проектировании позволяет выбрать наилучший вариант размещения дренажной сети в плане. В указанных условиях одним из вариантов такого решения задачи может быть применение горизонтального дренажа с дополнительным введением скважин-усилителей (комбинированный дренаж) на участках покровных отложений низкой и средней проницаемости. Необходимо отметить, что в таких условиях система комбинированного дренажа оказывается не только эффективной, но и достаточно мобильной с точки зрения корректирования принятых решений в процессе рабочего проектирования и мелиоративного строительства. В этом случае более эффективное действие системы достигается устройством дополнительных скважин-усилителей при неизменных расстояниях между линиями дренажа.
В ряде осушительных систем (особенно в тяжелых минеральных почвах или при близком водоупоре) или при реконструкции мелиоративных систем методика геофильтрационной схематизации может быть существенно модифицирована в зависимости от конкретных природных условий и поставленных задач.
Тем не менее рассмотренные общие принципы и методика схематизации природных условий позволяют подойти целенаправленно как к обоснованию инженерных решений, так и к постановке изысканий на объектах мелиорации.