![点击分享此内容可以赚币 分享](/master/images/share_but.png)
微动信号模拟及其在微动勘探中的应用.pdf
《微动信号模拟及其在微动勘探中的应用.pdf》由会员分享,可在线阅读,更多相关《微动信号模拟及其在微动勘探中的应用.pdf(8页珍藏版)》请在咨信网上搜索。
1、 第 47 卷第 4 期物 探 与 化 探Vol.47,No.4 2023 年 8 月GEOPHYSICAL&GEOCHEMICAL EXPLORATION Aug.,2023doi:10.11720/wtyht.2023.1303李传金,王强,渐翔,等.微动信号模拟及其在微动勘探中的应用J.物探与化探,2023,47(4):1040-1047.http:/doi.org/10.11720/wtyht.2023.1303Li C J,Wang Q,Jian X,et al.Microtremor signal simulation and its application in microtre
2、mor explorationJ.Geophysical and Geochemical Exploration,2023,47(4):1040-1047.http:/doi.org/10.11720/wtyht.2023.1303微动信号模拟及其在微动勘探中的应用李传金1,2,3,王强1,渐翔1,郑涛1,詹素华1,2,陈绍伟1(1.福建理工大学 土木工程学院,福建 福州 350118;2.福建理工大学 岩土工程研究所,福建 福州 350118;3.福建理工大学 地下工程福建省高校重点实验室,福建 福州 350118)摘 要:对于简正振型叠加法合成的微动信号,使用空间自相关法计算的频散曲线与理
3、论频散曲线对比,发现两者具有较好的一致性。模拟微动信号时涉及许多参数,这些参数对微动勘探开展具有重要意义。对这些参数进行了数值试验,试验结果表明:震源距离和台阵尺寸对结果有明显影响,在复杂地层情况下应考虑高阶面波,在安静的观测场地应用小尺寸台阵往往能获得更高频率的频散信息从而提高浅层分辨率。因此,在实际开展微动勘探时,应考虑震源分布、台阵尺寸和高阶面波等因素对勘探结果的影响。关键词:微动信号;简正振型叠加;微动勘探;频散曲线;空间自相关法中图分类号:P631 文献标识码:A 文章编号:1000-8918(2023)04-1040-08收稿日期:2022-06-16;修回日期:2023-02-0
4、1基金项目:中国科学院重点部署项目(ZDRW-ZS-2021-3-1);福建理工大学科研发展基金(GY-Z17157)第一作者:李传金(1982-),男,博士,副教授,主要从事微动勘探研究工作。Email:licj04 0 引言微动(microtremor)是指地球表面的微弱震动,来源于天气、海浪和潮汐变化等自然现象和车辆行驶、机械以及生产活动等人类活动,由体波和面波组成的一种复杂振动。微动勘探通过阵列观测,从微动信号中提取面波频散曲线,再反演频散曲线,最终获得地下介质的横波速度结构。近些年,微动勘探在工程物探中的应用越来越广泛1-5。在微动勘探理论方面,学者们对阵列形式6、高阶频散曲线提取7
5、-10等方面进行了深入研究。在微动勘探研究中,微动信号模拟在算法验证、阵列适用性评价、高阶频散曲线提取等方面起重要作用7,11-12。微动信号模拟方法有广义射线法8、简正振型叠加法11、有限差分法12等。每种方法有各自特点。由于在微动信号模拟过程中,需要计算大量的理论地震图,所需时间一般比较长,而简正振型叠加法具有计算速度快、时间短的优点13。因此,本文应用简正振型叠加法合成微动信号,并研究其在微动勘探中的应用。1 理论1.1 简正振型叠加法简正振型叠加法是一种理论地震图计算的常用方法。该方法通过地球自由振荡的各个振型叠加的方式表示位移场。Gilbert14首先提出用振型叠加的形式表示位移,给
6、出了激发函数表达式。Wood-house15提出了径向本征函数的计算方法,使得通过简正振型叠加计算理论地震图成为可能。Dahlen 等16在前人工作基础上,系统总结了简正振型理论、简正振型叠加法等内容。Yang 等17推导了更简洁实用的位移表达式。对于给定点源激发的位移场,可以通过简正振型线性叠加计算得出。假设在 x 处有一点源激发,按简正振型叠加法,在 x 处的位移表示为16-17g x,t;x()=kAkx,t()skx(),(1)式中:sk为第 k 个振型的本征函数;Ak为点源在 x处的第 k 个振型的激发函数;下标 k 封装了 4 个指标 q,n,l,m(),其中 q 只有两个取值(球
7、状振型或 4 期李传金等:微动信号模拟及其在微动勘探中的应用环状振型),n 为径向阶,l 为角阶,m 为方位角阶。本征函数 sk由球谐函数和径向本征函数构成,具体表达式参考 Gilbert 等14和 Yang 等17的文章。振型激发函数 Ak的表达式为16-17 Akx,t()=M:krs()()1-exp ikt()coskt2k,(2)式中:M 为震源矩张量;k为模态应变;k和 k分别为简正振型本征频率的实部和虚部。1.2 瑞利波频散曲线在层状介质中瑞利波具有频散特性,不同频率的瑞利波具有不同的相速度。瑞利波频散曲线的计算方 法 有 Haskell 法18、Knopoff 法19、矩 阵法
8、20、Abo-Zena 法21、广义反射透射系数法22和快速矢量传递算法23等。这些方法都通过求解瑞利波频散方程获得相速度,从而得到理论瑞利波频散曲线。求解瑞利波频散方程时,同一个频率可能有多个解,即同一频率的瑞利波可以以不同相速度传播。因此,理论频散曲线不只有 1 条相速度曲线,而是有多条相速度曲线,称这些相速度曲线为多阶瑞利波频散曲线。在多阶瑞利波频散曲线中,具有最低相速度值的频散曲线称为基阶瑞利波频散曲线,具有较高相速度值的频散曲线称为高阶瑞利波频散曲线11。根据相速度的大小,高阶瑞利波频散曲线进一步细分为一阶波、二阶波等。不同阶的频散曲线可以按能量占比复合而成的一条曲线,称为视频散曲线
9、。在实际观测中,采用空间自相关法从实测微动数据中提取的频散曲线不是多阶瑞利波频散曲线,而是视频散曲线7,11。通过介质响应函数可以计算理论上的视频散曲线,其关键是计算各阶瑞利波能量在总能量中的占比。各阶瑞利波能量在总能量中的占比可表示为24pif()p f()=cif()A2if()icif()A2if(),(3)式中:f 为频率;i 为第 i 阶瑞利波;ci(f)为相速度;pi(f)为各瑞利波的能量;pi(f)/p(f)为各瑞利波能量在总能量中的占比;Ai(f)为介质响应函数。介质响应函数 Ai(f)是地层介质对瑞利波传播的响应函数,可表示为Ai(f)=NRFR/k()i,(4)式中:i 表
10、示瑞利波的阶数,k 是波数,NR是与瑞利波的传递矩阵相关的表达式,FR是瑞利波的频散函数25。2 微动信号模拟本文采用的微动信号模拟方法是先通过简正振型叠加法计算理论地震图,再由多个理论地震图合成微动信号。微动信号模拟涉及地层模型、阵列布置、震源分布、面波阶数等多个参数。2.1 地层模型选取 3 层结构的速度递增型地层模型,参数如表 1 所示。表 1 地层模型参数Table 1 Parameters of geological model序号层厚度/m纵波速度/(ms-1)横波速度/(ms-1)密度/(gcm-3)1101 3002001.92501 8005002.232 5001 0002
11、.52.2 台阵布置按图 1 所示布置观测台阵,采用微动勘探中常用的双重圆形阵列,共布置 7 个台站,观测半径分别为 10 m 和 20 m。2.3 震源分布在距离台阵中心5001 000 m 的圆环范围内随机布置 1 000 个震源(如图 2 所示)。震源位于地表,方位和强度随机选取,在观测时长内随机出现。图 1 台阵布置Fig.1 Array configuration图 2 震源分布Fig.2 Sources distribution1401物 探 与 化 探47 卷 2.4 合成微动信号每个震源生成 1 个地震记录,所有的地震记录按时间先后顺序合并成 1 个微动信号。合成时,根据需要选
12、择不同的面波阶数。当地质模型存在高低速夹层时,高阶模式起重要作用,必须考虑高阶模式。与基阶瑞利波频散曲线对比时,只选择基阶模式,不考虑高阶模式。合成微动信号的时长为 30 min,采样率为 100 Hz。合成微动信号有 3 个分量,其中 2 个为水平分量,1 个为垂直分量。空间自相关法主要使用垂直分量。图 3 为合成微动信号垂直分量的波形。图 3 合成微动信号垂直分量的波形Fig.3 The vertical component wave of the synthetic microtremor signal3 合成微动信号在微动勘探中的应用微动勘探方法的关键技术是从微动信号中提取频散曲线。从
13、微动信号中提取频散曲线方法有频率波数法(FK 法)和空间自相关法(SPAC 法)。本文采用空间自相关法计算频散曲线。3.1 空间自相关系数将时长 30 min 的合成微动信号按 40.96 s 分段。对于每段数据,计算两个台站间的空间自相关函数,叠加平均所有台站距离相等的空间自相关函数得到空间自相关系数。计算完所有数据段后,计算空间自相关系数的均值和标准差。按图 1 所示的台站布置,计算两个台站之间的距离并进行分组排序,获得 5 个不同的台站距离,每 1 个距离对应 1 个空间自相关系数。图 4 显示了这 5 个不同台站距离的空间自相关系数。3.2 频散曲线采用扩展空间自相关法(ESPAC)2
14、6计算频散曲线。相速度选取范围为 101 000 m/s,间隔为 10 m/s。频率选取范围为0.5 20 Hz,间隔为0.1 Hz。图 4 空间自相关系数Fig.4 Spatial autocorrelation coefficients对于某一频率 f,由不同距离 d 和对应空间自相关系数组成一序列。用不同相速度生成零阶贝塞尔函数,并计算零阶贝塞尔函数和该序列之间残差。在选取的频谱范围内计算残差,对残差进行归一化。如图 5 所示,以横坐标为频率,纵坐标为相速度,按残差大小生成频谱(图 5)。在频谱图中用不同颜色表示残差的大小,用红色表示最小的残差。为突显残差较小值区域的特征,将频谱图中残差
15、大于 0.3的部分统一用蓝色表示。拾取残差最小值处的相速度并连成线,获得频散曲线。图 5 中不同频率范围的频谱特征有明显差异。频率在 412 Hz 时,频谱能量集中,容易拾取频散曲线。频率小于 2 Hz 时,频谱颜色几乎都为红色,结果有明显偏差,没有拾取频散曲线。频率在 24 Hz 范围,在较大相速度范围内,归一化后残差都小于 0.05,频谱颜色都接近红色。对于这一范围的频谱,虽然颜色差异较小,但依然可以通过残差的大小提取频散曲线。频率大于 12 Hz 时,频谱能量分散,对频散曲线的拾取造成一定的干扰。图 5 合成微动信号的频谱Fig.5 Frequency spectrogram of sy
16、nthetic microtremor signal2401 4 期李传金等:微动信号模拟及其在微动勘探中的应用3.3 和理论频散曲线对比使用表 1 所示的地层模型,采用 Haskell 传递矩阵法18计算理论频散曲线。在图 6 中,对比合成微动信号计算的频散曲线与理论频散曲线,发现两者在整体上具有较好的一致性。合成频散曲线和理论频散曲线具有相似特征,都可以划分出 3 个近似水平段,频率范围分别为 2.53 Hz、57.5 Hz 和大于10 Hz,对应的相速度范围分别为 800850 m/s、450500 m/s 和 190200 m/s。频散曲线的这种特征与地层模型的 3 层结构有关,和横波
17、速度由深到浅分别为 1 000 m/s、500 m/s 和 200 m/s 有一一对应关系。通过对比还发现,在不同频率段合成频散曲线和理论频散曲线略有差异。在 3.518 Hz 部分几乎完全重合;小于 3.5 Hz 的低频部分略有差异;大于 18 Hz 的高频部分合成频散曲线不连续,出现跳变,和理论频散曲线有较大差异。导致这种差异的原因是台阵尺寸有限,无法准确计算高低频部分的合成频散曲线。图 6 理论频散曲线与合成频散曲线对比Fig.6 Comparison of theoretical andsynthetic dispersion curves4 参数分析与讨论微动信号模拟涉及许多参数,这
18、些参数对微动勘探开展具有重要意义。本文对这些参数进行了数值试验,不同的参数设置得到不同的结果。试验结果分析如下。4.1 震源距离对不同距离的震源进行数值试验。图 2 中所示的震源距离为 5001 000 m。图 7 中补充了震源距离为 150 m 和 50200 m 的试验结果。由图 7 可知,在不同震源距离情况下,频散曲线可以拾取的最低 频率有明显差别。当震源距离为1 50 m时,拾a震源距离 150 m;b震源距离 50200 m;c频散曲线对比asource distance in the range of 1 to 50 m;bsource distance in the range
19、of 50 to 200 m;ccomparison of dispersion curves图 7 使用不同震源距离计算的频散曲线Fig.7 Dispersion curves with different sources distances取的最低频率为 6.2 Hz;当震源距离为 50200 m时,拾取的最低频率为 3.3 Hz;当震源距离为 5001 000 m 时,拾取的最低频率为 2.3 Hz。这说明,模拟微动信号时,如果震源距离设置太小,缺乏远距离震源,将无法准确获得低频部分的频散曲线。因此,模拟微动信号时应合理设置震源距离。4.2 台阵尺寸图 6 对比了合成频散曲线与理论频散曲
20、线,发现两者总体上具有较好的一致性,但在小于 3 Hz 的低频部分和大于 18 Hz 的高频部分略有差异。这种差异可能与台阵尺寸有关。对不同尺寸的阵列进行数值试验。台阵的类型还是图 1 所示双重圆形阵3401物 探 与 化 探47 卷 列,台阵尺寸发生变化,半径 R(指内圆的半径)由10 m 变成 1 m 和 100 m。图 8 为不同台阵尺寸得到的频散曲线。由图 8可知,表明不同尺寸的台阵得到的频散曲线是有所差别的。在图 8a 中,频率小于 10 Hz 时,在较大相速度范围频谱颜色几乎都为红色,说明较小尺寸(半径 R=1 m)的台阵难以获得低频的频散曲线。在图 8b 中,较大尺寸(半径 R=
21、100 m)的台阵可以准确估算低频的频散曲线,但在大于 8 Hz 的高频部分频散曲线有明显偏差。在图 8c 中,合成频散曲线和理论频散曲线对比结果表明:较大尺寸的台阵得到的频散曲线在低频部分较准确,在高频部分误差a台阵半径 R=1 m;b台阵半径 R=100 m;c频散曲线对比aR=1 m;bR=100 m;ccomparison of dispersion curves图 8 使用不同台阵尺寸计算的频散曲线Fig.8 Dispersion curves with different array sizes较大;较小尺寸的台阵得到的频散曲线在高频部分较准确,低频部分误差较大。因此,实际观测时应
22、根据勘探深度的要求,合理选择台阵尺寸。4.3 基阶和高阶图 6 中所示的合成频散曲线,采用表 1 所示的速度递增型地层模型,数值试验时只用面波的基阶模式,没用考虑高阶模式。除了速度递增型地层模型外,常见的地层模型还有含低速软夹层型和含高速硬夹层型7,11-12。表 2 显示了含高速硬夹层的地层模型,其中第 2 层的横波速度为 500 m/s,比相邻上下两层的速度都高。表 2 含高速硬夹层的地层模型参数Table 2 Parameters of geological model with high velocity interlayer序号层厚度/m纵波速度/(ms-1)横波速度/(ms-1)密
23、度/(gcm-3)11013002001.923018005002.232020003002.04250010002.5 图 9 为用表 2 的地质模型参数进行数值试验得到的结果。图 9a 为频谱,从中可以提取频散曲线。图 9b 为不同阶瑞利波的能量占比。从图 9b 中可知,在特定频率范围高阶瑞利波能量占主导地位。如在 2 Hz 附近和 58 Hz 范围,一阶瑞利波占主导地位,在 1420 Hz 部分,二阶瑞利波的能量占比随频率的增大而增大。图 9c 为多阶瑞利波频散曲线、理论视频散曲线和合成频散曲线对比图。由图 9c 可知,频散曲线呈现“之”字形。在 58 Hz 范围,合成频散曲线和一阶频散
24、曲线具有较好的一致性。在 2 Hz 附近和 45 Hz 范围,合成频散曲线主要由基阶和一阶的频散曲线按能量占比复合而成。在其他频率范围,合成频散曲线和基阶频散曲线有较好的一致性。在图9c 中,理论视频散曲线和合成频散曲线整体上具有较好的一致性,在高频部分略有差别。频率大于 16 Hz 时,理论视频散曲线呈现上升趋势,这与二阶瑞利波的能量占比不断提高有关。而合成频散曲线不存在明显的上升趋势,并且出现不连续现象。理论视频散曲线和合成频散曲线在高频处的差异可能与台阵尺寸有关,试验所用的台阵难以准确估算大于 18 Hz 部分的合成频散曲线,只有更小尺寸的台阵才能准确地估算高频部分的频散曲线。徐佩芬、杜
25、亚楠等7,11指出用空间自相关法计算的频散曲线是由不同阶的频散曲线按能量占比复合而成。当地下介质存在低速软弱夹层或高速硬夹4401 4 期李传金等:微动信号模拟及其在微动勘探中的应用a合成微动信号的频散曲线;b各阶瑞利波能量占比;c频散曲线对比adispersion curve calculated by synthetic microtremor signal;bpower fractions to different modes of Rayleigh waves;ccomparison of disper-sion curves图 9 含高速硬夹层的地层模型计算的频散曲线Fig.9 Di
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 微动 信号 模拟 及其 勘探 中的 应用
![提示](https://www.zixin.com.cn/images/bang_tan.gif)
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【自信****多点】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【自信****多点】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。