Файл: Шарапанов, Н. Н. Методика геофизических исследований при гидрогеологических съемках с целью мелиорации земель.pdf

ность корреляции результатов количественной интерпретации соседних кривых рк.

После установления пределов изменения величины Kh, целе­ сообразно, исходя из конкретных геологических условий с учетом абсолютной величины пределов изменения Kh, выбрать либо одно (средневзвешенное) значение Kh для всего интервала глубин, либо принять несколько значений Kh для различных полуразносов ABI2. В обоих случаях необходимо оценить возможную ошибку определения глубины, получаемую за счет осреднения величины Kh.

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

Количественная пнтерпретацпя кривых р,( является основой для определения глубины залегания геоэлектрпческпх горизонтов. Кривые ВЭЗ ВП (кривые комплексного параметра А к и его отно­ сительной величины /1,*) интерпретируются лишь совместно с кри­ выми рк способами, наложенными в работах В. А. Комарова [19], [39], и, в частности, методом, предложенным А. М. Ярослав­ цевым.

Рассмотрим некоторые особенности количественной интерпре­ тации данных сейсморазведки, полученных при определении положения уровня грунтовых вод в рыхлых породах [9, 22, 23].

В случае плоской поверхности наблюдений годографы про­ дольной, преломленной на уровне грунтовых вод (УГВ) волны iy, отличаются гладкостью и прямолинейностью; кажущиеся ско­ рости, определенные по встречным годографам, близки, нагоня­ ющие годографы практически параллельны нагоняемым. Ско­

н затухает тем быстрее, чем меньше глубина Н у и мощность водо­ носного горизонта. В целом, однако, затухание волны ty меньше по сравнению с преломленными волнами, связанными с другими границами в зоне аэрации.

100

где жт "п — точка пересечения годографа прямой (рефрагированной) волны и преломленной волны ty, связанной с УГВ.

Формулу (38) используют в том случае, когда при наблюдениях бывают получены встречные годографы (простые или сводные), имеющие участок перекрытия. Она позволяет определять глу­ бину Ну во всех точках, где известна величина toy. Формула (39) дает возможность определять глубину Н у только в точках воз­ буждения колебаний и поэтому применяется главным образом при интерпретации встречных годографов, не имеющих перекры­ того участка, или в случае одиночного годографа. Условием при­ менения формулы (39) является также отсутствие ясно выражен­ ных границ в зоне аэрации.

Величину vcp находят обычно по точке пересечения годографов волны ty и волны, предшествующей ей. При наличии заметной рефракции значение vcp рекомендуется находить при помощи приемов интерпретации, разработанных для градиентных сред [9].

Скорости иГ' у определяют обычно по разностному годографу. В случае, если встречные годографы не имеют перекрытого участка или используется одиночный годограф, значение vr принимают равным кажущейся скорости волны ty.

При наличии перекрытого участка у встречных годографов линию t0y находят во всех интервалах перекрытия обычным спо­ собом [9, 11]. При использовании одиночных годографов вели­ чину t0y определяют лишь в точках возбуждения колебаний, для чего годографы волны ty продолжают до пересечения с осью времен.

Вслучае, когда в зоне аэрации имеется ясно выраженная граница, глубину Ну рекомендуется вычислять по формуле пла­ стовых скоростей [9, И].

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

соотношение:

Данное соотношение позволяет нередко получить достаточно хорошие результаты. Рекомендуется использовать его при наблю­ дениях в поле для выбора оптимальной методики исследований. Кроме того, именно соотношение (40) применяется для оценки положения УГВ в случае его глубокого залегания и невозмож­ ности четкого прослеживания волны ty.

После выполнения количественной интерпретации материалов

•сейсморазведки, ВЭЗ и ВЭЗ ВП оказывается возможным перейти непосредственно к построению соответствующих геолого-гидро­ геологических карт и разрезов.

Исследования з различных геолого-гидрогеологических усло­ виях позволили установить, что геолого-гидрогеологическая ин­ терпретация результатов, несмотря на сохранение общего прин­

101

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

К а р т а г л у б и н ы з а л е г а н и я у р о в н я п е р ­ в о г о от п о в е р х н о с т и в о д о н о с н о г о г о р и ­ з о н т а . Эту карту строят по результатам определения глубины

Рпс. 26. Зависимость разности Д7/глубин до уровня грунтовых вод от степени глинистости пород Г (а)

и от параметра А* (6)

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

102

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

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

На рис. 26 представлена зависимость разности АН глубин до УГВ по данным бурения IIб и сейсморазведки Нср от содержа­ ния пылеватых и глинистых частиц Г для Сарпинской низменности в Калмыкии: АН — IIб — Н ср. Сопоставляя эту зависимость с за­ висимостью для относительной величины параметра Л* = / (Г), можно составить зависимость Лугв = / (Л*), которую необхо­ димо использовать для введения поправки Д//угв в получаемую по данным сейсморазведки глубину до УГВ. Изложенное позво­ ляет представить следующую общую схему построения карты глу­ бины залегания первого от поверхности водоносного горизонта.

1.Карту строят на основе уже построенных литологических разрезов.

2.В зависимости от гидродинамических особенностей водо­

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

ДЯугв - / (Л*).

3. На основе данных, полученных в тех же параметрических скважинах, оценивают ошибку определения глубины до уровня по результатам сейсморазведки и выбирают сечение изолиний на карте, которое должно минимум в 2 раза превышать ошибку опре­ деления УГВ.

Л и т о л о г и ч е с к а я к а р т а и р а з р е з ы . Расчле­ нение разреза по литологическому составу и дальнейшее по­ строение карт выполняют по результатам интерпретации электроразведочных данных. Геолого-гидрогеологическую интерпретацию начинают с разработки наиболее рациональной схемы интерпре­ тации, которая определяется гидрогеологическими условиями района работ, а следовательно, геоэлектрическим разрезом. Пре­ жде всего целесообразно провести качественный анализ параме­ трических кривых ВЭЗ ВП, выполненных вблизи скважин.

На рис. 27 приведены кривые ВЭЗ ВП, выполненных вблизи скважины, разрез которой является достаточно полным и харак­ теризует геологическое строение площади съемки (Сарпинская низменность).

103

Рис. 27. Параметрические кривые ВЭЗ ВП, выполненного у скв. 7к (Калмыкия).

} суглинки; 2 — супеси; 3 — пески; 4 — глины. Пунктирными линиями показаны ориентировочные асимптотические значения параметроц

