弥河下游浅埋古河道的勘探识别

图1 研究区地理位置和电测深测点分布

Fig.1 Illustration map showing location of the study area and VES measured points

20世纪70年代由于气候连续干旱,区内地下水超采严重,先后出现咸水入侵和海水入侵^[11]。浅埋古河道带是机井分布最密集的区域,随着地下水资源与用水需求间矛盾的日益尖锐,含水层被疏干,形成区域性地下水降落漏斗,诱导北部浅层咸水横向入侵,咸水南侵现象严重^[11-12]。疏干的古河道是咸水入侵的主要通道,也是防治咸水入侵的天然工程^[12-14]。为了缓解当地的用水危机,防止海(咸)水入侵持续加剧,该项勘探工作刻不容缓。

在第四纪松散层中探测古河道是基于古河道砂层与黏土层有明显的电阻率对比度特性^[15]。电阻率对电阻结构比较敏感,可以很好地定义浅层电阻区域^[16]。研究区地处莱州湾南岸全淡水区,通过分析物探、钻探和水文地质资料,统计出该区第四纪岩层的电性参数(图2)^[17]。

图2

图2 研究区第四纪岩层电性参数统计

Fig.2 Electrical parameters cartogram of Quaternary rocks in the study area

2 勘探设计与方法

对称四极电测深法,又称垂向电测深法(vertical electrical sounding,简称VES),是以岩层的电性差异为基础的一种体积勘探方法。电测深采集的数据根据下式计算视电阻率ρ_s:

$ρ_{s} = K \frac{Δ U_{MN}}{I}$ ,

$K = \frac{2 π}{\frac{1}{AM} - \frac{1}{BM} - \frac{1}{AN} + \frac{1}{NB}}$ 。

式中:ΔU_MN为测量电极M、N处电位差,I为供电电极A、B间电流,K为装置系数,AM、BM、AN、NB分别为各电极间的距离。

根据勘探前实地考察信息,结合历史记载资料,本次勘探剖面呈EW走向,在空间上垂直于弥河冲洪积扇主轴方向。电测深剖面全长23.2 km,横跨化龙、文家街道和古城街道三个城镇,共布设75个测深点(见图1),以初步确定古河道砂层分布区。使用重庆地质仪器厂生产的DZD-6A多功能直流电法仪进行数据采集,应用对称四极装置,最大供电极距AB=180 m。极距布置方式见表1。

表1 电测深数据采集极距

Table 1 Data acquisition list of VESm

序号	1	2	3	4	5	6	7	8	9	10	11	12	13
AB/2	4	6	9	2	16	20	25	32	40	50	60	74	90
MN/2	1.5	1.5	1.5	1.5	1.5	1.5	4	4	4	4	6	6	6

新窗口打开| 下载CSV

3 结果分析与讨论

3.1 电测深结果分析

在第四纪古河道富集区,由细砂、中粗砂和砂砾石构成的古河道沉积物的电阻率高于上覆及基底的黏土层。根据研究区的地质条件,筛选电测深法的勘探数据,初步判定20#~28#测深点处属异常区,预期在一定深度上存在古河道。

图3给出异常区9个测深点的视电阻率曲线,可以看出区内电测深曲线大致分为三个电性层,基本属ρ₁<ρ₂>ρ₃的“K”型曲线。

图3

图3 电测深法视电阻率曲线及其一维反演电阻率

Fig.3 The VES apparent resistivity curves and One-dimensional inversion resistivity

为了宏观体现异常区岩层的电阻率分布特征,根据9个测点的分布位置,绘制20#~28#测深点视电阻率等值线拟断面图(图4)。整个异常区域长约2 km,测点平均间距250 m;视电阻率等值线在水平方向的连续性体现了该区第四纪地层厚度大且连续,探测范围内未见基岩高阻异常;而视电阻率值在垂向上分布的不均匀性反映地层岩性在空间上的不均匀性。从数值上看,20#和28#测点处视电阻率值小于28 Ω·m,是整个异常区的低阻区;21#和27#整个探测范围及22#~26#测点的浅地表的视电阻率值在28~37 Ω·m;22#~26#测点属于整个区域的高阻分布区,尤其23#、24#和25#测点在AB/2=30~50 m范围内视电阻率值高达70 Ω·m以上。整个异常区高阻是低阻的3.5倍,明显的电性差异反映探测范围内存在不同电性的岩层,体现该区存在古河道砂层富集的可靠性。

