CN106291725A - 一种快速反演地下地质体空间位置的方法 - Google Patents
一种快速反演地下地质体空间位置的方法 Download PDFInfo
- Publication number
- CN106291725A CN106291725A CN201510241172.2A CN201510241172A CN106291725A CN 106291725 A CN106291725 A CN 106291725A CN 201510241172 A CN201510241172 A CN 201510241172A CN 106291725 A CN106291725 A CN 106291725A
- Authority
- CN
- China
- Prior art keywords
- grid
- euler
- value
- magnetic
- locus
- 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
Abstract
本发明属于地球物理重磁数据处理领域,具体涉及一种快速反演地下地质体空间位置的方法。该方法具体包括以下步骤:步骤一、获得磁测数据;步骤二、得到测点异常△T值;步骤三、对磁异常△T值进行网格化处理;步骤四、得到化极磁异常△T值网格和X、Y、Z三个方向导数网格;步骤五、对网格进行分割;步骤六、得到小网格内的欧拉解;步骤七、小网格欧拉解集;步骤八、得出最后解译成果图件。利用本发明技术方案后,对地下地质体进行反演推测能够反映实际地质特征,有效地评价了地下地质体的分布规律。从而达到快速、精确地探测地下矿产资源,对地质找矿工作具有重要的实际意义。
Description
技术领域
本发明属于地球物理重磁数据处理领域,具体涉及一种快速反演地下地质体空间位置的方法。
背景技术
在地球物理重磁勘探领域,位场资料处理和解释中如何准确推断地下地质体形态、空间位置和分布规律往往是地球物理勘查的难点。常规的重磁异常自动反演技术如欧拉反褶积法可以自动或半自动化确定场源位置,但构造指数的选取对最终反演结果影响较大,且在地质构造复杂区域,选取固定的构造指数来进行反演计算不能准确地计算和模拟各个子区域场源形态与空间位置,导致直接影响解译结果精度。
因此,亟需研制一种快速反演地下地质体埋深的方法,以避免不同构造地区按同一构造指数进行运算的影响,准确对地下地质信息进行有效的探测,同时反映出不同区域地质体形态和空间位置的变化。
发明内容
本发明要解决的技术问题是提供一种快速反演地下地质体空间位置的方法,能够有效、真实地探测地下地质体形态和空间位置。
为解决上述技术问题,本发明一种快速反演地下地质体空间位置的方法,依次包括以下步骤:
步骤一、在勘查区开展高精度磁法测量,获得磁测数据;
步骤二、对磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正,得到测点异常△T值;
步骤三、将步骤二所得到的磁异常△T值进行数据平滑处理,消除高频干扰,得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理;
步骤四、对步骤三中的处理后的磁异常△T值进行化极处理得到化极磁异常△T值网格,对化极磁异常△T值网格进行X、Y、Z三个方向求取导数并得到X、Y、Z三个方向导数网格;
步骤五、步骤四中化极磁异常△T值网格和步骤四中X、Y、Z三方向导数网格进行分割,把以上图件剖分成小网格;
步骤六、对每个小网格构造指数进行赋值,并且利用欧拉三维反褶积反演计算,得出所有小网格内的欧拉解集,设置解集误差容许范围指数为15%并筛选出符合条件的小网格内的欧拉解;
步骤七、提取出小网格欧拉解集,把所有欧拉解的坐标、深度数据编号并放在同一数据库;
步骤八、对所有解集数据主要是解的坐标和深度进行成图处理,对步骤六所述小网格求出的欧拉解的小地质单元汇聚起来并且在三维空间展示则可形成模拟的地质体,最终形成欧拉解集三维立体图,返回步骤六对小网格单元构造指数进行重新赋值,重复步骤七、步骤八,得出最后解译成果图件。
所述的步骤一中磁测总精度应高于2nT。
所述的步骤五中,所述的小网格个数为n个,大小为10m×10m。
所述的步骤六中,对每个小网格构造指数进行赋值,对为岩体接触带的网格赋值为0,对岩墙或岩席的网格赋值为1,对水平管状体的网格赋值2,对近似球体形状的岩体的网格赋值为3。
欧拉搜索窗口大小为5m×5m。
所述的欧拉三维反褶积运算采用下式进行计算:
式中N为构造指数值,B为背景磁场值,f为地下地质体场源(x0,y0,z0)在观测点(x,y,z)处的磁场函数表达式。
本发明的有益技术效果在于:利用本发明技术方案后,对地下地质体进行反演推测能够反映实际地质特征,有效地评价了地下地质体的分布规律。从而达到快速、精确地探测地下矿产资源,对地质找矿工作具有重要的实际意义。
具体实施方式
下面结合实施例对本发明作进一步详细说明。
本发明一种快速反演地下地质体空间位置的方法,依次包括以下步骤:
步骤一、在勘查区开展高精度磁法测量,获得磁测数据,为了达到更好的效果,磁测总精度应高于2nT,测区选在地形平坦的区域;
步骤二、对磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正,得到测点异常△T值,在对磁数据进行日变校正时,需查找地球基本磁场参数,通过利用http://www.ngdc.noaa.gov/seg/geomag/jsp/IGRF.jsp网站提供的球谐模型计算出该勘查区背景磁场强度值、磁倾角和磁偏角;
步骤三、将步骤二所得到的磁异常△T值进行数据平滑处理,消除高频干扰,得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理,网格化插值方法优选最小曲率法,对磁异常△T值进行多点数据平滑方法优选七点滑动平均方法进行计算;
步骤四、对步骤三中的处理后的磁异常△T值进行化极处理得到化极磁异常△T值网格,对化极磁异常△T值网格进行X、Y、Z三个方向求取导数并得到X、Y、Z三个方向导数网格;
步骤五、利用Matlab、c语言等程序对步骤四中化极磁异常△T值网格和步骤四中X、Y、Z三方向导数网格进行分割,把以上图件剖分成n个10m×10m的小网格,每一套网格包括10m×10m的化极磁异常△T网格和10m×10m的X、Y、Z三方向导数网格;
步骤六、根据实际地质及已有钻探资料对每个小网格构造指数进行赋值,对为岩体接触带的网格赋值为0,对岩墙或岩席的网格赋值为1,对水平管状体的网格赋值2,对近似球体形状的岩体的网格赋值为3,欧拉搜索窗口大小为5m×5m,并且利用欧拉三维反褶积反演计算,得出所有n个小网格内的欧拉解集,设置解集误差容许范围指数为15%并筛选出符合条件的n个小网格内的欧拉解;
所述的欧拉三维反褶积运算采用下式进行计算:
式中N为构造指数值,B为背景磁场值,f为地下地质体场源(x0,y0,z0)在观测点(x,y,z)处的磁场函数表达式;
步骤七、提取出n个小网格欧拉解集,包括每个解集坐标、深度、搜索窗口大小、窗口中心到两边的距离,把所有欧拉解的坐标、深度数据编号并放在同一数据库;
步骤八、利用三维地学成图软件Surpac或Discover对所有解集数据主要是解的坐标和深度进行成图处理,每一个解实际上代表模拟地下地质体的一个地质单元,对步骤六所述n个小网格求出的欧拉解有几十万个,几十万个小地质单元汇聚起来并且在三维空间展示则可形成模拟的地质体,最终形成欧拉解集三维立体图,将最终成果图件与实际地质资料如钻探、测井等资料进行对比,返回步骤六对n个小网格单元构造指数进行重新赋值,重复步骤七、步骤八,得出最后解译成果图件。
Claims (6)
1.一种快速反演地下地质体空间位置的方法,依次包括以下步骤:
步骤一、在勘查区开展高精度磁法测量,获得磁测数据;
步骤二、对磁测数据进行日变校正、正常场梯度改正、高度改正及基点改正,得到测点异常△T值;
步骤三、将步骤二所得到的磁异常△T值进行数据平滑处理,消除高频干扰,得到光滑的磁异常曲线并对所有光滑处理后的磁异常△T值进行网格化处理;
步骤四、对步骤三中的处理后的磁异常△T值进行化极处理得到化极磁异常△T值网格,对化极磁异常△T值网格进行X、Y、Z三个方向求取导数并得到X、Y、Z三个方向导数网格;
步骤五、步骤四中化极磁异常△T值网格和步骤四中X、Y、Z三方向导数网格进行分割,把以上图件剖分成小网格;
步骤六、对每个小网格构造指数进行赋值,并且利用欧拉三维反褶积反演计算,得出所有小网格内的欧拉解集,设置解集误差容许范围指数为15%并筛选出符合条件的小网格内的欧拉解;
步骤七、提取出小网格欧拉解集,把所有欧拉解的坐标、深度数据编号并放在同一数据库;
步骤八、对所有解集数据主要是解的坐标和深度进行成图处理,对步骤六所述小网格求出的欧拉解的小地质单元汇聚起来并且在三维空间展示则可形成模拟的地质体,最终形成欧拉解集三维立体图,返回步骤六对小网格单元构造指数进行重新赋值,重复步骤七、步骤八,得出最后解译成果图件。
2.根据权利要求1所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的步骤一中磁测总精度应高于2nT。
3.根据权利要求2所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的步骤五中,所述的小网格个数为n个,大小为10m×10m。
4.根据权利要求3所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的步骤六中,对每个小网格构造指数进行赋值,对为岩体接触带的网格赋值为0,对岩墙或岩席的网格赋值为1,对水平管状体的网格赋值2,对近似球体形状的岩体的网格赋值为3。
5.根据权利要求4所述的一种快速反演地下地质体空间位置的方法,其特征在于:欧拉搜索窗口大小为5m×5m。
6.根据权利要求5所述的一种快速反演地下地质体空间位置的方法,其特征在于:所述的欧拉三维反褶积运算采用下式进行计算:
式中N为构造指数值,B为背景磁场值,f为地下地质体场源(x0,y0,z0)在观测点(x,y,z)处的磁场函数表达式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510241172.2A CN106291725B (zh) | 2015-05-13 | 2015-05-13 | 一种快速反演地下地质体空间位置的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510241172.2A CN106291725B (zh) | 2015-05-13 | 2015-05-13 | 一种快速反演地下地质体空间位置的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291725A true CN106291725A (zh) | 2017-01-04 |
CN106291725B CN106291725B (zh) | 2018-11-30 |
Family
ID=57631447
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510241172.2A Active CN106291725B (zh) | 2015-05-13 | 2015-05-13 | 一种快速反演地下地质体空间位置的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291725B (zh) |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108492204A (zh) * | 2018-03-07 | 2018-09-04 | 吉林大学 | 一种高精度空地井不同空间重磁数据变换方法 |
CN108508479A (zh) * | 2018-03-07 | 2018-09-07 | 吉林大学 | 一种空地井立体重磁数据协同目标位置反演方法 |
CN108732622A (zh) * | 2018-05-18 | 2018-11-02 | 吉林大学 | 一种不同高度数据融合联合反演地质体几何形态的方法 |
CN109307847A (zh) * | 2018-11-06 | 2019-02-05 | 哈尔滨工程大学 | 磁体的二阶欧拉反演法及一种标量磁强计阵列 |
CN109752767A (zh) * | 2018-12-12 | 2019-05-14 | 核工业北京地质研究院 | 一种热液型铀矿床地球物理场源边界识别方法 |
CN111708088A (zh) * | 2020-06-28 | 2020-09-25 | 中国矿业大学 | 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统 |
CN112379448A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 磁异常磁源深度反演方法及装置 |
CN112462442A (zh) * | 2020-11-30 | 2021-03-09 | 山东大学 | 重磁位场场源位置估计方法、系统、介质及电子设备 |
CN112668146A (zh) * | 2020-12-03 | 2021-04-16 | 重庆科技学院 | 一种基于欧拉反褶积法实用性改进的场源位置估算方法 |
CN116047617A (zh) * | 2023-03-10 | 2023-05-02 | 中国地质科学院地球物理地球化学勘查研究所 | 一种井间地质特征识别方法及装置 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2369566A1 (en) * | 1999-04-02 | 2000-10-12 | Alan Royce Huffman | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
EA201070952A1 (ru) * | 2008-02-11 | 2011-02-28 | Шеврон Ю.Эс.Эй. Инк. | Способ определения геологической информации, относящейся к интересующему объему недр |
US8055446B2 (en) * | 2006-10-11 | 2011-11-08 | Byerly Kent A | Methods of processing magnetotelluric signals |
CN104570070A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油天然气集团公司 | 一种建立二维近地表地质模型的方法及设备 |
CN104597494A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 地震地层体分析方法及装置 |
-
2015
- 2015-05-13 CN CN201510241172.2A patent/CN106291725B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2369566A1 (en) * | 1999-04-02 | 2000-10-12 | Alan Royce Huffman | A method for gravity and magnetic data inversion using vector and tensor data with seismic imaging and geopressure prediction for oil, gas and mineral exploration and production |
US6615139B1 (en) * | 2002-03-28 | 2003-09-02 | Council Of Scientific & Industrial Research | Digitally implemented method for automatic optimization of gravity fields obtained from three-dimensional density interfaces using depth dependent density |
US8055446B2 (en) * | 2006-10-11 | 2011-11-08 | Byerly Kent A | Methods of processing magnetotelluric signals |
EA201070952A1 (ru) * | 2008-02-11 | 2011-02-28 | Шеврон Ю.Эс.Эй. Инк. | Способ определения геологической информации, относящейся к интересующему объему недр |
CN104570070A (zh) * | 2013-10-12 | 2015-04-29 | 中国石油天然气集团公司 | 一种建立二维近地表地质模型的方法及设备 |
CN104597494A (zh) * | 2013-10-31 | 2015-05-06 | 中国石油天然气集团公司 | 地震地层体分析方法及装置 |
Non-Patent Citations (4)
Title |
---|
G. LEBLANC ET AL.: "Denoising of Aeromagnetic Data via the Wavelet Transform", 《GEOPHYSICS》 * |
VALERIA C. F. BARBOSA ET AL.: "Stability analysis and improvement of structural index estimation in Euler-deconvolution", 《GEOPHYSICS》 * |
喻翔等: "河南栾川罗村钼多金属矿区的磁法、激电异常特征", 《东华理工大学学报(自然科学版)》 * |
张濡亮等: "音频大地电磁测深法在尼日尔阿泽里克铀成矿区的应用研究", 《世界核地质科学》 * |
Cited By (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108492204B (zh) * | 2018-03-07 | 2021-06-01 | 吉林大学 | 一种高精度空地井不同空间重磁数据变换方法 |
CN108508479A (zh) * | 2018-03-07 | 2018-09-07 | 吉林大学 | 一种空地井立体重磁数据协同目标位置反演方法 |
CN108492204A (zh) * | 2018-03-07 | 2018-09-04 | 吉林大学 | 一种高精度空地井不同空间重磁数据变换方法 |
CN108732622A (zh) * | 2018-05-18 | 2018-11-02 | 吉林大学 | 一种不同高度数据融合联合反演地质体几何形态的方法 |
CN109307847A (zh) * | 2018-11-06 | 2019-02-05 | 哈尔滨工程大学 | 磁体的二阶欧拉反演法及一种标量磁强计阵列 |
CN109752767A (zh) * | 2018-12-12 | 2019-05-14 | 核工业北京地质研究院 | 一种热液型铀矿床地球物理场源边界识别方法 |
CN111708088A (zh) * | 2020-06-28 | 2020-09-25 | 中国矿业大学 | 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统 |
CN111708088B (zh) * | 2020-06-28 | 2022-04-15 | 中国矿业大学 | 一种基于磁梯度张量的瞬变电磁实时动态超前探测方法及系统 |
CN112379448A (zh) * | 2020-10-30 | 2021-02-19 | 中国石油天然气集团有限公司 | 磁异常磁源深度反演方法及装置 |
CN112462442A (zh) * | 2020-11-30 | 2021-03-09 | 山东大学 | 重磁位场场源位置估计方法、系统、介质及电子设备 |
CN112462442B (zh) * | 2020-11-30 | 2022-04-08 | 山东大学 | 重磁位场场源位置估计方法、系统、介质及电子设备 |
CN112668146A (zh) * | 2020-12-03 | 2021-04-16 | 重庆科技学院 | 一种基于欧拉反褶积法实用性改进的场源位置估算方法 |
CN116047617A (zh) * | 2023-03-10 | 2023-05-02 | 中国地质科学院地球物理地球化学勘查研究所 | 一种井间地质特征识别方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106291725B (zh) | 2018-11-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106291725A (zh) | 一种快速反演地下地质体空间位置的方法 | |
CN105589108B (zh) | 基于不同约束条件的瞬变电磁快速三维反演方法 | |
CN104133245B (zh) | 一种地震资料的静校正方法及系统 | |
CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
CN106650192B (zh) | 一种火山岩型铀矿床磁性界面反演方法 | |
CN104537714A (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
Hannemann et al. | Three‐dimensional shallow structure from high‐frequency ambient noise tomography: New results for the Mygdonia basin‐Euroseistest area, northern Greece | |
CN106483559A (zh) | 一种地下速度模型的构建方法 | |
CN105277978A (zh) | 一种确定近地表速度模型的方法及装置 | |
CN102116869A (zh) | 高精度叠前域最小二乘偏移地震成像技术 | |
CN105242328B (zh) | 古热岩石圈厚度的确定方法及装置 | |
CN102901985A (zh) | 一种适用于起伏地表的深度域层速度修正方法 | |
CN104181593B (zh) | 一种三维无射线追踪回折波层析成像方法及装置 | |
CN103728667A (zh) | 一种视三维高密度电法的浅表层地质结构建模方法 | |
CN102877828A (zh) | 一种三维多井联合井地ct成像方法 | |
Hu et al. | Ray-illumination compensation for adjoint-state first-arrival traveltime tomography | |
Tang et al. | Topographic effects on long offset transient electromagnetic response | |
CN109902315A (zh) | 一种圈定隐伏花岗岩岩体深部边界的方法 | |
CN105573963B (zh) | 一种电离层水平不均匀结构重构方法 | |
CN103217715B (zh) | 多尺度规则网格层析反演静校正方法 | |
CN114488327B (zh) | 基于地面基点的水平磁场与井中垂直磁场联合测量方法 | |
CN107918147B (zh) | 绕射波成像方法和装置 | |
CN102841374A (zh) | 基于扫描面正演的伪三维快速微地震正演方法 | |
CN113608262B (zh) | 利用平动分量计算旋转分量的地震数据处理方法及装置 | |
CN107765341A (zh) | 一种确定地层剩余密度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |