SV 冲击和振动 1875 - 9203 1070 - 9622 Hindawi 10.1155 / 2021/4133402 4133402 研究文章 讨论交叉方法的理论基础应用于井下波速测试 https://orcid.org/0000 - 0002 - 4291 - 9141 越南盾 https://orcid.org/0000 - 0001 - 6906 - 8335 Zheng-hua https://orcid.org/0000 - 0002 - 9011 - 6992 必应 Yuan-dong Bi 开明 交通科学与工程学院 南京理工大学 南京 中国 njtech.edu.cn 2021年 15 6 2021年 2021年 30. 4 2021年 3 6 2021年 15 6 2021年 2021年 版权©2021清盾等。 这是一个开放的文章在知识共享归属许可下发布的,它允许无限制的使用,分布和繁殖在任何媒介,提供最初的工作是正确的引用。

在工程实践,理论依据交叉方法,用于获得横波到达时间在井下的波速测试方法表面向前和向后罢工,是p波的极性保持不变,而横波的极性转换的方向倒。然而,信号记录在测试的特点往往是发现与这个理论依据交叉冲突的方法,也即p波的极性转换的作用下表面向前和向后罢工。因此,进行了三维有限元数值模拟研究交叉的理论基础方法的有效性。结果表明,剪切和压缩波都观察到在180°水平之间的相位差信号痕迹,符合激发所产生的扭转方向的冲动。此外,分析数值模拟结果证明是可靠的解决方案;这表明交叉的理论基础应用于井下波速测试方法不当。与此同时,在数值模拟揭示了因素(倾斜激发、检波器偏转倾向,和背景噪音),可能会导致纵波不扭转下表面的极性向前和向后罢工。减少影响因素,我们提出一个方法为井下波速测试下表面罢工,到达时差之间基于源峰和响应峰,横波速度和数值仿真结果表明,用这种方法接近土壤的横波速度理论。

中国国家自然科学基金 U2039208 U1839202 江苏省研究生创新研究与实践项目
1。介绍

横波速度 V 年代 站点的土壤是一个关键的动态参数在地震工程中,主要反映网站动态特性测定,场地地震反应分析和网站分类( 1- - - - - - 5]。等原位地震测试的无创性方法(单或多通道光谱分析表面波)和侵入性方法(金属或井下方法和PS日志)现在通常用于获得 V 年代 小应变水平下土壤的层( 3, 6- - - - - - 10]。一般认为,通过简单的算术,侵入性方法可以获得更可靠 V 年代 在厚覆盖物网站相比,计算的无创方法 V 年代 使用一个反演过程和多个解决方案。在侵入性方法中,井下方法下表面向前和向后罢工其少等优点被广泛用于水井井间比方法,操作方便。

在井下方法,表面激发后,土壤中的身体波浪传播层和由接收器内的钻孔记录。基于信号的传播体波的旅行时间在土壤层可以计算,将最终用于 V 年代 概要文件的网站。因此,精确地决定了横波的传播时间。一般来说,估算横波的旅行时间从一个波形是困难的到达时间以来横波并不总是明确的。克服缺点的单一的波形,交叉的识别方法是常用的在实践中通过前后表面激发横波到达时间( 11, 12]。交叉方法的理论基础是p波的极性保持不变,而横波的极性改变时,引人注目的是倒的方向。然而,在原位测试, V 年代 获得直接的交叉方法显著的误差。一方面,横波的到达时间两组信号记录向前和向后的罢工在同一深度不显眼的由于环境噪声或罢工的弱能量 13, 14]。另一方面,发病纵波分量与反极性和小振幅出现在一些信号横波组件(一个清晰的 6, 15),从而导致误判的交叉S-waveform点。所有这些问题会导致显著的横波速度和实际之间的错误评价。可以消除环境噪声,提高信噪比或更少的干扰测试的选择。然而,逆转的一部分信号到来之前横波往往忽略由于其小的振幅,很少是调查及其特征。因此,它具有重要的工程意义可靠测定陆上或海上工程站点属性的特征信号记录在井下方法深入研究。与此同时,这项研究提供了一个理论依据的判断横波发病和更好的交叉方法适用于陆上或海上工程的横波速度测试网站。

在此,一个三维(3 d)建立了集中质量有限元模型来模拟井下的方法 V 年代 测试。的时域显式动态分析方法逐步集成,得到表面激励下的振动响应,用于记录的信号的特点,分析在不同观察深度和交叉方法的合理性。此外,数值模拟结果的可靠性和有效性的理论基础的交叉验证方法应用于井下方法的解析解。最后,根据数值模拟结果,提出了横波传播时间的一种计算方法,并能获得更准确的 V 年代 用更少的影响通过原位测试条件。这项工作有利于理解交叉方法的理论依据,准确地定横波的传播时间。

2。讨论交叉方法的理论基础

在井下方法中,探针由三个检波器的元素 X- - - - - - Y- - - - - - Z正交配置密封在一个圆柱形包放置在井下。年底水平的双方的板(见图 1(一)),一组荷载产生相反的方向。交叉方法,理论基础,只有横波的极性逆转,用于接收到信号的解释水平检波器获得横波的旅行时间。

井下的方法。(a)井下速度测试方法的原理与板材表面源。(b)信号跟踪记录在不同的井下测试深度的方法在北京一个网站。(c)大小16米的三维数值分析模型 × 16米 × 25米。

V 年代现场测试获得的交叉方法并不总是令人满意和很容易受到“选择时间。“图 1 (b)给井下方法在一个浅的部分测试结果在密云区表土网站,北京。通过交叉方法,绿点在图 1 (b)可以认为是第一个到达的横波的考虑下初始偏振部分。然而,蓝点会选为横波最终的到来,加上横波和纵波振幅的特征。在两种情况下,剪切波速计算的蓝点约为489.1 m / s,这是接近于经验 V 年代的价值比绿点计算的鹅卵石。它显示了一个情况,如果蓝色点被选为横波的到来,蓝点前的逆转的信号应该是纵波,冲突和交叉方法的理论基础。从这个,在井下信号接收方法的特点以及是否交叉方法的理论基础是理性的吗?更好的解释上述问题,建立了一个3 d集中质量有限元模型来模拟井下的剪切波速测试方法和解决的时域显式动态分析方法逐步集成。

3所示。井下的有限元数值模拟方法

模拟信号产生的前后表面在井下激发方法,建立三维分析模型。计算区域的大小16 m × 16米 × 25米被切断从地球均匀各向同性半无限空间中通过引入一个人工边界条件。上边界是免费的;底部边界和外侧边界人工边界。据坎帕内拉和斯图尔特( 16和石原 17),土剪切波速的菌株在井下方法测试微小应变水平,因此土壤介质是线性弹性,介质参数 V 年代= 200 m / s,密度 ρ = 1900公斤/米<年代up>3和泊松比 μ = 0.45。锤子罢工被水平模拟 X方向脉冲励磁模型表面的中心,和截止频率的激约80赫兹(见图 1 (c))。根据动态有限元计算精度,分析模型是由一个3 d离散eight-node六面体的固体有限元的尺寸0.2米 × 0.2米 × 0.2米。时间步决心0.0002年代基于稳定条件的时域显式动力学有限元方法 18]。Multitransmitting公式(MTF) [ 18- - - - - - 20.)是利用模型来模拟人工边界的弹性半空间内辐射阻尼的影响,为确保散射波的nonreflection当通过人工计算区域的边界。等参数单元的大小和时间步中设置的数值计算模型的选择根据稳定性判据数值积分的 18),单位规格满足以下方程: (1) Δ x 1 6 1 10 V 年代 f , 在哪里 Δ X 是单位的大小, V 年代 是模型的横波速度, f是输入的最大频率激励与工程意义。