图4

图4 勘探异常区的视电阻率拟断面

Fig.4 The apparent resistivity pseudosection in exploration anomaly area

3.2 基于反演模型分析古河道特征

为了获得地电体更精确的电性特征,应用最小二乘法进行一维模拟反演,结果见图3。

从图3中可以看出,电阻率反演模型直观展现了异常区在勘探深度范围内的3个不同电性的电性层排列。该区地层都呈现明显的中间高阻、上下层低阻的“K”型断面地电特征,反映地下岩层在埋深范围内的相似性,而不同测点的电阻率和埋深位置的不同,则反映不同测点岩层的岩性、粒径和空间分布的差异。为了定性描述岩层的电性特征,根据反演后的模型提取各测点的电性分层数据,详见表2。结合研究区地质、地球物理概况分析如下。

表2 异常区一维反演的岩层电阻率—深度统计

Table 2 The table of One-dimensional inversion resistivity-depth in anomaly area

测深点编号	表层		中层		底层
测深点编号	ρ/(Ω·m)	埋深/m	ρ/(Ω·m)	厚度/m	ρ/(Ω·m)	埋深/m
20#	19	6	39	24	19	30
21#	17	6	56	21	21	27
22#	15	7	73	27	34	34
23#	18	7	102	27	47	34
24#	17	6	104	24	44	30
25#	19	6	102	21	39	27
26#	19	5	82	19	21	24
27#	15	6	60	21	23	27
28#	20	6	40	24	20	30
均值	18	6	73	23	30	29

新窗口打开| 下载CSV

异常区近地表层的电阻率变化范围较小,最大值位于28#测点处,最小值位于22#和27#测点处,9个测点的平均表层电阻率为18 Ω·m,该层整体属黏土层,是古河道的顶板结构。

中部岩层的电阻率值在39~104 Ω·m,其中20#、28#的电阻率值小于40 Ω·m,属黏质砂土分布,21#、22#和27#测点的电阻率值在50~80 Ω·m范围内,是细砂层富集,23#~26#测深处的电阻率值大于80 Ω·m,属中粗砂分布区,尤其23#、24#和 25#高达100 Ω·m以上,反映中粗砂富集,含少量砂砾石。该层电阻率的最大值(24#测点处)是最小值(20#测点处)的2.7倍,反映砂的粒径及富集程度的区别,粒径越大,富集程度越好,阻值越大,富水性越高。该层9个测点的电阻率均值为73 Ω·m,是上覆黏土层阻值的4倍,明显反映古河道富集区的高阻分布特征。

底层的电阻率在19~47 Ω·m,最大、最小值分别位于23#和20#测点处,其中,20#、21#、26#、27#和28#测点处的电阻率值小于25 Ω·m,属黏土层沉积区,22#和25#测点处的电阻率值在30~40 Ω·m之间,属粉质砂土富集区,23#和24#处的电阻率在40~50 Ω·m范围内,岩性上属黏质砂土或细砂或黏质砂土和细砂互层区。整个底层的电阻率均值为29 Ω·m,属黏土和粉质砂土互层,夹少量细砂,构成古河道砂层的隔水底板结构。

从岩层的埋深数据上分析,异常区整个勘探线推断的古河道砂层的顶部平均埋深为6 m,最大埋深为7 m,位于22#和23#测深点,最小埋深5 m,位于26#测深点;底部平均埋深为29 m,最大埋深为34 m,位于22#和23#位置,最小埋深24 m,位于26#测深点。根据砂层的顶底板埋深,推断该区古河道砂层的平均厚度为23 m,其中最大厚度27 m,最小厚度19 m。

