Characteristics of shallow underground ancient channels in urban areas: A case study of Quaternary underground palaeochannel in Tonghu Basin, Guangdong Province
-
摘要:
人类文明与城市发展往往与河流息息相关。在地质历史时期继承性形成的古河道被埋藏在城市地下。由于其自身垂向和横向上的演化特性，造成城市地下不良岩土体分布较复杂，进一步造成古河道分布区地下空间开发过程中工程地质问题突出，因此，需要摸清古河道的空间特征。以地处大湾区东翼的广东省惠州市潼湖盆地第四系为例，综合利用钻探及地球物理资料，结合地表继承性沉积，从精细地层划分与沉积相识别入手，通过"垂向分期、侧向划界、平剖互动"的流程，逐级对古河道及沉积内部构型进行深入的解剖分析，形成了一套城市地下古河道分析描述方法，，对大湾区其他区块和具有相似沉积特征的其他城市地下空间表征提供借鉴。
Abstract:Human civilization and urban development are closely related to rivers. Due to its longitudinal and horizontal evolution, the complex distribution of palaeochannel formed in the geological history period and buried underground, the distribution in urban underground is relatively complex. However, in view of the prominent engineering geological problems such as the sensitivity and vulnerability of the distribution area of palaeochannel in the process of space development, it is urgent to find out the spatial characteristics and engineering characteristics of palaeochannel course. The quaternary system of Tonghu basin in Huizhou City, located in the east wing large bay area of Guangdong Province, was taken as an example. Comprehensive use of drilling and geophysical data, combined with surface inherited deposits, this work started from fine stratigraphic division and sedimentary facies identification, through the process of "vertical staging, lateral demarcation, flat profile interaction". A set of methods for analyzing and describing ancient rivers were establish by analyzing the structure of ancient river and sedimentary interior in depth step by step, which can provide reference for the characterization of underground space in other areas of the Greater Bay Area and other cities with similar sedimentary characteristics.
-
Keywords:
- urban geology /
- palaeochanne /
- deposition configuration /
- the Greater Bay Area
-
地下水氮污染是世界性难题[1-2]，而中国90%的城市地下水不同程度遭受N等元素污染，地下水污染已呈现由点向面、由城市向农村扩展的趋势[3]。近年来，随着煤炭资源开发促进能源化工基地经济的发展，煤炭开采过程中工业污染物的排放极有可能对地下水水质安全造成严重的威胁[4-12]。90%的煤炭含氮量为0.52%~1.41%[13]。前人探讨了煤矸石回填对地下水氮的影响[14]、采煤沉陷积水区地表水与浅层地下水氮的影响[6]，对能源化工基地煤炭资源开发区地下水三氮污染影响因素尚未开展研究。
煤炭资源开发过程中，改变原有包气带中所含的潜在污染物改变了原有的循环模式，而与地表水联系最紧密的是浅层地下水，开采过程在地表、浅层地下水中形成了特有的转化模式[15]，尤其是在采动过程中可能形成的一系列波状地貌，如大量发育的裂隙、裂缝、阶地等。研究区为沙地区，会形成大量的波状起伏的丘陵地貌，对地表、包气带的结构及水、污染物的运移产生了潜在的影响，区域物源补给直接影响氮的含量。故本文立足于地下水径流系统和物源条件探讨转化机理。
宁东煤炭基地鸳鸯湖矿区是国家煤炭规划建设的重要矿区，梅花井井田为鸳鸯湖矿区的5个大型煤矿之一。宁东煤炭基地属2个含煤时代，石炭纪—二叠纪和侏罗纪煤田8个矿区共含煤38层，梅花井矿区主要产侏罗纪煤炭。煤炭开采过程中工业污染物可能含有大量的氮污染物，导致含水层氮污染，生态环境脆弱，水资源相对匮乏，煤炭开采区地下水资源的保护及矿区可持续发展具有重要意义。前人研究了梅花井矿区水文地质条件[16]、采矿对含水层水量及结构的破坏类型和程度[17-19]，而对煤矿区乃至周边的影响区水质尚未探讨。梅花井井田位于宁东夹山梁地貌单元，本文基于地下水补径排条件选择宁东夹山梁地貌为基本研究单元（包括了梅花井井田和周边环境），研究风积沙滩地区煤矿井开采地下水氮污染的程度和影响因素。
1. 研究区水文地质条件
夹山梁位于宁夏回族自治区灵武市以东33km处。梅花井井田位于夹山梁的中东部，鸳鸯湖背斜东翼中部，矿区面积78.96km2，开采深度1300~ 200m。属半沙漠低山丘陵地形。地表为沙丘掩盖，多系风成新月形和垄状流动沙丘。区内无常年地表径流，雨季降水及矿井排水在井田北部低洼地区形成盐碱湖，冬季干枯为碱滩，基岩被第四系风积沙覆盖。
夹山梁地区含水层按岩性组合特征及地下水水力性质、埋藏条件等，主要含水层由上而下划分为第四系孔隙潜水含水层和侏罗系直罗组、延安组砂岩裂隙-孔隙承压水含水层（图 1）。
![]()
图 1 水文地质图及样点分布Q—第四系；K—白垩系；J—侏罗系；T—三叠系Figure 1. The hydrological geological map and sample distribution（1）本区第四系厚2.75~17m，平均厚5.15m①，地下水主要赋存于风积沙、小型洼地及沟谷冲洪积层中。按地下水赋存条件，可分为风积沙潜水层、风积-冲洪积潜水层。地下水主要赋存于风积沙、小型洼地及沟谷冲洪积层中。含沙漠凝结水，地下水位埋深1.00~3.48m，多随地形起伏而异，水位、水量随季节变化，矿化度为0.8~8.03g/L，水化学类型主要为Na- HCO3、SO4 · HCO3- Na · Mg、Cl·SO4-Na·Mg。
（2）侏罗系含水层组是影响矿区煤矿开采的主要含水层，包括上侏罗统直罗组含水层、中统延安组含水层，厚度为13.01~618.62m。煤层厚度主要分布在6.85~337.45m之间，垂向上上部直罗组含水层富水性较强，对梅花井煤层开采影响较大，下部延安组含水层结构较致密，裂隙不发育，富水性较差，对煤层开采影响较小，水化学类型为SO4· Cl-Na、HCO3-Mg·Na·Ca或SO4·HCO3-Na·Mg。
（3）隔水层岩性多为中细砂岩与粉砂岩、泥岩互层，岩性致密，与煤层共同形成良好的隔水层。
2. 研究方法
采样点位于鸳鸯湖矿区夹山梁，采集第四系含水层水样10个和侏罗系含水层样品3个（图 1），利用GPS精确定位，潜水水位埋深一般为6~8m，侏罗系含水层水位埋深一般为65~100m。采样时间为2016年9月12日—9月17日，统测1次。水样采集、固定及保存均按照《水质采样样品的保存和管理技术规定》（HJ 493—2009）②进行，并及时运往实验室检测。依据水和废水监测分析方法[20]对浅层地下水中的氮分布及限制性特征进行分析；现场测试指标为水温，室内测试分析TDS、NO3-、NO2–、NH4+共4项。水化学测试由中国地质调查局西安地质调查中心实验测试中心完成。评价方法参照《矿山地质环境调查评价规范DD2014— 05》 8.3.2.2条款，地下水污染评价采用单项超标倍数法。
3. 结果
3.1 三氮空间污染特征
对夹山梁地区浅层地下水中主要离子及部分水化学指标的数据统计显示（表 1），本区浅层地下水pH值范围为7.13~9.63，浅层地下水普遍处于偏碱性环境，且酸碱度空间差异不大。矿化度变化范围为185.2~3278.22mg/L，表明矿化度指标空间分布差异较大。水质为淡水-咸水之间。其中一处水样中NO3–毫克当量百分数超过25%（表 1），对水化学类型产生较大影响。
表 1 地下水测试数据Table 1. The list of groundwater test data样号 地貌类型 水温/℃ pH 井水埋深/m 含水层类型 总矿化度mg/L NH4+mg/L NO3–mg/L NO2–mg/L 水化学类型 J1 风积沙滩 16 7.73 65 承压水 491.3 0.06 15.7 < 0.01 HCO3-SO4-Na-Mg-Ca J2 风积沙滩 15 8.53 70 承压水 569.32 0.07 20.2 < 0.01 HCO3-SO4-Na-Ca-Mg J3 风积沙滩 14 7.96 100 承压水 630.02 0.06 6.71 < 0.01 HCO3-Mg-Na-Ca J4 风积沙滩 18 7.71 0 泉水 826.62 0.08 80.8 0.74 HCO3-SO4-Na-Mg-Ca J5 风积沙滩 17 7.89 5 潜水 1368.28 0.12 200 < 0.01 SO4-HCO3-Cl-Na-Mg J6 风积沙滩 17 7.95 2 潜水 1038.91 0.11 234 0.92 HCO3-NO3-Na-Mg J7 风积沙滩 16 8.01 4 潜水 1538.07 0.11 67.1 0.49 HCO3-Cl-SO4-Na-Mg J8 风积沙滩 12 7.94 4.3 潜水 1668.22 0.1 171 2.01 Na-HCO3 J9 冲积平原 7.13 1 潜水 185.2 0.06 7.97 < 0.01 HCO3-SO4-Ca-Mg J10 冲积平原 17 8.51 潜水 1622.24 0.08 43.2 0.84 Cl-SO4-HCO3-Na-Mg J11 冲洪积平原 16 9.63 2 潜水 1549.54 0.1 6.3 0.69 Cl-SO4-Na-Mg J12 冲洪积平原 8.19 3.5 潜水 3278.22 0.1 79.6 < 0.01 Cl-SO4-Na-Mg 由表 1、图 2可以看出，研究区13组水样中三氮NH4+、NO3–、NO2–水化学特征，含量范围分别为0.06~0.12mg/L、4.67~234mg/L、 < 0.01~2.01mg/L，与国家地下水质量标准Ⅲ类水质限值对比，氨氮均检出，但未污染；NO2–均检出，其中6组含量为重度或极严重污染；NO3–的7组样品未检出，其余5组NO3–污染级别为中度、轻度污染。空间上无论矿权范围还是矿权外，污染样点均存在，不同的是NO3–污染主要在第四系潜水含水层中，而NO3–污染已经迁移至承压水中。超标样点占调查样点的75%。
![]()
图 2 宁东梅花井地下水三氮污染评价Figure 2. Evaluation graph of groundwater three nitrogen pollutants in Meihuajing field, Ningdong3.2 矿区地下水氮物源
本次采集了煤矸石和区域土壤，pH为碱性。煤矸石产生量约占煤炭开采量的5%~25%[20]。从表 2可知，区域土壤氮含量平均值为346mg/kg，氮含量为5级，氮养分缺乏；煤矸石氮含量平均值为1213mg/kg，氮含量为3级，氮养分足量。矸石中含有丰富的N元素，即矿业开发为地下水氮污染提供了丰富的氮源。刘钦甫等[21]在煤矸石中氮溶出的动态淋滤实验中，得出煤矿区地表矸石中含有吸附状态的硝态氮，而这部分氮较容易被水所溶出，特别是中性水比酸性水更容易使其溶出。在3种形态的氮中，硝态氮明显占优势，且其变化趋势与总氮一致。
表 2 煤矸石、表土中N元素含量Table 2. The content of nitrogen in gangue and soil物源类型 样品编号 N pH 附近布设井孔 GS1 301 8.7 GS2 1556 7.58 研究区内 GS3 1406 8.02 GS4 257 6.78 GS5 519 8.31 GS6 424 8.70 GS7 899 8.15 GS8 859 8.23 GS9 1523 8.12 煤矸石 GS10 173 9.22 GS11 2382 8.15 研究区外 GS12 1668 7.78 GS13 703 8.39 GS14 2331 7.89 GS15 3123 7.52 GS16 2131 7.13 GS17 1595 7.14 GS18 770 10.97 GS19 418 8.36 煤矸石平均值 1213 8.17 土壤 T1 472 9.21 J3 T2 290 9.34 J5 T3 440 9.26 J2 T4 386 9.41 J1 T5 170 8.97 J7 T6 320 8.95 J11 土壤平均值 346 9.19 中国土壤普查技术含量分级[22] 6级 < 50 5级 50~750 4级 750~1000 3级 1000~1500 2级 1500~2000 1级 >2000 粉煤灰场旁水样取自电厂灰渣贮灰场旁边的水塘，废水排放标准暂无对硝酸盐的规定，本文参照国家地表水环境质量标准（GB3838—2002）和国家地下水环境质量标准（GB/T 14848—93）Ⅲ类水限值。从表 3可知，淋滤液、河水、湖水、降雨中含有硝酸盐和氨氮成分，但相对国家标准其水质是安全的。
表 3 粉煤灰淋滤液、地表水、雨水中N元素含量Table 3. The content of nitrogen in fly ash filter, surface water and rain water水样类型 样品编号 NH4+ NO3– NO2– 粉煤灰场旁水样 F1 0.06 10.8 < 0.01 F2 0.07 8.09 < 0.01 F3 0.07 11.4 < 0.01 河水 H1 < 0.02 1.36 < 0.01 H2 < 0.02 4.09 < 0.01 湖水 H3 0.07 1.9 < 0.01 雨水 Y1 0.04 2 0 国家地表水环境质量标准（GB38382002）Ⅲ类水限值 1 10 研究区为牧区，大量的牧羊粪便随着地表或雨水进入地下水，施用动物粪便等的地区，地表水回补后的地下水NO3–值较高；地下水中N元素的浓度逐年增加[22-24]。
4. 讨论
地下水中三种主要的无机氮形态NH4+、NO3–、NO2–，以NO3–、NO2–为主，其中地下水中的NO3-占无机氮的比例最高，氨氮未出现超标，超标成分为NO3–、NO2–。那么地下水中NH4+不超标而NO3–、NO2–超标的影响因素是什么？
首先，从地下水补给-径流-排泄条件考虑：①垂向补给条件，由前述可知，目前淋滤液、河水、湖水中三氮相对国家标准水质是安全的，引起地下水硝酸盐和亚硝酸盐超标的可能性较小。但是由于有丰富的物源，在地表水或降雨长期淋滤作用下，沿垂向可能引起地下水中硝酸盐和亚硝酸盐的超标。②受地形地貌影响明显，地下水补给区位于夹山梁和布朗山，研究区只有J9井孔位于丘陵高地，水质良好，而处于经过长期的蒸发作用和溶滤作用的排泄区的样点，如J8、J10、J11井点中，NO3和NO2–含量增加，出现超标现象，J8井点中的NO3–、NO2–超标倍数分别为99.5倍、7.55倍，J10井点中的NO3–、NO2–超标倍数分别为41、1.16倍，J11井点中的NO3–超标倍数为33.5倍。③煤矿采掘扰动改变含水层围岩和径流条件，也会促进N的释放。如J1与J2井点具相似的地形条件，不同的是J2位于矿权范围，井孔中NO3–含量出现超标现象，而在矿权范围外的J1井水质良好，煤炭资源在开采过程中，煤层、围岩中的氮化物矿物与氧气和水接触，在微生物的催化作用下，经过一系列复杂的地球化学反应，可促进N元素的释放。由于矿井水在地下与围岩裂隙水存在着一定的水力联系，这些含氮物质的释放将会对水环境造成严重的威胁[25-28]。
其次，地下水水文地球化学特性的众多研究表明，地下水氮污染形式主要是NO3–的污染，其是国内外最普遍、污染面积最大的地下水污染问题[29-31]，在NH4+、NO3–、NO2–系统中三氮变化关系通常总是向NO3–转化，而NO2–浓度增高的罕见情况也只是具有暂时性意义[6]。因此长期以来，NO2–在潜水中出现增高甚至超标，取决于复杂的地下水中氮的转化过程，主要包括有机氮的矿化作用、NH4+吸附作用、硝化作用、反硝化作用、异化还原等[31]。同时温度、pH、Eh、土层介质厚度、土壤透气性等也是重要的影响因素[25, 32-35]。Sitaula等[35]认为，硝化反应的最佳pH值范围为8~8.4，pH值低于7时，硝化速率明显降低，低于6或高于9.6时，硝化反应几乎停止，研究区pH值范围为7.73~9.63，是硝化反应的最佳区域。王晓娟等[36]认为，在氮的转化过程中，细菌起着重要的作用，几乎所有微生物在好氧环境转变为厌氧环境后均可参与硝化过程；有研究表明，某些反硝化细菌在好氧条件下也可以进行反硝化[37]，一定温度范围内，氮矿化（氮由有机态转化为无机态NH4+或NH3的过程）随温度的升高而升高，随土层深度增加而降低，随土壤通透性的降低而降低，这也可能是潜水中NO3–、NO2–含量高而承压水中只有NO3–含量高的原因。
5. 结论
（1）NO3–、NO2–与国家地下水质量标准Ⅲ类水质限值对比污染严重，NO2–污染主要在第四系潜水含水层中，而NO3–在潜水和承压水中均有污染。研究区超标样点占调查样点的76.92%。根据物源和氮污染空间分布情况推测，煤矿区氮的污染很可能是研究区高地球化学背景引起的。
（2）NO3–在研究区浅层地下水中的相对含量较高，NO3–毫克当量百分数超过25%的水样，对水化学类型产生影响。
（3）NO3–、NO2–的影响因素很大程度受限于煤矿开采、地形地貌条件、垂向补给及水文地球化学条件。低山丘陵高地及煤矿开采未影响到的地方，水质较好。温度、pH、Eh、土层介质厚度、土壤透气性等也是重要的影响因素。
-
![]()
图 5 潼湖盆地解剖区第一期(a)和第二期(b)古河道分布
Figure 5. Distribution of the first(a) and the second(b) paleochannel in the anatomic area of Tonghu Basin
-
Miall A D. Architectural-element analysis: A new method of facies analysis applied to fluvial deposits[J]. Earth Science Reviews, 1985, 22(4): 261-308. doi: 10.1016/0012-8252(85)90001-7
裘亦楠, 肖敬修, 薛培华. 湖盆三角洲分类的探讨[J]. 石油勘探与开发, 1982, (1): 1-11. https://www.cnki.com.cn/Article/CJFDTOTAL-SKYK198201000.htm 吴崇筠. 构造湖盆三角洲与油气分布[J]. 沉积学报, 1983, (1): 5-23. https://www.cnki.com.cn/Article/CJFDTOTAL-CJXB198301002.htm Purkait B, Majumdar D. Distinguishing different sedimentary facies in a deltaic system[J]. Sedimentary Geology, 2014, 308(7): 53-62.
Deveugle P E K, Jackson M D, Hampson G J, et al. Characterization of stratigraphic architecture and its impact on fluid flow in a fluvial-dominated deltaic reservoir analog: Upper Cretaceous Ferron Sandstone Member, Utah[J]. AAPG Bulletin, 2011, 95(5): 693-727. doi: 10.1306/09271010025
裘怿楠. 中国陆相碎屑岩储层沉积学的进展[J]. 沉积学报, 1992, 10(3): 16-23. https://www.cnki.com.cn/Article/CJFDTOTAL-CJXB199203003.htm Bates C D. Rational theory of delta formation[J]. American Association of Petroleum Geologists Bulletin, 1953, 37: 2119-2162.
Halfar J, Riegel W, Walther H. Facies architecture and sedimentology of a meandering fluvial system: a palaeogene example from the Weisselster Basin, Germany[J]. Sedimentology, 1998, 45(6): 1-17.
Wu C L, Bhattacharya J P, Ullah M S. Paleohydrology and 3D facies architecture of ancient point bars, Ferron sandstone, Notom Delta, South-central Utah, U.S. A[J]. Journal of Sedimentary Research, 2015, 85(4): 399-418. doi: 10.2110/jsr.2015.29
吴胜和, 纪友亮, 岳大力, 等. 碎屑沉积地质体构型分级方案探讨[J]. 高校地质学报, 2013, 19(1): 12-22. doi: 10.3969/j.issn.1006-7493.2013.01.004 梁宏伟, 穆龙新, 范子菲, 等. 基于正演模拟的曲流河构型层次研究[J]. 中国矿业大学学报, 2014, 6: 8. https://www.cnki.com.cn/Article/CJFDTOTAL-ZGKD201406017.htm Coleman J M, Wright L D, Modern river deltas: variability of processes and sand bodies[C]//Broussard M L. Deltas, Models for Exploration: Houston, Houston Geological Society, 1975: 99-149.
Orton G. A spectrum of Middle Ordovician fan deltas and braid plain deltas, North Wales: a consequence of varying fluvial clastic input[C]//Nemec W, Steel R. Fan-deltas: Sedimentology and tectonic settings. London: Blackie and Son, 1988: 23-49.
Orton G J, Reading H G. Variability of Deltaic Processes in Terms of Sediment Supply, with Particular Emphasis on Grain Size[J]. Sedimentology, 1993, 40(3): 475-512. doi: 10.1111/j.1365-3091.1993.tb01347.x
Bhattacharya J P, Davies R K. Sedimentology and structure of growth faults at the base of the Ferron Member along Muddy Creek, Utah[C]//Chidsey T C, Morris T H. Regional to wellbore analog for fluvial-deltaic reservoir modeling: The Ferron Sandstone of Utah: American Association of Petroleum Geologists, Studies in Geology, 2004, 50: 279-304.
Olariu C, Bhattacharya J P. Terminal distributary channels and delta front architecture of river-dominated delta systems[J]. Journal of Sedimentary Research, 2006, 76(2): 212-233. doi: 10.2110/jsr.2006.026
Gani M R, Bhattacharya J P. Bedding correlation vs. facies correlation in deltas: Lessons for Quaternary stratigraphy[C]//Giosan L, Bhattacharya J P. River Deltas: Concepts, Models, and Examples. SEPM, 2005, 83: 31-47.
李伟荣. 潼湖湿地地质环境演变特征[C]//中国地质学会2007年生态环境脆弱区工程地质学术论坛, 2007. 王涛, 徐怀民, 江同文, 等. 滨岸相砂质输导层非均质性对油气运聚影响研究: 以塔里木盆地哈得逊油田东河砂岩为例[J]. 高校地质学报, 2018, 24(1): 128-138. https://www.cnki.com.cn/Article/CJFDTOTAL-GXDX201801013.htm Leeder M R. Fluviatile fining-upwards cycles and the magnitude of palaeochannels[J]. Geological Magazine, 1973, 110(3): 265-276. doi: 10.1017/S0016756800036098
下载:
计量
- 文章访问数: 3097
- HTML全文浏览量: 388
- PDF下载量: 1804