时间一步∆ t还应符合下列条件:稳定 (2) Δ t 最小值 Δ x V , 在哪里 V 纵波的速度在3 d模型。

2给出了信号跟踪观察点不同深度水平 X向前和向后激下方向。所有信号的极性反转时相反的励磁的应用,如图 2(一)- 2(c))。第一个交叉点的时候交叉方法来判断信号痕迹,标记为 T 一个,被认为是横波的到达时间,但这并非如此。

信号在不同深度痕迹观测点的水平 X方向向前和向后下激发;(一)——(c)原始信号痕迹在不同选择深度;(d) - (f)放大部分信号的波形在8米,16米,24 m观察点。

众所周知的理论分析,在无限均匀各向同性介质,没有纯横波因为之间的内在耦合产生的纵波和横波点源( 21, 22]。根据介质的波速模型,准确到达时间的纵波和横波在观察点不同深度可以被标记为计算 T P T 年代分别在信号痕迹(如图 2)。相比之下, T P T 年代,它是发现 T 一个之间的是 T P T 年代, V 年代将高估了 T 一个。除此之外,它是信号放大痕迹的从前面部分8 m, 16米,和24 m观察点,如图 2(d) - 2(f),最初的纵波也显示了反极性向前和向后激。因此,判断使用交叉的横波到达时间的方法是不准确的。

观察图 2当观察点在较浅的深度,横波发生干扰的纵波早到,很难被识别。放大波形也不帮助(如图 2(d));的识别 T P可以还是错误的 T 年代很容易。观测点的深度增加,以来逐渐纵波与横波分离纵波速度大于 V 年代(见图 2(f))。虽然p波随深度的增加急剧衰减,其能源仍然横波发生影响,这样的判断 T 年代仍然是困难的。横波发生在这种情况下,只有近似会发现即使放大信号,并与不同尺度不同的值。此外,难以识别p波发作,因为相对较小幅度的纵波与横波相比,产生了错误的工程纵波没有极化特征的认知。

理论结果表明,极性的纵波、横波时不能保持相同的逆向激励下的极性逆转。然而,它可以观察到在原位测试纵波向前和向后激发下具有相同的极性,这可能是因为测试状态和环境干扰的影响。

4所示。分析验证理论依据交叉方法

基于精确和封闭解的均匀各向同性弹性半空间内由于应用程序的集中力切向平面的边界和不同时间的海维塞单位函数曹国伟( 23),借助拉普拉斯和汉克尔变换,方向的位移与载荷方向一致直接应用于切向力以下行为表面上,相应的泊松比= 1/4,可以表示为方程( 3)的半逆方法: (3) u = 0 , τ < 1 3 , u = F 4 π μ Z τ τ 2 1 / 3 τ 2 + 2 / 3 1 / 2 2 τ 2 + 1 / 3 2 4 τ τ 2 1 / 3 τ 2 + 2 / 3 1 / 2 , 1 / 3 < τ < 1 , u = F 4 π μ Z τ 2 + 1 2 τ τ 2 1 / 3 τ 2 + 2 / 3 1 / 2 2 τ 2 + 1 / 3 2 4 τ τ 2 1 / 3 τ 2 + 2 / 3 1 / 2 + 2 τ 2 τ 2 1 2 τ τ 2 2 / 3 1 / 2 2 τ 2 1 2 τ 2 1 2 4 τ τ 2 1 τ 2 2 / 3 1 / 2 , 1 < τ , 在哪里 (4) τ = V 年代 t Z , F 海维塞单位函数, Z 加载点之间的距离, μ 是剪切模, V 年代剪切波速度。

