Mechanisms of water-heat dynamics and earth-heat dynamics of well water temperature micro-behavior
-
摘要: 以井-含水层热系统分析为基础， 结合井水温度(水温)动态观测到的实际结果， 提出了水温微动态形成的两类基本机制， 即水热动力学机制与地热动力学机制. 水热动力学机制指井水温度的动态是由于水流动产生的热对流引起的变化； 地热动力学机制指井区岩土中大地热流作用或热传导引起的井水温度变化. 结果表明， 多数同震响应是水热动力学机制下生成的水温微动态； 水温的震后效应或变化则多是热传导或深部热流变化引起的， 属于地热动力学机制下生成的水温微动态. 少数水温微动态特征较为复杂， 多是源于井-含水层系统中包含有多个水文地质特征差异明显的含水层， 尚不能简单地用水热或地热动力学机制予以解释.Abstract: On the basis of analysis on heat system of well-aquifer, combining with practical observation results of well water temperature behaviors, two mechanisms of water-heat dynamics and earth-heat dynamics of well water temperature micro-behavior are presented. The mechanism of water-heat dynamics is that the well water temperature micro-behavior originates from heat transport related to moving of water flow. The mechanism of earth-heat dynamics is that the well water temperature micro-behavior originates from heat flow or heat transfer in the earth. Results show that most coseismic responses of water temperature in a well are micro-behaviors produced under water-heat dynamics mechanism, whereas post-seismic effects or changes of water temperature in a well are micro-behaviors caused by heat conduction or change of deep heat flow, and belong to earth-heat dynamics mechanism. The characteristics of a few micro-behaviors of water temperature in a well are more complicated, which is due to existence of multilayer aquifers in a well-aquifer system. And the hydrogeological characteristics of these aquifers are more obviously different, which cannot be yet explained under the mechanisms of water-heat or earth-heat dynamics.
-
引言
流动数字地震台网作为“‘十五’中国数字地震观测网络”项目建成的新一代中国数字地震观测系统的一部分，在我国地震观测事业中发挥着重要作用(刘瑞丰等，2008； 郑秀芬等，2009)． 流动数字地震观测设备必须满足流动性要求(郭建，刘光鼎，2009; 陆其鹄等，2009)，因此组成流动数字地震观测系统的各个设备必须具有体积小、 重量轻、 功耗低等特点(游庆瑜等，2003； 李江等，2010)． 另外，流动观测现场一般都不具备交流电源，而是采用蓄电池或太阳能电池供电. 在这种情况下，观测设备的功耗越低，则耗电越少，采用容量不太大的蓄电池或太阳能电池即可工作，从而有助于降低观测成本． 当前，我国流动地震观测中使用的观测仪器，有进口的，也有国内自制的，基本都是将独自存在的地震计与地震数据采集器通过电缆连接而成(陈祖斌等，2006； 王超等，2007). 在这种情况下，进行流动观测时，需要同时运输地震计和地震数据采集器，增加了观测成本和安装成本．
目前，国际上一体化低功耗宽频带数字地震仪的种类并不多，成为正常销售产品的只有英国Guralp公司生产的CMG-40TDE． 在本项设备完成研制之前，我国还没有自行研制的用于野外流动地震观测的一体化数字地震仪． 为此，在中国地震局统一组织下，中国地震局地震预测研究所开发出了新型的集地震信号提取、 数据采集、 记录和服务为一体的低功耗宽频带数字地震仪，并于2010年5月将该仪器投入到青海玉树MS7.1地震活动现场进行流动观测．
1. 整机系统总体设计与结构
一体化低功耗宽频带数字地震仪研制项目组在经过充分的预研、 调查和研制准备后，确定了总体技术方案． 图 1是该地震仪的外观和内部实体图． 从图 1可以看出，整机系统外观结构由以下部分组成： 底座、 三分向宽频带数字机械摆、 4层电路板和机壳． 底座上安装航空插座，向外提供各种输入/输出接口． 为了减少接口数量，我们将网络通信接口和串口通信接口合二为一，只向外提供电源、 通信和授时模块3个接口. 另外，底座上还添加了3个底脚螺丝和一个水平气泡用于系统安装时调平． 宽频带数字地震仪的机械摆通过扁平电缆与上层电路板之间进行连接． 从底座伸出3根立柱，用于支撑4块电路板： 前级放大电路板(AMP板)、 反馈及采集电路板(ADC板)、 电源及控制电路板(PWR板)和CPU板． AMP板生成驱动地震计电容极板的振荡信号，并完成地震计输出小信号的放大和解调； ADC板负责将地震计输出的模拟电平转换成数字量，并完成地震计反馈控制； PWR板提供仪器所有的数字电路供电(+3.3 V)、 模拟电路供电(±12 V，±3.5 V)以及部分控制电路供电； CPU板上运行Linux操作系统，提供数据采集管理、 数据存贮、 数据网络服务等功能． 4层电路板之间分别用一定高度的有机玻璃套环进行定位，下方通过立柱上的小台阶进行定位，上方通过金属套环和顶丝固定在立柱上． 为了便于携带安装，机壳上设计有拎手，同时在拎手上增加“N”方向指示． 机壳上还设计了一个CF卡盖，方便用户在野外观测时随时更换CF卡． 下面从硬件模块、 低功耗处理和仪器校正等3个方面进行详细说明．
1.1 硬件功能设计
该地震仪对应的硬件结构如图 2所示，由4大部分组成： 宽频带地震计、 数据采集、 FPGA控制和微控制单元．
1.1.1 宽频带地震计
宽频带地震计是由一个垂直分向和两个水平分向组成的三位一体式结构的地震计. 每个分向均为一体式独立结构，可以单独拆卸下来，其工作原理如图 3所示． 当地面运动时，由电容换能器产生的电压信号经放大后进入反馈网络产生反馈电流，流经反馈线圈产生反馈力．
1.1.2 高精度数据采集
一体化地震仪数据采集部分由前置放大、 24位A/D转换、 辅助采集通道、 高稳定度晶体振荡电路等单元组成． 为了将地震计输出的微小地震信号调整成24位A/D转换器所需要的信号大小，采集器三通道前置放大由3级基本电路组成． 第一级小信号输入的放大电路采用场效应管组成的单端输入、 双端输出的差分放大电路． 使用场效应管可以提高输入阻抗，达到与电容换能器极板匹配的目的． 第二级解调电路部分通过控制信号对高速模拟开关MAX353进行通道切换，从而达到调制目的． 第三级放大电路由运算放大器LT6011和电容、 电阻构成，形成差分平衡输出． 实际使用中，该电路仅仅保留了直流增益．
24位A/D转换器采用TI公司生产的ADS1281，其内部集成了可配置的数字滤波器． 在电路连接上，采用引脚设置方式. 在4.096 MHz的工作时钟下，第一级SINC滤波器的抽取比为16; 第二级为最小相位FIR数字滤波，抽取比为32. ADS1281的输出采样率为2 kHz．
辅助通道采样率固定为10 Hz，用于电源电压量监测和宽频带地震计三分向零位监测. 所采用的ADS7822是一种12位串行高速、 采集速率可达75 kHz的微功耗ADC芯片． ADS7822的输入端为“多选一”模拟开关，输入端采用高阻值电阻分压网络实现±10 V的采集量程. 由于ADS7822为单电源供电设计，模拟输入端电压范围不能超过供电电压范围，因此需要通过运算放大器对采集到的模拟量进行电平转化及平移． 采集到的数字信号通过SPI接口送到CPU板．
晶体振荡器主要提供数据采集ADS1281的采样时钟和可编程门阵列FPGA的时钟计数脉冲． 选择稳定度优于10-6的TCXO模块，并带有压控输入端，在对钟的状态下，对晶振进行频率调节． TCXO晶振采用T11A(8.192 MHz)，其工作电压为3.3 V DC±5%，频率稳定度≤±1.0×10-6，控制电压变化范围为0.5—2.5 V(中心电压1.5 V DC)，频率牵引范围≥±10×10-6，斜率为正． 由12位DAC7513提供TCXO的控制电压，通过CPU的SPI接口进行控制．
1.1.3 FPGA控制逻辑
一体化地震仪对数据采集通道、 GPS码和SSC数据收发进行逻辑控制的所有功能仅由1枚FPGA芯片完成(Yu et al，2010). 该芯片采用Altera公司出产的具有20060个逻辑门的EP1C3T144C8芯片． FPGA控制逻辑采用模块化的设计方法，其中包括： ① 数据采集模块. 每个ADS1281对应一个采集模块，当检测到DRDY有效时，FPGA向ADS1281连续发送32个SCLK信号，同时读取输出的32位数据到32位移位寄存器，等32位数据读取完成时，移动该数据到锁存寄存器，并设置数据标识为“有效”，以便SSC数据上行模块获取该数据； ② IRIG-B码接收模块. 根据IRIG-B码格式对GPS输入信号进行解析并生成6个16位数据缓存到6个寄存器中； ③ 秒沿处理模块. 当PPS秒沿到时，获取4.096 MHz晶振计数值并缓存； ④ SSC数据上行模块. 定时检测每个锁存的寄存器是否有新的数据，如果有，则根据类型对数据增加数据类型标识，生成32位数据并通过SSC将数据发送到CPU； ⑤ SSC数据下行模块. 当进行地震计标定处理时，将标定数据从CPU通过SSC传输到1024字节的FIFO缓冲区，并按照每1 ms 1个数据的方式发送到DA控制器进行标定控制．
1.1.4 微控制单元
微控制单元选用Atmel公司的AT91SAM9263 CPU. 该CPU采用ARM926EJ-S核心，最高时钟为240 MHz． CPU与FPGA之间采用同步串行控制器SSC进行数据接收与命令发送操作. 由于SSC接口配置了专用外设DMA控制器PDC，显著降低了外设与存储器之间数据传输所需的时钟周期数，因此提高了微控制器的性能，数据传输效率高，CPU负荷小，不会出现采集数据丢失的现象．
微控制单元内置嵌入式Linux操作系统，同时在操作系统核心内实现了数据采集管理、 FIR数字滤波运算(同时进行1，10，20，50，100，200和500 Hz线性相位和最小相位滤波计算)、 地震计控制等功能，并通过网络接口提供实时数据服务和参数设置． 软件系统的主要功能包括： 参数设置与管理、 实时数据记录(存储于8 GB CF卡上)、 地震事件检测与记录、 网络数据服务(提供低延时网络数据传输功能，用于“预警”处理)、 主动发送服务、 串口数据服务、 标定处控制(提供脉冲、 正弦、 伪随机二进制码等3种标定信号)、 GPS时间服务、 地震计调零、 存储空间监测与管理和系统运行监控等(王洪体等，2006)．
1.2 低功耗处理
为了降低整机系统功耗，使其能够适合流动观测的要求，在仪器设计时作了以下4种处理： ① 通过修改U-Boot和Linux内核程序，将ARM CPU主频频率由原来的200 MHz降低到143 MHz； ② 增加了网络电源开断控制按钮，使得在不进行网络数据传输的情况下，用户可以关闭网络电源； ③ 授时服务由原来的连续GPS授时模式，修改为间隔式GPS授时，即当一次GPS授时完成后，自动关闭GPS电源，并在下一次授时时，重新打开GPS电源进行服务； ④ LED显示15分钟后，自动关闭其背光电源．
1.3 仪器校正
传统的地震计与地震数据采集器分开设计的方式极少考虑对仪器的输出数据进行各种校正处理，主要是因为地震计与地震数据采集器不进行成套销售，或者由于某些公司只生产地震计或者地震数据采集器，所以用户方进行数据校正难度很大． 采用一体化的设计，可以大大降低该项工作的难度，直接在仪器生产时就可以进行数据校正．
1.3.1 灵敏度校正
由于各台仪器的灵敏度往往都会有一些偏差，为了将3个通道的灵敏度都调节到500 counts/μm/s，需要对灵敏度因子进行修正． 具体方法为： 首先通过对比法(专业设备性能指标测试检测方法标准研究课题组，2011)得到地震仪的三分向灵敏度S，然后分别计算三分向灵敏度调节因子：
将上式得到的调节因子K存入仪器的存储器中，最后在获取数据时自动乘以该因子就可以得到灵敏度校正后的输出结果．
1.3.2 标准方位和正交性校正
由于地震计内部装配误差等因素，地震计实际传感方向与标示方向并不一致. 实际传感方向用xyz来表示，它们与ENU不重合，而是有一个小角度的误差，可采用如下公式进行消除：
式中，vx，vy和vz是地震计原始记录的EW，NS和UD 3个地动速度分量； v<EW，v<NS和v<UD是通过校正后输出的数据； α<NS，α<EW为地震计NS分向和EW分向的方位角； α<UD为地震计垂直分向在水平面投影的方位角; β<NS和βEW分别为地震计NS分向和EW分向相对于水平面的仰角，β<UD为地震计垂直分向偏离铅锤线的角度． 6个角度量可通过对比法或者振动台法测得的结果经过拟合后得出(专业设备性能指标测试检测方法标准研究课题组，2011)．
2. 系统测试
2.1 简介
样机完成后，我们分别于2010年3月和9月在中国地震局地震预测研究所测震实验室和河北沽源九连城地震计测试基地依据相关测试标准(Charles et al，2010； 专业设备性能指标测试检测方法标准研究课题组，2011)对其进行了各项参数指标测试． 在测试中使用的具体测试设备有：
1)超低失真度信号发生器(DS360)： 失真度为-110 dB.
2)温湿度试验箱(SETH-020L)： 温度范围为-50—100 ℃，温控精度为0.1 ℃.
3)标准时钟： GPS授时，具有时、 分、 秒脉冲输出，准确度为1 μs.
4)甚低频振动测试系统： 频率范围为0.0002—250 Hz(有中国计量科学院出具的检定证书和校准证书)．
2.2 主要参数指标测试及其结果
1)灵敏度测试. 测试时，将样机平稳地放在振动台台面中心，调整灵敏轴使其与振动方向平行，记录振动台输出信号的速度值. 从计算机记录文件读取地震仪的输出数字值，获取到的测试结果如表 1所示．
表 1 灵敏度测试结果Table 1. Test result of sensitivity测试频率/Hz 响应灵敏度/(count·μm-1·s-1) UD EW NS 5.0 503.774 502.226 506.015 2)线性度测试. 由于地震计的频率特性为带通滤波器，不能使用静态直流信号输入进行测试，只能使用正弦信号进行交流测试． 具体测试时，测试频率选择5 Hz，然后通过调整标准信号源输出幅度x，同时记录样机的输出幅值y，最后进行线性拟合并求出线性度误差，测试结果见表 2．
表 2 线性度测试结果Table 2. Test result of linearity频率/Hz 信号源输入/V UD EW NS 被测仪器输出/count 线性拟合结果/count 被测仪器输出/count 线性拟合结果/count 被测仪器输出/count 线性拟合结果/count 1 1321265 1318186 1046320 1044447 1015009 1013827 2 2640482 2637392 2091207 2088685 2028991 2027645 3 3957207 3956597 3134036 3132922 3041847 3041463 4 5274885 5275803 4176722 4177160 4055317 4055280 5 6591787 6595009 5218334 5221397 5067499 5069098 5.0 6 7910342 7914214 6261688 6265635 6080917 6082916 7 9228467 9233417 7306788 7309872 7094492 7096734 8 10549850 10552625 8352665 8354109 8109529 8110552 9 11872650 11871831 9399388 9298347 9124821 9124369 10 13199180 13191037 10448010 10442584 10141650 10138187 线性度 0.234％ 0.179％ 0.117％ 3)幅频特性测试. 测试频点选择16个，分别为0.0083，0.01，0.0167，0.025，0.05，0.1，1，5，9，19，33，39，44，49，52和58 Hz． 系统实际测得的经过归一化后的幅频特性曲线如图 4所示．
4)动态范围测试. 使用正弦波测量地震计有效观测频带内限幅电平与频率的关系． 限幅电平描述了地震计观测大信号的能力，具体测试结果见图 5． 根据测试结果可得系统的动态范围大于140 dB．
5)系统功耗测试. 按照1.2节中描述的方式对4种处理方法分别进行了功耗测试，具体测试结果见表 3． 从表 3可以看出，4种处理方式均可明显降低整机系统功耗. 在省电模式下可节省1 W左右，而整机系统功耗也从原来的1.6 W降低到了0.6 W．
表 3 4种处理方式可节省的功耗Table 3. Saved power consumption for four processing modes处理方式 可节省功耗/W 降低CPU主频 0.2 关闭网络电源 0.3 关闭GPS电源 0.2 LED显示关闭 0.3 6)多采样率同时输出和低延时输出测试. 整机系统在上电启动后，通过在台式机上同时运行多个客户端软件，并在每个软件界面上针对每个分向数据设置不同的采样率，验证了多采样率同时输出功能； 通过客户端软件设置不同的输出间隔(分别为100，200，500 ms和1 s)，验证了低延时输出功能．
2.3 与同类型仪器的比较
通过与Guralp公司生产的一体化宽频带地震仪CMG-40TDE对比，可以得出，该地震仪在频带范围、 动态范围、 系统功耗、 灵敏度等方面与CMG-40TDE非常接近． 从噪声水平方面来讲，CMG-40TDE是一款中等自噪声(Ringler，Hutt，2010)的地震仪，而该地震仪的噪声相对来说要低一些． 另外，该地震仪在软件功能方面还具有CMG-40TDE和其它地震数据采集器不具备的一些特点： 多采样率同时输出和低延时数据输出． 多采样率同时输出可以大大提高地震仪产出数据的利用率，使得地震仪在同一时刻既可以输出50/100/200 Hz数据用于常规的台站观测或者流动观测，又能输出1 Hz数据用于全球地震学研究; 而低延时数据输出方式则可以将该地震仪用于地震预警数据传输．
3. 试验研究
2010年4月14日7时49分，在青海省玉树县(33.23°N，96.61°E)发生MS7.1地震，震源深度18 km． 玉树地震是继汶川地震后的又一次震灾严重的地震．
玉树地震后，中国地震局科学技术司立即组织了“玉树地震综合科学考查”. 其中，野外流动地震观测由中国地震局地球物理研究所与青海省地震局共同组成的流动地震观测科考队负责． 根据中国地震局的统一部署和要求，整个科考队于2010年5月2日出发，共建立了两个观测项目: 一个是在震源区架设由14个流动地震台站组成的密集台阵，分布在玉树地震断层两侧并覆盖震源区； 另一个是北起青海花石峡，南至西藏内乌齐架设的由16个(其中4个与密集台阵共用)流动地震台站组成的600 km地震大剖面(图 6)． 一体化宽频带地震仪共计12台，全部应用于第二个观测项目上，包括QSX(清水乡)、 MAD(玛多台)、 HHX(黑河乡)、 YNG(野牛沟)、 ZMT(扎马台)、 ZDT(扎朵镇)、 ZQT(珍秦台)、 NQT(囊谦)、 BZT(白扎乡)、 JQT(吉曲台)、 JSK(甲桑卡)和JDX(吉多乡)．
![]()
图 6 南北剖面台站(三角形)位置Figure 6. Deployment of stations (denoted by triangles) for the NS observation profile流动地震观测持续了半年，记录了大量的区域地震事件和远震事件. 其中3级以上地震有30余次，4级以上地震有10余次，5级以上地震有2次. 记录到的最大地震为MS5.9，发生于2010年5月29日10点29分，位于(33.3°N，96.3°E)，震源深度为10 km(波形见图 7)． 根据这些记录到的余震资料，可以对地震序列的空间分布进行精确定位，反演地震的震源过程，并根据地震学方法得到震源区的深部结构图像，获得该震源区的构造异常分布和区域深部构造环境．
![]()
图 7 JDX台站记录到的2010年5月29日青海省玉树MS5.9地震Figure 7. The 29 May 2010 MS5.9 earthquake in Yushu, Qinghai of China, recorded by the station JDX4. 结论
针对流动观测中传统的由于地震计与地震数据采集器单独存在而出现的运输不方便、 功耗大的问题，本文详细介绍了一款将地震信号提取、 数据采集、 记录和服务集为一体的低功耗宽频带数字地震仪． 该地震仪具有频带宽(60 s—80 Hz)、 动态范围高(>140 dB)、 功耗低(0.6 W左右)、 携带方便(整机系统在15 kg左右，包括供电系统、 GPS天线和包装箱)等特点． 另外，该仪器还具有多种采样率同时输出和低延时数据输出等其它同类仪器不具备的特点，进一步扩大了本仪器的应用范围． 该仪器在研制过程中进行了大量的实验室和台站测试． 仪器研制完成后，投入到了青海玉树MS7.1地震震后实际流动观测试验中，并取得了详尽的第一手观测资料． 通过在青海玉树地震震后流动观测试验中的具体应用，不管从仪器的运输和携带上，还是从仪器的稳定运行上都可以得出，该仪器是一款非常适合流动观测的设备． 与国外同类仪器相比，该仪器还具有很好的性价比． 下一步将从更低的功耗、 更轻的整机系统重量和更稳定的系统运行等3个方面着手，进一步提高该系统的各项性能指标．
-
![]()
![]()
![]()
图 4 井水位与水温同震响应组合关系多样性的水热动力学机制解释
虚线为井水温度的原梯度线； 实线为同震响应后井水位变化引起的水温梯度线变化(a) 井水位振荡时，井筒上半部水温上升，下半部水温下降； (b) 水温梯度值为正时，井水位上升，井水温上升； (c)水温梯度值为正时，井水位下降，井水温下降；(d) 水温梯度值为负时，井水位上升，井水温下降； (e)水温梯度值为负时，井水位下降， 井水温上升
Figure 4. Explanation of several combination relations of coseismic response of water level and water temperature by using water-heat dynamic mechanism where dashed line denotes original gradient of water temperature,solid line denotes the gradient of water temperature resulted from change in water level after coseismic response
(a) If water level oscillates,water temperature rises in upper half of a well and drops in lower half of the well； (b) If the gradient of water temperature in a well is positive,water level rises, water temperature also rises； (c) If the gradient of water temperature in a well is positive,water level drops,water temperature also drops； (d) If the gradient of water temperature in a well is negative,water level rises,water temperature drops; (e) If the gradient of water temperature in a well is negative,while water level drops,water temperature rises
![]()
图 5 汶川MS8.0地震的3个典型的水温同震响应及震后效应记录图
(a) 同震下降-震后恢复型； (b) 同震上升-震后高值型； (c) 同震下降-震后低值型
Figure 5. Three recordings of typical coseismic response to the Wenchuan MS8.0 earthquake and post-seismic change of water temperature in some wells
(a) Coseismic response drops,post-seismic effect is resumable (Mengyin well,Shandong Province);(b) Coseismic response rises,post-seismic effect keeps on a high value (Nanxi well,Sichuan Province); (c) Coseismic response drops,post-seismic effect keeps on a lower value (Zhouzhi well,Shanxi Province)
![]()
图 6 青海省德令哈干井中地温对汶川MS8.0地震的同震响应及震后效应曲线
Figure 6. Recordings of coseismic response to the Wenchuan MS8.0 earthquake and post-seismic change of earth temperature in a dry well,Delingha,Qinghai Province
![]()
图 7 四川西昌川03井3个不同深度上水温微动态(a)、 井孔水文地质特征(b)和井水温度梯度变化(c)图
Figure 7. Micro-behavior of water temperature at three depths (a),hydrogeological characteristic of well hole (b) and change in gradient of water temperature (c) in the well Chuan03,Xichang,Sichuan Province
表 1 水热动力学机制下井水温度变化特征
Table 1 Characteristics of water temperature variation in a well under water-heat dynamics mechanism
含水层内的变化 井-含水层间水流运动 井筒内水流
运动方向井水位
变化井水温变化 正梯度 负梯度 地下水补排
关系变化补给量增多或排泄量减小
排泄量增多或补给量减少含水层地下水→井水
井水→含水层地下水向上
向下上升
下降上升
下降下降
上升变形破坏与孔
隙压力变化被压缩使孔压增大
被拉张或破裂使孔压减小含水层地下水→井水
井水→含水层地下水向上
向下上升
下降上升
下降下降
上升
下载: 导出CSV
表 2 塔院井6个深度上水温潮汐基本特征比较
Table 2 Basic features of water temperature tide at six depths of the well Tayuan
观测深度/m 日潮差/℃ 观测日期 月相 公历 农历 48 0 2008-04-04—08 二月廿八至三月初三 朔 85 0 2008-08-12—16 七月十二至十六 望 130 0.0005 2008-01-11—15 十二月初四至初八 上弦 178 0.0007 2007-12-21—25 十一月十二至十六 望 184 0.0035 2008-11-12—16 十月十五至十九 望 187 0.0020 2008-11-26—30 十月廿九至十一月初三 朔
下载: 导出CSV
表 3 汶川MS8.0地震水温与水位同震响应的组合类型
Table 3 Combination types of coseismic response of water temperature and water level of Wenchuan MS8.0 earthquake
响应形态 井数 井数百分比 组合类型 水位 水温 A 振荡 下降 19 19.6% B 振荡 上升 10 10.3% C 上升 上升 23 23.7% D 上升 下降 17 17.5% E 下降 下降 18 18.6% F 下降 上升 10 10.3%
下载: 导出CSV
表 4 常见岩石的导热系数
Table 4 Heat conductivity of common rocks
岩浆岩导热系数/(×4.23 mW·m-1·K-1) 变质岩导热系数/(×4.23 mW·m-1·K-1) 花岗岩 闪长岩 安山岩 玄武岩 石英岩 大理岩 千枚岩 页岩 2.09—3.10 1.85—2.10 1.10 1.87 1.60—4.80 2.60—3.20 1.94—2.66 2.16 沉积岩导热系数/(×4.23 mW·m-1·K-1) 砂岩 砾岩 石灰岩 白云岩 凝灰岩 1.1—2.6 1.8 1.88 0.94—4.3 0.61—1.37
下载: 导出CSV
表 5 川03井3个不同深度水温的汶川地震同震响应与震后效应特征比较
Table 5 Characteristics of coseismic response to Wenchuan MS8.0 earthquake and post-seismic change of water temperature at three depths of the well Chuan03
传感器置深/m 同震响应 震后效应 说明 形态 幅度/(10-2 ℃) 形态 持续时间/d 395 振荡 1.96 缓降 >3.5 水温持续下降，水位持续上升 595 急升 2.36 缓降—转平 1.5 水温在高值上稳定 765 缓升 0.76 缓升 >3.5 水温持续上升
下载: 导出CSV
表 6 川03井3个不同深度水温的潮汐效应特征比较
Table 6 Tide response of water temperature at three depths of the well Chuan3
传感器置深/m 形态显示 日潮差/(10-2 ℃) 水位潮汐关系 395 很明显 2.64 相位一致，关系密切 595 较明显 0.83 相位基本一致，有关系 765 不明显 ？ 关系不明显
下载: 导出CSV
-
车用太, 刘成龙, 鱼金子. 2008. 井水温度微动态及其形成机制[J]. 地震, 28 (4): 20-28 Che Y T, Liu C L, Yu J Z. 2008. Micro-behavior of well temperature and its mechanism[J]. Earthquake, 28 (4): 20-28 (in Chinese)
陈大庆, 刘耀炜, 杨选辉, 刘永铭. 2007. 远场大震的水位、 水温同震响应及其机理研究[J]. 地震地质, 29 (1): 122-130. Chen D Q, Liu Y W, Yang X H, Liu Y M. 2007. Co-seismic water level, temperature responses of some wells to far-field strong earthquakes and their mechanisms[J]. Seismology and Geology, 29 (1): 122-130 (in Chinese)
程建, 王多义, 李得力, 庞河清, 王帅成, 李军. 2010. 汶川大地震"远端效应"[J]. 成都理工大学学报: 自然科学版, 37 (2): 155-159 Cheng J, Wang D Y, Li D L, Pang H Q, Wang S C, Li J. 2010. Research on the "endpoint effect": On the geothermal anomaly in Longquanyi, Chengdu, China[J]. Journal of Chengdu University of Technology: Science & Technology Edition, 37 (2): 155-159 (in Chinese)
付子忠. 1988. 地热动态观测与地热前兆[G]//地质构造与地壳应力(1). 北京: 地震出版社: 1-8. Fu Z Z. 1988. The observation and earth-thermal precursor[G]//Thesis on Geological Tectonic and Stress in Earth Crust (1). Beijing: Seismological Press: 1-8 (in Chinese).
高文学, 蒋凤亮, 高庆华, 张业成, 朱克文. 2000. 地球化学异常[M]. 北京: 地震出版社: 163-174. Gao W X, Jiang F L, Gao Q H, Zhang Y C, Zhu K W. 2000. The Anomalies of Geochemistry[M]. Beijing: Seismological Press: 163-174 (in Chinese).
李美, 康春丽, 李志雄, 荆风, 薛艳, 闫伟. 2010. 汶川MS8.0地震前地表潜热通量异常[J]. 地震, 30 (3): 64-71 Li M, Kang C L, Li Z X, Jing F, Xue Y, Yan W. 2010. Abnormal surface latent heat flux prior to the Wenchuan MS8.0 earthquake[J]. Earthquake, 30 (3): 64-71 (in Chinese)
刘耀炜. 2009. 动力加载作用与地下水物理动态过程研究[D]. 北京: 中国地质大学(北京): 68-73 Liu Y W. 2009. Dynamic Loading and Physical Dynamics Process of Groundwater[D]. Beijing: China University of Geosciences (Beijing): 68-73 (in Chinese)
马玉川. 2010. 井水温度潮汐效应及其应变响应能力研究[D]. 北京: 中国地震局地壳应力研究所: 26-28. Ma Y C. 2010. Tidally Induced Water Temperature Changes in Boreholes[D]. Beijing: Institute of Crustal Dynamics, China Earthquake Administration: 26-28 (in Chinese).
任雅琼, 陈顺云, 马瑾. 2012. 依兰-伊通断裂带地表温度变化分析[J]. 地震学报, 34 (5): 698-705 Ren Y Q, Chen S Y, Ma J. 2012. Variation of land surface temperature in Yilan-Yitong fault zone of northeastern China[J]. Acta Seismologica Sinica, 34 (5): 698-705 (in Chinese)
石耀霖, 曹建玲, 马丽, 尹宝军. 2007. 唐山井水温的同震变化及其物理解释[J]. 地震学报, 29 (3): 265-273 Shi Y L, Cao J L, Ma L, Yin B J. 2007. Tele-seismic coseismic well temperature changes and their interpretation[J]. Acta Seismologica Sinica, 29 (3): 265-273 (in Chinese)
史浙明, 王广才, 刘成龙, 梅建昌, 王金维, 方慧娜. 2013. 三峡井网地下水位的同震响应及其含水层参数关系[J]. 科学通报, 58 (25): 3080-3087. Shi Z M, Wang G C, Liu C L, Mei J C, Wang J W, Fang H N. 2013. Co-seismic response of groundwater level in the Three Gorges well network and its relationship to aquifer parameters[J]. Chinese Sci Bull, 58 (25): 3080-3087 (in Chinese).
王成喜, 邓斌, 朱利东. 2009. 四川省青川县东河口地震遗址公园发现温泉及天然气溢出[J]. 地质通报, 28 (7): 991-994. Wang C X, Deng B, Zhu L D. 2009. Discovery of hot spring and nature gas exposure spots in Donghekou, Qingchuan County, Sichuan Province, China[J]. Geological Bulletin of China, 28 (7): 991-994 (in Chinese)
杨竹转. 2011. 地震引起的井水位水温变化及其机理研究[D]. 北京: 中国地震局地质研究所: 70-73 Yang Z Z. 2011. Variations of Well Water Level and Temperature Caused by Earthquakes and Their Generating Mechanism[D]. Beijing: Institute of Geology, China Earthquake Administration: 70-73 (in Chinese)
尹宝军. 2010. 唐山井地下水动态特征的研究[D]. 北京: 中国地震局地球物理研究所: 145-150 Yin B J. 2010. Study on the Fluctuant Characteristics of Underground Water at Tangshan Well[D]. Beijing: Institute of Geophysics, China Earthquake Administration: 145-150 (in Chinese)
鱼金子, 车用太, 刘五洲. 1997. 井水温度微动态形成的水动力学机制研究[J]. 地震, 17 (4): 389-395 Yu J Z, Che Y T, Liu W Z. 1997. Preliminary study on hydrodynamic mechanism of microbehavior of water temperature in well[J]. Earthquake, 17 (4): 389-395 (in Chinese)
鱼金子, 车用太, 刘成龙, 李万明. 2012. 金沙江水网对日本9.0级地震的同震响应及其特征与机制[J]. 地震, 32 (1): 59-69 Yu J Z, Che Y T, Liu C L, Li W M. 2012. Co-seismic responses of water level and temperature in well of the Jinshajiang groundwater observation network to Japan MS9.0 earthquake: Characteristics and mechanism[J]. Earthquake, 32 (1): 59-69 (in Chinese)
张永仙, 石耀霖, 张国民. 1991. 流量与水温关系的模型研究及地震水温前兆机制的探讨[J]. 中国地震, 7 (3): 88-94 Zhang Y X, Shi Y L, Zhang G M. 1991. Model studies on relationship between underflow and its temperature and discussions about mechanism of groundwater temperature precursor[J]. Earthquake Research in China, 7 (3): 88-94 (in Chinese)
张元生, 郭晓, 钟美娇, 沈文荣, 李稳, 何斌. 2010. 汶川地震卫星热红外亮温异常[J]. 科学通报, 55 (10): 904-910 Zhang Y S, Guo X, Zhong M J, Shen W R, Li W, He B. 2010. Wenchuan earthquake: Brightness change from satellite infrared information[J]. Chinese Sci Bull, 55 (10): 904-910 (in Chinese)
-
期刊类型引用(37)
1. 孙强，王鹏，张正帅. 基于正态分布的双差波速比方法研究-以山东长岛震群为例. 震灾防御技术. 2024(01): 119-129 . 
百度学术
2. 帕肉合·卡依提，南芳芳，潘振生. 基于和达法分析2012年新源、和静交界M_S 6.6地震波速比变化特征. 地震地磁观测与研究. 2024(03): 23-30 . 百度学术
3. 邢康，贾漯昭，刘晓锋，吴亚锋，孙宗强. 2023年12月18日积石山M_S6.2地震前单台波速比变化特征. 地震地磁观测与研究. 2024(05): 28-33 . 百度学术
4. 陈丽娟，陈学忠，李艳娥. 2021-09-16泸县M_S6.0地震前波速比变化研究. 大地测量与地球动力学. 2023(05): 529-533 . 百度学术
5. 李文君，曾宪伟，任雪梅，许英才，朱鹏涛，李自芮，马小军. 吴忠——灵武地区M_S≥5.0地震前波速比变化特征分析. 地震地磁观测与研究. 2022(02): 1-9 . 百度学术
6. 陈丽娟，吴绍宇，董蕾，魏红梅. 2013-12-16巴东M_S5.1地震前后波速比异常特征. 大地测量与地球动力学. 2022(10): 1080-1084 . 百度学术
7. 张天宇，王光明，刘自凤，黑贺堂. 2022年云南宁蒗5.5级地震前后波速比变化特征. 大地测量与地球动力学. 2022(11): 1156-1160+1183 . 百度学术
8. 李冬梅，郑建常，张正帅，赵金花，崔华伟. 基于双差法研究长岛震群震源区介质波速比. 地震地磁观测与研究. 2022(S1): 129-132 . 百度学术
9. 王亚茹，刘晓丹，郭蕾，王妍，周依，王想. 河北涿鹿M_L4.8地震前后单台波速比特征. 华北地震科学. 2021(02): 87-95 . 百度学术
10. 李文君，曾宪伟，罗国富，李国斌，罗恒之，司学芸. 2016年1月21日青海门源M_S 6.4地震的前兆异常研究. 地震地磁观测与研究. 2021(05): 68-75 . 百度学术
11. 张小涛，苑争一，姚丽，田雷，薛艳，杨文，姜祥华，史海霞，解孟雨，邓世广. 2021年5月21日云南漾濞M_S6.4地震总结. 地震地磁观测与研究. 2021(06): 165-180 . 百度学术
12. 郑建常，李冬梅. 基于误差分布的震源区波速比反演及其应用:乳山震群源区介质性质变化研究. 地球物理学报. 2019(05): 1693-1703 . 百度学术
13. LIU Zifeng，YE Beng，CHEN Jia，ZHOU Qingyun，WANG Guangming，PENG Guanling，LI Yongli. Study of the Change in Wave Velocity Ratio before and after Two Strong Earthquakes Using Airgun Source Data. Earthquake Research in China. 2019(02): 305-319 . 必应学术
14. 信世民，郭学增，刘静，张蕾，李瑞卿. 唐山震区上地壳平均波速比分析. 地震地磁观测与研究. 2019(03): 15-20 . 百度学术
15. 李婷婷，刘利，胡光武，陈飞. 江苏地区波速比与泊松比分布特征. 地震. 2019(03): 149-157 . 百度学术
16. 李冬圣，王莉婵，王宁，蔡玲玲，贾炯，李小军. 2012年河北唐山4.8级地震前后波速比变化特征. 华北地震科学. 2019(03): 62-66 . 百度学术
17. 何应文，杨晓琳，王建芳，杨建文，杞应平. 滇西至滇西南地区近期中强地震前单台波速比变化特征研究. 大地测量与地球动力学. 2019(08): 798-803+868 . 百度学术
18. 刘春，邱玉荣，王莹，呼楠. 2013年甘肃岷县漳县6.6级地震前后波速比变化特征分析. 高原地震. 2019(S1): 1-4 . 百度学术
19. 邢康，贾漯昭，程楠，吴亚锋. 九寨沟M_S 7.0地震前单台波速比变化特征. 地震地磁观测与研究. 2018(03): 39-44 . 百度学术
20. 张慧，李志雄，李盛. 海南岛及近海波速比特征初步研究. 地震工程学报. 2018(03): 524-534 . 百度学术
21. 侯跃伟，杨云，郑建华. 南京及邻区数字地震波速比变化特征分析. 防灾科技学院学报. 2018(03): 17-22 . 百度学术
22. 刘自凤，张天继，付虹，李智蓉，王光明，赵小艳，余达远. 2018年云南通海5.0级地震前波速比变化特征. 地震研究. 2018(04): 494-502 . 百度学术
23. 方禹心，顾强强，荆涛，包秀敏. 基于沈阳台的辽宁省地区波速比变化研究. 防灾减灾学报. 2018(04): 59-62 . 百度学术
24. 王小娜，刘锦. 广东新丰江库区波速比变化特征研究. 防灾减灾学报. 2017(04): 14-20 . 百度学术
25. 黄金刚，薛艳，宋美琴，周龙泉，易桂喜，吕坚，张致伟，任枭，韩立波，郑建常. 基于震例的测震学异常辅助判定系统. 中国地震. 2017(01): 154-165 . 百度学术
26. 李艳娥，王林瑛，宋美卿，陈丽娟. 从波速比变化看汶川与芦山地震的孕震过程. 大地测量与地球动力学. 2016(11): 991-997 . 百度学术
27. 李丽，宋美琴，吴昊昱，李宏伟. 2010年河津M_S4.8地震前后单台波速比变化特征分析. 地震地磁观测与研究. 2016(01): 8-14 . 百度学术
28. 刘文邦. 青海东北部地区地震直达波波速比及视速度测定与研究. 地震地磁观测与研究. 2016(02): 21-25 . 百度学术
29. 刘文邦，万玉杰，李玮杰. 2016年青海门源6.4级地震前后波速比变化研究. 地震研究. 2016(S1): 49-54+133 . 百度学术
30. 赵敏，龙锋，宫悦. 锦屏一级水电站库区蓄水前后波速比变化研究. 四川地震. 2016(02): 22-26 . 百度学术
31. 王亚茹，王想，宫猛，郭蕾，刘晓丹. 河北及邻区平均波速比变化特征分析. 中国地震. 2016(04): 747-755 . 百度学术
32. 张洪艳，张广伟，王晓山，张晨侠，康建红. 吉林地区波速比分布特征及构造意义. 地震地质. 2015(03): 829-839 . 百度学术
33. 韩晓明，陈文凯，张文韬. 河套地震带波速比分区特征. 地震地磁观测与研究. 2015(04): 6-13 . 百度学术
34. Weng Zhaoqiang，Liang Xiangjun，Wu Haoyu，Liu Linfei，Li Li. Study on V_P/ V_S Variation Behavior in the Shanxi Region by the Wadati Single Station and Multi-earthquake Method. Earthquake Research in China. 2015(04): 514-526 . 必应学术
35. 翁钊强，梁向军，吴昊昱，刘林飞，李丽. 单台多震和达法研究山西地区的波速比变化特征. 中国地震. 2015(01): 89-100 . 百度学术
36. 张洪艳，苗春兰，卢燕红，康建红，张宇. 吉林省前郭5.8级震群前后波速比变化特征. 震灾防御技术. 2015(03): 605-614 . 百度学术
37. 刘文邦，王培玲，马玉虎. 青海东北部地区多台波速比研究. 地震研究. 2014(S1): 45-49 . 百度学术
其他类型引用(4)
计量
- 文章访问数: 581
- HTML全文浏览量: 271
- PDF下载量: 23
- 被引次数: 41