Кривая рк в целом отражает все литологические горизонты в разрезе и может быть использована для количественной интер­ претации, в результате которой должны быть определены глубина залегания и мощность выделяемых горизонтов. Геологическая интерпретация по р затруднена пз-за слабой дифференциации кривой рк. Характерно некоторое увеличение рк для глин по сравнению с песками, насыщенными водами с минерализацией

20—40 г/л.

Кривые поляризуемости рк и постоянной времени спада тк более дифференцированы и могут быть использованы для геолого­ гидрогеологической интерпретации. Еще более дифференцирована кривая комплексного параметра А к\ это вызвано тем, что увеличе­ ние поляризуемости и одновременное увеличение скорости спада для суглинков и глин влечет за собой более резкий рост величины комплексного параметра. Обратная картина наблюдается для песков.

Кривая относительной величины параметра Л,* по степени дифференциации близка к кривой А к, однако большая ее неста­ бильность в области малых сопротивлений позволяет отдать предпочтение использованию параметра А для литологического расчленения разреза.

После того как намечена общая схема интерпретации, необхо­ дим детальный анализ составленных корреляционных зависи­ мостей и выбор литологических показателей для построения карт и разрезов. В'качестве корреляционных зависимостей следует выбрать те, которые характеризуются достаточно высокими коэф­ фициентами корреляции. При этом необходимо проанализировать зависимости между геофизическими параметрами и всеми показа­ телями гранулометрического состава песчано-глинистых пород (содержание глинистых, пылеватых, песчаных частиц■, число пластичности). В том случае, если получены зависимости от ряда показателей гранулометрического состава (и особенно от числа пластичности), можно, пользуясь той пли иной классификацией грунтов, составить таблицу геолого-гидрогеологической интер­ претации.

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

глины (частиц <( 0,01 мм — более 70%); тяжелые суглинки (частиц <0,01 мы — 50—70%);

средние суглинки (частиц <0,01 мм — 40—50%); легкие суглинки (частиц < 0,01 мм — 20—40%); супесь (частиц < 0,01 мм — 10—20%); песок (частиц < 0,01 мм — менее 10 %).

105