获得方程的数值解( 3),一个近似海维塞单位函数对应的海维塞单位函数方程( 3介绍了),见以下方程: (5) F t ¯ = 16 t ¯ 3 , 0.0 t ¯ < 0.25 , 1 48 t ¯ t ¯ 0.5 2 , 0.25 t ¯ < 0.5 , 1 , t ¯ 0.5 , 在哪里 t ¯ 是无量纲时间, t ¯ = 0.5 t / t 0 , t是时间因素, t0的上升时间单位亥维赛函数近似。当 t0= 0.05,近似海维塞单位函数图所示 3(一个)。不难证明( 5), F t ¯ 是连续可微的,它在该地区的一阶衍生品是连续的吗 t ¯ 0

分析方法的结果。(一)近似海维塞单位函数用于分析方法。(b)位移弹性半空间内2米和8米的直接在表面切向点源。泊松比是0.25。

通过改变极性 F t ¯ ,两组对应位移不同 Z2米和8米的距离直接加载点的计算是通过以下的方程 μ = 68 MPa和 V 年代= 200 m / s,如图 3 (b)。的到达时间 tnondimensionalized由 V 年代 t / Z 。这样,横坐标的值1.0对应于横波的到来,而0.577的值对应于纵波的到来,如图 3 (b)。可以观察到两个痕迹的第一个交叉点的到来是纵波,纵波是完全相反的极性 F t ¯ 是相反的。

尽管分析解决方案无法模拟井下的复杂几何方法,它有助于理解p波初始运动的特点,这与上述模拟结果是一致的。因此,当脉冲的方向逆转,P -和s波应该反向极性。因此,它可以达到的理论基础交叉方法应用于井下方法不当。

5。剪切波速在井下方法的计算方法

上述结果表明,p波的极性逆转时激发的方向是反向;因此,横波的到达时间不能准确、有效地识别根据交叉方法。除了交叉方法,另外两个横波旅行时间计算方法可以用于井下的方法,即。、峰法(PPM)和互相关法(CCM) [ 16, 24- - - - - - 26]。这两种方法可以获得可靠的 V 年代当横波波形的信号是明显的峰值。然而,相邻观测点的信号之间的相关性很差,当只有一个接收器接收信号在每个罢工。此外,由于兴奋倾向的影响,偏转,倾向峰值检波器的位置,和振幅的信号,PPM的旅行时间计算的准确性受到影响。

倾向的影响励磁、偏转和检波器的倾向,横波将三个检波器元素中发现的 X- - - - - - Y- - - - - - Z正交配置。因此,一个新的横波旅行时间计算方法。首先,均方根两个水平分量和垂直分量的信号在时域方程(收购 6)。指的是常见的最大峰值的方法,选择复合信号的波形的峰值非常横波运动的代表。通过计算时移的触发信号的峰值和峰值之间的复合信号或相邻观测点的复合信号的峰值,可以获得横波传播时间,如图所示 4: (6) 信号 混合 3 = 标记 l H 1 2 + 标记 l H 2 2 + 标记 l V 2 3

示意图的NPPM横波旅行时间计算。

旅行时间的横波传播源和观测点和相邻之间间隔1米观察点计算通过PPM, CCM和NPPM分别。根据速度方程,横波速度计算。图 5显示之间的相对误差 V 年代介质的模型和计算 V 年代 倾向于激发的条件下,检波器偏转和倾向。为 V 年代错误计算PPM和CCM,信号在两个水平方向(NH签署<年代ub>1和NH<年代ub>2)使用。

相对误差的 V 年代 PPM、CCM和NPPM在不同的测试深度:(a)激发20°倾角,检波器偏转20°,检波器10°倾角,使用触发信号;(b)激发20°倾角,检波器偏转20°,检波器10°倾角,基于邻接收信号;(c)激发20°倾角,检波器偏转50°,检波器10°倾角,使用触发信号;(d)激发20°倾角,检波器偏转50°,检波器10°倾角,基于邻接收信号;(e)激发倾斜45°,检波器偏转20°,检波器10°倾角,使用触发信号;(f)激发倾斜45°,检波器偏转20°,基于相邻检波器10°倾角,收到信号。

