CN114137609B - 线性微动数据校正方法及装置 - Google Patents
线性微动数据校正方法及装置 Download PDFInfo
- Publication number
- CN114137609B CN114137609B CN202111320785.7A CN202111320785A CN114137609B CN 114137609 B CN114137609 B CN 114137609B CN 202111320785 A CN202111320785 A CN 202111320785A CN 114137609 B CN114137609 B CN 114137609B
- Authority
- CN
- China
- Prior art keywords
- data component
- axis data
- data
- axis
- component
- 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
- 238000012937 correction Methods 0.000 title claims abstract description 27
- 238000012545 processing Methods 0.000 claims abstract description 19
- 238000004458 analytical method Methods 0.000 claims description 21
- 230000010287 polarization Effects 0.000 claims description 12
- 238000001228 spectrum Methods 0.000 claims description 12
- 238000004891 communication Methods 0.000 claims description 7
- 230000008878 coupling Effects 0.000 claims description 2
- 238000010168 coupling process Methods 0.000 claims description 2
- 238000005859 coupling reaction Methods 0.000 claims description 2
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000011218 segmentation Effects 0.000 claims description 2
- 238000010586 diagram Methods 0.000 claims 1
- 238000005259 measurement Methods 0.000 abstract description 5
- 238000000516 activation analysis Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 230000005284 excitation Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 206010037844 rash Diseases 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 230000001131 transforming effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/53—Statics correction, e.g. weathering layer or transformation to a datum
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
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
本发明涉及地球物理勘探技术领域,提供一种线性微动数据校正方法及装置。本发明提供的线性微动数据校正方法,包括:将多个检波器按照编号顺序依次等间距线性排布,并将每个所述检波器的指北针指向北方设置;测量起点检波器和终点检波器的位置坐标;获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面。本发明提供的线性微动数据校正方法,采用线性布设方式,不受场地噪声源的方向影响,在检波器线性排列的方向与天然源噪声的优势方向不同时,也能获得正确的横波速度值,从而增强了线性微动测量的适用性。
Description
技术领域
本发明涉及地球物理勘探技术领域,尤其涉及一种线性微动数据校正方法及装置。
背景技术
微动法是近几年发展起来的一种新兴的勘探方法,常用于水利、交通、市政等行业的地质勘探,其原理是利用自然界存在的震动,通过专业的仪器设备采集震动信号,采用专业的数据软件处理得到地下地层的面波速度结构,从而为规划、设计和地质等专业提供数据支撑。目前,微动法的数据采集方式台阵布设方式主要有嵌套三角形、圆形、L形和线性等,其中线性的排列方式采集效率最高、布设方式最简单。
目前,微动法所采用的线性采集方式有以下几点不足:①微动采集的是天然源的噪声,在自然条件下,天然源的噪声来自四面八方,但在某些情况下,由于人为的干扰或突发的自然现场(地震,火山爆发等),天然源的噪声具有优势的来源方向,因此理论上来说采用线性排列时排列的方向应指向噪声的优势来源方向,从而保证采集的资料信噪比最高、反演的速度最真实,但是在实际生产过程中,往往很难将排列指向噪声的优势来源方向,原因是:一、天然源的噪声很微弱,工作现场很难有条件进行噪声源调查;二、线性排列的方向一般与设计好的工作测线方向重叠,不允许偏离设计好的工作测线。②当线性排列的方向不指向天然源的优势来源方向时,采集的数据经过处理以后,得到的横波速度值与真实的横波速度值不符,从而影响资料解释的精度。
发明内容
本发明提供一种线性微动数据校正方法及装置,用以解决现有技术中线性排列的方向不指向天然源的优势方向时,所得到的横波速度值不准确的缺陷。
本发明提供一种线性微动数据校正方法,包括:将多个检波器按照编号顺序依次等间距线性排布,并将每个所述检波器的指北针指向北方设置;测量起点检波器和终点检波器的位置坐标;获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面。
根据本发明提供的一种线性微动数据校正方法,所述将多个检波器按照编号顺序依次等间距线性排布,并将每个所述检波器的指北针指向北方设置的步骤还包括:将每个所述检波器与地面垂直摆放,并使每个所述检波器均与地面耦合。
根据本发明提供的一种线性微动数据校正方法,所述测量起点检波器和终点检波器的位置坐标的步骤还包括:当多个所述检波器线性排布存在拐点时,测量处于拐点位置的检波器的位置坐标。
根据本发明提供的一种线性微动数据校正方法,所述获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤进一步包括:利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角。
根据本发明提供的一种线性微动数据校正方法,所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤进一步包括:计算所述x轴数据分量和所述y轴数据分量的激化分析图,确定方位角的象限。
根据本发明提供的一种线性微动数据校正方法,所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤还包括:计算所述x轴数据分量和所述y轴数据分量的振幅谱的激化分析图,得到天然源噪声来源的优势方向与正东方向的夹角。
根据本发明提供的一种线性微动数据校正方法,所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤还包括:当所述x轴数据分量和所述y轴数据分量的振幅谱散点数据存在两个旁支时,对所述x轴数据分量和所述y轴数据分量的振幅谱数据按照不同的频率进行分段切割,得到方位角数据。
根据本发明提供的一种线性微动数据校正方法,所述获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤还包括:利用所述z轴数据分量和所述方位角数据,得到所述z轴数据分量校正后的数据sTz,对所述z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面;其中,sTz为z轴数据分量校正后的数据。
根据本发明提供的一种线性微动数据校正方法,所述利用所述z轴数据分量和所述方位角数据,得到所述z轴数据分量校正后的数据sTz,对所述z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面的步骤进一步包括:对所述z轴数据分量进行傅里叶变换,得到频率域数据Fsz,利用所述方位角数据对频率域数据Fsz进行方位角校正得到校正后的频率域数据FsTz,对所述校正后的频率域数据FsTz进行傅里叶逆变换得到所述z轴数据分量校正后的数据sTz;其中,Fsz为频率域数据,FsTz为校正后的频率域数据。
本发明还提供一种执行如上所述的线性微动数据校正方法的装置,包括:多个检波器,多个所述检波器按照编号顺序等间距线性排布,多个所述检波器之间通讯连接,每个所述检波器与地面垂直设置,每个所述检波器的指北针指向北方设置;全站仪,所述全站仪用于测量首个检波器和末个检波器的位置坐标;数据终端,所述数据终端与每个所述检波器通讯连接;其中,所述检波器为三分量检波器,所述检波器的数量大于7个。
本发明提供的线性微动数据校正方法,采用线性布设方式,不受场地噪声源的方向影响,在检波器线性排列的方向与天然源噪声的优势方向不同时,也能获得正确的横波速度值,从而增强了线性微动测量的适用性。
附图说明
为了更清楚地说明本发明或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1是本发明提供的线性微动数据校正方法的流程图;
图2是原始数据的激化分析图;
图3是原始数据的振幅谱的激化分析图;
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合本发明中的附图,对本发明中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的说明书和权利要求书中的术语“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本发明的描述中,除非另有说明,“多个”的含义是两个或两个以上。
下面结合图1-图3描述本发明的线性微动数据校正方法及装置。
如图1所示,在本发明的一个实施例中,线性微动数据校正方法,具体包括以下步骤:
步骤101:将多个检波器按照编号顺序依次等间距线性排布,并将每个检波器的指北针指向北方设置。
具体来说,在数据采集前,先对多个检波器进行编号,然后将多个检波器按照编号顺序进行线性排布,并保证相邻两个检波器之间的间距相等。在将检波器摆放时,多个检波器之间无线连接。可选地,在本发明的实施例中,检波器的数量应大于7个,检波器为三分量检波器。
步骤102:测量起点检波器和终点检波器的位置坐标。
具体来说,在本实施例中,多个检波器沿地面等间距线性排布,其排布方式为一条直线,全站仪仅需要测量首个检波器和末个检波器的位置坐标,中间各个检波器的位置坐标可根据间隔排布的数值推算出来,降低了现场数据的采集难度,也简化了测量工序。
步骤103:获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面。
具体来说,检波器摆放完毕以后,开始数据采集。当探测目的层较浅时,数据采集时间一般为10-15分钟;当探测目的层较深时,数据采集时间一般为30-60分钟,采集的数据均来自于自然界存在的噪声,通过数据终端可实时查看采集到的数据的质量。
数据采集完毕以后,从检波器中导出三个分量数据,分别为x轴数据分量、y轴数据分量和z轴数据分量,设这三个方向采集到的数据分量分别为(sx,sy,sz),利用sx和sy做极化分析,找出天然源噪声来源的优势方向,计算得到优势方向的方位角α。利用z轴数据分量和方位角,得到z轴数据分量校正后的数据sTz,对z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面。
本发明实施例提供的线性微动数据校正方法,采用线性布设方式,不受场地噪声源的方向影响,在检波器线性排列的方向与天然源噪声的优势方向不同时,也能获得正确的横波速度值,从而增强了线性微动测量的适用性。
进一步地,在本发明的一个实施例中,将多个检波器按照编号顺序依次等间距线性排布,并将每个检波器的指北针指向北方设置的步骤还包括:将每个检波器与地面垂直摆放,并使每个检波器均与地面耦合。
具体来说,在检波器摆放完毕以后,开始数据采集之前,应确保检波器与地面完全耦合,每个检波器均与地面垂直摆放,保证调平气泡在居中位置,同时,每个检波器的指北针需要完全指北。
在本发明的一个实施例中,测量起点检波器和终点检波器的位置坐标的步骤还包括:当多个检波器线性排布存在拐点时,测量处于拐点位置的检波器的位置坐标。
具体来说,当多个检波器的排列方式不是一条直线,而是一条曲线时,还需要检测处于拐点处的检波器的位置坐标。
在本发明的一个实施例中,获取x轴数据分量、y轴数据分量和z轴数据分量,并对x轴数据分量、y轴数据分量和z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤进一步包括:利用x轴数据分量和y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角。
具体来说,计算x轴数据分量和y轴数据分量的激化分析图,确定优势方向的方位角的象限。如图2所示,图2中确定出的优势方向的方向角的象限为第二象限和第四象限。在图2中,方位角分布在两个象限中,说明噪声源的位置随着时间的推移发生了变化,在这种情况下,需要对原始的x轴数据分量和y轴数据分量进行分时分段切割,将在同一象限的数据进行单独处理。
进一步地,通过sx和sy的振幅谱散点数据线性拟合得到y轴数据分量关于x轴数据分量的最佳拟合曲线Y=f(x,b),并得到参数b的最佳估计值,参数b包含了噪声源来源的优势方向与正东方向的锐角的夹角β。将上一步计算得到的象限值带入夹角β可得到噪声来源的优势方向的方位角α。
如图3所示,图3所示的sx和sy的振幅谱散点数据存在2个旁支(旁支1和旁支2),这说明不同频率的噪声的优势方向不一样,因此需对sx和sy的振幅谱数据按照不同的频率进行分段切割,将相同来源方向的频率段进行单独的处理,从而计算得到方位角数据数据A(α1,α2,α3...)。
在本发明的一个实施例中,获取x轴数据分量、y轴数据分量和z轴数据分量,并对x轴数据分量、y轴数据分量和z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤还包括:利用z轴数据分量和方位角数据,得到z轴数据分量校正后的数据sTz,对z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面;其中,sTz为z轴数据分量校正后的数据。
具体来说,对sz进行傅里叶变换,得到频率域数据Fsz,利用方位角数据数据A对频率域数据Fsz进行方位角校正得到FsTz=Fsz·(CAD),其中C和D为变换矩阵,FsTz为校正后的频率域数据。对校正后的频率域数据FsTz进行傅里叶逆变换得到z轴数据分量校正后的数据sTz。利用z轴数据分量校正后的数据sTz做微动数据处理分析,消除角度偏差,从而处理得到正确的横波速度剖面。
本发明实施例还提供了一种执行线性微动数据校正方法的装置,包括:多个检波器、全站仪和数据终端。多个检波器按照编号顺序等间距线性排布,多个检波器之间通讯连接,每个检波器与地面垂直设置,每个检波器的指北针指向北方设置。全站仪用于测量首个检波器和末个检波器的位置坐标,数据终端与每个检波器通讯连接,其中,检波器为三分量检波器,检波器的数量大于7个。
具体来说,在数据采集前,先对多个检波器进行编号,然后将多个检波器按照编号顺序进行线性排布,并保证相邻两个检波器之间的间距相等。在将检波器摆放时,多个检波器之间无线连接。每个检波器将采集到的数据发送至数据终端,用户可在数据终端查看各采集数据的质量。当检波器与数据终端的距离较远,无法通讯时,距离较近的检波器可作为中继器,所有的检波器均具有中继功能,远端的检波器可通过无线自动搜索到近端的检波器的中继信号,最终与数据终端连接,在中继信号搜索过程中,能够自动匹配最佳的中继信号,从而保证远端检波器、中继检波器以及数据终端之间的无线信号强度最优,最终保证数据传输的速度满足数据采集的要求。
进一步地,在本实施例中,在检波器摆放完毕以后,开始数据采集之前,应确保检波器与地面完全耦合,每个检波器均与地面垂直摆放,保证调平气泡在居中位置,同时,每个检波器的指北针需要完全指北,以保证计算出的噪声源的优势方向的方位角α的准确性。
进一步地,检波器为三分量检波器,在数据采集完毕以后,可从检波器中导出x轴数据分量、y轴数据分量和z轴数据分量。
进一步地,全站仪用于测试首个检波器和末个检波器的位置坐标,当检波器线性排列存在拐点时,全站仪还用于测量处于拐点处的检波器的位置坐标。
本发明实施例提供的执行线性微动数据校正方法的装置,通过将多个检波器线性排布,并将每个检波器的指北针指向北方,使线性排布方向不受场地噪声源的方向影响,在检波器线性排列的方向与天然源噪声的优势方向不同时,也能获得正确的横波速度值,从而增强了线性微动测量的适用性。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。
Claims (7)
1.一种线性微动数据校正方法,其特征在于,包括:
将多个检波器按照编号顺序依次等间距线性排布,并将每个所述检波器的指北针指向北方设置;
测量起点检波器和终点检波器的位置坐标;
获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面;
所述获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤进一步包括:利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角;
所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤进一步包括:计算所述x轴数据分量和所述y轴数据分量的极化分析图,确定方位角的象限;
所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤还包括:计算所述x轴数据分量和所述y轴数据分量的振幅谱的极化分析图,得到天然源噪声来源的优势方向与正东方向的夹角;
当方位角分布于多个象限时,对原始的所述x轴数据分量和所述y轴数据分量进行分时分段切割,将在同一象限的数据进行单独处理。
2.根据权利要求1所述的线性微动数据校正方法,其特征在于,所述将多个检波器按照编号顺序依次等间距线性排布,并将每个所述检波器的指北针指向北方设置的步骤还包括:
将每个所述检波器与地面垂直摆放,并使每个所述检波器均与地面耦合。
3.根据权利要求1所述的线性微动数据校正方法,其特征在于,所述测量起点检波器和终点检波器的位置坐标的步骤还包括:
当多个所述检波器线性排布存在拐点时,测量处于拐点位置的检波器的位置坐标。
4.根据权利要求1所述的线性微动数据校正方法,其特征在于,所述利用所述x轴数据分量和所述y轴数据分量进行极化分析,得到天然源噪声来源的优势方向的方位角的步骤还包括:
当所述x轴数据分量和所述y轴数据分量的振幅谱散点数据存在两个旁支时,对所述x轴数据分量和所述y轴数据分量的振幅谱数据按照不同的频率进行分段切割,得到方位角数据。
5.根据权利要求4所述的线性微动数据校正方法,其特征在于,所述获取x轴数据分量、y轴数据分量和z轴数据分量,并对所述x轴数据分量、所述y轴数据分量和所述z轴数据分量进行数据处理,得到正确的横波速度剖面的步骤还包括:
利用所述z轴数据分量和所述方位角数据,得到所述z轴数据分量校正后的数据sTz,对所述z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面;
其中,sTz为z轴数据分量校正后的数据。
6.根据权利要求5所述的线性微动数据校正方法,其特征在于,所述利用所述z轴数据分量和所述方位角数据,得到所述z轴数据分量校正后的数据sTz,对所述z轴数据分量校正后的数据sTz进行微动数据处理分析,得到正确的横波速度剖面的步骤进一步包括:
对所述z轴数据分量进行傅里叶变换,得到频率域数据Fsz,利用所述方位角数据对频率域数据Fsz进行方位角校正得到校正后的频率域数据FsTz,对所述校正后的频率域数据FsTz进行傅里叶逆变换得到所述z轴数据分量校正后的数据sTz;
其中,Fsz为频率域数据,FsTz为校正后的频率域数据。
7.一种执行权利要求1-6中任一项所述的线性微动数据校正方法的装置,其特征在于,包括:
多个检波器,多个所述检波器按照编号顺序等间距线性排布,多个所述检波器之间通讯连接,每个所述检波器与地面垂直设置,每个所述检波器的指北针指向北方设置;
全站仪,所述全站仪用于测量首个检波器和末个检波器的位置坐标;
数据终端,所述数据终端与每个所述检波器通讯连接;
其中,所述检波器为三分量检波器,所述检波器的数量大于7个。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111320785.7A CN114137609B (zh) | 2021-11-09 | 2021-11-09 | 线性微动数据校正方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111320785.7A CN114137609B (zh) | 2021-11-09 | 2021-11-09 | 线性微动数据校正方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114137609A CN114137609A (zh) | 2022-03-04 |
CN114137609B true CN114137609B (zh) | 2023-12-01 |
Family
ID=80392534
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111320785.7A Active CN114137609B (zh) | 2021-11-09 | 2021-11-09 | 线性微动数据校正方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114137609B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105929445A (zh) * | 2016-04-18 | 2016-09-07 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
WO2017071438A1 (zh) * | 2015-10-28 | 2017-05-04 | 中国石油天然气股份有限公司 | 全向矢量地震数据处理方法、装置、计算机可读存储介质及设备 |
RU2620785C1 (ru) * | 2016-06-14 | 2017-05-29 | Общество с ограниченной ответственностью "Макросейс" | Способ определения местоположения очага микросейсмического события |
CN107290722A (zh) * | 2017-06-29 | 2017-10-24 | 中国石油大学(北京) | 微震源的定位方法和装置 |
CN107807381A (zh) * | 2017-12-01 | 2018-03-16 | 招商局重庆交通科研设计院有限公司 | 基于岩体破裂微震波活动规律的边坡失稳风险的动态监测方法及装置 |
CN107918156A (zh) * | 2017-10-30 | 2018-04-17 | 中国石油天然气集团公司 | 检测海底节点采集地震数据极性的方法及装置 |
CN110488349A (zh) * | 2019-08-20 | 2019-11-22 | 福建省建筑设计研究院有限公司 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
CN112526597A (zh) * | 2020-11-09 | 2021-03-19 | 中铁第四勘察设计院集团有限公司 | 一种微动技术的勘探方法、装置和存储介质 |
CN113267806A (zh) * | 2021-05-28 | 2021-08-17 | 长江水利委员会长江科学院 | Tbm刀盘破岩噪声源的多波采集系统及超前探测方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9001619B2 (en) * | 2011-10-19 | 2015-04-07 | Global Microseismic Services, Inc. | Method for imaging microseismic events using an azimuthally-dependent focal mechanism |
US9784866B2 (en) * | 2013-07-28 | 2017-10-10 | Geokinetics Usa, Inc. | Method and apparatus for enhanced monitoring of induced seismicity and vibration using linear low frequency and rotational sensors |
-
2021
- 2021-11-09 CN CN202111320785.7A patent/CN114137609B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2017071438A1 (zh) * | 2015-10-28 | 2017-05-04 | 中国石油天然气股份有限公司 | 全向矢量地震数据处理方法、装置、计算机可读存储介质及设备 |
CN105929445A (zh) * | 2016-04-18 | 2016-09-07 | 北京京援伟达技术有限公司 | 微破裂向量扫描方法 |
RU2620785C1 (ru) * | 2016-06-14 | 2017-05-29 | Общество с ограниченной ответственностью "Макросейс" | Способ определения местоположения очага микросейсмического события |
CN107290722A (zh) * | 2017-06-29 | 2017-10-24 | 中国石油大学(北京) | 微震源的定位方法和装置 |
CN107918156A (zh) * | 2017-10-30 | 2018-04-17 | 中国石油天然气集团公司 | 检测海底节点采集地震数据极性的方法及装置 |
CN107807381A (zh) * | 2017-12-01 | 2018-03-16 | 招商局重庆交通科研设计院有限公司 | 基于岩体破裂微震波活动规律的边坡失稳风险的动态监测方法及装置 |
CN110488349A (zh) * | 2019-08-20 | 2019-11-22 | 福建省建筑设计研究院有限公司 | 基于微动三分量谱比vhsr的无损探测方法及应用 |
CN112526597A (zh) * | 2020-11-09 | 2021-03-19 | 中铁第四勘察设计院集团有限公司 | 一种微动技术的勘探方法、装置和存储介质 |
CN113267806A (zh) * | 2021-05-28 | 2021-08-17 | 长江水利委员会长江科学院 | Tbm刀盘破岩噪声源的多波采集系统及超前探测方法 |
Non-Patent Citations (4)
Title |
---|
Multichannel analysis of passive surface waves based on crosscorrelations;Feng Cheng 等;《GEOPHYSICS》;第81卷(第15期);57-66 * |
微动勘探技术在城市轨道交通勘察中的应用;张建清 等;人民长江;47(1);45-48 * |
微动勘查技术在城市地下空洞探测应用研究;董耀 等;《河南省地质学会2020年学术年会论文集》;209-214 * |
董耀 等.微动勘查技术在城市地下空洞探测应用研究.《河南省地质学会2020年学术年会论文集》.2020,209-214. * |
Also Published As
Publication number | Publication date |
---|---|
CN114137609A (zh) | 2022-03-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103114848B (zh) | 一种基于岩心测量的地层裂缝空间重构方法 | |
WO2016141598A1 (zh) | 共姿态道集地表一致性振幅补偿方法及装置 | |
CN106525976B (zh) | 基于声发射层析成像定量分析混凝土结构损伤部位的方法 | |
CN112965135B (zh) | 一种石窟崖体裂隙空间异质分布的无损探测综合方法 | |
CN106646375A (zh) | 一种基于all‑fft相位差的岩石声发射源定位方法 | |
CN106291542A (zh) | 一种隧道三维成像方法 | |
CN112666554A (zh) | 一种沥青路面雷达振幅特征裂缝宽度识别方法 | |
CN106705830A (zh) | 一种基于北斗卫星的超高层建筑高精度变形监测系统及监测方法 | |
CN113156427A (zh) | 探地雷达数据的反演方法和装置 | |
CN109541695B (zh) | 人工场源频率域电场梯度远区视电阻率快速成像方法 | |
CN107063597B (zh) | 基于pos系统的桥梁检测方法 | |
CN112326785A (zh) | 同步注浆充填效果冲击映像法检测与评价方法 | |
CN109541689B (zh) | 一种基于反射波能量特征的介质密实度评价方法 | |
CN114137609B (zh) | 线性微动数据校正方法及装置 | |
CN111751259A (zh) | 确定不规则岩石样品有效孔隙度的方法及装置 | |
CN112325808B (zh) | 一种基于多psd的平面度实时校准补偿测量方法 | |
CN113777654B (zh) | 一种基于伴随状态法初至波走时层析的海水速度建模方法 | |
CN114415245A (zh) | 一种适用于中深部蚀变岩型金矿找矿方法 | |
CN113900145B (zh) | 一种基于微动信号的地下横波速度分布的推测方法和装置 | |
CN113075732B (zh) | 一种消除高分辨率小多道地震地层异常起伏的方法 | |
KR101799813B1 (ko) | 공통반사표면법을 이용한 지하투과레이더 탐사의 자료처리 결과 신호대잡음비 향상과 브이알엠에스 역산과정을 이용한 씽크홀 탐사장치 및 그 방법 | |
CN109521464B (zh) | 一种燃气pe管线定位方法 | |
CN109143405B (zh) | 一种观测系统有效采样均匀性定量评价方法 | |
CN109211151A (zh) | 型材的检测装置、方法、设备和介质 | |
EP4253931A1 (en) | Non-invasive diagnostic system for real time detection of near-surface water leakage in urban and suburban contexts |
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 |