CN110146924A - 基于水波初至偏振方位的海底地震仪位置及方位反演方法 - Google Patents
基于水波初至偏振方位的海底地震仪位置及方位反演方法 Download PDFInfo
- Publication number
- CN110146924A CN110146924A CN201910598077.6A CN201910598077A CN110146924A CN 110146924 A CN110146924 A CN 110146924A CN 201910598077 A CN201910598077 A CN 201910598077A CN 110146924 A CN110146924 A CN 110146924A
- Authority
- CN
- China
- Prior art keywords
- orientation
- submarine seismograph
- modified
- ripples
- arrival
- 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.)
- Granted
Links
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/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
- G01V1/3843—Deployment of seismic devices, e.g. of streamers
- G01V1/3852—Deployment of seismic devices, e.g. of streamers to the seabed
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/14—Signal detection
- G01V2210/142—Receiver location
- G01V2210/1427—Sea bed
Abstract
本发明公开了一种基于水波初至偏振方位的海底地震仪位置及方位反演方法,采用非线性反演方法,对于一个待修正的海底地震仪,利用水波初至偏振方位残差之差构造目标函数,并求取令目标函数最小的地震仪位置,在求得精确的地震仪位置之后,再以炮点方位与偏振方位之差求得地震仪的实际方位指向。尤其是以多对水波初至偏振方位残差之差的平方和作为目标函数,应用DE最优化方法求取使目标函数取最小值的地震仪位置的方案,在目标函数中消去了地震仪方位的影响,具有较快的收敛速度和较好的稳定性,相比线性方法具有较好的精度。
Description
技术领域
本发明涉及海底地震仪的位置及方法反演方法,具体涉及一种基于水波初至偏振方位的海底地震仪位置及方位反演方法。
背景技术
OBS海底地震仪在布放到海底之后,其方位指向通常是随机的,由于海流等因素的影响,其实际位置通常也会偏离预定位置,因此需要对其位置及方位进行反演处理。传统方法采用线性反演的方式,在一些特殊情况下,存在不收敛的情况。
发明内容
针对现有技术的不足,本发明旨在提供一种基于水波初至偏振方位的海底地震仪位置及方位反演方法,采用非线性反演方法,对于一个待修正的海底地震仪,以多对水波初至偏振方位残差之差的平方和作为目标函数,应用DE最优化方法求取使目标函数取最小值的地震仪位置,在目标函数中消去了地震仪方位的影响,具有较快的收敛速度和较好的稳定性。
为了实现上述技术目的,本发明采用如下技术方案:
基于水波初至偏振方位的海底地震仪位置及方位反演方法,具体为:
S1、根据待修正的海底地震仪的设计布放位置,利用该地震仪接收N个不同位置的炮点的三分量地震波形信号;
S2、对于一个待修正的海底地震仪,利用信号的水波初至偏振方位残差构造目标函数,并求取使目标函数取最小值的待修正的海底地震仪位置;
S3、利用步骤S1求得使目标函数取最小值的待修正的海底地震仪位置计算待修正的海底地震仪位置指向炮点方向的方位角,再以该待修正的海底地震仪位置指向炮点的方位角与炮点信号的水波初至偏振方位角之差求得待修正的海底地震仪的实际方位指向,并据此得到海底地震仪方位的修正量。
进一步地,步骤S2中,以多对信号的水波初至偏振方位残差之差的平方和作为目标函数,应用DE最优化方法求取使目标函数取最小值的待修正的海底地震仪位置。
进一步地,步骤S2的具体过程为:
设待修正的海底地震仪在海底的设计布放位置为(x,y),实际位置为(x+δx,y+δy),(δx,δy)为偏移量;
用多个炮点对的方位残差之差的平方和构造目标函数S2:
和别为待修正的海底地震仪位置指向第i号和第j号炮点的方位角的观测值,i≠j,表达式如下:
和分别为待修正的海底地震仪位置指向第i号和第j号炮点方向的方位角的理论值,i≠j,表达式如下:
φ0表示待修正的海底地震仪的方位的修正量;和分别为根据第i号和第j号炮点的信号的三分量水波初至波形计算得到的偏振方位角,和分别表示待修正的海底地震仪在初始位置偏离(δx,δy)时指向第i号和第j号炮点的理论方位角;
由此可得:
采用DE非线性最优化方法求取使S2取最小值的偏移量(δx,δy)。
作为另一种实施方式,步骤S2中,设待修正的海底地震仪在海底的设计布放位置为(x,y),实际位置为(x+δx,y+δy),(δx,δy)为偏移量;
构造目标函数S1,求取使S1取最小值的(φ0,δx,δy):
其中,为待修正的海底地震仪位置指向第i号炮点的方位角的观测值,表达式如下:
为待修正的海底地震仪位置指向第i号炮点方向的方位角的理论值,表达式如下:
S1是所有炮点的方位残差的平方和,表示在初始位置偏离(δx,δy)时指向i号炮点的理论方位角,N为炮点的数量;φ0表示待修正的海底地震仪的方位的修正量;为根据第i号炮点的信号的三分量水波初至波形计算得到的偏振方位角。
更进一步地,步骤S3具体为:
针对第i号炮点计算的方位修正量φ0为:
对于一个待修正的海底地震仪,针对每个炮点i,1≤i≤N都可计算得到一个方位修正量φ0,对所有结果求平均值,即得到最终的海底地震仪方位修正量φ'0:
本发明的有益效果在于:采用非线性反演方法,对于一个待修正的海底地震仪,利用水波初至偏振方位残差之差构造目标函数,并求取令目标函数最小的地震仪位置,在求得精确的地震仪位置之后,再以炮点方位与偏振方位之差求得地震仪的实际方位指向,准确性较好。尤其是以多对水波初至偏振方位残差之差的平方和作为目标函数,应用DE最优化方法求取使目标函数取最小值的地震仪位置的方案,在目标函数中消去了地震仪方位的影响,具有较快的收敛速度和较好的稳定性。
附图说明
图1为本发明实施例3中实际数据采用第二种方法的迭代收敛情况示意图;
图2为本发明实施例3中实际数据采用第一种方法的迭代收敛情况示意图。
具体实施方式
以下将结合附图对本发明作进一步的描述,需要说明的是,本实施例以本技术方案为前提,给出了详细的实施方式和具体的操作过程,但本发明的保护范围并不限于本实施例。
实施例1
本实施例提供一种基于水波初至偏振方位的海底地震仪位置及方位反演方法,具体为:
S1、根据待修正的海底地震仪的设计布放位置,利用该地震仪接收N个不同位置的炮点的三分量地震波形信号
S2、设待修正的海底地震仪在海底的设计布放位置为(x,y),实际位置为(x+δx,y+δy),(δx,δy)为偏移量。
根据纵波的性质,纵波的偏振方向与传播方向相同,因此炮点放出的信号的水波初至偏振方位与待修正的海底地震仪位置指向炮点的方向一致,待修正的海底地震仪位置指向第i号炮点方向的观测值可以表示为:
其中为根据第i号炮点的三分量水波初至波形计算得到的偏振方位角,φ0表示待修正的海底地震仪的方位的修正量;待修正的海底地震仪位置指向第i号炮点方向的理论值可由待修正的海底地震仪位置与炮点位置计算得到:
本实施例提供两种方式计算待修正的海底地震仪的方位的修正量φ0。
1、作为一种实施方式,直接的方法是,构造目标函数S1,求取使S1取最小值的(φ0,δx,δy):
S1是所有炮点方位残差的平方和,表示在初始位置偏离(δx,δy)时指向i号炮点的理论方位角,N为炮点的数量。
2、作为另一种实施方式,可以用多个炮点对的方位残差之差的平方和构造目标函数S2:
将前述各式代入上式,得到:
φ0表示待修正的海底地震仪的方位的修正量;和分别为根据第i号和第j号炮点的信号的三分量水波初至波形计算得到的偏振方位角,和分别表示待修正的海底地震仪在初始位置偏离(δx,δy)时指向第i号和第j号炮点的理论方位角;
式(5)中,对于同一待修正的海底地震仪,方位修正量相同,因而被消掉,因此只含有两个未知数。
采用DE(deferential evolution)非线性最优化方法求取使S2取最小值的偏移量(δx,δy),针对第i号炮点计算的方位修正量φ0为:
对于一个待修正的海底地震仪,针对每个炮点i,1≤i≤N都可计算得到一个方位修正量φ0,对所有结果求平均值,即得到最终的海底地震仪方位修正量φ仩0:
实施例2
本实施例旨在提供一种如实施例1所述方法所得反演结果的误差计算方法。
本实施例采用bootstrap方法给出反演结果的误差估计,具体做法是根据给定的实际初至偏振数据和方位角误差信息,按实际的炮点/台站位置生成一组理论观测数据,并将原有数据的残差随机加到理论数据中,然后对这些理论观测数据进行反演计算,以结果的统计方差统作为误差。
对于一个海底地震仪,假定有N个炮点数据参与计算,所涉及数据包括台站水平位置及方位(lat,lon,φ0),炮点坐标、方位角及残差(lati,loni,θi,ri),根据这组数据生成一组“新的理论数据”的具体步骤如下:
1.生成随机序列Ri,i=1,N。先对Ri序列赋初值Ri=i;
2.对Ri中的每个元素进行循环,每次循环生成[1,N]区间内的一个随机数j,然后Ri和Rj的值对调。循环完成,则随机序列Ri确定。实际程序中采用C语言中的系统调用函数srand()/rand()生成随机数,并以调用时的系统时间作为seed。
3.则新的一组观测数据为
重复上述步骤生成M组“观测数据”,对每组数据进行反演得到M个结果,这些结果的统计方差即是反演结果的误差。
实施例3
本实施例旨在对实施例1所述的方法进行进一步的检验。
为检验上述方法,本实施例分别进行了理论模型试验和实际数据的处理,所用数据来自2003年在我国南海某区域进行的OBS观测,这次观测有5台OBS台站。理论模型采用其中的obs4的实际观测系统,在理论计算得到的方位角中加上均值为0并且符合高斯分布的随机扰动,随机扰动的方差为6°,(φ0,δx,δy)各值的真实值均为0。表1是进行了1000次这种计算的统计结果,从表1中的结果可以看出,非线性反演方法的结果好于线性反演的结果,各反演量的平均值更接近真实值,方差也较小。
表2是5个台站实际数据的反演结果。从结果可以看出,两种方法得到的地震仪方位φ0非常接近,而地震仪的位置有些台站两种方法偏差较大,达到30米左右,但结合表1理论模型的结果,这个偏差也在正常范围内。
表1
表2实际数据的反演结果
理论模型实验和实际数据的结果表明,采用公式(3)(5)的两种方法得到的结果几乎完全相同,但迭代收敛情况公式(5)明显好于公式(3)的方法。如图1-2所示,图1-2是一个实际的obs台站数据分别采用公式(5)(3)的迭代收敛情况,公式(5)在迭代到600多次时即达到最优稳定状态,而采用公式(3)的方法,在迭代1000多次之后才达到稳定最优状态。
相关参考文献:
[1]Christeson,G.L.,K.D.McIntosh,T.H.Shipley,E.R.Flueh,and H.Goedde(1999),Structure of the costa rica convergent margin,offshore nicoyapeninsula,JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH,104(B11),25,443–25,468.
[2]JURKEVICS,A.(1988),Polarization analysis of 3-component arraydata,B.SEISMOL.SOC.AM.,78(5),1725–1743.
[3]Li,X.Y.,andJ.X.Yuan(1999),Geophone orientation and coupling inthree-component sea-floor data:a case study,GEOPHYS.PROSPECT.,47(6),995–1013.
[4]NAKAMURA,Y.,P.L.DONOHO,P.H.ROPER,and P.M.MCPHERSON(1987),Large-offset seismic surveying using ocean-bottom seismographs and air guns-instrumentation and field technique,GEOPHYSICS,52(12),1601–1611.
[5]PARK,J.,F.L.VERNON,andC.R.LINDBERG(1987),Frequency-dependentpolarization analysis of high-frequency seismograms,J GEOPHYS RES-SOLID,92(B12),12,664–12,674.
[6]Ruzek,B.,and M.Kvasnicka(2001),Differential evolution algorithm inthe earthquake hypocenter location,PURE APPL GEOPHYS,158(4),667–693.
[7]Ruzek,B.,P.Hrubcova,M.Novotny,A.Spicak,and O.Karousova(2007),Inversion of travel times obtained during active seismic refractionexperiments celebration 2000,alp 2002 and sudetes 2003,STUD GEOPHYS GEOD,51(1),141–164.
对于本领域的技术人员来说,可以根据以上的技术方案和构思,给出各种相应的改变和变形,而所有的这些改变和变形,都应该包括在本发明权利要求的保护范围之内。
Claims (5)
1.基于水波初至偏振方位的海底地震仪位置及方位反演方法,其特征在于,具体为:
S1、根据待修正的海底地震仪的设计布放位置,利用该地震仪接收N个不同位置的炮点的三分量地震波形信号;
S2、对于一个待修正的海底地震仪,利用信号的水波初至偏振方位残差构造目标函数,并求取使目标函数取最小值的待修正的海底地震仪位置;
S3、利用步骤S1求得使目标函数取最小值的待修正的海底地震仪位置计算待修正的海底地震仪位置指向炮点方向的方位角,再以该待修正的海底地震仪位置指向炮点的方位角与炮点信号的水波初至偏振方位角之差求得待修正的海底地震仪的实际方位指向,并据此得到海底地震仪方位的修正量。
2.根据权利要求1所述的基于水波初至偏振方位的海底地震仪位置及方位反演方法,其特征在于,步骤S2中,以多对信号的水波初至偏振方位残差之差的平方和作为目标函数,应用DE最优化方法求取使目标函数取最小值的待修正的海底地震仪位置。
3.根据权利要求2所述的基于水波初至偏振方位的海底地震仪位置及方位反演方法,其特征在于,步骤S2的具体过程为:
设待修正的海底地震仪在海底的设计布放位置为(x,y),实际位置为(x+δx,y+δy),(δx,δy)为偏移量;
用多个炮点对的方位残差之差的平方和构造目标函数S2:
和别为待修正的海底地震仪位置指向第i号和第j号炮点的方位角的观测值,i≠j,表达式如下:
和分别为待修正的海底地震仪位置指向第i号和第j号炮点方向的方位角的理论值,i≠j,表达式如下:
φ0表示待修正的海底地震仪的方位的修正量;和分别为根据第i号和第j号炮点的信号的三分量水波初至波形计算得到的偏振方位角,和分别表示待修正的海底地震仪在初始位置偏离(δx,δy)时指向第i号和第j号炮点的理论方位角;
由此可得:
采用DE非线性最优化方法求取使S2取最小值的偏移量(δx,δy)。
4.根据权利要求1所述的基于水波初至偏振方位的海底地震仪位置及方位反演方法,其特征在于,步骤S2中,设待修正的海底地震仪在海底的设计布放位置为(x,y),实际位置为(x+δx,y+δy),(δx,δy)为偏移量;
构造目标函数S1,求取使S1取最小值的(φ0,δx,δy):
其中,为待修正的海底地震仪位置指向第i号炮点的方位角的观测值,表达式如下:
为待修正的海底地震仪位置指向第i号炮点方向的方位角的理论值,表达式如下:
S1是所有炮点的方位残差的平方和,表示在初始位置偏离(δx,δy)时指向i号炮点的理论方位角,N为炮点的数量;φ0表示待修正的海底地震仪的方位的修正量;为根据第i号炮点的信号的三分量水波初至波形计算得到的偏振方位角。
5.根据权利要求3或4所述的基于水波初至偏振方位的海底地震仪位置及方位反演方法,其特征在于,步骤S3具体为:
针对第i号炮点计算的方位修正量φ0为:
对于一个待修正的海底地震仪,针对每个炮点i,1≤i≤N都可计算得到一个方位修正量φ0,对所有结果求平均值,即得到最终的海底地震仪方位修正量φ'0:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910598077.6A CN110146924B (zh) | 2019-07-03 | 2019-07-03 | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910598077.6A CN110146924B (zh) | 2019-07-03 | 2019-07-03 | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110146924A true CN110146924A (zh) | 2019-08-20 |
CN110146924B CN110146924B (zh) | 2020-05-26 |
Family
ID=67596916
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910598077.6A Active CN110146924B (zh) | 2019-07-03 | 2019-07-03 | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110146924B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824551A (zh) * | 2019-11-14 | 2020-02-21 | 自然资源部第二海洋研究所 | 一种基于声学测距和多波束地形的obs精确定位方法 |
CN111257941A (zh) * | 2020-02-18 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111352160A (zh) * | 2020-03-19 | 2020-06-30 | 中国科学院地质与地球物理研究所 | 一种海底地震仪自动重定位装置及方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154334A (zh) * | 2015-04-13 | 2016-11-23 | 中石化石油工程地球物理有限公司胜利分公司 | 基于网格搜索的井下微地震事件实时反演定位方法 |
CN106249295A (zh) * | 2015-06-15 | 2016-12-21 | 中国石油化工股份有限公司 | 一种井中微地震p、s波联合快速定位方法及系统 |
WO2017062322A1 (en) * | 2015-10-07 | 2017-04-13 | Schlumberger Technology Corporation | Seismic sensor orientation |
-
2019
- 2019-07-03 CN CN201910598077.6A patent/CN110146924B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106154334A (zh) * | 2015-04-13 | 2016-11-23 | 中石化石油工程地球物理有限公司胜利分公司 | 基于网格搜索的井下微地震事件实时反演定位方法 |
CN106249295A (zh) * | 2015-06-15 | 2016-12-21 | 中国石油化工股份有限公司 | 一种井中微地震p、s波联合快速定位方法及系统 |
WO2017062322A1 (en) * | 2015-10-07 | 2017-04-13 | Schlumberger Technology Corporation | Seismic sensor orientation |
Non-Patent Citations (2)
Title |
---|
ADRIAN K. DORAN ET AL.: "Ocean-Bottom Seismometer Instrument Orientations via Automated Rayleigh-Wave Arrival-Angle Measurements", 《BULLETIN OF THE SEISMOLOGICAL SOCIETY OF AMERICA》 * |
J. C. STACHNIK ET AL.: "Determination of New Zealand Ocean Bottom Seismometer Orientation via Rayleigh-Wave Polarization", 《SEISMOLOGICAL RESEARCH LETTERS》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110824551A (zh) * | 2019-11-14 | 2020-02-21 | 自然资源部第二海洋研究所 | 一种基于声学测距和多波束地形的obs精确定位方法 |
CN111257941A (zh) * | 2020-02-18 | 2020-06-09 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111257941B (zh) * | 2020-02-18 | 2020-09-04 | 中国科学院地质与地球物理研究所 | 一种组合式海底地震仪方位角自动识别装置及方法 |
CN111352160A (zh) * | 2020-03-19 | 2020-06-30 | 中国科学院地质与地球物理研究所 | 一种海底地震仪自动重定位装置及方法 |
Also Published As
Publication number | Publication date |
---|---|
CN110146924B (zh) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Kent et al. | Distribution of magma beneath the East Pacific Rise near the 9 03′ N overlapping spreading center from forward modeling of common depth point data | |
CN110146924A (zh) | 基于水波初至偏振方位的海底地震仪位置及方位反演方法 | |
Kissling et al. | Three-dimensional structure of the Long Valley Caldera, California, region by geotomography | |
CN106814378B (zh) | 一种gnss位置时间序列周期特性挖掘方法 | |
CN108254780A (zh) | 一种微地震定位及各向异性速度结构层析成像方法 | |
CN106680869B (zh) | 微地震事件的检测和定位方法与装置 | |
CN106896403B (zh) | 弹性高斯束偏移成像方法和系统 | |
Nakatsuka et al. | Reduction of magnetic anomaly observations from helicopter surveys at varying elevations | |
CN113740915B (zh) | 一种球坐标系重力和接收函数联合反演地壳结构参数的方法 | |
CN112099082B (zh) | 一种共面元共方位角道集的地震回折波走时反演方法 | |
CN113960532A (zh) | 一种基于假想源的二次定位计算的微地震定位方法 | |
CN109100798B (zh) | 实现折射多次波层析反演的方法、装置及处理终端 | |
CN108508479B (zh) | 一种空地井立体重磁数据协同目标位置反演方法 | |
CN111142170B (zh) | 一种基于重力梯度极值点的潜艇位置探测方法 | |
WO2018130272A1 (en) | A method for fast determination of the location of an acoustic emission event in a vti/tti medium | |
CN107589446A (zh) | 利用高斯束计算波路径的层析成像速度建模方法 | |
Abubakirov et al. | The mechanism of the deep-focus, Sea of Okhotsk earthquake of May 24, 2013 as inferred from static displacements and broadband seismograms | |
Whitcombe et al. | The response of the time-term method to simulated crustal structures | |
CN106646527A (zh) | 基于三测站数据特征点的电离层扰动传播测量方法及系统 | |
Wang et al. | Ray-based crosshole radar traveltime tomography using MSFM method | |
Qingpo et al. | First Arrival Traveltime Tomography with Near Surface Model Constraints | |
Xiong‐Wei et al. | Comparison of Inversion Method of Wide Angle Ocean Bottom Seismometer Profile: A Case Study of Profile OBS973‐2 Across Liyue Bank in the South China Sea | |
CN111045078A (zh) | 一种复杂近地表条件下的初至波旅行时层析反演方法 | |
Guandong et al. | Method and application of offset vector summation in node repositioning based on refractions | |
CN117368849A (zh) | 一种基于炮点差分的全息二次定位方法及其应用 |
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 |