观察到,在各种组合条件下,相对误差 V 年代 NPPM小于5%时,深度超过3米,甚至不到2.5%的触发信号。至于结果通过PPM和CCRM,横波在9米的相对误差图 5逐渐增加,深度逐步下降,10%在5米的深度。此外,从数据可以看出 5(一个) 5 (c)PPM的计算精度和CCRM影响检波器的偏转和水平的选择信号。当励磁的倾角改变从20°、45°(如图 5(一个)- - - - - - 5 (f)),旅行时间计算PPM和CCRM显著的影响和相应的相对误差 V 年代 增加。结果表明,提出的峰间值的新方法(NPPM)可以减少倾斜激励的影响,检波器偏转,并倾向剪切波传播时间的计算,和确定 V 年代是真正的一个关闭。

6。结论

交叉的可靠性研究方法应用到井下的离岸和在岸网站波速测试方法,我们采用了一个3 d时域显式动力学有限元模型结合MTF调查网站的P - / S波速度提出了基于交叉方法。结论和结果如下:

本文中的结果表明,该理论基础、p波的极性相同,而横波与表面极性逆向前和向后罢工,是无效的

本文的数值模拟表明,该条件交叉方法的理论基础是由倾斜的因素激发、检波器偏转,倾向,和背景噪音

基于我们的数值模拟和分析解决方案,基于到达时间差的源峰和响应峰与横波速度的计算方法提出了

数据可用性

使用的数据来支持本研究的发现可以从相应的作者。

的利益冲突

作者宣称没有利益冲突有关的出版。

作者的贡献

周董,Zhenghua,杰苏导致的概念和设计研究。东清导致数据的分析和解释。东清写手稿的草案。周Zhenghua做了重要修改手稿的重要知识内容。郝Bing和李Yuan-Dong执行统计分析。周Zhenghua获得资金。周Zhenghua监督这项研究。

确认

这项工作得到了国家自然科学基金(U2039208和U1839202)和江苏省研究生创新研究与实践项目。

