Зададим себе элементарный вопрос: что нужно, чтобы открытая наблюдательная скважина зафиксировала понижение S в пласте? Ответ прост: необходимо, чтобы из
скважины вытек объем воды V — nr^S. Если поступление
этой воды в пласт затруднено сравнительно невысокой проницаемостью пород, то уровень воды в скважине будет снижаться с некоторым отставанием от пьезометрической поверхности пласта, т.е. наша измерительная система - скважина оказывается инерционной. Очевидно, значение инерционности максимально на ранних этапах понижения напора в данной точке пласта, когда скорости изменения уровня максимальны, причем время существенного ее проявления tn будет расти с увеличением поперечного сечения (ос скважины и с падением проницаемости водоносных пород. Так, для совершенных скважин погрешность, обусловленная инерционностью, оказывается меньше 10% при времени [23]
0)
f. ~ (15+20) (5 14)
ВОПРОС. Как влияет на инерционность длина фильтра несовершенного пьезометра?
Отсюда нетрудно видеть, что в типичных супесчаных и суглинистых грунтах, например, измерения понижений оказываются сильно искаженными инерционностью скважин в течение нескольких часов. Тот же порядок имеет время tn и для скважин в песчаных или трещиноватых породах, прифильтровая зона которых существенно затронута процессами кольматации (см. раздел 5.4.4).
Для учета влияния инерционности необходима предварительная тарировка наблюдательных скважин, для чего проще всего использовать кратковременные наливы воды в скважину [23 ]. В слабопроницаемых суглинистых и глинистых грунтах правильнее стремиться к
ный в координатах IgS, lg(t/r2), то эти графики должны иметь одинаковую форму. Параллельным перемещением по осям координат можно добиться совпадения графиков, причем они будут сдвинуты по вертикали на величину
lg -г^-Ьр , а по горизонтали — на величину lg(4a*), как
это1 иллюстрируется рис. 5.17. Снимая с совмещенного графика эти величины, определяем параметры Гид*.
Рис. 5.17. Совмещение опытных точек ( показаны кружками ) с эталонной кривой
Способ прямой линии, ввиду его особой важности, уже освещен в разделе 5.1 применительно к временному прослеживанию. Это, в частности, основной способ обработки данных одиночных откачек (как правило, обрабатывается кривая восстановления уровня — см. раздел
. Подчеркнем, что он годится только для времени- наступления квазистационарного режима, т.е. при выполнении условия (4.29). При этом, следовательно, упускается информация о начальном участке кривой прослеживания. Используя структуру формулы (4.30), нетрудно показать возможности способа прямой линии и для пространственного прослеживания — в координатах S — Igr, что для двух наблюдательных скважин, находящихся в зоне квазистационарного режима (см. условие (4.29)), равносильно применению формулы (4.31). Наконец, при комбинированном прослеживании график спрямляется в координатах S — lg(t/г2).
Точечные способы основаны на использовании отдельных замеров, либо никак не объединяемых друг с другом в интерпретационной схеме, либо увязанных, согласно той или иной аналитической зависимости, еще с одной замерной точкой (так, только что упомянутый расчет по формуле Дюпюи может рассматриваться и как двухточечный способ).
Способы, основанные на интегральных преобразованиях,чаще всего используют операционный метод (см. раздел 4.2). Для этого по замеренным значениям понижений в наблюдательной скважине вычисляются их изображения У для ряда значений параметра преобразования tp (см. формулы (4.53) и (4.53а)), после чего для интерпретации используется соответствующее решение в изображениях, рассматриваемое как график связи $=f( tp). Если при этом справедлива, например, логарифмическая связь 5 и t (см., в частности, формулу (4.57а), то используется соответствующая модификация способа прямой. При данном способе обработки легко учитывается непостоянство расхода откачки (см. формулы (4.56) и (4.57)).
Способ совмещения кривых прослеживания предполагает последовательный подбор параметров, исходя из требования совмещения фактического и расчетного графиков временного прослеживания. По этому способу, определив предварительно параметры, строят расчетную кривую прослеживания и сопоставляют ее с фактической, после чего, учитывая характер зависимости хода процесса от параметров, их изменяют, добиваясь наилучшего совмещения расчетной и фактической кривых; при этом могут эффективно использоваться специально ориентированные алгоритмы для «автоматического» подбора параметров на ЭВМ.
Сопоставительная оценка способов обработки. Все упомянутые здесь способы обработки являются равноценными, если принятая расчетная схема точно отражает природные и технические условия эксперимента, а замеры понижений и расходов при откачке сделаны без погрешностей. Поскольку на практике всегда бывает иначе, разные способы оказываются существенно различающимися с точки зрения их чувствительности к погрешностям расчетной схемы и замеров, а следовательно, и надежность получаемых по ним параметров может оказаться существенно различной.
С этих позиций ясно, что способы, объединяющие информацию, оказываются, при прочих равных условиях, эффективней, чем точечные способы, если, конечно, при такой свертке не допускается больших дополнительных погрешностей (усреднения, интерполяции, экстраполяции и т.д.). Последнее требование обычно справедливо для временного прослеживания и относительно реже выполняется при площадном прослеживании - ввиду ограниченного числа замерных точек (наблюдательных скважин).
При относительно хорошем соответствии расчетной схемы условиям эксперимента наиболее удобным из способов, свертывающих информацию, является способ прямой линии для временного прослеживания (в обычной или операционной модификации), который широко и успешно используется на практике. В частности, график S = / (In t) служит важнейшим основанием для проведения диагностики откачек (см. раздел 5.5.2). При большом числе наблюдательных скважин способ прямой линии может эффективно применяться при комбинированном (площадном и временном) прослеживании.
Широкое распространение в практических расчетах получили также различные модификации способа эталонной кривой, позволяющего провести обработку и диагностику данных по всей имеющейся информации временного или комбинированного прослеживания (включая начальный этап нестационарного режима). Однако совмещение опытной кривой с эталонной всегда привносит в расчет повышенную долю субъективизма.
Способ совмещения кривых прослеживания можно применять для обработки крупномасшатбных откачек в сложных геофильтрационных условиях.
Таким образом, выбор эффективного способа обработки ОФР — задача не из легких. Поэтому попытаемся дать хоть какие-то общие ориентиры для ее решения. Будем исходить здесь из того, что при отсутствии существенного влияния плановых границ пласта или плановой неоднородности (см. раздел 5.3) и несовершенства опытных скважин (см. раздел 5.4) важнейшим показателем качества интерпретации является совпадение значений проводимости, полученных временным и площадным прослеживанием. Отсюда можно рекомендовать такую последовательность расчетных операций при обработке данных кустовой откачки:
|~Г| используя одновременные замеры по двум ближним наблюдательным скважинам на заключительном этапе откачки1, определяют проводимость из формулы (4.31), предполагающей квазистационарный режим:
т - 0,37 Qc IsТг ■ ы S(r,) - s(ri) ri ' (5.16)
при откачках вблизи реки, имеющих конечный стационарный этап, вместо (5.16) используется формула (3.47);
[Т] по графикам временного прослеживания способом прямой линии определяют проводимость Твр. Если ~ Твр, то выполняются следующие этапы расчета (в противном случае приходится искать какие-то важные просчеты в принятой интерпретационной схеме или — для планово-однородных пластов — отнестись с большим до- верием к величине Т^);
оценивают коэффициент пьезопроводнсоти
(уровнепроводности) способом прямой (если результаты
расчетов проводимости этим способом признаны удовлетворительными) или точечным способом;
ПГ| проверяют выполнение условия квазистационарности (4.29) для расчетных диапазонов, принятых в операциях 1 и 2 (если оно не выполнено, то проводят дополнительные уточнения);
~5~| проводят более детальные расчеты, максимально использующие расчетные точки графиков и временного, и площадного прослеживания; при этом могут широко присменяться способы эталонной кривой или способы, основанные на интегральных преобразованиях;
[~6 | достоверность найденных расчетных параметров оценивают путем сопоставления получаемых по ним расчетных графиков с фактическими. Для этого, как уже отмечено, можно эффективно использовать ЭВМ.
В то же время следует очень осторожно относиться к широкой передаче ЭВМ всех функций по интерпретации опыта, особенно с учетом многочисленных сложностей диагностики.
Важнейшим моментом при оценке надежности полученных параметров является анализ чувствительности расчетной модели. Для этого поочередно меняют значения параметров в тех или иных диапазонах и сравнивают полученные таким образом новые расчетные кривые с ранее построенной расчетной (или опытной) кривой. Если большие изменения параметра вызывают лишь малые расхождения в кривых, то это свидетельствует о малой надежности расчетного значения параметра, о слабой чувствительности к нему выбранного расчетного алгоритма. Такой анализ позволяет часто выделить те (более чувствительные) участки опытных кривых, по которым предпочтительнее оценивать данный параметр. В противном случае желательно строить или изменять расчетный алгоритм таким образом, чтобы мало надежные параметры не использовались (как промежуточные данные) для расчета других параметров.
ВОПРОС. В каких точках области влияния откачки — ближних или дальних — одинаковые изменения проводимости окажут более заметный эффект на понижение в центральной скважине? Почему?
ЗАДАНИЕ. Попытайтесь объяснить, исходя из смысла и структуры формул (5.16) и (4.30), почему результаты площадного прослеживания наиболее чувствительны к изменениям проводимости в зоне, охваченной наблюдательными скважинами ; почему, наоборот, результаты временного прослеживания по способу прямой на эти изменения почти не реагируют? (для ответа на эти вопросы еще раз уясните физический смысл понятия квазистационарного режима (см. раздел 4.1.2), которому только и отвечают упомянутые зависимости, лежащие в основе площадного (4.31) и временного (4.30) прослеживания).
В целом анализ чувствительности подтверждает практические представления о том, что даже ОФР, проведенные на высоком качественном уровне, могут давать весьма ощутимые погрешности в определяемых параметрах. Причина этому — различные неучтенные отклонения от расчетной схемы (например, фоновые колебания напоров или площадная изменчивость фильтрационных свойств). Так, для водопроводимости, определенной качественными кустовыми откачками, характерны погрешности порядка 10%, возрастая до нескольких десятков процентов для качественных одиночных опробований. В этой связи нет смысла настаивать на излишне высокой точности алгоритмов интерпретации ОФР, считая во всех случаях вполне допустимыми для них погрешности порядка первых процентов для кустовых откачек и 5-10% — дли одиночных.
Конечной целью диагностики является выбор расчетной интерпретационной схемы, учитывающей с оптимальной полнотой и надежностью природные и технические условия проведения эксперимента. Понятие оптимальности в данном случае одновременно предполагает, что расчетная схема, во-первых, является, по возможно-
Точнее, теми из них, что использованы при прослеживании.
сти, простой, исключающей влияние малозначащих факторов, а во-вторых, находится в соответствии с объемом полученной при эксперименте информации. К сожалению, на практике эти требования к расчетной схеме всегда выполняются лишь с некоторым приближением - как в силу сложности реальных условий, так и ввиду ограниченности исходной информации. Поэтому диагностика опытной откачки может быть формализована лишь в ограниченной степени, решающую роль здесь играют опыт и интуиция специалиста. Вместе с тем, даже опытные специалисты часто допускают просчеты в интерпретации откачек, обусловленные недостаточным пониманием всех тонкостей и многообразных проявлений фильтрационного процесса при откачке, которые вытекают из его детального теоретического анализа. В связи с этим чрезвычайно важно, чтобы основные принципы диагностики опытных откачек строились на сочетании качественного гидрогеологического анализа и количественных методов интерпретации. В то же время, диагностика откачек, как и любая схематизация, должна вестись поэтапно, путем постепенного уточнения интерпретационной схемы, причем условия, обеспечивающие надежность диагностики, следует принимать во внимание уже при постановке и проведении эксперимента (см. раздел 7.1).
Диагностика откачки начинается с качественного анализа условий ее проведения, после чего, исходя из принципа обратной связи, все шире привлекаются и различные элементы количественного анализа. Естественно, исходным элементом диагностики является тщательный геологический анализ. Предполагается, что такой анализ, по возможности, проводится уже при постановке Опытной откачки, существенно предопределяя тем самым и круг возможных интерпретационных схем. Дальнейшая детализация геофильтрационного строения опробуемых пластов базируется на полученных при бурении данных по литологическому и механическому составу пород, по геофизическому и гидрогеологическому опробованию в процессе проходки и оборудования скважин. Все эти данные в сочетании с результатами замеров напоров в процессе и по окончании бурения позволяют построить гео- фильтрационный разрез и в первом приближении наметить типовую расчетную схему. Как правило, это должна быть одна (иногда две) из основных схем, рассмотренных в разделе 5.2.
Вспомогательная информация для предварительных качественных построений может быть получена при анализе и лабораторном исследовании ненарушенных образцов водоносных и относительно водоупорных пород. Так, для выделения относительных водоупоров и предварительной оценки возможного перетекания полезно оценить расход его по известным из лабораторных экспериментов параметрам проницаемости (хотя подобные оценки обычно дают заниженные показатели). Для тонкослоистых фильтрующих толщ результаты испытаний на фильтрацию вдоль й вкрест напластования дают представление о степени профильной анизотропии. Особое значение имеет использование лабораторных данных компрессионных испытаний для представлений о коэффициентах упругоемкости водоносных и водоупорных пород.
Полезную информацию дают также сведения о напорах подземных вод, полученные в процессе и по окончании бурения. Так, разрыв естественных напоров отдельных пластов вдали от койтуров питания (стока) чаще всего безусловно свидетельствует о пренебрежимо малом их взаимодействии и в процессе откачки.
ВОПРОСЫ. 1. Обязательно ли факт противоположного свойства (совпадение значений естественных напоров) является доводом в пользу возможности существенного проявления перетекания при откачке? (будем иметь в виду, что в естественных условиях перетекание могло идти длительное время и на больших площадях. 2. Кай должны меняться с глубиной напоры в слоистой водоносной тблщё вблизи ее контура стока, несовершенного по степени вскрытия толщи?
Подчеркнем еще раз, что большинство из упомянутых здесь вспомогательных диагностических приемов и методов могут и должны использоваться уже при проектировании откачки, способствуя тем самым ее большей целенаправленности и выработке исходных представлений о возможной интерпретационной схеме. Последующее уточнение этой схемы ведется непосредстсвенно по материалам откачки, в первую очередь по составляемым в процессе эксперимента индикаторным графикам площадного — S =f, (lg г), временного — S =/ (lg t) и комбинированног о прослеживания.
Важнейшим элементом при диагностике откачки являются индикаторные графики временного прослеживания. При достаточной длительности опы- а та и полноте измерений
сама форма этих гравис. 5.18. Характерные индикатор- фиков часто может сви- ные графики откачки: детельствовать О тех
1 - в изолированном напорном пласте; 2 - в ” ^ пгпбеттостях
гетерогенном пласте; 3 - вблизи водонепрони- и т " ‘ uluuchhuwijia
цаемого контура; 4 - вблизи контура питания фИЛЬТраЦИОННОГО ПрО- или при наличии перетекания из весьма водо- цесСЗ. Для Примера На обильного пласта рИС 5 j g ПрИведены ха
рактерные индикаторные графики некоторых типовых расчетных схем. С помощью подобных графиков идентифицируется наличие непроницаемых или слабопроницаемых контуров (кривая 3), дополнительное питание из поверхностных водоемов или мощных водоносных пластов (кривая 4) и т.п. Вместе с тем, многие факторы проявляются на индикаторных графиках качественно идентично (см. разделы 5.3 и 5.4), что резко снижает возможности дагностики.
ЗАДАНИЕ. Пользуясь материалом разделов 5.3 и 5.4, приведите возможные объяснения характера графиков 2, 3 и 4 на рис. 5.18, применяя их к различным типовым схемам. Не забудьте при этом, что графики вида 3 и 4 могут являться «усеченными» вариантами графика вида 2.
С учетом упомянутых обстоятельств, диагностика индикаторных графиков одиночных откачек оказывается весьма субъективной. С одной стороны, эти графики в максимальной степени подвержены влиянию различных аномальных факторов, а с другой — практически отсутствует контрольная информация, позволяющая выявить и устранить их влияние.
С этой точки зрения кустовые откачки обладают неоспоримыми преимуществами: возможность взаимного дополнения и сопоставления информации по различным наблюдательным скважинам, а также устранение погрешностей, связанных с измерениями уровней по центральной скважине, позволяют рассчитывать на сравнительно надежную диагностику гораздо чаще, чем в одиночных экспериментах. Это, в частности, касается выбора представительных участков индикаторных графиков для последующего расчета: в разделах 5.3 и 5.4 показано, что разные участки графиков чаще всего отвечают различным аналитическим зависимостям. Например, при обработке способом прямой линии (см. раздел 5.5.1) необходимо иметь гарантию, что выбранный прямолинейный участок не искажен влиянием разного рода технических факторов (неравномерность работы насоса, скин-эффект централь-1 ной скважины, несовершенство по степени вскрытия и инерционность пьезометров — см. раздел 5.4). В этом плане важно подчеркнуть следующие элементы, связанные; с диагностикой:
1J возможность проверки интерпретационной схемы по параллельности конечных участков графиков временного прослеживания или по разбросу точек на графике комбинированного прослеживания (в координатах S—■
2 возможность проверки — для планово-однородных пластов — представительности выбранного участка графика по близости расчетных значений проводимости, полученных временным и площадным прослеживанием (см. раздел 5.5.1). Подобные возможности полностью исключаются для одиночных откачек, что делает выделение расчетного участка графика крайне субъективным.
ПРИМЕР [23 ]. Кустовая откачка продолжительностью 7 сут проводилась с устойчивым расходом 1200 м /сут в безнапорном водоносном горизонте, приуроченном к среднезернистым пескам мощностью около 40 м. Проводимость2определенная по наблюдательным скважинам, составляла — 750 м /сут. На графике восстановления (рис. 5.19) для центральной скважины диаметром 300 мм, пройденной вращательным способом на глинистом растворе, выделяются три (1-3) прямолинейных участка (измерения, сделанные в течение первой минуты и отражающие доминирующее влияние емкости скважины, на графике опущены). Расчетные значения проводимости, полученные по уклонам разны£ участков, составляют, м /сут: Т' = 90, Т" = 800 и Т'" = 3700. Величина Т', рассчитанная по участку 1 (для t < 10мин), отражает влияние скин-эффекта (с емкостью скважины) и отвечает проводимости прифильтровой зоны скважины. Величина Т'' примерно соответствует истинной проводимости пласта (участок 2, время t от 10 мин до 2 ч). Наконец, значение Т,п явно завышено, что обусловлено выполаживанием графика на участке 3 вследствие эффектов упругогравитационного режима (см. раздел 5.3.4), а также влияния истории откачки.