CN109541690B - 一种浅层介质结构面松散程度评价方法 - Google Patents
一种浅层介质结构面松散程度评价方法 Download PDFInfo
- Publication number
- CN109541690B CN109541690B CN201811458091.8A CN201811458091A CN109541690B CN 109541690 B CN109541690 B CN 109541690B CN 201811458091 A CN201811458091 A CN 201811458091A CN 109541690 B CN109541690 B CN 109541690B
- Authority
- CN
- China
- Prior art keywords
- medium
- velocity
- measuring point
- speed
- structural surface
- 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 38
- 238000012545 processing Methods 0.000 claims abstract description 5
- 230000005284 excitation Effects 0.000 claims description 4
- 238000005259 measurement Methods 0.000 claims description 4
- 230000001419 dependent effect Effects 0.000 claims description 2
- 238000003672 processing method Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 abstract description 15
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000009659 non-destructive testing Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 238000001028 reflection method Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
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/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6224—Density
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
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)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
Abstract
Description
技术领域
本发明涉及地质检测方法技术领域,具体地指一种浅层介质结构面松散程度评价方法。
背景技术
结构面检测(structure surface detection),工程建设中,通常出现大量固定构建的安装,随着工程项目的运营使用,结构会出现不同的变形、脱落等病害,所以需要定期对结构面进行病害检测。结构面检测的内容主要有脱空检测、变形检测等,采用手段主要采用地质雷达法进行检测。
地震反射波法(Seismic Reflection Method),利用地震反射波进行人工地震勘探的方法,是资源勘探的首要手段,勘探结果能较准确地确定界面的深度和形态,圈定局部构造,判断地层岩性。主要采用单道连续剖面法和多道连续剖面法两种,而为了提高信噪比,在多道连续剖面法中还广泛采用共深点反射技术。获得采集数据后,通过特殊处理流程获得地层反射波剖面,综合相关资料进行解释最终获得地层属性特征。
结构面检测技术中地质松散度是反映地质结构情况的重要参数,当前针对松散度检测主要是依靠地质雷达法进行检测,如专利号为“CN10564188A”的名为“基于探地雷达成套设备的路面结构内部质量状况评价方法”的中国发明专利,该专利介绍了一种通过地质雷达法探测路面结构内部对其进行状况评价的方法,该方法中针对地质松散度的检测是通过雷达反射波振的大小、波形散乱现象程度来进行判断,这种方法比较简单,但是实际上这种方法需要采用专利的设备,而且这种方法只能大致判断疏松程度,并不能对其量化,缺少评价参数标准,难以进行大范围的推广。
发明内容
本发明的目的就是要解决上述背景技术中现有松散度测量方法存在缺少量化手段、无法精确评判的问题,提供一种浅层介质结构面松散程度评价方法。
本发明的技术方案为:一种浅层介质结构面松散程度评价方法,其特征在于:包括以下步骤:1)、沿检测或勘察的介质表面布置测点;
2)、逐点激发采集每个测点处的直达波和反射波信号;
3)、对信号进行数据处理,获得测点处的直达波速度和反射波速度;
4)、按照如下公式计算该测点处介质结构面速度系数;
其中:Cv——结构面速度系数;
Vf——反射波速度;
Vz——直达波速度;
当结构面速度系数Cv小于a时,认为该测点处的介质为松散结构,当结构面速度系数Cv大于或等于a时,认为该测点处的介质为密实结构。
进一步的所述的步骤1)中,按照地震反射波多次覆盖原理,在待检测或待勘察的介质表面沿线性布置测点。
进一步的所述的步骤2)中,每个测点在采集数据时,每个测点上激发一次机械波信号,每个测点上设置有多个传感器接受激发的信号,逐个测点依次进行,直至完成所有测点激发。
进一步的所述的步骤3)中,数据处理的方法,通过带通滤波器对获取数据进行处理,去除数据中的零飘和野值。
进一步的所述的步骤3)中,获得测点处的直达波速度的方法为,拾取直达波起跳时间(直达波到达时间),对拾取的直达波时间散点(每一个传感器记录的振动波形记录,拾取直达波到达时间。这个到达时间,在以距离为横轴,时间为纵轴的坐标系里,就是一系列的近似直线的离散点),采用直线方程进行拟合,以激发点到接收点的距离即炮检距为自变量,直达波起跳时间为因变量,拟合的直线斜率即为直达波速度。
进一步的所述的步骤3)中,获得测点处的反射波速度的方法为,对已知深度的浅层介质界面,根据直达波速度设定速度范围,在该范围选取速度值计算该速度值下的反射波时距曲线(时距曲线(time distance curve)是指地震波到达地面各检波点的时间与震源点到各检波点的距离之间的关系曲线),计算时距曲线与炮集记录(一点激发,多道传感器接收的振动记录)的最佳拟合即为反射波速度。
进一步的所述的速度范围的设定方法为:根据直达波速度Vz划定速度范围为[1/2Vz,2/3Vz]。
进一步的反射波时距曲线与炮集记录的最佳拟合方法为:以时距曲线为基础,设定时间窗口,将该时间窗口内的炮集记录振幅值积分求和,最大振幅能量对应的速度即为最佳拟合反射速度。
进一步的所述的时间窗口为15~30us。
进一步的所述的a的取值为0.5≤a≤0.9。
本发明的评价方法在无需对现有设备进行改变的情况下,通过对地震反射波法获取的数据进行特定流程数据处理,获得浅层地质介质已知界面的松散程度分布情况,对工程应用有直接指导意义,填补了已知介质结构面松散程度无损检测该领域的空白。
附图说明
图1:高铁无碴轨道板超声波检测测线布置图;
图2:超声波共炮点道集原始记录;
图3:超声波共炮点道集预处理记录;
图4:超声波共炮点道集直达反射拾取记录。
具体实施方式
下面结合附图和具体实施例对本发明作进一步的详细说明。
对某高铁轨道板结构离缝进行检测,使用超声波反射法进行检测,沿着轨道横向方向间隔布置多条测线,每条测线沿轨道纵向布置,每条测线包括多个测点,如图1所示,为本实施例的测点布置结构。每条测线均匀布置测点,相邻测点之间的间距为0.2m,每个测点都进行数据采集测量时,将带有超声换能器的装置放置在测点上,激发超声换能器,超声波的频率为30KHz,装置上设置有11个接受传感器,相邻传感器之间间隔0.03m,每个传感器都接受到激发的超声波信号,然后收集处理。
每个测点采集得到一个共炮点道集,采集获得原始数据如图2所示。对原始记录进行带通滤波处理,去零飘处理,去野值处理,得到如图3所示的记录。对预处理记录拾取波形起跳点得到初至曲线,如图4中的曲线1,通过线性拟合得到该曲线的斜率为0.000336,求倒数即可得到直达波速度Vz为2980m/s。
以Vz为中心,划定速度范围[Vz*1/2,Vz*3/2](即[1490,4470]),以dv=Vz/100=29.8m/s为间隔,计算Vi=Vz*1/2+dv*i=1490+29.8*i(i为该速度范围内按固定间隔取值序号,0≤i<100,Vi为i对应的速度值)。根据已知结构面的层厚度D,计算速度Vi对应的时距曲线Si,以Si为中心,划定长度为20us的时间窗口,如图4中曲线2和曲线3所围成的窗口范围。
对该窗口范围内炮集记录振幅绝对值进行积分求和,得到该速度Vi所对应的能量Ei。搜索所有能量值最大的Emax所对应的速度Vi,即为该层结构对应的最佳反射速度Vf=2920.4m/s。
计算k=Vf/Vz=0.98,并进一步计算结构面速度系数Cv=1-|k-1|=0.02。该测点处的结构面速度系数小于0.9,本实施例中a根据经验取值0.9,则证明该测点处的介质为松散状态,通过观察轨道板上的离缝分布情况,该测点处的轨道板出现离缝情况在,证明因为轨道板下方介质为松散状态,轨道板由于下方受力不均匀,松散介质难以提供足够的支撑作用力,因此出现离缝,证明本实施例的测量评价方法是准确的。
对轨道板检测范围内,所有的测点计算结构面速度系数Cv,当结构面速度系数Cv小于0.9时,认为该测点处的介质为松散结构,当结构面速度系数Cv大于或等于0.9时,认为该测点处的介质为密实结构,并综合成图,可以清晰的得到轨道板检测情况的松散度示意图。通过松散度示意图能够清晰的获得轨道板下方介质的松散情况,为轨道板后期维修提供良好的指导意义。
以上显示和描述了本发明的基本原理、主要特征和本发明的优点。本行业的技术人员应该了解,本发明不受上述实施例的限制,上述实施例和说明书中描述的只是说明本发明的原理,在不脱离本发明精神和范围的前提下本发明还会有各种变化和改进,这些变化和改进都落入要求保护的本发明范围内。本发明要求保护范围由所附的权利要求书及其等同物界定。
Claims (7)
1.一种浅层介质结构面松散程度评价方法,其特征在于:包括以下步骤:1)、沿检测或勘察的介质表面布置测点;
2)、逐点激发采集每个测点处的直达波和反射波信号;
3)、对信号进行数据处理,获得测点处的直达波速度和反射波速度;
4)、按照如下公式计算该测点处介质结构面速度系数;
其中:Cv——结构面速度系数;
Vf——反射波速度;
Vz——直达波速度;
当结构面速度系数Cv小于a时,认为该测点处的介质为松散结构,当结构面速度系数Cv大于或等于a时,认为该测点处的介质为密实结构;
所述的步骤3)中,获得测点处的直达波速度的方法为,拾取直达波起跳时间,对拾取的直达波时间散点,采用直线方程进行拟合,以炮检距为自变量,直达波起跳时间为因变量,拟合的直线斜率即为直达波速度;
所述的步骤3)中,获得测点处的反射波速度的方法为,对已知深度的浅层介质界面,根据直达波速度设定速度范围,在该范围选取速度值计算该速度值下的反射波时距曲线,计算反射波时距曲线与炮集记录的最佳拟合即为反射波速度;所述的速度范围的设定方法为,根据直达波速度Vz划定速度范围为[1/2Vz,2/3Vz]。
2.如权利要求1所述的一种浅层介质结构面松散程度评价方法,其特征在于:所述的步骤1)中,按照地震反射波多次覆盖原理,在待检测或待勘察的介质表面沿线性布置测点。
3.如权利要求1所述的一种浅层介质结构面松散程度评价方法,其特征在于:所述的步骤2)中,每个测点在采集数据时,每个测点上激发一次机械波信号,每个测点上设置有多个传感器接收激发的信号,逐个测点依次进行,直至完成所有测点激发。
4.如权利要求1所述的一种浅层介质结构面松散程度评价方法,其特征在于:所述的步骤3)中,数据处理的方法,通过带通滤波器对获取数据进行处理,去除数据中的零飘和野值。
5.如权利要求1所述的一种浅层介质结构面松散程度评价方法,其特征在于:反射波时距曲线与炮集记录的最佳拟合方法为:以时距曲线为基础,设定时间窗口,将该时间窗口内的炮集记录振幅值积分求和,最大振幅能量对应的速度即为最佳拟合反射速度。
6.如权利要求5所述的一种浅层介质结构面松散程度评价方法,其特征在于:所述的时间窗口为15~30us。
7.如权利要求1所述的一种浅层介质结构面松散程度评价方法,其特征在于:所述的a的取值为0.5≤a≤0.9。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811458091.8A CN109541690B (zh) | 2018-11-30 | 2018-11-30 | 一种浅层介质结构面松散程度评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811458091.8A CN109541690B (zh) | 2018-11-30 | 2018-11-30 | 一种浅层介质结构面松散程度评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109541690A CN109541690A (zh) | 2019-03-29 |
CN109541690B true CN109541690B (zh) | 2020-05-19 |
Family
ID=65851911
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811458091.8A Active CN109541690B (zh) | 2018-11-30 | 2018-11-30 | 一种浅层介质结构面松散程度评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109541690B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110426739B (zh) * | 2019-08-02 | 2021-07-16 | 中铁第四勘察设计院集团有限公司 | 一种地质勘探检测方法、装置和存储介质 |
CN110632667B (zh) * | 2019-10-18 | 2020-12-18 | 徐州工程学院 | 一种基于冲击波激震条件下的隐伏陷落柱超前探测方法 |
CN111622737B (zh) * | 2020-05-20 | 2023-03-10 | 陕西省引汉济渭工程建设有限公司 | 一种基于钻孔声波分布快速确定地层岩体rqd的方法 |
CN112379418B (zh) * | 2020-10-16 | 2022-06-17 | 山东大学 | 一种地震直达波波速计算方法及*** |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6925031B2 (en) * | 2001-12-13 | 2005-08-02 | Baker Hughes Incorporated | Method of using electrical and acoustic anisotropy measurements for fracture identification |
CN1245637C (zh) * | 2003-11-25 | 2006-03-15 | 饶其荣 | 管波探测法 |
JP4277219B2 (ja) * | 2005-02-28 | 2009-06-10 | 清水建設株式会社 | 岩盤の亀裂密度と飽和度の推定方法 |
CN102087367B (zh) * | 2010-11-19 | 2012-10-03 | 中冶集团武汉勘察研究院有限公司 | 纵波-面波速度比值法划分基岩风化带的方法 |
CN102288990B (zh) * | 2011-08-03 | 2012-12-19 | 浙江大学 | 超浅层三维地震反射波勘探方法 |
CN103217703B (zh) * | 2013-04-11 | 2016-03-23 | 中国科学院武汉岩土力学研究所 | 一种以掌子面放炮为震源的隧道超前地质预报方法及装置 |
CN104264722B (zh) * | 2014-09-30 | 2016-08-24 | 云南航天工程物探检测股份有限公司 | 成桥基桩无损检测方法 |
-
2018
- 2018-11-30 CN CN201811458091.8A patent/CN109541690B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109541690A (zh) | 2019-03-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109541690B (zh) | 一种浅层介质结构面松散程度评价方法 | |
CN104007176B (zh) | 一种复杂岩土工程介质的全波场检测***及方法 | |
US10451475B2 (en) | Gauge length optimization in distributed vibration sensing | |
WO2017024523A1 (zh) | 一种射线弹性参数的反演方法 | |
CN102680575B (zh) | 一种复杂岩土介质的冲击映像方法及*** | |
DK177865B1 (da) | Fremgangsmåde til detektering eller monitorering af en subsurface-struktur af carbonhydridreservoirstørrelse | |
Park | MASW horizontal resolution in 2D shear-velocity (Vs) mapping | |
CN104264718A (zh) | 一种桩基完整性的检测方法 | |
Shirole et al. | Damage monitoring in rock specimens with pre-existing flaws by non-linear ultrasonic waves and digital image correlation | |
CN109884709B (zh) | 一种基于面波旅行时层析的转换波静校正方法 | |
CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
CN113156427A (zh) | 探地雷达数据的反演方法和装置 | |
CN112857698B (zh) | 一种基于声表面波检测墙体渗漏的方法 | |
CN110988143A (zh) | 一种隧道混凝土管片缺陷检测方法及设备 | |
JP6531934B2 (ja) | ハイブリッド表面波探査方法及びハイブリッド表面波探査システム | |
CN108508093A (zh) | 一种工件缺陷高度的检测方法及*** | |
CN114994752A (zh) | 一种弹性波地震频率谐振勘探方法 | |
CN109541688B (zh) | 一种基于能量特征的结构面缺陷评价方法 | |
CN112257241B (zh) | 一种三角网菲涅尔带时差层析反演方法 | |
CN110133721B (zh) | 一种水力压裂过程监测方法及*** | |
Muste et al. | Acoustic Mapping Velocimetry proof-of-concept experiment | |
CN112526611A (zh) | 表层地震波品质因子的提取方法及装置 | |
CN109024198A (zh) | 一种基于被动源随机激励的道面内部损伤检测方法 | |
CN105866394A (zh) | 岩土工程介质的信号处理方法、装置及*** | |
CN111175379B (zh) | 一种基于同步压缩小波变换的Lamb波板结构健康监测方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |