CN101937106A - 一种海底大起伏测线磁测数据的处理方法 - Google Patents
一种海底大起伏测线磁测数据的处理方法 Download PDFInfo
- Publication number
- CN101937106A CN101937106A CN 201010238099 CN201010238099A CN101937106A CN 101937106 A CN101937106 A CN 101937106A CN 201010238099 CN201010238099 CN 201010238099 CN 201010238099 A CN201010238099 A CN 201010238099A CN 101937106 A CN101937106 A CN 101937106A
- Authority
- CN
- China
- Prior art keywords
- magnetic
- survey
- data
- seafloor
- parameter
- 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
- 238000000034 method Methods 0.000 title abstract description 26
- 238000012545 processing Methods 0.000 title abstract description 4
- 238000012876 topography Methods 0.000 claims description 18
- 238000003672 processing method Methods 0.000 claims description 8
- 238000001914 filtration Methods 0.000 abstract description 2
- 238000006243 chemical reaction Methods 0.000 abstract 1
- 238000011160 research Methods 0.000 description 12
- 239000011159 matrix material Substances 0.000 description 10
- 230000000694 effects Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 230000018109 developmental process Effects 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- UCKMPCXJQFINFW-UHFFFAOYSA-N Sulphide Chemical compound [S-2] UCKMPCXJQFINFW-UHFFFAOYSA-N 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000013499 data model Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000007598 dipping method Methods 0.000 description 1
- 238000005290 field theory Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 230000005415 magnetization Effects 0.000 description 1
- 238000012821 model calculation Methods 0.000 description 1
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种海底大起伏测线磁测数据的处理方法,首先获得海底目标区的磁测数据、磁测点位置数据和海底地形数据;基于所得到的海底地形数据对海底目标区进行三维网格划分,建立磁参数模型;建立反映磁异常理论值与观测值之间误差大小的目标反演函数,以目标反演函数最小值为拟合标准,修改磁参数;最终由反演得到的磁参数计算目标区域的磁异常参量。本发明方法不需向上延拓、高截滤波、将起伏面上的观测数据转化到水平面上,避免了引起分辨率降低的中间过程,利用三维空间域带地形直接反演法对磁测数据反演,不仅提高了分辨率,而且对测点位置的起伏量没有限制,适用于大起伏测线的海底磁测。
Description
技术领域
本发明涉及海底地球物理勘探技术领域,尤其涉及一种海底大起伏测线磁测数据的处理方法。
背景技术
随着全球资源的日趋紧张,从上世纪80年代开始,世界各国都逐渐把目光转移到海洋上来,开始了对海洋领域资源的大力勘探。磁法勘探是地球物理勘探中一种不可替代的重要探测手段,主要有航空、地面、海面和海底磁测。
从理论上来看,几种磁测方式并无本质区别。但从勘探技术、方法和勘探效果上讲却各有特点。由于磁源场满足平方反比定律,磁测分辨率和勘探效果取决于源收距(目标体和观测点间的距离)。航空磁测效率高、不受地域限制,但由于源收距比地面磁测和海底磁测的源收距要大,分辨率相对地面和海底磁测要低,适合大规模普查。地面磁测最简易、勘探效果最好,但受地域限制。海面磁测方法和地面磁测相似,其勘探目标体在海底下面,海面磁测的源收距远大于航空和地面磁测,所以分辨率和勘探效果相比航空磁测和地面磁测较差,适合海域调查和深部地质研究。海底磁测的源收距比海面磁测和海域上的航空磁测的源收距要小得多,海底磁测相比海面磁测将大大提高勘探效果和分辨率,是航空磁测和海面磁测无法替代的,可提供海底火山、热液活动、硫化物、金属矿床等独特的成像技术,具有广泛的应用前景。
但是海底磁测相对于航空磁测和海面磁测在观测技术上非常复杂,难度很大。如图1所示,比较了几种磁测方式的磁测线。其中,航空磁测由搭载在飞机上的磁力仪测量飞行区域的磁数据。由于磁测时飞行起伏一般不超几十米,其磁测线的波动范围很小。同样,海面磁测由搭载在考察船上的磁力仪测量航行区域的磁数据,由于海面一般认为处于同一平面,其磁测线较之航空磁测线更为平滑。而海底磁测中磁力仪一般封装在耐压舱内,耐压舱放置在与考察船之间用缆绳连接的海底拖体内。这样磁测点位置就由拖***置确定,但由于深海海况一般比较复杂,拖***置不像航空磁测的飞机或海面的考察船那样能稳定地控制,造成测点位置起伏非常大,测点间起伏可达到几百上千米甚至几公里。导致海底磁测线的波动剧烈。同时洋脊海底的地形复杂,其轮廓线起伏也很剧烈。对于具有这种特征的磁测数据不能利用航空磁测或陆地磁测现有的常规处理,如延拓、曲化平、反演解释等技术来处理海底磁测数据,否则可能获得错误的结果。
这为海底磁测数据处理反演解释技术提出了新的研究课题,迫切需要发展有关海底测数据处理方法、技术等基础理论研究。因此开展海底磁测工作具有重要的现实意义并对磁法勘探理论发展具有重要意义。该项同类技术目前国外主要在频率域进行二、三维反演,它们需要滤波、向上延拓、要求将观测数据转化到水平面上,而这些会损失高频信号而降低分辨率。
发明内容
本发明旨在针对深海海底磁测数据具有测点位置起伏特别大、海底地形复杂等特征,不能用现有曲化平方法将观测数据转化到某一水平面上,从而无法用现有常规方法进行数据处理和反演解释这一问题,提供一种适合于这种特征的海底磁测数据的空间域带地形3D直接反演方法。
一种海底大起伏测线磁测数据的处理方法,包括:
(1)获得磁测数据、磁测点位置数据和海底地形数据;
在海底磁测中,磁测数据由拖体内的磁力仪采集获得,拖体被放置在深海中,与考察船之间用缆绳相连接,磁测点位置即为拖体的位置,海底地形数据由考察船上的水下定位装置(如声纳等)采集获得;
(2)对海底目标区域进行三维网格划分,建立磁参数模型;
根据获得的海底地形数据将海底地形轮廓划分为若干个三维网格,如每个单位网格为一个磁数据的测点。则海底地形可由三维网格的组合近似,划分的网格数越多,对海底地形轮廓的近似越精确。
针对划分后的三维网格建立磁参数的理论模型,模型由数学方程表示,方程可写为:
其中,m、n分别表示三维网格的行列数,gij表示第(i,j)个网格的磁异常,它是模型方程的理论值,aij表示磁异常gij对第(i,j)个单元网格磁参数σij的偏导数。
(3)建立海底磁测数据的目标反演函数,以目标反演函数最小值为拟合标准,修改磁参数;
目标反演函数表示为:
式中,gij表示第(i,j)个网格上的磁异常的观测值,目标反演函数实际上反映了磁异常的理论值与观测值之间的误差。首先根据经验给出每个单元网格的磁参数初始值,计算磁异常初始值,然后对磁参数进行不断修改,修改量一般采用广义逆法求解,直到目标反演函数的值最小,表示磁异常的理论值与观测值之间的误差最小,此时的磁参数模型方程达到反演理论拟合标准。
(4)由反演得到的磁参数计算目标区域的磁异常。
本发明方法不需对海底测线进行向上延拓、高截滤波或将起伏面上的观测数据转化到水平面上,避免了引起分辨率降低的中间过程,利用三维空间域带地形直接反演法对磁测数据反演,不仅提高了分辨率,而且对测点位置的起伏量没有限制,适用于大起伏测线的海底磁测。
附图说明
图1为几种磁测方式所得磁测线的特性比较图;
图2为本发明方法的流程图;
图3为将海底地形划分为三维单元网格示意图。
具体实施方式
以下结合附图,详细说明本发明方法的具体实施过程。
如图2所示,本发明一种海底大起伏测线磁测数据的处理方法包括:
(1)获得海底目标区域的磁测数据、磁测点位置数据和海底地形数据;
在海底磁测中,磁测数据由拖体内的磁力仪采集获得,拖体被放置在深海中,与考察船之间用缆绳相连接,磁测点位置即为拖体的位置,由于大海洋流的随机波动性,使拖体会在缆绳的牵动下上下起伏,而使得拖体内的磁力仪所获得的磁测线波动较为剧烈,如图1所示。磁测数据构成了海底磁测大起伏测线,它是一个较复杂的曲面。海底地形数据由考察船上的水下定位装置(如声纳等)采集获得。
(2)对海底目标区域进行三维网格划分,建立磁参数初始模型;
根据海底地形数据将海底地形轮廓划分为若干个规则的长方体网格,如图3所示,长方体网格顶面的长度和宽度根据研究目标区域范围及研究要求等具体情况确定为已知量。若海底地形起伏较大,可将网格顶面的长度和宽度设置为较小的值,这样对目标区域划分为更多的网格,提高对复杂地形的拟合精度。若海底地形较平坦,则可将网格顶面的长度和宽度设置为较大的值。长方体顶面的z轴坐标取此处对应的海底深度,网格高度通过网格顶面的z轴坐标与底面埋深之差确定,而长方体网格的底面埋深根据地质、地球物理及其它资料确定。各单元网格中心作为一个测点位置则海底地形可由长方体单元网格的顶面近似表示,划分的网格数越多,对海底地形轮廓的近似越精确。
根据位场理论,海底磁测大起伏测线构成的任意复杂曲面上的磁异常可由各单元网格的测点磁异常叠加得到,海底磁测线曲面上的磁异常可用以下的模型方程表示:
其中,m、n表示所划分的网格有m行n列,σij表示第(i,j)个网格的磁参数,如磁化强度或磁化率等,gij表示第(i,j)个测点的磁异常,它是模型方程的理论值,aij表示磁异常gij对第(i,j)个单元网格磁参数σij的偏导数。式(1)可写成矩阵形式:
G=A·B (2)
其中,G表示在每个测点由式(1)得到的磁异常模型方程理论值gij组成的磁异常矩阵,A是磁异常矩阵G对第(i,j)个单元网格磁参数的偏导数aij组成的偏导数矩阵,B为每个网格磁参数σij组成的磁参数矩阵,分别有:
GT=[g11 g12 … g1n g21 g22 … g2n … gm1 gm2 … gmn]
(3)
BT=[σ11 σ12 … σ1n σ21 σ22 … σ2n … σm1 σm2 … σmn]
其中,T表示转置。
(3)建立海底磁测数据的目标反演函数,以目标反演函数最小值为拟合标准,不断修改磁参数;
反演目标函数F可写为:
式中,G′表示在每个测点观测得到的磁异常g′ij组成的矩阵。按照反演理论,在大起伏任意观测曲面上的3D磁异常反演就是确定各单元网格的磁参数,由每个单元网格的磁参数经过叠加求和所得到的磁异常与观测得到的磁异常数据相吻合,使理论计算结果与观测结果之间的拟合差最小,使得反演目标函数F的值最小。
由式(3)的极小值为拟合标准可得反演方程如下:
ΔG=A·ΔB (5)
其中,ΔG=G′-G,ΔB为各单元网格磁参数的修改量Δσij组成的矩阵:
ΔBT=[Δσ11 Δσ12 … Δσ1n Δσ21 Δσ22 … Δσ2n … Δσm1 Δσm2 … Δσmn](6)
首先,根据经验给出磁参数矩阵B的初始猜测值B0,由广义逆法求解式(3)可得第一次磁参数模型修改量ΔB0,由B1=B0+ΔB0得到第一次迭代修改后的磁参数矩阵。然后将B1作为新的初始模型参数,同样由广义逆法得到第二次修改量ΔB1。重复以上过程,直到第k次模型参数修改量ΔBk-1,然后按照以下迭代公式对磁参数矩阵进行修改,直到反演目标函数F的值最小。
Bk=Bk-1+ΔBk-1(k=1,2,3...)(7)
(4)由反演后得到的磁参数模型计算某一高度水平面或曲面上的磁异常。
由本发明方法计算磁测数据的过程无需经过延拓或曲化平的操作,但是由本发明方法计算得到最终的磁测数据模型方程后,同样可以将磁测数据演化到任意延拓面(延拓面一般为某一水平面,如图1所示),只需在模型方程中输入延拓面的位置坐标,可以很方便的计算出延拓面的磁响应值。
Claims (3)
1.一种海底大起伏测线磁测数据的处理方法,其特征在于,包括:
(1)获得海底目标区域的磁测数据、磁测点位置数据和海底地形数据;
(2)对海底目标区域进行三维网格划分,建立磁参数模型;
(3)建立海底磁测数据的目标反演函数,以目标反演函数最小值为拟合标准不断修改磁参数;
(4)由反演得到的磁参数计算目标区域的磁异常。
2.根据权利要求1所述的海底大起伏测线磁测数据的处理方法,其特征在于,所述的磁参数模型为:
其中,m、n分别表示网格的行、列数,gij表示第(i,j)个网格的磁异常,σij表示第(i,j)个网格上的磁参数,aij表示磁异常gij对磁参数σij的偏导数。
3.根据权利要求2所述的海底大起伏测线磁测数据的处理方法,其特征在于,所述的目标反演函数为:
其中,g′ij表示第(i,j)个网格磁异常的观测值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102380990A CN101937106B (zh) | 2010-07-27 | 2010-07-27 | 一种海底大起伏测线磁测数据的处理方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010102380990A CN101937106B (zh) | 2010-07-27 | 2010-07-27 | 一种海底大起伏测线磁测数据的处理方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN101937106A true CN101937106A (zh) | 2011-01-05 |
CN101937106B CN101937106B (zh) | 2012-07-25 |
Family
ID=43390532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010102380990A Expired - Fee Related CN101937106B (zh) | 2010-07-27 | 2010-07-27 | 一种海底大起伏测线磁测数据的处理方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101937106B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316972A (zh) * | 2014-10-16 | 2015-01-28 | 中国海洋石油总公司 | 一种磁源重力视强度反演成像方法 |
CN106291723A (zh) * | 2016-07-25 | 2017-01-04 | 中国石油大学(北京) | 基于双参数正则化的核磁共振回波数据反演方法及装置 |
CN107797148A (zh) * | 2017-10-18 | 2018-03-13 | 中国国土资源航空物探遥感中心 | 一种基于三维地质建模的航磁异常场分离方法及*** |
CN107942399A (zh) * | 2017-11-23 | 2018-04-20 | 桂林理工大学 | 一种大距离位场向上延拓计算方法 |
CN113189559A (zh) * | 2021-05-10 | 2021-07-30 | 中国人民解放军海军潜艇学院 | 一种星载成像高度计遥感数据海底地形反演方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998027444A2 (en) * | 1996-12-19 | 1998-06-25 | Schlumberger Holdings Limited | Performing geoscience interpretation with simulated data |
US6603313B1 (en) * | 1999-09-15 | 2003-08-05 | Exxonmobil Upstream Research Company | Remote reservoir resistivity mapping |
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN101382070A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种油藏注采动态监测的电磁方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN101676745A (zh) * | 2008-09-18 | 2010-03-24 | 核工业北京地质研究院 | 一种dct域内的重磁场分离方法 |
-
2010
- 2010-07-27 CN CN2010102380990A patent/CN101937106B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO1998027444A2 (en) * | 1996-12-19 | 1998-06-25 | Schlumberger Holdings Limited | Performing geoscience interpretation with simulated data |
US6603313B1 (en) * | 1999-09-15 | 2003-08-05 | Exxonmobil Upstream Research Company | Remote reservoir resistivity mapping |
CN101004454A (zh) * | 2006-01-20 | 2007-07-25 | 中国石油天然气集团公司 | 目标最小化的三维电磁快速反演方法 |
CN101382070A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种油藏注采动态监测的电磁方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN101676745A (zh) * | 2008-09-18 | 2010-03-24 | 核工业北京地质研究院 | 一种dct域内的重磁场分离方法 |
Non-Patent Citations (2)
Title |
---|
《中国优秀博硕士学位论文全文数据库(硕士)基础科技辑》 20040915 刘百红 高精度磁测资料处理与反演方法研究 全文 1-3 , 第03期 * |
《中国博士学位论文全文数据库 基础科技辑》 20061115 郭志宏 航磁及梯度数据正反演解释方法技术实用化改进及应用 全文 1-3 , 第11期 * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104316972A (zh) * | 2014-10-16 | 2015-01-28 | 中国海洋石油总公司 | 一种磁源重力视强度反演成像方法 |
CN106291723A (zh) * | 2016-07-25 | 2017-01-04 | 中国石油大学(北京) | 基于双参数正则化的核磁共振回波数据反演方法及装置 |
CN106291723B (zh) * | 2016-07-25 | 2018-04-17 | 中国石油大学(北京) | 基于双参数正则化的核磁共振回波数据反演方法及装置 |
CN107797148A (zh) * | 2017-10-18 | 2018-03-13 | 中国国土资源航空物探遥感中心 | 一种基于三维地质建模的航磁异常场分离方法及*** |
CN107797148B (zh) * | 2017-10-18 | 2018-08-21 | 中国国土资源航空物探遥感中心 | 一种基于三维地质建模的航磁异常场分离方法及*** |
CN107942399A (zh) * | 2017-11-23 | 2018-04-20 | 桂林理工大学 | 一种大距离位场向上延拓计算方法 |
CN113189559A (zh) * | 2021-05-10 | 2021-07-30 | 中国人民解放军海军潜艇学院 | 一种星载成像高度计遥感数据海底地形反演方法 |
Also Published As
Publication number | Publication date |
---|---|
CN101937106B (zh) | 2012-07-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Blondel | The handbook of sidescan sonar | |
Lehmann et al. | Semiautomated georadar data acquisition in three dimensions | |
Searle | Mid-ocean ridges | |
US7113868B2 (en) | Method and system for processing geophysical survey data | |
CN103345759B (zh) | 一种海底大型复杂沙波地貌的精确探测方法 | |
CN114137624B (zh) | 一种基于卫星高度计反演海底地形的方法和*** | |
CN106405589A (zh) | 一种全球电离层格网模型的确定方法及装置 | |
CN103389077B (zh) | 一种基于mbes的海底沙波地貌运动探测方法 | |
CN108896040B (zh) | 天空海一体化水下潜器惯性/重力组合导航方法和*** | |
CN101937106B (zh) | 一种海底大起伏测线磁测数据的处理方法 | |
CN106443776B (zh) | 一种基于时间切片法的海底地震仪重定位方法 | |
CN103926625A (zh) | 一种利用地磁总场对水下磁目标高精度远距离的定位方法 | |
CN113189559B (zh) | 一种星载成像高度计遥感数据海底地形反演方法 | |
CN104268848A (zh) | 一种海洋内波波速监测的方法 | |
CN111220146B (zh) | 一种基于高斯过程回归学习的水下地形匹配定位方法 | |
CN102176012B (zh) | 一种二维浅海水下地形合成孔径雷达影像仿真的方法 | |
CN103837140B (zh) | 一种初始弧后盆地大陆架划界的测绘方法 | |
Muliati et al. | Validation of Altimeter Significant Wave Height Using Wave Gauge Measurement in Pacitan Coastal Waters, East Java, Indonesia | |
Asakawa et al. | Development of vertical cable seismic (VCS) system for seafloor massive sulfide (SMS) | |
Zhu et al. | Monitoring of Surface Subsidence of the Mining Area Based on SBAS. | |
Asakawa et al. | Vertical cable seismic surveys for SMS exploration in Izena Cauldron, Okinawa-Trough | |
Melvold et al. | Subglacial topography of Jutulstraumen outlet glacier, East Antarctica, mapped from ground-penetrating radar, optical and interferometric synthetic aperture radar satellite data | |
Gizlenmistir | Production of airborne laser scanner skilled advanced unmanned air vehicle and the potential of preliminary data | |
Liu et al. | Integrated navigation system based on correlation between gravity gradient and terrain | |
Mourad et al. | Interaction of marine geodesy, satellite technology and ocean physics |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
C17 | Cessation of patent right | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20120725 Termination date: 20130727 |