根据一维电阻率反演模型的解释结果绘制异常区古河道剖面图,初步重建该区的地球物理模型。如图5所示,20#和28#测深范围内主要是黏土分布,只在中部是粉质黏土和细砂互层;21#~27#除顶底板由黏土、粉质黏土或两者互层构成外,中层体现为细砂、中粗砂、砂砾石分布,是该区古河道砂层的主要富集区。该区地层从上到下整体呈现“黏土—砂层—黏土”的岩性规律。图中连接整条勘探线的细砂、中粗砂和砂砾石分布区,初步得到异常区的古河道剖面,可见整体呈河流相沉积,该处砂层富集明显,连续性较好,是良好的含水层结构,顶板为沉积黏土层,底板是黏土、粉质砂土或细砂沉积,反映稳定的沉积环境。

图5

图5 古河道电测深解释成果及钻孔岩芯柱状图

Fig.5 The VES interpretation profile of paleochannel and column diagram of drilling cores

3.3 与水文地质钻探结果对比分析

为了进一步揭露古河道特征,验证电法勘探结果,分别在22#和26#测深点附近布置K1、K2两个水文地质钻孔,根据取芯样本绘制钻孔柱状图,详见图5。

从图中可以看出,钻孔揭露的异常区地层岩性与物探解释结果大致相同。根据物探、钻探揭露的古河道砂层的埋深范围作误差分析,其中,测深点22#附近K1钻孔揭露的古河道顶板埋深7.0 m、底板埋深33.5 m,该点电测深反演的古河道顶板埋深7 m、底板埋深34 m,两者的相对误差分别为0%、1.5%;测深点26#附近K2钻孔揭露的古河道顶板埋深4.9 m、底板埋深23.8 m,该点电测深反演的古河道顶板埋深5 m、底板埋深24 m,两者的相对误差分别为2.0%、0.8%。可以看出,电法反演与钻孔揭露两种结果的相对误差都在2%以内,体现两者的勘探结果吻合度高,证明应用电测深法寻找古河道含水层的可行性和可靠性。

4 结论

1) 勘探区域内富集古河道地层的电性特征,在以地质条件和第四纪松散层的电性参数为基础的VES勘探方法所获得的电测深曲线中,总体上呈现出“K”型曲线特征,其曲线的上升段反映了古河道中砂砾层信息。根据反演后的电阻率曲线的特征,可以较为明确地推断出地层中的古河道砂砾富集信息。

2) 本次电法勘探结果表明,研究区古河道砂层反映为56~104 Ω·m的高阻区,埋深在5~7 m不等,平均厚度23 m,整体以河流相沉积在近地表黏土层下部,岩性为细砂、中粗砂和砂砾石,是区内主要含水层。近地表和下覆岩层属低阻区,其中近地表黏土层平均电阻率18 Ω·m,下覆岩层平均电阻率30 Ω·m,主要为稳定黏土层,局部粉质黏土富集,分别属古河道的顶、底板结构。

3) 根据水文地质钻探结果,电测深法所反映出的古河道特征与钻孔所揭示的古河道特征基本一致。在以黏土和砂砾为主要岩性的第四系松散地层中,电测深法可以有效地定位古河道信息。但由于物探结果往往具有多解性,对于地层成因复杂的区域中古河道勘探可采用两种或多种物探方法组合,相互校正以提高勘探结果的准确性。

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

吴忱, 许清海, 赵明轩.

世界所有大河都有埋藏古河道

[J]. 地理学与国土研究, 1992,8(2):29-34.

, Xu

Q H

, Zhao

M X

All rivers in the world have buried palaeochannel

[J]. Geography and Geo-Information Science, 1992,8(2):29-34.

[2]

吴忱, 朱宣清, 何乃华, 等. 华北平原古河道研究[M]. 北京: 中国科学技术出版社, 1991:41-132.

[本文引用: 4]

, Zhu

X Q

, He

N H

, et al. The research of palaeochannel in the North China Plain [M]. Beijing: China Science and Technology Press , 1991:41-132.

[本文引用: 4]

[3]

赵艳霞, 徐全洪, 刘方圆, 等.

近20年来中国古河道研究进展

[J]. 地理科学进展, 2013,32(1):3-19.

