Определение миграционных параметров в полевых условиях проводятся посредством запуска индикатора в пласт через опытные скважины с последующим его улавливанием по соседним скважинам. Опыты можно проводить как в естественном, так и в нарушенном (путем откачки ил инагнетания) потоке, причем наиболее достоверны эксперименты в возмущенном потоке. Они отличаются более четкой фильтрационной картиной, позволяют увеличить размеры зоны опробования и тем самым снизить влияние масштабных эффектов. Кроме того, важно, что они позволяют одновременно изучать фильтрационные характеристики участка, необходимые для интерпретации индикаторного опробования. Однако и эти эксперименты не свободны от серьезных недостатков, нередко сильно ограничивающих их ценность.
Говоря ранее об определении фильтрационных параметров, мы подчеркивали, что в подавляющем большинстве случаев их надежная оценка возможна лишь на основе полевых опытных работ. Существенно иначе обстоит дело с миграционными параметрами, когда речь идет об их определении для комплексов пористых песчано-глинистых пород. Здесь полевые работы часто либо вообще неэффективны, либо не имеют ощутимых преимуществ в сравнении с лабораторными методами. В частности, опыт показывает, что активная пористость песков, определяемая опытными запусками индикатора в пласт, оказывается, как правило, резко заниженной (подчас — в несколько раз). Объясняется это в первую очередь фильтрационной неоднородностью реальных водоносных пластов.
Рассмотрим, например, двухслойный (слои 1 и 2) пласт, изображенный на рис. 6.19. Пусть проницаемости к1 и к2 различаются в пять раз, а мощности находятся в обратном отношении. Значения активной пористости (п) будем считать равными. Индикатор подается в виде «пакета» вместе с водой, нагнетаемой в скважину А, и улавливается в скважине Б. Первая «волна» индикатора будет принесена в скважину Б по более проницаемому слою 1 через время tl * Ln/v{, откуда сразу определяется активная пористость п, если задана скорость фильтрации этого слоя v[. На деле, однако, обычно известна усредненная по всему пласту скорость фильтрации, которую определяют исходя из решения соответствующей плановой задачи. Поэтому вместо действительного значения п, пропорционального величине будет получено расчетной значение п , пропорциональное усредненному коэффициенту фильтрации к - (kjirij + к2т2)/(ntj + т2). Следовательно,
ki l1i+^i ml+m2 mx+m2 k{
n
в нашей задаче составит 1/3. Итак, получаемое значение активной пористости занижено в три раза. Из формулы (6.64) также видно, что наличие тонких хорошо проницаемых прослоев (mt « т?) может приводить к занижению активной пористости примерно в ку/к2 раз (при близких значениях п, и п2)-
|
А |
/////// |
/////////// |
|
|
О о о о |
о о о о о о с |
ООО |
о о о о о о |
Рис. 6.19. Схема опробования двухслойного водоносного пласта
о
Б
' ‘
-7^с
Таким образом, для расшифровки результатов проведенного эксперимента необходима послойная оценка проницаемости, что, однако, практически исключается при общепринятой методике Проведения опытно-фильтрационных работ.
Важно, что профильная фильтрациоынная неоднородность — отнюдь не единственная причина больших погрешностей в определяемых значениях миграционных параметров. С такими погрешностями приходится сталкиваться и при индикаторном опробовании даже однородных песков. Объясняется это тем, что за время опробования, при относительно высоких скоростях фильтрации, индикатор проникает лишь в ограниченную часть порово- го пространства (см. раздел 6.1.1).
С другой стороны, напомним, что в песчано-глинистых грунтах определение коэффициента микродисперсии обычно лишено практического смысла (см. раздел 6.3), а основной прогнозный параметр — активная пористость, близкая к общей пористости, — вполне уверенно определяется лабораторными экспериментами и даже по аналогии. Все это подтверждает весьма ограниченную ценность полевых опытно-миграционных работ (ОМР) в толщах пористых песчано-глинистых пород.
Совсем другое дело — трещиноватые породы. С одной стороны, два основных параметра — активная трещиноватость и коэффициент микродисперсии — практически неопределяемы в лабораторных условиях. С другой стороны, именно трещиноватые породы являются наиболее быстрыми транспортерами вещества в подземных водах. Все это делает полевые опыты с индикаторами необходимым элементом изучения миграционных свойств трещиноватых горных пород.
Опытные полевые скважины лучше всего приурочивать к участкам проведения кустовых откачек, с тем чтобы считать фильтрационные параметры пласта известными. Кроме того, при этом можно использовать часть скважин фильтрационного куста. Нужно, однако, иметь в виду, что требования к расположению и оборудованию скважин для миграционного опробования могут существенно отличаться от требований к опытно-фильтрационным работам. Прежде всего это относится к масштабам опробования: даже в трещиноватых породах, характеризующихся относительно большими действительными скоростями фильтрации, область влияния эксперимента обычно ограничивается первыми десятками метров (при продолжительности опыта, измеряемой первыми сутками). С другой стороны, ОМР требуют очень детальной расшифровки профильной фильтрационной неоднородности опробуемого разреза (см. раздел 6.4). Для этого можно с успехом использовать гидрогеофизические методы, особенно расходометрию и резистивиметрию. Обнаруживая с их помощью наиболее фильтрующие зоны, можно затем участь реальную фильтрационную неоднородность и при выделении интервалов опробования, и при интерпретации опыта.
Рассмотрим теперь возможные схемы собственно индикаторного поробования в нарушенном фильтрационном потоке, когда индикатор запускают в условиях режима фильтрации, близкого к стационарному (квазистацио-
нарному).
ЗАДАНИЕ. Покажите, что последнее условие нетрудно выполнить на практике, ибо скорость распространения фильтрационных возмущений существенно выше, чем у индикаторных возмущений. Используйте для этого формулу (4.18) в сопоставлении с характерными масштабами конвективного переноса (за одно и то же расчетное время t).
При проведении ОМР в точке ввода индикатора в водоносный пласт создается один из следующих концентрационных режимов:
поддержание ее неизменного значения в период всего
опыта;
L2J пакетный ввод — поддержание постоянной концентрации только в течение определенного времени запуска;
От типа входного индикаторного сигнала в значительной степени зависят информативность опыта по отношению к определяемым миграционным параметрам, а также простота и доступность обработки кривых «отклика». В частности, при применении пакетного и, особенно, импульсного режимов ввода в условиях больших разбавлений индикатора, а также в средах с высокими рассеивающими (например, в трещиновато-пористых породах) или поглощающими свойствами, надежность опыта резко падает и существенно ограничивается чувствительностью способов индикации.
По типу фильтрационного возмущения возможны различные схемы миграционного опробования. В некоторых случаях опыт проводят в фильтрационном потоке, создаваемом откачивающей скважиной, когда запуск индикатора осуществляется в близрасположенную наблюдательную скважину, а регистрация — в откачиваемой воде. При этом можно использовать для опыта эксплуатационные дренажные скважины рудников, водозаборные скважины хозяйственного назначения и т.п. Однако такая экспериментальная схема обладает рядом недостатков, из которых основными являются сильное разбавление индикатора в откачиваемой воде и весьма малое время для его регистрации (при обычно применяемых пакетном и импульсном запусках).
ЗАДАЧА. Покажите, что в этом опыте коэффициент разбавления имеет порядок —г- (г — удаление наблюдательной скважины от
с
центральной, dc — даметр наблюдательной скважины). Не забудьте при этом учесть коэффициент искажения потока вблизи наблюдательной скважины (см. раздел 6.1.4).
Учитывая отмеченные недостатки опытов с откачивающей скважиной, наиболее широко стали проводить эксперименты по схеме нагнетания (налива) индикаторного раствора в центральную скважину с последующим прослеживанием его распространения в пласте по наблюдательным скважинам. Интерпретация таких опытов ведется по аналитическим зависимостям для плоскорадиальной миграции, которые получают аналогично приведенным ранее зависимостям для конвективно-дисперсионного массопереноса (см. разделы 6.3 и 6.4) или конвективно- кондуктивного теплопереноса (см. раздел 6.5) — в случае использования теплового индикатора. При этом можно применять различные физические модели, отвечающие конкретным условиям эксперимента. Так, запуск солевого индикатора в пласт трещиноватых пород описывается расчетной схемой микродисперсии, т.е. решением, подобным формуле (6.27). Более того, можно использовать это решение и непосредственно, аналогично интерпрета-
ции лабораторных экспериментов (см. раздел 6.6), если заменить в нем скорость фильтрации v(г) условной средней величиной
v =2 v(r) = . _
р w Ttrm (6.65)
где г — расстояние от центральной скважины (работающей с постоянным расходом Q) до наблюдательной;
т— мощность пласта.
При использовании солевого индикатора в трещиновато-пористых породах опробование обычно протекает в рамках расчетной схемы неограниченной емкости (см. раздел 6.4), в чем нетрудно убедиться с помощью критерия (6.49). Для интерпретации поэтому можно использовать решение (6.47), заменяя в нем х на г, a v — на vp согласно формуле (6.65).
При всех достоинствах у рассмотренной схемы опробования имеются и свои дефекты - резкая неравномерность поля скоростей фильтрации с быстрым их падением по мере удаления от центральной скважины, а также искажающая роль естественного потока, влияющего на перенос индикатора в краевых частях опробуемой зоны. Приходится, кроме того, считаться и с техническим несовершенством систем слежения за индикаторной меткой в пласте по наблюдательным скважинам (см. разд. 7.3.2).
Улучшение ситуации с этих позиций может быть достигнуто проведением опыта по двухскважинной (дуплет- ной) схеме, когда индикатор запускается в нагнетательную скважину , а его регистрация производится в воде откачивающей скважины. Тем самым резко сокращается требуемый для опыта объем бурения. Дуплетная схема весьма эффективна при опробовании глубоко залегающих водоносных горизонтов, когда оперативный отбор гидрохимических проб связан со значительными техниче- скими трудностями.
* Благодаря «замыканию* всех меченых индикатором линий тока на ('укачивающей скважине, при дуплетном опробовании обеспечивается наилучшее, в определенном смысле, осреднение характеристик изучаемой среды в пределах
области влияния опыта.
Более подробно вопросы постановки и интерпретации ОМР рассмотрены в разделе 7.3.
На этом, собственно, изложение гидродинамических основ теории тепло- и массопереноса можно считать законченным. Напомним, однако, что мы договорились включить в эту главу и материал, касающийся еще одной специальной проблемы, — влагопереноса в ненасыщенных грунтах. Ей будут посвящены заключительные разделы главы.
Рассмотрение терии влагопереноса имеет первостепенное значение для анализа динамики влаги в зоне аэрации при гидрогеологическом изучении ряда важных практических проблем. Среди них:
1
изучение условии питания подземных вод;
прогнозная оценка связи подземных вод с вновь создаваемыми водоемами или водотоками;
|з] обоснование ч:хем искусственного пополнения подземных вод;
изучение «защитных» свойств пород зоны аэра-
ции с позиции возможного загрязнения подземных вод;
гидрогеологические исследования на мелиорируемых территориях;
[б~] изучение режима влажности грунтов в связи с решением различных задач инженерной геологии (имея в виду прежде всего, что от этого режима во многом зависят изменения прочностных и деформационных характеристик грунтов).
Значимость этих проблем и специфика изучаемых процессов требуют, вообще говоря, выделения теоретических основ влагопереноса в специальную главу нашего курса. Мы, однако, не будем этого делать, и прежде всего потому, что теория здесь по ряду объективных причин пока слишком слабо увязывается с практикой современных гидрогеологических исследований. Соответственно, мы ограничимся в данной книге лишь самыми оОщими положениями этой теории.
В отличие от фильтрации в водонасыщенных породах, для влагопереноса в зоне аэрации характерно наличие в поровом пространстве горных пород двух фаз — воды и воздуха1, причем последний, находясь в свободном состоянии *, образует непрерывную фазу и обычно имеет прямую связь с атмосферным воздухом. При этом соотношение фаз, т.е. степень водонасьпцения, является тем определяющим фактором в процессе влагопереноса, от которого зависят характер и интенсивность внутренних сил, действующих на частицы жидкости в порах. Чем меньше степень насыщения данного грунта, тем сложнее в целом комплекс этих сил и, соответственно, механизм переноса влаги: все большую роль начинает играть перемещение парообразной воды, — наряду с другими факторами, взаимосвязанными с массовым потоком (температурными, электрокинетическими, электромагнитными, химическими). Однако в условиях пород со сравнительно высокой влажностью (большей максимальной молекулярной вла- гоемкости), обычно характерной для зоны аэрации, главное значение сохраняется за силами гравитации, а также капиллярными и сорбционными силами . Так как в пределах каждого равномерно насыщенного элементарного объема грунта (рассматриваемого изолированно от смежных объемов) последние две силы направлены на удержание воды в порах, но их суммарное влияние принято характеризовать сорбционно-капиллярным потенциалом поровой влаги — работой необходимой для извлечения из фунта некоторой (например, единичной) массы
Ц1ШШПШГЕ ' приходящийся на еди-
J ницу объема воды, —
•j жидкости. Практически измеряется не ра- о бота, а ее эквивалент,
о
давление, которое принято называть всасывающим давлением: если отнести его к единице веса воды и, по аналогии, с
пьезометрическои высотой, ввести высоту всасывания
Ч*, то она будет выражаться в метрах водя-
Рис. 6.20. Схематическое представле- ного столба. Схема- ние измерения всасывающего давления: тически такой ВОДЯ- 1 - образец; 2 - тонкопористая мембрана „ * *