CN105676281B - 利用瑞雷面波波速确定地层力学参数的方法 - Google Patents
利用瑞雷面波波速确定地层力学参数的方法 Download PDFInfo
- Publication number
- CN105676281B CN105676281B CN201610044666.6A CN201610044666A CN105676281B CN 105676281 B CN105676281 B CN 105676281B CN 201610044666 A CN201610044666 A CN 201610044666A CN 105676281 B CN105676281 B CN 105676281B
- Authority
- CN
- China
- Prior art keywords
- wave
- rayleigh
- wave velocity
- surface wave
- velocity
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 26
- 238000012360 testing method Methods 0.000 claims abstract description 26
- 239000002689 soil Substances 0.000 claims abstract description 24
- 238000005553 drilling Methods 0.000 claims abstract description 12
- 239000006185 dispersion Substances 0.000 claims description 15
- 230000035515 penetration Effects 0.000 claims description 11
- 238000010008 shearing Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 6
- 230000005284 excitation Effects 0.000 claims description 5
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000000926 separation method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000032258 transport Effects 0.000 claims description 3
- 238000007689 inspection Methods 0.000 claims 1
- 239000011435 rock Substances 0.000 abstract description 7
- 238000012795 verification Methods 0.000 abstract description 3
- 239000011800 void material Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000010219 correlation analysis Methods 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/284—Application of the shear wave component and/or several components of the seismic signal
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种利用瑞雷面波波速确定地层力学参数的方法。该方法充分利用了瑞雷面波勘测带来的信息,根据瑞雷面波波速大小,对地层的岩土工程力学参数进行反演,结合钻探资料、剪切波速及室内土工试验进行对比验证,最终达到采用瑞雷波波速确定地层的各种力学参数,为没有钻探、坑探资料的岩土工程勘测提供地层力学指标参考值和地层划分。本发明高效、经济、实用,减少了大量的钻探、坑探工作,尤其适用于钻机难以到达的地形地质条件,结果合理有效,大大提升了工作效率,降低了成本,为地层的岩土工程力学性质参数的确定提供了一个有效的解决方法。
Description
技术领域
本发明属于岩土工程勘测领域,具体涉及一种利用瑞雷面波波速反推地层各种岩土工程力学参数的方法。
背景技术
瑞雷面波在研究地球内部结构、近地表各种工程勘查问题中,尤其在近地表地层结构与属性特征获取方面,是一种非常有效的地震勘探方法。由于瑞雷面波与横波速度密切相关,因此,通过提取瑞雷面波的频散曲线来反演横波速度,进而获取地层结构、构造等信息以实现勘探目的。其实,瑞雷面波中隐含的信息是非常丰富的,而目前利用瑞雷面波进行勘探的方法提取的信息有限,体现的地质信息过于单一、不够直观,影响了瑞雷面波在地质勘测领域的应用。
发明内容
本发明需要解决的技术问题是提供一种利用瑞雷面波勘探,从而减少钻探、坑探工作量,扩展岩土工程勘测适用性,提高勘测速率,降低勘测成本的地层岩土工程力学参数确定方法。
为解决上述技术问题,本发明所采用的技术方案是:
一种利用瑞雷面波波速确定地层力学参数的方法,包括以下步骤:
a、设置检波器:首先,在地面上沿瑞雷面波的传播方向,以一定的道间距Δx设置M个检波器;
b、测量瑞雷面波波速VR:震源激发后,在相邻检波器之间接收到的瑞雷波的时间差为Δt或瑞雷波的相位差则测量范围MΔx内瑞雷面波波速为:或其中i是自然数,在同一测点测量出一系列频率fi的VRi值,得到一条VR-f曲线,即频散曲线,频散曲线也可转换为VR-λR曲线,λR为瑞雷面波波长,λR=VR/f;
c、求取瑞雷面波波速VR与剪切波波速VS相关性:通过频散曲线进行反演解释,可得到地下不同深度的瑞雷面波波速VR值;并得出VR与VS的相关式:VR=VS·(0.87+1.12μ)/(1+μ),其中μ为泊松比;
d、根据VR与VS相关式求取地基力学参数,并适当进行钻探、剪切波法测试、土工试验,对瑞雷波勘探结果进行验证。
其中地基的各种力学参数如下:
动剪切模量Gd=ρVS 2(KPa),其中ρ是密度,
泊松比μ=(VP 2-2VS 2)/2(VP 2-VS 2),其中VP是纵波波速
静止土压力系数K0=(lgVS-1.875)/0.398,
先期固结压力PC=100(lgVS-1.998)/0.51,单位KPa,
无侧限抗压强度Q=100(lgVS-2.127)/0.443,单位KPa,
内摩擦角单位°,
内摩擦力C=Q/2,单位KPa,
确定地基承载力fk=12·N,N是土层的标准贯入击数,N与地基土剪切波速VS之间存在着相关关系:A和B为经验系数,不同的地区和地质情况,取值不同。
进一步的,步骤d中适当进行钻探、剪切波法测试、土工试验,对瑞雷波勘探结果进行验证的过程如下:
1)在工区选择3个以上钻孔,进行标准贯入实验获得不同深度标准贯入击数,然后在对应钻孔位置分别进行瑞雷波法勘探,获取该勘探点不同深度瑞雷波速,依据公式计算不同深度标准贯入击数,验证是否吻合。
2)随后,将该钻孔的土样运回土工实验室,进行土工实验,测试土样内摩擦角内摩擦力C,承载力fk物理力学指标,再依据公式计算对应力学指标,验证是否吻合。
由于采用了上述技术方案,本发明取得的技术进步是:
本发明根据瑞雷面波波速的大小与介质的物理特性之间的相关关系,根据岩土工程勘测实例,找出地层物性和瑞雷面波波速之间的函数关系,结合钻探资料、剪切波速及室内土工试验进行对比验证,能准确的计算出地层的岩土工程力学性质参数,为没有钻探、坑探资料的岩土工程勘测提供地层力学指标参考值和地层划分。本发明高效、经济、实用,减少了大量的钻探、坑探工作,尤其适用于钻机难以到达的地形地质条件,结果合理有效,大大提升了工作效率,降低了成本,为地层的岩土工程力学性质参数的确定提供了一个有效的解决方法。
附图说明
图1标准贯入击数N等值线剖面图,图中●为标准贯入试验击数;N为瑞雷横波换算贯入击数;
图2是本发明瑞雷面波频散曲线图,图中H为等波速层层厚,Vs为瑞雷波波速;
图3是本发明实验工区钻孔内剪切波波速数据表;
图4是本发明中心测点面波频散曲线;
图5是本发明瑞雷波法计算剪切波波速成果表;
图6是本发明瑞雷面波深度—速度剖面图。
具体实施方式
下面结合实施例对本发明做进一步详细说明:
一种利用瑞雷面波勘探地层岩土工程力学参数的方法,具体实现包括如下步骤:
a、首先,在地面上沿瑞雷面波的传播方向,以一定的道间距Δx设置M个检波器;
b、震源激发后,在相邻检波器之间接收到的瑞雷波的时间差为Δt或相位差则相邻道Δx长度内瑞雷波的传播速度为:
VRi=Δx/Δti或
在满足空间采样定理的条件下,则测量范围MΔx内瑞雷面波波速为:
或在同一测点测量出一系列频率fi的VRi值,就可以得到一条VR-f曲线,即所谓的频散曲线,当然,频散曲线还可以转换为VR-λR曲线,其中λR为波长,λR=VR/f;
c、通过对频散曲线(VR-f曲线或VR-λR曲线)进行反演解释,即可得到地下不同深度的瑞雷面波波速VR值;
d、VR值的大小与介质的物理特性有关,这是其作为岩土的物理性质评价指标的根本原因;
e、VR与VS(剪切波、横波)具有很好的相关性,其相关式为:
VR=VS·(0.87+1.12μ)/(1+μ),
式中μ为泊松比,在场地随机选取几个位置,进行钻探,可测得不同深度的剪切波波速数值,与使用本式计算得到的剪切波速进行对比,若存在误差则须进行校正;
f、国内外的大量试验研究成果表明,地基土剪切波速VS、土层的标准贯入击数N之间存在着相关关系:
式中A和B为经验系数,不同的地区和地质情况,取值不同,根据e步骤中的钻探实测的标准贯入击数,对式中经验参数进行确定,以便通过瑞雷面波波速计算整个工区的标准贯入击数;
g、确定地基的各种力学参数:
动剪切模量Gd=ρVS 2(KPa),其中ρ是密度,
泊松比μ=(VP 2-2VS 2)/2(VP 2-VS 2),其中VP是纵波波速,
静止土压力系数K0=(lgVS-1.875)/0.398,
先期固结压力PC=100(lgVS-1.998)/0.51,单位KPa,
无侧限抗压强度Q=100(lgVS-2.127)/0.443,单位KPa,
内摩擦角单位°,
内摩擦力C=Q/2,单位KPa,
h、利用以上公式,即可根据瑞雷面波波速计算得到各种参数,进而得到各参数的剖面等值线图;
i、确定地基承载力fk,其与剪切波速VS之间服从指数相关关系:fk=12·N;
j、由于各地区地质情况的不同,经验参数选取则不尽相同,在工区适当进行钻探、剪切波法测试、土工试验,以确定新工区经验参数,保证瑞雷波勘探结果的正确性,同时也能辅证勘探效果。过程如下:
1)在工区选择3个以上钻孔,进行标准贯入实验获得不同深度标准贯入击数,然后在对应钻孔位置分别进行瑞雷波法勘探,获取该勘探点不同深度瑞雷波速,依据公式计算不同深度标准贯入击数,验证是否吻合。
2)随后,将该钻孔的土样运回土工实验室,进行土工实验,测试土样内摩擦角内摩擦力C,承载力fk等物理力学指标,再依据公式计算对应力学指标,验证是否吻合。
实施例1
以下是国电吴起周湾一期49.5MW风电场工程风机位F17的测试过程。陕北黄土地区,地质情况较为简单,但地形起伏高差较大,很多地点钻机难以达到。
首先按照勘探目的布设瑞雷面波勘探测线,完全包含探测区,这样能采集到足够的信息。然后激发震源,一般采用锤击震源激发瑞雷面波,多激发几次以确保检波器采集到的信息具备高信噪比。分析检波器采集到的信息,可以得到瑞雷面波的频散曲线(即VR-f曲线)。
根据瑞雷波波速与剪切波波速的关系,计算剪切波波速。因为剪切波波速与瑞雷面波波速的对应关系式存在经验系数。因此需要在场地随机选取一些位置进行钻探,以校正经验系数,并辅证瑞雷面波勘探效果。钻探测量得到的剪切波波速如图3所示,其对应的瑞雷面波频散曲线如图4所示,而依据瑞雷面波计算的剪切波波速如图5所示。对比可知,由瑞雷面波波速计算得到的剪切波波速与钻孔实测波速非常接近,验证了本方法的有效性。
图1是整条瑞雷面波测线测得的瑞雷面波波速和深度直接的关系图。可以根据各个测点的频散曲线,将地下介质分层,如图2所示。根据各层剪切波速度的不同,可以计算不同深度的标准贯入值,即这里的经验系数的求取则是根据钻探已知钻孔的标准贯入值与剪切波波速的数值做相关分析,即可得到。
另外的地层岩土工程力学参数,也都可以由剪切波速度进行求取。我们将瑞雷面波勘探得到的这些参数与实测参数进行列表对比,见表1。
表1瑞雷面波勘探的标准贯入值与实验室实测值对比
由上表对比可知:根据瑞雷面波勘探所得到的地层岩土工程力学参数与实际数值非常接近,因此瑞雷波勘探是一种快速、经济地确定地基各种力学参数的有效方法。
经过实测钻探资料的校正之后,修正后的经验参数可以应用于同一工区未经钻探测量的地点的地层岩土工程力学参数的求取,可得不同深度详细参数如表2所示。
表2风机位H4瑞雷面波资料解释
Claims (1)
1.一种利用瑞雷面波波速确定地层力学参数的方法,包括以下步骤:
a、设置检波器:首先,在地面上沿瑞雷面波的传播方向,以一定的道间距Δx设置M个检波器;
b、测量瑞雷面波波速VR:震源激发后,在相邻检波器之间接收到的瑞雷波的时间差为Δt或瑞雷波的相位差则测量范围MΔx内瑞雷面波波速为:或其中i是自然数,在同一测点测量出一系列频率fi的VRi值,得到一条VR-f曲线,即频散曲线,频散曲线也可转换为VR-λR曲线,λR为瑞雷面波波长,λR=VR/f;
c、求取瑞雷面波波速VR与剪切波波速VS相关性:通过频散曲线进行反演解释,可得到地下不同深度的瑞雷面波波速VR值;并得出VR与VS的相关式:VR=VSg(0.87+1.12μ)/(1+μ),其中μ为泊松比;
d、根据VR与VS相关式求取地基力学参数,并适当进行钻探、剪切波法测试、土工试验,对瑞雷波勘探结果进行验证;
地基的各种力学参数如下:
动剪切模量Gd=ρVS 2(KPa),其中ρ是密度,
泊松比μ=(VP 2-2VS 2)/2(VP 2-VS 2),其中VP是纵波波速,
静止土压力系数K0=(lgVS-1.875)/0.398,
先期固结压力PC=100(lgVS-1.998)/0.51,单位KPa,
无侧限抗压强度Q=100(lgVS-2.127)/0.443,单位KPa,
内摩擦角单位°,
内摩擦力C=Q/2,单位KPa,
确定地基承载力fk=12·N,N是土层的标准贯入击数,N与地基土剪切波速VS之间存在着相关关系:A和B为经验系数,不同的地区和地质情况,取值不同;
步骤d中适当进行钻探、剪切波法测试、土工试验,对瑞雷波勘探结果进行验证的过程如下:
1)在工区选择3个以上钻孔,进行标准贯入实验获得不同深度标准贯入击数,然后在对应钻孔位置分别进行瑞雷波法勘探,获取该勘探点不同深度瑞雷波速,依据公式计算不同深度标准贯入击数,验证是否吻合;
2)随后,将该钻孔的土样运回土工实验室,进行土工实验,测试土样内摩擦角内摩擦力C,承载力fk物理力学指标,再依据公式计算对应力学指标,验证是否吻合。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044666.6A CN105676281B (zh) | 2016-01-22 | 2016-01-22 | 利用瑞雷面波波速确定地层力学参数的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610044666.6A CN105676281B (zh) | 2016-01-22 | 2016-01-22 | 利用瑞雷面波波速确定地层力学参数的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105676281A CN105676281A (zh) | 2016-06-15 |
CN105676281B true CN105676281B (zh) | 2018-11-06 |
Family
ID=56302205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610044666.6A Active CN105676281B (zh) | 2016-01-22 | 2016-01-22 | 利用瑞雷面波波速确定地层力学参数的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105676281B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107092030B (zh) * | 2017-05-12 | 2019-04-12 | 武汉港湾工程质量检测有限公司 | 通过瑞雷波探测软黏土包的方法 |
CN107677197B (zh) * | 2017-09-19 | 2019-08-30 | 贵州电网有限责任公司 | 一种地下目标体的体积测量装置及测量方法 |
CN109101684B (zh) * | 2018-07-02 | 2022-11-22 | 武汉工程大学 | 一种规则层状半无限体中瑞利波基阶模态频散曲线的解析算法 |
CN109752262B (zh) * | 2019-01-18 | 2020-10-27 | 中国水利水电科学研究院 | 一种基于原位相对密度确定覆盖层土体动剪模量参数的方法 |
US11754744B2 (en) | 2020-03-04 | 2023-09-12 | Southern University Of Science And Technology | Surface wave prospecting method for jointly extracting Rayleigh wave frequency dispersion characteristics by seismoelectric field |
CN111722281B (zh) * | 2020-06-10 | 2022-07-12 | 中国铁路设计集团有限公司 | 基于面波勘探技术的基础沉降计算方法 |
CN112765774B (zh) * | 2020-12-25 | 2022-07-01 | 青岛黄海学院 | 一种铁路震源瑞雷面波力学模型及其数值模拟方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102943461B (zh) * | 2012-12-06 | 2014-08-13 | 天津市市政工程设计研究院 | 利用瑞雷面波评价地基强夯加固效果的量化分析方法 |
-
2016
- 2016-01-22 CN CN201610044666.6A patent/CN105676281B/zh active Active
Non-Patent Citations (1)
Title |
---|
瑞利波法评价地基土特性的应用;石中明;《中国优秀硕士学位论文全文数据库 工程科技II辑》;20031231;第5-55页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105676281A (zh) | 2016-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105676281B (zh) | 利用瑞雷面波波速确定地层力学参数的方法 | |
CN103258091B (zh) | 非常规油气藏水平井段三维岩体力学模型建立的方法及装置 | |
CN105277982B (zh) | 一种泥页岩总有机碳含量地震预测方法 | |
CN105134189A (zh) | 一种LogGMIR测井地应力识别优质储层的方法 | |
CN101231346A (zh) | 应用地震波速度预测煤、岩体物理力学参数方法 | |
CN103256046A (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
CN111273375B (zh) | 一种应用于缺水地区浅埋地下工程的地质勘察方法 | |
Stolte et al. | Towards consideration of epistemic uncertainty in shear-wave velocity measurements obtained via seismic cone penetration testing (SCPT) | |
CN110058323A (zh) | 一种致密砂岩地层脆性指数计算方法 | |
CN104374827A (zh) | 横观各向同性岩体原位动弹性模量的各向异性系数的测量方法 | |
CN102182437B (zh) | 煤矿井下钻孔水力压裂应力边界确定及消除方法 | |
CN108037184A (zh) | 隧道仰拱施工质量检测方法 | |
Ma et al. | Geophysical monitoring technology for CO2 sequestration | |
Yousefzadeh et al. | Microseismic 101: monitoring and evaluating hydraulic fracturing to improve the efficiency of oil and gas recovery from unconventional reservoirs | |
Pérez-Santisteban et al. | Comparison among SASW, ReMi and PS-logging techniques: Application to a railway embankment | |
Kurtuluş | Istanbul geotechnical downhole arrays | |
Lo et al. | MEASUREMENT OF UNKNOWN BRIDGE FOUNDATION DEPTH BY PARALLEL SEISMIC METHOD. | |
Ovando-Shelley et al. | Microtremor measurements to identify zones of potential fissuring in the basin of Mexico | |
Lidén | Ground vibrations due to vibratory sheet pile driving | |
Yordkayhun et al. | Joint analysis of shear wave velocity from SH-wave refraction and MASW techniques for SPT-N estimation. | |
CN104179164A (zh) | 用标贯击数N探测长江河口、砂嘴、砂岛地层砂性土相对密度Dr的方法 | |
CN104516013A (zh) | 一种利用垂直地震数据进行时频分析与油气检测的方法 | |
Kamalb et al. | Comparison of Shear Wave Velocity Derived from PS Logging and MASW–A Case Study of Mymensingh Pourashava, Bangladesh. | |
Theland et al. | Assessment of small-strain characteristics for vibration predictions in a Swedish clay deposit | |
Park et al. | Horizontal stresses at shallow depths in Seoul (Korea) gneissic region |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 050031 Changan District, Hebei, China, North Street, building No. 6, No. Applicant after: China Electric Power Construction Group Hebei electric survey and Design Research Institute Co., Ltd. Address before: 050031 Changan District, Hebei, China, North Street, building No. 6, No. Applicant before: Hebei Province Power Surveying Design & Research Institute |
|
GR01 | Patent grant | ||
GR01 | Patent grant |