DOI:10.11820/dlkxjz.2013.01.001 URL [本文引用: 1]

近20年来,我国在嫩江大安古河道、莱州湾南岸平原古河道、长江南京以下段古河道及古深槽、海底古河道等区域的研究,获得了大量研究成果,深化了理论认识,拓展了应用领域。这些新成果与新认识进一步表明,中国各外流大河均有末次盛冰期埋藏古河道,该期古河流深切于末次冰期间冰阶的陆地及浅海大陆架,形成切割谷和谷底部的深槽,谷内堆积了末次盛冰期和早全新世的河流相滞留物质和粗颗粒物质,上部被中全新世海相、海陆交互相或陆相细粒物质覆盖,构成了埋深约20~50 m的浅埋古河道带。古河道中蕴藏着比较丰富的地下淡水,是洪冲积平原及浅海陆架地区的重要水源。

Zhao

Y X

, Xu

Q H

, Liu

F Y

, et al.

Progresses of palaeochannel studies in China in the past 20 years

[J]. Progress in Geography, 2013,32(1):3-19.

DOI:10.11820/dlkxjz.2013.01.001 URL [本文引用: 1]

The last 20 years have seen new progresses in palaeochannel studies in China, including the investigations on palaeochannels in Nenjiang Daan, palaeochannels in the plain south of Laizhou bay, and palaeochannels, deep ancient grooves and submarine palaeochannels in the Yangtze river downstream of Nanjing, and so on. As a result, a large number of achievements have been accomplished, and new theories proposed; the research achievements have found new applications as well. The new knowledge and new theories indicate that buried palaeochannels exist in all of large outflowing rivers in China since the LGM(Last Glacial Maximum). The palaeochannel runs deep into the stadial terrestrial and the sea bed of the continental shelf from the Last Glacial Epoch, forming a cutting valley with deep grooves at the bottom, stacked with fluvial carryover and coarse particulate material from the LGM and early Holocene. The upper part of the valley is covered with marine, paralic, and continental fine-grained material from the Holocene, which frames the shallow-buried palaeochannel with 20 m to 50 m depth. Palaeochannels are rich in underground water, and an important water source for the areas of alluvial plain and shallow continental shelf.

[4]

Henk

J A

, Berendsen

, Esther

Late Weichselian and Holocene palaeogeography of the Rhine-Meuse delta, the Netherlands

[J]. Palaeogeography, Palaeoclimatology, Palaeoecology, 2000,(161):311-335.

[5]

金永念, 季克其, 仇慎平, 等.

江苏阴平、华冲地区第四纪古河道电性特征研究

[J]. 水文地质工程地质, 2004(3):99-101.

URL [本文引用: 1]

从电性、电测深曲线类型和地电结构研究入手,利用电测深技术对江苏东北部阴平、华冲地区第四纪古河道进行了推断解译,指出区内第四纪古河道发育在朱墩-桑墟以东地区,埋深25～50m,尤其是50m左右古河道内砂层颗粒粗、层厚且富水.

Jin

Y N

, Ji

K Q

, Qiu

S P

, et al.

Electrical property of Quaternary ancient watercourse of Yinping and Huachong area, Jiangsu Province

[J]. Hydrogeology & Engineesing Geology, 2004(3):99-101.

[6]

付新建, 刘克宇, 李改梅.

对称四极电阻率测深曲线在古河道带及河间带的电性特征

[J]. 地下水, 2008,30(5):107.

X J

, Liu

K N

, Li

G M

The electrical characteristics of symmetrical four resistivity sounding curves in the ancient river belt and interchannel areas

[J]. Groundwater, 2008,30(5):107.

[7]

郭高轩, 刘文臣, 辛宝东, 等.

利用电测深法探测泃河与错河古河道

[J]. 工程勘察, 2010,(4):87-90.

Guo

G X

, Liu

W C

, Xin

B D

, et al.

Old course detection of Juhe and Cuohe river by electric sounding method

[J]. Geotechnical Investigation & Surveying, 2010,(4):87-90.

[8]

郭龙凤, 黄少文, 李亮亮, 等.

