基于分布式声波传感阵列的地震动事件定位可行性研究 Feasibility study of seismic events positioning based on distributed acoustic sensing array 1 first-author 1 2 3 corresp 朱鸿鹄，教授，博士生导师，从事地质工程与岩石力学研究；E-mail：zhh@nju.edu.cn 朱鸿鹄，教授，博士生导师，从事地质工程与岩石力学研究；E-mail：zhh@nju.edu.cn zhh@nju.edu.cn 4 1 3 1 1 1 南京大学 地球科学与工程学院，江苏 南京，210023 南京大学 地球科学与工程学院，江苏 南京，210023 School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China School of Earth Sciences and Engineering, Nanjing University, Nanjing 210023, China 南京大学(苏州) 高新技术研究院，江苏 苏州，215123 南京大学(苏州) 高新技术研究院，江苏 苏州，215123 Nanjing University(Suzhou) High-Tech Institute, Suzhou 215123, China Nanjing University(Suzhou) High-Tech Institute, Suzhou 215123, China 南京大学 大地探测与感知研究院，江苏 南京，210023 南京大学 大地探测与感知研究院，江苏 南京，210023 Institute of Earth Exploration and Sensing, Nanjing University, Nanjing 210023, China Institute of Earth Exploration and Sensing, Nanjing University, Nanjing 210023, China 中国科学院 深海科学与工程研究所，海南 三亚，572000 中国科学院 深海科学与工程研究所，海南 三亚，572000 Institute of Deep-Sea Science and Engineering, Chinese Academy of Sciences, Sanya 572000, China Institute of Deep-Sea Science and Engineering, Chinese Academy of Sciences, Sanya 572000, China 赵俊 基于在云南省宾川县开展的现场试验，探究频率-波数分析技术应用于分布式光纤声波传感（DAS）数据处理的可行性，分析该技术在判断地震动事件所在方向的应用效果，在此基础上，探讨频率-波数分析法的偏差来源，并提出减小定位误差、提高分辨能力的措施。研究结果表明，DAS阵列的方向偏差为2.37°，相较于检波器阵列，其判断精确度更高。使用特殊设计的光缆提升应变的传递效率、使用锚固装置增强缆-土之间耦合效果以及改变DAS阵列的几何形态捕获多方向的震动信息等措施可以有效减少定位误差，提高分辨能力。 Based on the field test carried out in Binchuan county, Yunnan Province, the feasibility using F-K analysis method to process distributed acoustic sensing(DAS) data was studied, and the application effect in judging the direction of seismic events was analyzed. On this basis, the source of deviation in positioning using F-K analysis method was discussed, and three methods to reduce the positioning error and improve the positioning resolution capability were proposed. The results show that the direction deviation of the DAS array is 2.37°. Compared with the dense seismic array, DAS has higher judgment accuracy. The methods such as using specially designed optical cables to improve the efficiency of strain transfer, using anchoring devices to enhance the coupling effect between the cable and soil, and changing the geometry of the DAS array to capture vibration information from multiple directions can reduce the positioning errors and improve the resolution. 分布式声波传感(DAS) 地球物理探测 分布式光纤传感(DFOS) 地震动事件定位 distributed acoustic sensing(DAS) geophysical exploration distributed fiber optic sensing(DFOS) positioning of seismic events 密集地震台阵监测技术因其记录的资料丰富、信噪比高，在地震震源定位、近地表结构成像及地球深部结构探测等领域扮演着重要角色 [ 1 1 - 3 3 1-3 ] 。但由于密集地震台阵的部署和维护成本很高，因此，该技术通常只用于短期探测。近年来，蓬勃发展的分布式光纤声波传感(distributed acoustic sensing，DAS)技术为长时间、高密度的地震监测提供了新的技术方法。该技术以传感光缆为探测元件，通过解调沿光缆背向传播的瑞利散射光，实现光缆沿线动态应变(震动、声波)的长距离、分布式、实时定量监测 [ 4 4 ] 。与传统的节点式地震仪相比，传感光缆价格低廉、体积小、后期维护成本低，并且具有极强的环境适应性(如抗电磁干扰、耐腐蚀、防水等) [ 5 5 - 7 7 5-7 ] ，因此，可以部署于海底 [ 8 8 - 9 9 8-9 ] 、冰川 [ 10 10 - 11 11 10-11 ] 等传统地震仪无法大规模布设的恶劣场景中。此外，DAS还具有长距离、密集空间采样的优势，即它能将1根数十公里的光缆变成一连串以一定间距排列的“地震仪”，其空间分辨率可根据需要在数米至数十厘米之间调整。DAS自20世纪90年代问世以来广受关注，并在VSP采集、天然地震观测、近地表结构探测等领域得到成功应用 [ 12 12 - 15 15 12-15 ] 。若将地震仪替换为DAS阵列，则可以捕获更加丰富的波场记录。然而，DAS和传统地震仪在测量参数、响应频带、仪器敏感度、方向敏感性等方面存在着较大差异，一些传统的地震数据处理方法难以直接用于处理DAS数据。这成为DAS应用于地球物理探测的一个巨大障碍。 作为一种地震数据处理技术，频率-波数分析(F-K分析)通过对各地震仪记录的震动数据进行聚束(beam forming)处理，能准确、快速地获得震动波的背方位角和慢度，进而较准确地判断出地震动事件所处的方位。这种数据处理技术为DAS应用于地震动事件定位提供了可能。近年来，基于DAS的震动事件定位研究逐渐受到研究者的关注 [ 16 16 - 17 17 16-17 ] 。LIANG等 [ 16 16 ] 利用开展的小型模型试验，成功定位多个具有相同频率的窄带信号源，展示出DAS震动事件定位的能力，但由于其试验的尺度较小，工况较为简单，因此，定位的效果还需在复杂的工况下通过大型现场试验来进一步验证。本文基于在云南省宾川县开展的DAS现场试验，探究频率-波数分析技术应用于DAS阵列数据处理的可行性，通过对比DAS阵列和附近检波器阵列的聚束结果，验证DAS F-K分析技术的应用可行性，并提出减小聚束定位误差的几种思路。 1 基本原理及方法 1 1.1 　 DAS 原理 2 DAS是一种利用相干瑞利散射光的相位信息测量应变率(应变)的技术，它由解调仪和传感光纤2部分组成。解调仪不断地向连接的光纤发射激光脉冲，由于光纤内部存在不均匀的散射体，部分入射的脉冲光会发生散射。其中，沿光纤背向传播的瑞利散射光是DAS测量的基础。当光纤某个位置受外界扰动而产生应变时，该处光纤的折射率会发生改变，导致背向瑞利散射光的相位也会发生改变。解调仪通过分析光纤各个位置背向瑞利散射光的相位信息，得到应变率(应变)信息，其传感原理如 图1 图1 所示。 1.2 　 F-K 分析技术原理 2 F-K分析是一种地震阵列处理技术，常被用于估算地震波入射阵列的方向(背方位角) [ 18 18 - 19 19 18-19 ] 。由于各监测台站的空间位置不一( 图2 图2 (a))，因此，当地震波到达阵列时，各监测台站记录的地震波到达时间会不同，即波形记录之间存在时移( 图2 图2 (b))。这种时移与地震波在台站下方传播的慢度矢量(速度矢量的倒数)和台站间的位移矢量密切相关。 F-K分析技术通过设定不同的地震波入射方向及传播速度，利用已知的台站空间间距计算时移，然后，对各监测台站的波形记录进行时移叠加( 图 2 图 2 (c))，最后，基于叠加后的能量估算地震波入射阵列的方向和在阵列中的传播速度。 (a) 地震波以一定夹角入射监测阵列；(b) 各监测台站原始的波形记录；(c) 经过时移后各监测台站的波形记录 假定阵列几何中心的台站为参考台站，其记录的地震信号为 <math id="M1"><mi>s</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704284&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704293&type= 4.65666676 3.38666677 ，则阵列中第 i 个台站记录的地震信号 <math id="M2"><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704278&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704257&type= 6.43466663 4.48733330 可以表示为 <math display="block" id="M3"><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>s</mi><mo stretchy="false">(</mo><mi>t</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>·</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704288&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704272&type= 24.63800049 4.48733330 式中： <math id="M4"><msub><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>i</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704303&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704279&type= 2.28600001 3.72533321 为连接参考台阵到第 i 个台站的位移矢量； <math id="M5"><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704280&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704289&type= 3.38666677 3.72533321 为地震波的水平慢度矢量。 若已知地震波的入射方向及传播速度，则基于 式(1) 式(1) 可对阵列中的 N 个台站记录进行时移叠加，其叠加结果 <math id="M6"><mi>b</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704304&type= 4.99533367 3.38666677 可以表示为 <math display="block" id="M7"><mi>b</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>x</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">(</mo><mi>t</mi><mo>+</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>·</mo><msub><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704305&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704281&type= 33.61266708 8.72066593 由 式(2) 式(2) 可知，叠加结果与慢度矢量和台站间的位移矢量有关。当地震波的入射方向和传播速度越接近真实值时，其叠加的能量越强，聚束的效果越好。虽然地震波的入射方向和传播速度是未知量，但是可以假定入射角和传播速度的大致区间，即假定不同的地震波慢度矢量，通过计算叠加能量，搜索最佳的入射角和传播速度组合，该组合即为F-K分析的最终结果。其中，聚束能量 E 可用 式(3) 式(3) 来计算： <math display="block" id="M8"><mi>E</mi><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mo>+</mo><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo><mi mathvariant="normal">d</mi><mi>t</mi></mrow><mo>=</mo><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mo>+</mo><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mi>s</mi><mo stretchy="false">(</mo><mi>t</mi><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">u</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>-</mo><mi mathvariant="bold-italic">u</mi><mo stretchy="false">)</mo><mo>·</mo><msub><mrow><mi mathvariant="bold-italic">r</mi></mrow><mrow><mi>i</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mi mathvariant="normal">d</mi><mi>t</mi><mo>=</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704282&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704291&type= 71.62799835 10.49866676 <math display="block" id="M9"><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi></mrow></mfrac><msubsup><mo>∫</mo><mrow><mo>-</mo><mi mathvariant="normal">∞</mi></mrow><mrow><mo>+</mo><mi mathvariant="normal">∞</mi></mrow></msubsup><mrow><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mi>S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msup><mrow><mfenced open="|" close="|" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></mfrac><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msup><mrow><mi mathvariant="normal">e</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><mi>i</mi><mo>⋅</mo><mo stretchy="false">(</mo><mi mathvariant="bold-italic">k</mi><mo>-</mo><msub><mrow><mi mathvariant="bold-italic">k</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo stretchy="false">)</mo></mrow></msup></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi mathvariant="normal">d</mi><mi>ω</mi></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704292&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704306&type= 50.46133423 10.49866676 式中： <math id="M10"><mi>S</mi><mo stretchy="false">(</mo><mi>ω</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704307&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704283&type= 6.68866634 3.38666677 为地震信号 <math id="M11"><mi>s</mi><mo stretchy="false">(</mo><mi>t</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704284&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704293&type= 4.65666676 3.38666677 的傅里叶变换式； <math id="M12"><msub><mrow><mi mathvariant="bold-italic">k</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704294&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704308&type= 3.21733332 3.72533321 为真实地震波的波数矢量； k 为假定的波数矢量； u 0 为真实地震波的慢度矢量； u 为假定的慢度矢量。 1.3 　地震动事件定位方法 2 地震动事件定位主要基于DAS和F-K分析技术。对于F-K分析而言，阵列的孔径、台站的数目和台站的间距均会影响聚束的效果。一般而言，阵列孔径越大、台阵数目越多、台站间的距离越小，则聚束的分辨率越高、旁瓣数越少 [ 18 18 , 20 20 ] 。DAS技术具备的探测优势(如测量范围广、空间采样密度高等)能很好地提高聚束的分辨率，将DAS技术与F-K分析技术相结合能更准确地对地震动事件进行定位。 地震动事件定位的基本流程如 图3 图3 所示。首先，从DAS记录中拾取待分析的波形记录；其次，对波形记录进行预处理操作，过滤掉部分噪声记录；然后，进行F-K分析，根据聚束的能量可以大致判断出地震波的入射方向；最后，基于多个不同位置DAS阵列的聚束结果，对地震动事件进行定位。 2 现场试验 1 2.1 　场地概况 2 本试验场地位于云南省大理州宾川县东部郊外区域( 图4 图4 )，地理位置为100°35′E，26°00′N，海拔约1 600 m，属于高原季风气候区，干湿季明显。该区域位于滇中高原与横断山脉的交接地区，地质构造复杂。该区域总体地貌类型为山高谷深，地形起伏变化大，侵蚀切割强烈。地势呈现为东西高，中部低，南部地势略高于北部，山区面积占比约为85%。其中，山区出露的地层为碎屑岩，岩性为砂岩夹泥岩、砂岩、泥岩互层，平原区出露的地层为第四系残坡积或冲洪积层内的堆积土 [ 21 21 - 22 22 21-22 ] 。 2.2 　场地布设 2 本试验选择地势较为平坦的区域作为测试场地。沿着北东—南西向开挖长约600 m、宽约 15 cm、深约30 cm的沟槽。将通讯光缆埋置于沟槽底部，并用细砂混合泥土进行回填。通讯光缆的基本技术指标及结构如 表1 表1 所示。由前人的研究可知，光缆与周边岩土体的耦合效果会极大地影响记录数据的质量 [ 23 23 - 25 25 23-25 ] 。在本试验中，为了保证通讯光缆与周边土体之间具有较好的耦合效果，在回填时对土体进行了充分压实。为了验证DAS监测数据的准确性，在光缆埋设位置附近还布设了数台地震检波器。DAS解调仪与检波器的基本技术指标如 表2 表2 所示。 光纤类型 加强材料 纤芯数量/个 截面直径 φ /mm 单位质量/(kg·km -1 ) 基本结构 单模 钢绞线 1 5.0 3.8 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704334&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704316&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=44704324&type= 42.59999847 12.26666641 3.351 设备 采样频率/Hz 道间距/m 标距长度/m 监测分量 记录模式 授时模式 DAS解调仪 1 000 10 10 单分量 连续记录 GPS授时 检波器 500 — — 三分量 连续记录 GPS授时 2.3 　地震动事件模拟及数据采集 2 近年来，大容量气枪震源的发展引起了研究者的广泛关注 [ 26 26 - 27 27 26-27 ] 。与人工震源(炸药、落锤等)相比，气枪震源具有能量大、频率低、可重复性高、近/源场破坏小等优点，是一种理想的绿色环保震源 [ 28 28 - 29 29 28-29 ] 。一些科研人员利用气枪震源对地下介质的波速变化展开了研究 [ 30 30 - 31 31 30-31 ] 。本文利用气枪震源在测试场地进行了震动激发试验。气枪的激发位置及如 图4 图4 (c)所示，气枪的沉放深度为10 m，激发的压力稳定在7.5 MPa。为了捕捉完整的震动信号，DAS及检波器均设为连续记录模式。DAS上每个传感通道的空间位置已通过事先的敲击试验获知。 3 试验结果分析与讨论 1 3.1 　试验结果 2 图5 图5 所示为气枪震源激发前后DAS阵列及检波器阵列的记录。从 图5 图5 可见：DAS阵列和检波器阵列均能清晰地记录气枪震源激发的地震波信号(纵波到时、传播速度以及面波信号等)。与检波器记录的信号相比，单道DAS记录的信号信噪比较低。特别是在靠近气枪震源的部分传感通道，由于该处光缆与周边土体的耦合质量不佳，导致震动信号无法较清晰地被DAS捕获。 (a) DAS阵列；(b) 检波器阵列 为了进一步验证DAS记录数据的可靠性，对比了单道DAS与周边检波器的波形记录。因为两者记录的参数不同(应变率与速度)，为方便比较，将振幅信息进行了归一化处理。同时，为了展现信号的频率特征，对两者数据进行了时频分析，结果如 图6 图6 (a)和 图6 图6 (b)所示。从 图6 图6 (a)和 图6 图6 (b)可见：两者记录的频率信息较为一致，且纵波的到时也具有一致性，略微不同的是DAS记录的主频(1~80 Hz)要比检波器记录的主频(1~60 Hz)宽，且DAS对尾波的响应敏感性不如检波器。 图6 图6 (c)和 图6 图6 (d)所示为经过滤波处理(滤波频带1~20 Hz)后的波形记录。从 图6 图6 (c)和 图6 图6 (d)可见：DAS记录与检波器记录在相位上具有较好的一致性，两者的互相关系数为0.82，说明DAS能较准确地记录震动波的相位信息。 (a) DAS数据的时频分析图；(b) 检波器数据的时频分析；(c) 经过滤波处理后的DAS波形记录； (d) 经过滤波处理后的检波器波形记录 为了定位地震动事件，选取18:02:44至18:02:45的波形记录，基于 图3 图3 所示的处理流程分别对DAS及检波器阵列进行分析。F-K分析的参数设定如 表3 表3 所示。在给定慢度和背方位角的搜索范围下，不断地对慢度和背方位角进行网格搜索，基于 式(3) 式(3) 计算不同慢度和背方位角组合下的聚束能量。DAS阵列和检波器阵列聚束后的结果如 图7 图7 所示， 图7 图7 中，圆心位置表示传感阵列几何中心的位置，径向刻度为慢度，五角星表示真实的地震动事件所在方向。由 图7 图7 可见：DAS阵列聚束得到的背方位角为88°，真实的背方位角为90.37°，聚束误差为2.37°；检波器阵列聚束得到的背方位角为82°，真实背方位角为85.46°，聚束误差为3.46°。 滤波范围/Hz 背方位角搜索范围/( o ) 背方位角搜索步长/( o ) 慢度搜索范围/(s·km -1 ) 慢度搜索步长/(s·km -1 ) 1~20 0~360 2 0.1~2.0 0.01 (a) DAS阵列；(b) 检波器阵列 3.2 　分析与讨论 2 上述结果表明，无论是DAS阵列还是检波器阵列，均可以通过F-K分析技术判断地震动事件的大致方向。相较于检波器阵列，DAS阵列聚束的误差更小，这是因为DAS阵列能在空间上密集采样，捕获更丰富的地震波场信息。尽管DAS单道记录的信噪比要比检波器记录的低，但由于DAS具备密集空间采样的优势，因此，可以通过叠加多道邻近的DAS记录来提升数据的质量。 虽然DAS阵列采集了丰富的地震数据，但其聚束的结果与真实情况相比还存在一定偏差，其偏差可能来源于以下几个方面。 1) 缆-土耦合效果。光缆与周边岩土体之间的耦合性一直被视为影响DAS数据质量的关键因素之一。当光缆与周边介质之间耦合较差时，光缆无法准确地捕获介质震动信息，进而对数据分析造成影响 [ 32 32 - 34 34 32-34 ] 。本试验中，尽管事先已经通过压实回填土来增强光缆与周边土体的耦合效果，但由于光缆的铺设距离较长(约600 m)，部分光缆没有能很好地与周边岩土体耦合，其记录数据的信噪比明显低于其他的传感通道记录的数据的信噪比。这些低信噪比的数据对F-K分析产生了一定影响，使得聚束出现了一定误差；此外，本试验所选用的光缆为松套结构的通讯光缆，尽管松套结构的光缆能降低信号在传输过程中的光损，但也极大地影响了应变之间的传递，使得DAS无法准确地捕获周边岩土体的震动信息，这也是影响聚束效果的因素之一。 2) 阵列几何形状。聚束的分辨能力与阵列的几何形状密切相关。一般而言，阵列中需要存在相交的线性阵列才能较准确地约束地震波的入射方向(如“L”型)。当仅有1组线性的阵列时，阵列对平行于阵列方向传播的地震波约束较好，而对于其他方向入射的地震波，由于仅有单一方向的波形记录，阵列的聚束效果较差，旁瓣数较多。在本试验中，DAS阵列及检波器阵列均为近北东—南西向的布设，因此，在进行F-K分析时，对于垂直于阵列方向的约束较弱。从 图7 图7 可见，平行于阵列方向的分辨率要远远高于垂直于阵列方向的分辨率。阵列布设的单向性使得F-K分析的方向分辨率偏低，导致最终聚束的结果出现误差。 3) 方向敏感性。方向敏感性是DAS对不同方向入射地震波的一种响应特性。与传统的摆式地震仪不同，DAS只对沿光纤轴向方向传播的纵波以及与光纤成45°夹角的横波高度敏感，但对其他方向入射的地震波敏感性微弱 [ 35 35 ] 。因此，当用DAS采集震动数据时，对于部分方向入射的地震波，其波场信息可能会在DAS记录中缺失，这对聚束定位的结果会产生较大的影响。在本试验中，由于光缆各个位置的布设方向存在差异( 图4 图4 (c))，因此，各个传感通道对入射地震波的敏感性也会不一致。部分传感通道记录的震动数据会受到方向敏感性的干扰，当进行F-K分析时，这些受干扰的数据会影响聚束定位的分辨率，使得最终聚束的结果出现偏差。 基于上述分析，可采用以下3种方法来减小聚束的误差，提升定位的精度。 1) 在缆-土耦合效果方面，由前人的研究可知，光缆的种类、结构、锚固方式、土体围压、土体含水率等会影响缆-土之间的耦合效果 [ 23 23 , 32 32 , 36 36 - 37 37 36-37 ] 。因此，为了增强缆-土耦合效果，可以从以下几个方面入手：一是选用特殊设计的光缆进行探测。由于DAS的工作环境复杂，为了适应不同的工况、延长光缆的使用寿命，所用的一般是松套结构的光缆，但这种光缆会影响应变之间的传递效率。为了提升应变之间的传递效率，可以在光缆内部填充油膏，这样既能提升应变之间的传递效率，又能够防止水蒸气进入纤芯，损害光缆。二是可以在光缆上以一定间距增设锚固点(如布设十字锚、锚板等)，以扩大光缆与周边岩土体的接触面积。三是在回填时选用优质的回填料，使光缆与回填料之间充分接触，并且通过压实回填料进一步增强缆-土之间的耦合效果。 2) 在布设DAS阵列方面，应尽量避免布设线性光缆，应当在多个方向布设光缆来提升聚束的分辨能力。单点式的检波器布设简单、方便，阵列的布设形态可灵活多变(如小刀形、棋盘形、环形等)，而分布式的光缆布设较为复杂。布设光缆既要考虑阵列对地震波的响应能力，又要考虑施工的难度及布设的成本，其几何形态不能过于复杂。基于上述分析，列举了3种光缆的布设形态，具体如 图8 图8 所示。从 图8 图8 可以看到：这3种阵列的几何形状均较为简单，能降低光缆布设的难度，且各阵列至少在3个方向上具有传感通道，能对地震波的入射方向进行约束。 (a) 沙漏；(b) 直角形；(c) 雨伞形 3) 在光缆的方向敏感性方面，近年来不少研究者提出对光缆的内部结构进行设计(如将光缆螺旋缠绕)来满足多方向测量的需求 [ 38 38 - 39 39 38-39 ] 。虽然这种多分量传感技术还处于发展阶段，但其应用的可行性已得到初步验证。未来可以将这种螺旋缠绕光缆应用于DAS阵列中，提高聚束的分辨能力。 4 结论 1 1) DAS阵列以及检波器阵列均能很好地捕获地震波场信息。尽管单道DAS记录受到缆-土耦合等多方面因素的影响，其数据质量不如检波器记录，但由于DAS能在空间上密集采样，因此，可以叠加多道邻近的DAS记录来提升数据质量。此外，DAS阵列捕获的地震波场信息较检波器阵列而言更为丰富，其F-K分析产生的误差也更小。 2) 缆-土之间的耦合效果和阵列的几何布设形态是导致本试验聚束产生偏差的重要因素，光纤的方向敏感性也对聚束偏差有一定影响。 3) 使用特殊设计的光缆来提升应变的传递效率、使用锚固点来增强缆-土之间耦合效果、改变DAS阵列的几何形态来捕获多方向的震动信息等措施可减小聚束偏差。 4) 需要指出的是，单个DAS阵列捕获的波形记录只能大致判断出地震动事件所在的方向，如要定位地震动事件，还需结合其他DAS阵列的波形记录，通过进行F-K分析获得多组方向信息，进而判断出地震动事件的平面位置。 参考文献 朱子杰 , 王绪本 , 刘志强 , 等 . 基于密集台阵资料的背景噪声研究青藏高原东南缘地震各向异性 [J]. 地球物理学报 , 2021 , 64 ( 3 ): 823 - 837 . 基于密集台阵资料的背景噪声研究青藏高原东南缘地震各向异性 ZHU Zijie , WANG Xuben , LIU Zhiqiang , et al . Seismic anisotropy in the southeastern margin of the Tibetan Plateau revealed by ambient noise tomography based on high-density array [J]. Chinese Journal of Geophysics , 2021 , 64 ( 3 ): 823 - 837 . Seismic anisotropy in the southeastern margin of the Tibetan Plateau revealed by ambient noise tomography based on high-density array 李敏娟 , 沈旭章 , 张元生 , 等 . 基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析 [J]. 地球物理学报 , 2018 , 61 ( 5 ): 2075 - 2087 . 基于密集台阵的青藏高原东北缘地壳精细结构及九寨沟地震震源区结构特征分析 LI Minjuan , SHEN Xuzhang , ZHANG Yuansheng , et al . Fine crustal structures of northeast margin of the Tibetan Plateau and structural features of Jiuzhaigou earthquake focal area constrained by the data from a high-density seismic array [J]. Chinese Journal of Geophysics , 2018 , 61 ( 5 ): 2075 - 2087 . Fine crustal structures of northeast margin of the Tibetan Plateau and structural features of Jiuzhaigou earthquake focal area constrained by the data from a high-density seismic array 郝春月 , 郑重 . 地震台阵监测能力综述 [J]. 地震地磁观测与研究 , 2020 , 41 ( 6 ): 3 - 14 . 地震台阵监测能力综述 HAO Chunyue , ZHENG Zhong . An overview of seismic array monitoring capabilities [J]. Seismological and Geomagnetic Observation and Research , 2020 , 41 ( 6 ): 3 - 14 . An overview of seismic array monitoring capabilities 蔡海文 , 叶青 , 王照勇 , 等 . 分布式光纤声波传感技术研究进展 [J]. 应用科学学报 , 2018 , 36 ( 1 ): 41 - 58 . 分布式光纤声波传感技术研究进展 CAI Haiwen , YE Qing , WANG Zhaoyong , et al . Progress in research of distributed fiber acoustic sensing techniques [J]. Journal of Applied Sciences , 2018 , 36 ( 1 ): 41 - 58 . Progress in research of distributed fiber acoustic sensing techniques 崔允亮 , 王新 , 周联英 , 等 . 基于光纤监测的大直径变截面钢管复合桩承载性状 [J]. 中南大学学报(自然科学版) , 2020 , 51 ( 6 ): 1627 - 1636 . 基于光纤监测的大直径变截面钢管复合桩承载性状 CUI Yunliang , WANG Xin , ZHOU Lianying , et al . Bearing behavior of large-diameter variable cross-section steel composite pile based on optical fiber monitoring [J]. Journal of Central South University(Science and Technology) , 2020 , 51 ( 6 ): 1627 - 1636 . Bearing behavior of large-diameter variable cross-section steel composite pile based on optical fiber monitoring 朱鸿鹄 , 朱维申 , 殷建华 , 等 . 地下开挖模型试验的光纤监测 [J]. 中国矿业大学学报 , 2010 , 39 ( 6 ): 826 - 830 . 地下开挖模型试验的光纤监测 ZHU Honghu , ZHU Weishen , YIN Jianhua , et al . Fiber optic monitoring of an underground excavation model text [J]. Journal of China University of Mining & Technology , 2010 , 39 ( 6 ): 826 - 830 . Fiber optic monitoring of an underground excavation model text 蒋庆 , 王强 . 基于干涉型分布式光纤技术的水下管道泄漏检测实验分析 [J]. 中南大学学报(自然科学版) , 2012 , 43 ( 2 ): 576 - 580 . 基于干涉型分布式光纤技术的水下管道泄漏检测实验分析 JIANG Qing , WANG Qiang . Underwater gas pipeline leakage and location simulation in lab-scale based on optical fiber interferometer [J]. Journal of Central South University(Science and Technology) , 2012 , 43 ( 2 ): 576 - 580 . Underwater gas pipeline leakage and location simulation in lab-scale based on optical fiber interferometer LINDSEY N J , DAWE T C , AJO-FRANKLIN J B . Illuminating seafloor faults and ocean dynamics with dark fiber distributed acoustic sensing [J]. Science , 2019 , 366 ( 6469 ): 1103 - 1107 . Illuminating seafloor faults and ocean dynamics with dark fiber distributed acoustic sensing SLADEN A , RIVET D , AMPUERO J P , et al . Distributed sensing of earthquakes and ocean-solid Earth interactions on seafloor telecom cables [J]. Nature Communications , 2019 , 10 ( 1 ): 1 - 8 . Distributed sensing of earthquakes and ocean-solid Earth interactions on seafloor telecom cables WALTER F , GRÄFF D , LINDNER F , et al . Distributed acoustic sensing of microseismic sources and wave propagation in glaciated terrain [J]. Nature Communications , 2020 , 11 ( 1 ): 2436 . Distributed acoustic sensing of microseismic sources and wave propagation in glaciated terrain BOOTH A , CHRISTOFFERSEN P , SCHOONMAN C , et al . Distributed acoustic sensing of seismic properties in a borehole drilled on a fast‐flowing Greenlandic outlet glacier [J]. Geophysical Research Letters , 2020 , 47 ( 13 ): 1 - 10 . Distributed acoustic sensing of seismic properties in a borehole drilled on a fast‐flowing Greenlandic outlet glacier 李彦鹏 , 李飞 , 李建国 , 等 . DAS技术在井中地震勘探的应用 [J]. 石油物探 , 2020 , 59 ( 2 ): 242 - 249 . DAS技术在井中地震勘探的应用 LI Yanpeng , LI Fei , LI Jianguo , et al . Application of distributed acoustic sensing in borehole seismic exploration [J]. Geophysical Prospecting for Petroleum , 2020 , 59 ( 2 ): 242 - 249 . Application of distributed acoustic sensing in borehole seismic exploration LINDSEY N J , MARTIN E R , DREGER D S , et al . Fiber-optic network observations of earthquake wavefields [J]. Geophysical Research Letters , 2017 , 44 ( 23 ): 11792 - 11799 . Fiber-optic network observations of earthquake wavefields 林融冰 , 曾祥方 , 宋政宏 , 等 . 分布式光纤声波传感系统在近地表成像中的应用Ⅱ: 背景噪声成像 [J]. 地球物理学报 , 2020 , 63 ( 4 ): 1622 - 1629 . 分布式光纤声波传感系统在近地表成像中的应用Ⅱ: 背景噪声成像 LIN Rongbing , ZENG Xiangfang , SONG Zhenghong , et al . Distributed acoustic sensing for imaging shallow structure Ⅱ: ambient noise tomography [J]. Chinese Journal of Geophysics , 2020 , 63 ( 4 ): 1622 - 1629 . Distributed acoustic sensing for imaging shallow structure Ⅱ: ambient noise tomography 施斌 , 王宝善 , 张诚成 , 等 . 川西甲基卡锂矿3 211 m科学深钻多物理量分布式光纤观测 [J]. 科学通报 , 2022 , 67 ( 23 ): 2719 - 2726 . 川西甲基卡锂矿3 211 m科学深钻多物理量分布式光纤观测 SHI Bin , WANG Baoshan , ZHANG Chengcheng , et al . Multi-physical distributed fiber optic observation in a 3 211 m-deep scientific borehole at Jiajika lithium mine, western Sichuan [J]. Chinese Science Bulletin , 2022 , 67 ( 23 ): 2719 - 2726 . Multi-physical distributed fiber optic observation in a 3 211 m-deep scientific borehole at Jiajika lithium mine, western Sichuan LIANG Jiajing , WANG Zhaoyong , LU Bin , et al . Distributed acoustic sensing for 2D and 3D acoustic source localization [J]. Optics Letters , 2019 , 44 ( 7 ): 1690 - 1693 . Distributed acoustic sensing for 2D and 3D acoustic source localization PARKER T , SHATALIN S , FARHADIROUSHAN M . Distributed Acoustic Sensing: a new tool for seismic applications [J]. First Break , 2014 , 32 ( 2 ): 61 - 69 . Distributed Acoustic Sensing: a new tool for seismic applications ROST S , THOMAS C . Array seismology: methods and applications [J]. Reviews of Geophysics , 2002 , 40 ( 3 ): 2 - 1 . Array seismology: methods and applications 莫璧铭 , 李剑 , 孔慧华 , 等 . 基于FK能量聚束的传感器阵列优化布设方法 [J]. 计算机测量与控制 , 2019 , 27 ( 10 ): 285 - 288 . 基于FK能量聚束的传感器阵列优化布设方法 MO Biming , LI Jian , KONG Huihua , et al . Optimization method of shock sensor array layouts based on FK beam forming [J]. Computer Measurement & Control , 2019 , 27 ( 10 ): 285 - 288 . Optimization method of shock sensor array layouts based on FK beam forming 沈旭章 , 梅秀苹 , 张淑珍 , 等 . 兰州台阵响应函数及不同方位小地震事件FK分析结果 [J]. 西北地震学报 , 2010 , 32 ( 1 ): 59 - 64 . 兰州台阵响应函数及不同方位小地震事件FK分析结果 SHEN Xuzhang , MEI Xiuping , ZHANG Shuzhen , et al . The array response function of Lanzhou seismic array and results of FK analysis for small earthquakes in different azimuths [J]. Northwestern Seismological Journal , 2010 , 32 ( 1 ): 59 - 64 . The array response function of Lanzhou seismic array and results of FK analysis for small earthquakes in different azimuths 周鹏敏 , 薛鹏 , 范柱国 , 等 . 云南省宾川县地质灾害发育特征及形成条件 [J]. 安徽农业科学 , 2016 , 44 ( 1 ): 190 - 192 . 云南省宾川县地质灾害发育特征及形成条件 ZHOU Pengmin , XUE Peng , FAN Zhuguo , et al . Geological hazards characteristics and formation conditions in Binchuan County, Yunnan Province [J]. Journal of Anhui Agricultural Sciences , 2016 , 44 ( 1 ): 190 - 192 . Geological hazards characteristics and formation conditions in Binchuan County, Yunnan Province 范昆琨 , 范柱国 , 肖彬 , 等 . 宾川县地质灾害特征及防治措施 [J]. 中国水土保持 , 2011 ( 3 ): 17 - 19 . 宾川县地质灾害特征及防治措施 FAN Kunkun , FAN Zhuguo , XIAO Bin , et al . Characteristics and prevention measures of geological disasters in Binchuan County [J]. Soil and Water Conservation in China , 2011 ( 3 ): 17 - 19 . Characteristics and prevention measures of geological disasters in Binchuan County 施斌 , 张丹 , 朱鸿鹄 . 地质与岩土工程分布式光纤监测技术 [M]. 北京 : 科学出版社 , 2019 . 地质与岩土工程分布式光纤监测技术 SHI Bin , ZHANG Dan , ZHU Honghu . Distributed fiber optic sensing for Geoengineering monitoring [M]. Beijing : Science Press , 2019 . Distributed fiber optic sensing for Geoengineering monitoring 程刚 , 施斌 , 朱鸿鹄 , 等 . 光纤和砂土界面耦合性能的分布式感测试验研究 [J]. 高校地质学报 , 2019 , 25 ( 4 ): 487 - 494 . 光纤和砂土界面耦合性能的分布式感测试验研究 CHENG Gang , SHI Bin , ZHU Honghu , et al . Experimental study on coupling performance of fiber and sand interface based on distributed sensing [J]. Geological Journal of China Universities , 2019 , 25 ( 4 ): 487 - 494 . Experimental study on coupling performance of fiber and sand interface based on distributed sensing ZHU Honghu , HO A N , YIN Jianhua , et al . An optical fibre monitoring system for evaluating the performance of a soil nailed slope [J]. Smart Structures and Systems , 2012 , 9 ( 5 ): 393 - 410 . An optical fibre monitoring system for evaluating the performance of a soil nailed slope 李孝宾 , 宋政宏 , 杨军 , 等 . 利用分布式光纤声波传感器监测大容量气枪震源信号 [J]. 地震地质 , 2020 , 42 ( 5 ): 1255 - 1265 . 利用分布式光纤声波传感器监测大容量气枪震源信号 LI Xiaobin , SONG Zhenghong , YANG Jun , et al . Monitoring signal of airgun source with distributed acoustic sensing [J]. Seismology and Geology , 2020 , 42 ( 5 ): 1255 - 1265 . Monitoring signal of airgun source with distributed acoustic sensing 陈颙 , 王宝善 , 姚华建 . 大陆地壳结构的气枪震源探测及其应用 [J]. 中国科学: 地球科学 , 2017 , 47 ( 10 ): 1153 - 1165 . 大陆地壳结构的气枪震源探测及其应用 CHEN Yong , WANG Baoshan , YAO Huajian . Seismic airgun exploration of continental crust structures [J]. Scientia Sinica(Terrae) , 2017 , 47 ( 10 ): 1153 - 1165 . Seismic airgun exploration of continental crust structures 杨微 , 王宝善 , 王伟涛 , 等 . 陆地水体气枪震源重复性的影响因素及控制方法 [J]. 华南地震 , 2020 , 40 ( 4 ): 1 - 9 . 陆地水体气枪震源重复性的影响因素及控制方法 YANG Wei , WANG Baoshan , WANG Weitao , et al . Influence factors and control methods for the repeatability of air Gun source in land reservoirs [J]. South China Journal of Seismology , 2020 , 40 ( 4 ): 1 - 9 . Influence factors and control methods for the repeatability of air Gun source in land reservoirs 陈颙 , 张先康 , 丘学林 , 等 . 陆地人工激发地震波的一种新方法 [J]. 科学通报 , 2007 , 52 ( 11 ): 1317 - 1321 . 陆地人工激发地震波的一种新方法 CHEN Yong , ZHANG Xiankang , QIU Xuelin , et al . A new method of artificially exciting seismic waves on land [J]. Chinese Science Bulletin , 2007 , 52 ( 11 ): 1317 - 1321 . A new method of artificially exciting seismic waves on land 王宝善 , 曾祥方 , 宋政宏 , 等 . 利用城市通信光缆进行地震观测和地下结构探测 [J]. 科学通报 , 2021 , 66 ( 20 ): 2590 - 2595 . 利用城市通信光缆进行地震观测和地下结构探测 WANG Baoshan , ZENG Xiangfang , SONG Zhenghong , et al . Seismic observation and subsurface imaging using an urban telecommunication optic-fiber cable [J]. Chinese Science Bulletin , 2021 , 66 ( 20 ): 2590 - 2595 . Seismic observation and subsurface imaging using an urban telecommunication optic-fiber cable 魏芸芸 , 王海涛 , 苏金波 , 等 . 新疆2次中强地震前气枪震源反射波震相走时异常变化初步研究 [J]. 中国地震 , 2016 , 32 ( 2 ): 270 - 281 . 新疆2次中强地震前气枪震源反射波震相走时异常变化初步研究 WEI Yunyun , WANG Haitao , SU Jinbo , et al . The preliminary study on travel time abnormal variation of reflection wave phase of air-Gun in Xinjiang before two earthquakes with M S 5.0 [J]. Earthquake Research in China , 2016 , 32 ( 2 ): 270 - 281 . The preliminary study on travel time abnormal variation of reflection wave phase of air-Gun in Xinjiang before two earthquakes with MS5.0 ZHU Honghu , SHE Junkuan , ZHANG Chengcheng , et al . Experimental study on pullout performance of sensing optical fibers in compacted sand [J]. Measurement , 2015 , 73 : 284 - 294 . Experimental study on pullout performance of sensing optical fibers in compacted sand 张诚成 , 施斌 , 刘苏平 , 等 . 钻孔回填料与直埋式应变传感光缆耦合性研究 [J]. 岩土工程学报 , 2018 , 40 ( 11 ): 1959 - 1967 . 钻孔回填料与直埋式应变传感光缆耦合性研究 ZHANG Chengcheng , SHI Bin , LIU Suping , et al . Mechanical coupling between borehole backfill and fiber-optic strain-sensing cable [J]. Chinese Journal of Geotechnical Engineering , 2018 , 40 ( 11 ): 1959 - 1967 . Mechanical coupling between borehole backfill and fiber-optic strain-sensing cable ZHANG C C , ZHU H H , LIU S P , et al . Quantifying progressive failure of micro-anchored fiber optic cable-sand interface via high-resolution distributed strain sensing [J]. Canadian Geotechnical Journal , 2020 , 57 ( 6 ): 871 - 902 . Quantifying progressive failure of micro-anchored fiber optic cable-sand interface via high-resolution distributed strain sensing ZHAN Z . Distributed acoustic sensing turns fiber-optic cables into sensitive seismic antennas [J]. Seismological Research Letters , 2020 , 91 ( 1 ): 1 - 15 . Distributed acoustic sensing turns fiber-optic cables into sensitive seismic antennas 施斌 . 论大地感知系统与大地感知工程 [J]. 工程地质学报 , 2017 , 25 ( 3 ): 582 - 591 . 论大地感知系统与大地感知工程 SHI Bin . On the ground sensing system and ground sensing engineering [J]. Journal of Engineering Geology , 2017 , 25 ( 3 ): 582 - 591 . On the ground sensing system and ground sensing engineering 佘骏宽 , 朱鸿鹄 , 张诚成 , 等 . 传感光纤-砂土界面力学性质的试验研究 [J]. 工程地质学报 , 2014 , 22 ( 5 ): 855 - 860 . 传感光纤-砂土界面力学性质的试验研究 SHE Junkuan , ZHU Honghu , ZHANG Chengcheng , et al . Experiment study on mechanical properties of interface between sensing optical fiber and sand [J]. Journal of Engineering Geology , 2014 , 22 ( 5 ): 855 - 860 . Experiment study on mechanical properties of interface between sensing optical fiber and sand KUVSHINOV B N . Interaction of helically wound fibre-optic cables with plane seismic waves [J]. Geophysical Prospecting , 2016 ( 3 ): 671 - 688 . Interaction of helically wound fibre-optic cables with plane seismic waves NING I L C , SAVA P . High-resolution multicomponent distributed acoustic sensing [J]. Geophysical Prospecting , 2018 , 66 ( 6 ): 1111 - 1122 . High-resolution multicomponent distributed acoustic sensing 10.11817/j.issn.1672-7207.2023.05.015 P62 A 1672-7207（2023）05-1804-10 刘威, 朱鸿鹄, 张汉羽, 等. 基于分布式声波传感阵列的地震动事件定位可行性研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1804-1813. 刘威, 朱鸿鹄, 张汉羽, 等. 基于分布式声波传感阵列的地震动事件定位可行性研究[J]. 中南大学学报(自然科学版), 2023, 54(5): 1804-1813. LIU Wei, ZHU Honghu, ZHANG Hanyu, et al. Feasibility study of seismic events positioning based on distributed acoustic sensing array[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1804-1813. LIU Wei, ZHU Honghu, ZHANG Hanyu, et al. Feasibility study of seismic events positioning based on distributed acoustic sensing array[J]. Journal of Central South University(Science and Technology), 2023, 54(5): 1804-1813. 江苏省自然科学基金资助项目(BK20211086)；三亚崖州湾科技城科研项目(SKJC-2020-01-009)；国家自然科学基金资助项目(42077235)；中国科学院A类战略性先导科技专项子课题资助项目(XDA22040105) Project(BK20211086) supported by the Natural Science Foundation of Jiangsu Province; Project(SKJC-2020-01-009) supported by the Sanya Yazhou Bay Science and Technology City Scientific Research Project; Project(42077235) supported by the National Natural Science Foundation of China; Project(XDA22040105) supported by Chinese Academy of Sciences Class A Strategic Pilot Science and Technology Special Sub-Project Foundation 2022-05-14 2022-07-26 2023-05-26 2023-06-20 38909952 中南大学学报(自然科学版) 5 54