-
国内地震地下流体观测实践中,同井水温和水位的对比观测已有30多年历史。分析二者之间的影响关系,对认识井水位、水温的动态及其机理[1-2]、了解应力加载作用引起的地下介质变化[3],甚至捕捉地震前兆异常等[4]都具有重要的意义。根据前人研究,目前冷水下渗说、热弥散说是解释同井观测水位、水温同步动态变化的主要机理;水位、水温影响因素主要有:仪器观测系统的稳定性、观测环境的干扰(抽水和注水等)以及气象三要素(气温、气压、降雨)的影响,都会造成水位和水温反映构造应力的伪变化[5]。
根据中国台网中心震台网函2016(258)号文件要求,为了进一步提高井水温观测网运行质量与监测效能,规范井水温观测的技术关键环节,重新对观测效能较高的昌邑鲁02井进行水温加密梯度测量。2017年4月23日开展昌邑02井加强密度观测水温梯度测试验,同年8月进行了水位、水温仪器的更换;2018年8月上旬至9月上旬,该井水位和水温出现准同步的快速下降后,很快回升,接着又持续下降持续时间较长,该同步变化是观测环境造成,还是与井孔结构有关,目前尚不能确定。为了进一步判定其动态变化影响因素和性质,本文收集整理了2017年昌邑02井加强密度观测水温梯度测试验及区域水文地质条件、井孔基础资料、区域地下水开采和降雨等资料,分析了该井水位和水温动态变化的影响,从水位和水温同步动态变化机理上进行了讨论。
-
昌邑02井周边地势自南而北由高到低呈阶梯状,形成山丘、平原、滨海3个地貌单元[6]。该井位于山东潍河西岸,向北18 km临莱州湾,处于昌邑-大店断裂带西侧,海拔7.6 m;地下水富集形成构造单元,主要包括松散的岩石孔隙水和基岩裂隙水;昌邑区域地下水主要补给来源包括河流入渗补给、大气降水渗水、越流补给、侧向径流地下水补给以及地表灌溉的回渗补给(图1~2),排泄的方式主要是抽水,包括城市供水量、工业用水量、蒸发的地下水量、及农业用水开采地下水量等[6]。
-
昌邑井原属胜利油田钻探报废井,1978年在征得油管系统同意后,被山东省地震局改用为宏观观测井,编号为鲁02号。为防止杂物落入井孔,1978年5—8月对井孔进行改造,由井口焊接直径140 mm钢管,高出地面8 m,1984年9月进行SW-40型地下水观测。成井深度1 383.5 m,套管深1 171.90 m。1998年数字化改造,分别安装LN-3数字水位仪观测水位、SZW-1A型数字式温度计观测水温(水温探头置于井下400 m左右),2017年更换为地壳所DRSW-2水位水温观测仪。由于井水位持续下降,水位探头下放到9 m左右,水温探头置于井下350 m左右,采样率均为分钟值。2017年之前观测含水层400 m以下的全井段为第三系砂岩孔隙水(图1)。
-
2009年以来,该井水位至今仍保持趋势下降变化(图3a),通过井水位与地下水开采量动态对比统计结果(表1)看,地下水超采区地下水位平均埋深呈逐年增大的趋势,水位多年动态变化(趋势下降)受地下水开采有关[6]。该井水温在2009—2014年间一直趋势上升,平均水温22.54 ℃,属于上升漂移型动态(图3b)。根据前人研究[7],对于水温多年上升漂移型动态的成因一般有3种:①大地热流值的微弱改变在水温上的反映;②由于围岩在缓慢线性增加的应力场作用下,导致孔隙水上升引起观测点处温度的变化;③仪器“零漂” 引起的变化。结合全省观测井水温趋势变化(观测仪器系统事件)以及通过对该井水位和水温的动态变化曲线对比,认为第3种情况的可能性较大,即水温仪器“零漂”;同井观测水位、水温长趋势动态变化是水位趋势下降、水温趋势上升。
表 1 昌邑02井2009—2014年水位与地下水开采量[6]
年份 开采量/ m3 地下水位平均埋深/m 2009 51 791.3 1.767 2010 55 117.3 1.867 2011 58 496.0 2.062 2012 61 819.0 2.227 2013 65 019.3 2.522 2014 68 119.3 3.010 -
将昌邑水位整点值和水温日均值曲线,同时结合区域降雨量进行分析发现(图4),水位、水温季节性变化明显,在强降雨(夏季)时存在水位同步上升接着下降的突变形态;水位整体趋势向下、水温整体趋势向上,水温平均温度21.74 ℃。2017年更换新仪器后,水位水温趋势并没有改变,说明水位主要受地下水开采影响,水温长期向上不仅仅与仪器零点漂移相关。昌邑井周围观测环境周边有大片农田和村庄,农业和生活抽水多在150 m左右,可能会对观测井有一定的同层抽水干扰。水位在夏季降雨补给不明显(季节性小幅度上升),但水位年数据动态呈明显下降变化(图4)。当浅层冷水被开采,下层热水流入井筒内的比例就会增加,水温出现上升变化,可以用冷热水混合机理解释[8]。
2018年1月30日和3月30日水位下放探头造成水温突降,2019年3月30日、6月30日、9月30日、12月30日水位校测造成水温突降,12月6日氡背景值取水样造成水温干扰。对于2018年8月10日和9月20日水位下降、同步水温下降动态变化,经台站巡检未发现周围环境或现仪器异常。根据冷热水交换机理,水温迅速下降有2种可能:大量冷水下渗或深层抽水。由图4a~4b可见水温、水位同步下降台阶,昌邑井周围有油田勘探钻井较多,怀疑是周围油田勘探深层抽水有关。
-
根据中国台网中心震台网函2016(258)号文件要求,为了进一步提高井水温观测网运行质量与监测效能,规范井水温观测的技术关键环节,重新对观测效能较高的昌邑鲁02井进行水温加密梯度测量。2017年4月利用SZW-1A型数字式温度计进行昌邑井加密温度梯度试验,不同层位水温的测量方式是:井深365 m内,每20 m —个测点,共17个测点;365~401 m每4 m一个测点,共10个测点,同时对365~401 m深度位置按需要加密观测。每个测点持续观测30 min以上,对于365~401 m深度测点观测2 h以上。根据观测试验结果计算测点深度变化特征值见表2。昌邑试验分27层进行的测试根据水温分层试验测量分钟数据分析,绘制了昌邑水温梯度图、均方差图、以及幅度图(图5)。
表 2 昌邑02井水温试验各深度变化特征值
深度 / m 温度 / ℃ 一阶差分 变化幅度 最大日变差 均方差 深度差分/m 梯度 25 14.643 2 −0.031 150 −1.322 5 0.017 7 0.096 409 0 45 14.367 0 −0.007 500 −0.236 6 0.041 1 0.048 541 20 −0.037 500 65 14.974 1 0.020 293 0.600 2 −0.007 9 0.154 033 20 0.101 467 85 15.304 0 0.010 970 0.285 4 −0.043 0 0.068 536 20 0.054 850 105 15.601 1 0.010 234 0.278 1 −0.018 3 0.068 391 20 0.051 172 125 16.000 4 0.000 440 0.387 8 −0.043 5 0.018 562 20 0.002 201 145 16.404 9 0.012 252 0.403 0 −0.001 1 0.108 637 20 0.061 258 165 16.812 0 0.013 380 0.401 2 0.000 5 0.103 023 20 0.066 900 185 17.207 0 0.013 658 0.778 3 −0.027 3 0.103 023 20 0.068 288 205 17.618 1 0.014 616 0.383 8 −0.001 0 0.113 627 20 0.073 081 225 18.067 5 0.012 765 0.326 9 −0.002 5 0.097 126 20 0.063 823 245 18.465 8 0.017 272 0.522 6 −0.083 3 0.102 230 20 0.086 362 265 18.968 2 0.017 728 0.377 5 −0.021 2 0.121 254 20 0.088 638 285 19.481 6 0.012 952 0.257 7 −0.082 6 0.081 541 20 0.064 758 305 19.882 0 0.013 206 0.368 7 −0.054 3 0.088 804 20 0.066 031 325 20.306 6 0.018 343 0.491 3 −0.015 4 0.119 683 20 0.091 714 345 20.812 7 0.019 263 0.441 5 −0.006 3 0.110 261 20 0.096 315 365 21.339 3 0.001 400 0.100 5 −0.006 0 0.016 947 20 0.007 000 369 21.446 3 0.001 583 0.077 4 −0.009 6 0.014 144 4 0.039 569 373 21.533 4 0.001 238 0.075 5 −0.004 7 0.014 726 4 0.030 943 377 21.613 5 0.001 384 0.069 8 0.000 9 0.016 506 4 0.034 590 381 21.697 9 0.001 264 0.061 9 −0.013 1 0.012 627 4 0.031 610 385 21.773 0 0.001 315 0.078 7 −0.001 0 0.015 720 4 0.032 869 389 21.852 3 0.001 394 0.072 9 −0.014 7 0.013 696 4 0.034 844 393 21.940 7 0.001 964 0.089 1 −0.020 4 0.017 140 4 0.049 107 397 22.051 5 0.001 757 0.098 7 −0.003 1 0.017 891 4 0.043 934 401 22.153 6 0.000 161 0.016 6 −0.005 2 0.006 394 4 0.004 018 昌邑井下50 m深度段左右,水温由负梯度变为正梯度(图5d),推测井管可能长年受井水腐蚀,封闭性变差,该深度段内可能有潜水混入井孔所致[9];昌邑井水温日变幅从320 m深度向下总体上有逐步减小的变化特征,362 m深度及以下又逐步减小。在369~385 m深度间均在0.07 ℃ 及以内,变化较稳定。该井孔潜水段处变化幅度较大,体现了含水层潜水段冷热交换作用较强烈的特点。参考值日变差主要体现水温的稳定性(图5,表2),在365~401 m深度范围内,昌邑井水温日变差均小于0.009 ℃较为稳定。结合均方差图5c计算结果分析,在深度369~401 m,均方差均小于0.018 ℃,观测数据的离散程度最小、稳定性最好。由图5d看出,昌邑井水温从上至下梯度变化大致可分为3段:第1段是150 m以上;第2段是150~350 m;第3段则是350 m以下。第1段: 水温梯度变化幅度较大,其温度明显偏低,而其日变幅、日变差和均方差则明显突出偏大;第2段是150~350 m,变化明显比第1段趋于稳定,水温梯度变化具有线性化特征;第3段:350~400 m,水温梯度显示了与上面两层段完全不同的变化特征,随深度增加温度升高较快,水温梯度以及变化幅度、均方差变化趋于平稳,线性关系较强。
-
为了去除水温趋势和随机成分(毛刺)影响,选取2019年4月13—24日观测数据,利用最小二乘法拟合水温动态曲线(图6)。经数据处理后的水温动态曲线中,峰、谷出现时间与固体潮、水位潮汐的日变起伏在形态、幅度与相位上均表现出一致性,与井水位潮汐变化吻合,显示水温观测数据具有一定潮汐特征,水温潮差变幅0.005 56 ℃(图6)。
随着“十三五”仪器安装和改造的实施,高精度数字水温观测技术的不断升级,国内部分观测井记录到水温固体潮现象,即水温动态曲线中峰、谷出现时间与重力理论固体潮曲线中一致,认为水温潮汐效应[10]。
本文利用最小二乘法拟合处理后的水温数据,与水位计算相关系数为0.482 0,水位与固体潮的相关系数0.581 0。推测水温的峰谷变化是水位波动引起的次生效应,即在潮汐力作用下含水层发生变形,井含水层系统中水流状态或流量发生变化,由含水层带入井筒的热量也相应发生变化。
马玉川等[11]通过研究井水温潮汐认为:由潮汐力作用于含水层发生膨胀变形时,含水层上部水向下运动使水温下降(随季节性变化);当含水层受到潮汐引力作用而压缩变形时,下部水向上运动则使水温上升,认为水温变化是含水层受力变形结果,可称为岩体力学机制。
昌邑02井所处区内断裂构造发育,这些断裂可能沟通了不同地层中的含水层,使井区的地下水之间具有良好的水力联系。基于水温加密观测试验分析和水位、水温的多年动态趋势以及年变化影响因素分析,也佐证了这点。同井观测水位、水温同步变化影响主要体现在地下水开采引起水位大幅度下降和水温持续上升,以及强降雨引起水位、水温同步向上动态变化。昌邑水温主要受降雨和仪器干扰影响,相比水位,水温目前受地下水开采影响较小,水温动态曲线变化较稳定,结合加密水温试验,与探头放置位置属于梯度平稳区有关。
-
通过上述昌邑02井观测、水温加密试验等数据以及区域地下水开采、降雨等资料,探讨了该井水位和水温动态变化影响因素,获得了如下认识。
1)基于水温加密观测试验,结合区域水文地质条件、井孔基础资料、区域地下水开采和降雨等资料,发现水位易受浅层水干扰,水位水温同步变化影响主要是地下水开采引起水位大幅度下降和水温持续上升,以及强降雨引起水位、水温同步向上动态变化,同时例证该井地下水类型为混合水。
2)鉴于该井水温多年来一直趋势上升,一直被认为可能由仪器“零漂”引起的变化。2017年更换新仪器后,水位水温并没有改变向上的趋势,基于冷热水理论,考虑水温主要受地下水开采影响。
3)对水温实验资料分析的目的就是选择较好的水温观测层位,提高观测效能。对于这一问题,中国地震台网中心震台网函2016(258)号文件中印发的《观测井水位校测与水温梯度测量要求(试行)》提出了梯度大、背景噪声小、潮汐效应明显的技术要求,昌邑井应选择在350~400 m。综上来看,深度应选择在400 m上下比较理想。但2017年8月更换仪器水温探头后,水温显示平均温度22.54~21.74 ℃,提升50~350 m左右,数据质量效能受影响较小,数据稳定且能观测到潮汐效应。
4)利用最小二乘法拟合水温动态曲线和相关性分析,发现与井水位潮汐变化基本一致,显示水温观测数据具有一定潮汐特征,推测水温的峰谷变化是水位波动引起的次生效应。水温固体潮特征表明,该井水温受岩体动力学及地热动力学共同作用的现象,可能与该井特殊的水文地质条件有关。
Influential Factors of Changyi 02 Well Water Level and Water Temperature Dynamic Change
-
摘要: 基于昌邑02井观测、水温加密试验等数据以及区域地下水开采、降雨等资料,探讨了该井水位和水温动态变化影响因素;同时,利用最小二乘法拟合水温数据动态变化发现:拟合后的水温动态曲线在形态、幅度与相位上与水位动态变化均呈现出相关性,昌邑02井具有潮汐效应特征,并且与该井水位的潮汐动态同步;探讨02井水温潮汐变化机理,以及水温水位之间的水力联系。Abstract: Based on the data of observation and water temperature infill test of Changyi 02 well, as well as the data of regional groundwater exploitation and rainfall, the influential factors of dynamic change of water level and water temperature were discussed. At the same time, the least square method is used to fit the dynamic changes of water temperature data. It is found that the fitted water temperature dynamic curve is correlated with the dynamic changes of water level in terms of shape, amplitude and phase. Changyi 02 well has the characteristics of tidal effect, and synchronizes with the tidal dynamics of the water level of the well.
-
Key words:
- water level /
- water temperature /
- water temperature encryption test /
- tide force
-
表 1 昌邑02井2009—2014年水位与地下水开采量[6]
年份 开采量/ m3 地下水位平均埋深/m 2009 51 791.3 1.767 2010 55 117.3 1.867 2011 58 496.0 2.062 2012 61 819.0 2.227 2013 65 019.3 2.522 2014 68 119.3 3.010 表 2 昌邑02井水温试验各深度变化特征值
深度 / m 温度 / ℃ 一阶差分 变化幅度 最大日变差 均方差 深度差分/m 梯度 25 14.643 2 −0.031 150 −1.322 5 0.017 7 0.096 409 0 45 14.367 0 −0.007 500 −0.236 6 0.041 1 0.048 541 20 −0.037 500 65 14.974 1 0.020 293 0.600 2 −0.007 9 0.154 033 20 0.101 467 85 15.304 0 0.010 970 0.285 4 −0.043 0 0.068 536 20 0.054 850 105 15.601 1 0.010 234 0.278 1 −0.018 3 0.068 391 20 0.051 172 125 16.000 4 0.000 440 0.387 8 −0.043 5 0.018 562 20 0.002 201 145 16.404 9 0.012 252 0.403 0 −0.001 1 0.108 637 20 0.061 258 165 16.812 0 0.013 380 0.401 2 0.000 5 0.103 023 20 0.066 900 185 17.207 0 0.013 658 0.778 3 −0.027 3 0.103 023 20 0.068 288 205 17.618 1 0.014 616 0.383 8 −0.001 0 0.113 627 20 0.073 081 225 18.067 5 0.012 765 0.326 9 −0.002 5 0.097 126 20 0.063 823 245 18.465 8 0.017 272 0.522 6 −0.083 3 0.102 230 20 0.086 362 265 18.968 2 0.017 728 0.377 5 −0.021 2 0.121 254 20 0.088 638 285 19.481 6 0.012 952 0.257 7 −0.082 6 0.081 541 20 0.064 758 305 19.882 0 0.013 206 0.368 7 −0.054 3 0.088 804 20 0.066 031 325 20.306 6 0.018 343 0.491 3 −0.015 4 0.119 683 20 0.091 714 345 20.812 7 0.019 263 0.441 5 −0.006 3 0.110 261 20 0.096 315 365 21.339 3 0.001 400 0.100 5 −0.006 0 0.016 947 20 0.007 000 369 21.446 3 0.001 583 0.077 4 −0.009 6 0.014 144 4 0.039 569 373 21.533 4 0.001 238 0.075 5 −0.004 7 0.014 726 4 0.030 943 377 21.613 5 0.001 384 0.069 8 0.000 9 0.016 506 4 0.034 590 381 21.697 9 0.001 264 0.061 9 −0.013 1 0.012 627 4 0.031 610 385 21.773 0 0.001 315 0.078 7 −0.001 0 0.015 720 4 0.032 869 389 21.852 3 0.001 394 0.072 9 −0.014 7 0.013 696 4 0.034 844 393 21.940 7 0.001 964 0.089 1 −0.020 4 0.017 140 4 0.049 107 397 22.051 5 0.001 757 0.098 7 −0.003 1 0.017 891 4 0.043 934 401 22.153 6 0.000 161 0.016 6 −0.005 2 0.006 394 4 0.004 018 -
[1] 鱼金子, 车用太, 刘五洲. 井水温度微动态形成的水动力学机制研究[J]. 地震, 1997, 17(4): 389-396. [2] 车用太, 刘喜兰, 姚宝树, 等. 首都圈地区井水温度的动态类型及其成因分析[J]. 地震地质, 2003, 25(3): 403-420. doi: 10.3969/j.issn.0253-4967.2003.03.006 [3] 杨竹转. 地震波引起的井水位水温同震变化及其机理研究[D]. 北京: 中国地震局地质研究所, 2011. [4] 曹新来, 王贺生, 张素欣, 等. 昌黎井含水层系统的水位水温动态关系与地震活动[J]. 华北地震科学, 2001, 19(2): 55-61. doi: 10.3969/j.issn.1003-1375.2001.02.009 [5] 丁风和, 韩晓雷, 哈里白, 等. 通辽井水位和水温动态变化影响因素分析[J]. 华南地震, 2017, 37(1): 49-53. [6] 张辉, 董敏, 马娟, 等. 山东昌邑区深层地下水开采与地面沉降关系分析[J]. 地震工程学报, 2018, 40(S): 214-218. [7] 李仲巍, 梁国经, 郑双凤, 等. 抚松井水温动态特征分析[J]. 防灾减灾学报, 2012, 28(1): 75-79. doi: 10.3969/j.issn.1674-8565.2012.01.015 [8] 车用太, 王吉易, 李一兵, 等. 首都圈地下流体监测与地震预测[M]. 北京: 气象出版社, 2004. [9] 杨竹转, 邓志辉, 杨贤和, 等. 井孔水温动态变化的影响因素探讨[J]. 地震, 2010, 30(2): 71-79. doi: 10.3969/j.issn.1000-3274.2010.02.009 [10] 张子广, 张素欣, 李薇, 等. 昌黎井水温潮汐形成机理分析[J]. 地震, 2007, 27(3): 34-40. doi: 10.3969/j.issn.1000-3274.2007.03.005 [11] 马玉川、王博. 祁县井水温和水位的奇异潮汐关系探讨[J]. 中国地震, 2014, 30(1): 55-63. doi: 10.3969/j.issn.1001-4683.2014.01.006 -