综合物探方法在古河道型地下水库工程建设中的应用研究

[J]. 地球物理学进展, 2018,33(3):1205-1212.

Guo

L F

, Huang

S W

, Li

L L

, et al.

Application and research of integrated geophysical method in construction of palaeochannel groundwater reservoir

[J]. Progress in Geophysics, 2018,33(3):1205-1212.

[9]

孟庆海, 韩美, 赵明华, 等.

弥河冲洪积扇和古河道初步研究

[J]. 山东师大学报:自然科学版, 1999,14(1):46-50.

Meng

Q H

, Han

, Zhao

M H

, et al.

A preliminary Study of the Mihe River alluvial diluvial fan and the palaeochannels

[J]. Journal of Shandong Normal University:Natural Science, 1999,14(1):46-50.

[10]

韩美

莱州湾地区海水入侵与地貌的关系

[J]. 海洋与湖泊, 1996,27(4):414-420.

Han

Relationship between the seawater Intrusion and landforms in Laizhou Bay area

[J]. Oceanologia et Limnologia Sinica, 1996,27(4):414-420.

URL [本文引用: 2]

To study the relationship between seawater intrusion and landforms, field investigations, laboratory analyses, satellitic photos, drillings, and relief maps, were used in this May, 1993 to October, 1995 study on the relation between seawater intrusion and landforms. The results showed that sandy bays and argillaceous bays favor seawater intrusions. Because geomorphic types in Laizhou Bay's south coast are different from those in the east coast, seawater intrusions in the south coast were different from those in the east coast. Paleochannels and coastal landforms were the main geomorphic types that affected seawater intrusions. The modes and intensity of seawater intrusions were not alike in distinct landform units and types.

[11]

李道高, 赵明华, 韩美, 等.

莱州湾南岸平原浅埋古河道带研究

[J]. 海洋地质与第四纪地质, 2000,20(1):23-29.

[本文引用: 3]

D G

, Zhao

M H

, Han

, et al.

A study of the shallowly-buried palaeochannel zones in the south coast plain of the Laizhou Bay

[J]. Marine Geology & Quaternary Geology, 2000,20(1):23-29.

[本文引用: 3]

[12]

韩美, 赵明华, 李道高, 等.

莱州湾南岸平原古河道及其与海(咸)水入侵关系研究

[J]. 自然灾害学报, 1999,8(2):73-80.

Han

, Zhao

M H

, Li

D G

, et al.

Study on the ancient channels and the relationship between the ancient channels and the sea (salt) water intrusion of the south coastal plain of Laizhou Bay

[J]. Journal of Natural Disasters, 1999,8(2):73-80.

[13]

于得芹, 蒙永辉.

寿光市水文地质特征浅析

[J]. 环境地质, 2014,30(8):30-33.

D Q

, Meng

Y H

Analysis on hydrogeological characteristics of Shouguang city

[J]. Environmental Geology, 2014,30(8):30-33.

[14]

刘恩峰

莱州湾南岸滨海平原沉积环境变化与咸水入侵关系研究

[D]. 济南:山东师范大学, 2002.

Liu

E F

The relation research of the sedimentary environmental transformation and saline water encroachment in Laizhou Bay south sea shore plain

[D]. Jinan: Shandong Normal University, 2002.

[15]

Cardarelli

, Donno

G D

Multidimensional electrical resistivity survey for bedrock detection at the Rieti Plain (central Italy)

[J]. Journal of Applied Geophysics, 2017,141:77-87.

[16]

Martínez Moreno

F J

, Monteiro Santos

F A

, Bernardo

, et al.

Identifying seawater intrusion in coastal areas by means of 1D and Quasi-2D joint inversion of TDEM and VES data

[J]. Journal of Hydrology, 2017,552:609-619.

[17]

何玉海

高密度电法在莱州湾海水入侵调查中的研究与应用

[J]. 海洋环境科学, 2016,35(2):301-305.

Y H

Research and application of High-density resistance method in seawater invasion investigation of Laizhou Bay

[J]. Marine Environmental Science, 2016,35(2):301-305.