Kayabasi 一个。 Gokceoglu C。 液化势评价一个地区使用不同的技术(Tepebasi Eskişehir,土耳其) 工程地质 2018年 246年 139年 161年 10.1016 / j.enggeo.2018.09.029 2 - s2.0 - 85054293484 施耐德 j . A。 梅恩 p W。 一种音乐形式 g . J。 岩土网站描述大孟菲斯地区使用的针入度测试 工程地质 2001年 62年 1 - 3 169年 184年 10.1016 / s0013 - 7952 (01) 00060 - 6 2 - s2.0 - 0035477001 猎人 j . A。 Pullan s E。 伯恩斯 r。 井下地震测井高分辨率反射测量在疏松的表土 地球物理学 1998年 63年 4 1371年 1384年 10.1190/1.1444439 2 - s2.0 - 0032125357 猎人 j . A。 Benjumea B。 哈里斯 j·B。 表面和井下横波地震厚土现场调查的方法 土动力学和地震工程 2002年 22 9 - 12 931年 941年 10.1016 / s0267 - 7261 (02) 00117 - 3 2 - s2.0 - 0036817337 r D。 g·K。 横波速度结构下的上地幔台湾从数组中表面波的分析 《地球物理研究快报 2005年 32 7 10.1029/2004 gl021868 2 - s2.0 - 21844442650 Garofalo F。 Foti 年代。 产品的影响 F。 InterPACIFIC项目:比较的侵入性和非侵入性方法表征地震站点。第二部分:距离地波和钻孔方法之间进行对比 土动力学和地震工程 2016年 82年 1 241年 254年 10.1016 / j.soildyn.2015.12.009 2 - s2.0 - 84959199719 斯托科 k . H。 s . H。 森林 r D。 一些贡献的原位地球物理测量解决岩土工程问题 学报ISC-2: 2日国际岩土与地球物理会议网站描述 2004年9月 葡萄牙的波尔图街头 Viana塞卡 19 22 Foti 年代。 Parolai 年代。 Albarello D。 Picozzi M。 应用瞬态面波法的方法对地震遗址的表征 在地球物理调查 2011年 32 6 777年 825年 10.1007 / s10712 - 011 - 9134 - 2 2 - s2.0 - 80053915701 迪菲奥雷 V。 Cavuoto G。 Tarallo D。 Punzo M。 伊万格丽斯塔 l 多道面波分析和井下测试考古“腭山”地区(意大利罗马):评估和2 d的影响对剪切波速的影响 在地球物理调查 2016年 37 3 625年 642年 10.1007 / s10712 - 015 - 9350 - 2 2 - s2.0 - 84949495655 T.-F。 G.-X。 T.-Y。 W.-S。 F.-C。 Y.-M。 瑞利波非线性反演基于萤火虫算法 应用地球物理 2014年 11 2 167年 178年 10.1007 / s11770 - 014 - 0430 - 8 2 - s2.0 - 84904883474 ASTM D7400-4 2007 井下地震的标准测试方法测试:一年一度的ASTM标准 2007年 美国,西肯肖霍肯PA ASTM z Y。 岩土工程测试和监控手册 1994年 沈阳,中国 辽宁科学技术出版社出版 Baziw e . J。 数字滤波技术解释地震数据锥 岩土工程学报 1993年 119年 6 998年 1018年 10.1061 /(第3期)0733 - 9410 (1993)119:6 (998) 2 - s2.0 - 0027788359 x M。 y . X。 l Z。 b . J。 单孔法的优化方法基于相关性的剪切波速测量功能 岩石和土力学 2006年 27 1161年 1165年 在中国 Crice D。 井横波调查工程现场调查 2002年 美国纽约萨拉托加 Geostuff 坎帕内拉 r·G。 斯图尔特 w·P。 地震锥分析使用数字信号处理动态网站描述 加拿大岩土期刊 1992年 29日 3 477年 486年 10.1139 / t92 - 052 2 - s2.0 - 0027039162 石原 K。 土在地震土工技术行为 1996年 牛津大学,英国 牛津大学出版社 z . P。 介绍波动理论在工程 1996年 中国,北京 科学出版社 z . P。 h·L。 一个传输边界弹性波传播的数值模拟 土动力学和地震工程的国际期刊 1984年 3 4 174年 183年 10.1016 / 0261 - 7277 (84)90033 - 0 2 - s2.0 - 0021502791 Z.-P。 推断漫散射边界条件 波动 1996年 24 2 117年 138年 10.1016 / 0165 - 2125 (96)00010 - 8 2 - s2.0 - 0001028086 Sanchez-Salinero 我。 Roesset j . M。 斯托科 k . H。 分析研究体波的传播和衰减 1986年 美国奥斯汀,得克萨斯州 德州大学奥斯汀分校岩土工程中心 阿基 K。 理查兹 p·G。 定量地震学 1980年 美国弗吉尼亚州赫恩登 大学科学书 曹国伟 c c。 动态响应的弹性半空间表面切向载荷 应用力学学报 1960年 27 3 559年 567年 10.1115/1.3644041 2 - s2.0 - 85024258125 Teachavorasinskun 年代。 Lukkunaprasit P。 一个简单的相关性曼谷软粘土的剪切波速 岩土工程 2004年 54 5 323年 326年 10.1680 / geot.2004.54.5.323 2 - s2.0 - 3142726382 x M。 j·S。 x。 D。 在剪切速度相关函数的应用单孔测量方法 地震工程与工程动力学 2004年 24 60 65年 在中国 l J。 动态土壤结构与半空间嵌入在一个灵活的基础:封闭解析解对入射平面SH波 地震工程杂志 2021年 25 8 1565年 1589年 10.1080 / 13632469.2019.1586802 2 - s2.0 - 85063547944