CN112464521B - 走滑断裂确定方法及装置 - Google Patents
走滑断裂确定方法及装置 Download PDFInfo
- Publication number
- CN112464521B CN112464521B CN202011172082.XA CN202011172082A CN112464521B CN 112464521 B CN112464521 B CN 112464521B CN 202011172082 A CN202011172082 A CN 202011172082A CN 112464521 B CN112464521 B CN 112464521B
- Authority
- CN
- China
- Prior art keywords
- gravity
- data
- derivative
- development direction
- abnormal information
- 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 46
- 230000005484 gravity Effects 0.000 claims abstract description 179
- 238000011161 development Methods 0.000 claims abstract description 116
- 230000002159 abnormal effect Effects 0.000 claims abstract description 94
- 238000001914 filtration Methods 0.000 claims abstract description 58
- 238000012545 processing Methods 0.000 claims abstract description 49
- 230000005856 abnormality Effects 0.000 claims description 17
- 238000004590 computer program Methods 0.000 claims description 16
- 238000003860 storage Methods 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 4
- 230000003247 decreasing effect Effects 0.000 claims 2
- 238000010586 diagram Methods 0.000 description 16
- 238000005259 measurement Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007423 decrease Effects 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 238000006467 substitution reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000005336 cracking Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供了一种走滑断裂确定方法及装置,涉及地球物理勘探技术领域,该方法包括:获取布格重力异常网格化数据;根据布格重力异常网格化数据生成重力导数异常信息;获取第一发育方向数据和第二发育方向数据;第一发育方向数据和第二发育方向数据根据重力导数异常信息得到;根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据;根据重力导数异常信息和滤波处理数据生成叠加图,以利用叠加图确定走滑断裂数据。本发明对走滑断裂的重力异常进行信息提取,进而得到走滑断裂数据,解决了重力资料解释走滑断裂的多解性难题,可以更准确地确定走滑断裂的发育特征和断裂位置。
Description
技术领域
本发明涉及地球物理勘探技术领域,尤其是涉及一种走滑断裂确定方法及装置。
背景技术
重力勘探是油气及地质矿产勘探中的一种重要方法,在新区勘探、目标勘探、潜山、裂谷、断陷盆地、断裂特征研究中具有重要作用。现有重力资料处理解释对重力异常梯级带展示的断裂的处理应用极为广泛,重力梯级带只能反映出垂向存在断距的正断层、逆断层或逆掩断层,无法清晰刻画走滑断裂,在利用重力异常解释走滑断裂时则往往具有较大的人为因素,缺少必要的处理解释方法,致使其存在多解性、解释位置不准确等问题。
发明内容
本发明提供了一种走滑断裂确定方法及装置,可以降低重力资料解释的多解性和解释工作的难度。
第一方面,本发明实施例提供了一种走滑断裂确定方法,该方法包括:获取布格重力异常网格化数据;根据所述布格重力异常网格化数据生成重力导数异常信息;获取第一发育方向数据和第二发育方向数据;所述第一发育方向数据和所述第二发育方向数据根据所述重力导数异常信息得到;根据所述第一发育方向数据、所述第二发育方向数据和所述重力导数异常信息生成滤波处理数据;根据所述重力导数异常信息和所述滤波处理数据生成叠加图,以利用所述叠加图确定走滑断裂数据。
第二方面,本发明实施例还提供一种走滑断裂确定装置,该装置包括:第一获取模块,用于获取布格重力异常网格化数据;计算模块,用于根据所述布格重力异常网格化数据生成重力导数异常信息;第二获取模块,用于获取第一发育方向数据和第二发育方向数据;所述第一发育方向数据和所述第二发育方向数据根据所述重力导数异常信息得到;生成模块,用于根据所述第一发育方向数据、所述第二发育方向数据和所述重力导数异常信息生成滤波处理数据;确定模块,用于根据所述重力导数异常信息和所述滤波处理数据生成叠加图,以利用所述叠加图确定走滑断裂数据。
第三方面,本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述走滑断裂确定方法。
第四方面,本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述走滑断裂确定方法的计算机程序。
本发明实施例带来了以下有益效果:本发明实施例提供了一种走滑断裂确定方案,该方案首先获取布格重力异常网格化数据;之后根据布格重力异常网格化数据生成重力导数异常信息;再获取第一发育方向数据和第二发育方向数据;第一发育方向数据和第二发育方向数据根据重力导数异常信息得到;根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据;根据重力导数异常信息和滤波处理数据生成叠加图,以利用叠加图确定走滑断裂数据。本发明实施例通过对布格重力异常网格化数据进行处理转换,得到滤波处理数据,实现对走滑断裂的重力异常进行信息提取,进而得到走滑断裂数据,解决了重力资料解释走滑断裂的多解性难题,可以更准确地确定走滑断裂的发育特征和断裂位置。
本发明的其他特征和优点将在随后的说明书中阐述,并且,部分地从说明书中变得显而易见,或者通过实施本发明而了解。本发明的目的和其他优点在说明书、权利要求书以及附图中所特别指出的结构来实现和获得。
为使本发明的上述目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附附图,作详细说明如下。
附图说明
为了更清楚地说明本发明具体实施方式或现有技术中的技术方案,下面将对具体实施方式或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施方式,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的走滑断裂确定方法流程图;
图2为本发明实施例提供的走滑断裂确定方法实施步骤示意图;
图3为本发明实施例提供的某区布格重力异常图;
图4为本发明实施例提供的走滑断裂确定方法处理效果图;
图5为本发明实施例提供的走滑断裂确定装置结构框图;
图6为本发明实施例提供的生成模块结构框图;
图7为本发明实施例提供的计算机设备结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合附图对本发明的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
目前,走滑断裂的重力资料解释往往存在多解性,人为因素大,缺少实用性的技术方法,断裂***特征研究是重力资料解释中的重点和难题。
基于此,本发明实施例提供的一种走滑断裂确定方法及装置,是对重力资料处理解释技术的方法创新,为重力资料处理解释提供一种确定走滑断裂的重力处理方法,从而提高重力数据对走滑断裂位置、展布的定量描述能力,提供技术手段和资料依据,降低重力资料解释的多解性和解释工作的难度。
为便于对本实施例进行理解,首先对本发明实施例所公开的一种走滑断裂确定方法进行详细介绍。
本发明实施例提供了一种走滑断裂确定方法,参见图1所示的一种走滑断裂确定方法流程图,该方法包括以下步骤:
步骤S102,获取布格重力异常网格化数据。
在本发明实施例中,布格重力异常网格化数据中包括多个网格数据,每个网格点可以作为一个测点。布格重力异常网格化数据可以利用重力资料处理得到。
步骤S104,根据布格重力异常网格化数据生成重力导数异常信息。
在本发明实施例中,由布格重力异常网格化数据通过计算,求取重力导数异常信息。需要保证对网格中的每个网格点按照同样的方法进行计算,得到重力导数异常信息。
步骤S106,获取第一发育方向数据和第二发育方向数据。
在本发明实施例中,第一发育方向数据和第二发育方向数据可以由人工根据重力导数异常信息得到。第一发育方向数据和第二发育方向数据用于确定走滑断裂可能发育的方向,并且第一发育方向的垂直方向为第二发育方向。根据重力导数异常信息将在平面上呈现的走向作为第一发育方向,按第一发育方向走向的垂直方向确定第二发育方向。
步骤S108,根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据。
在本发明实施例中,计算目标测点的滤波处理数据,可以使用该点的重力导数异常信息,该点的第一发育方向数据和第二发育方向数据进行计算,对每个目标测点进行同样的处理。
步骤S110,根据重力导数异常信息和滤波处理数据生成叠加图,以利用叠加图确定走滑断裂数据。
在本发明实施例中,在得到重力导数异常信息之后,根据该重力导数异常信息生成重力导数异常图,将滤波处理数据叠加到重力导数异常图上,可以采用半透明彩色图上下图层叠加、不同颜色线条等值线叠加,也可以是,位于上层的重力导数异常图采用不充填颜色的等值线显示,位于下层的走滑断裂信息采用无等值线线条的彩色图,具体叠加步骤可以根据实际需求进行选取。
在得到叠加图后,根据叠加图,在重力导数异常沿走向追踪的扭断或错位处,信息叠加图上若有走滑断裂信息的条带状发育,追踪该条带状走滑断裂信息,以正负异常分界(0值线)为走滑断裂位置,进行连线,确定走滑断裂的存在和发育位置,得到走滑断裂数据。
本发明实施例提供了一种走滑断裂确定方案,该方案首先获取布格重力异常网格化数据;之后根据布格重力异常网格化数据生成重力导数异常信息;再获取第一发育方向数据和第二发育方向数据;第一发育方向数据和第二发育方向数据根据重力导数异常信息得到;根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据;根据重力导数异常信息和滤波处理数据生成叠加图,以利用叠加图确定走滑断裂数据。本发明实施例通过对布格重力异常网格化数据进行处理转换,得到滤波处理数据,实现对走滑断裂的重力异常进行信息提取,进而得到走滑断裂数据,解决了重力资料解释走滑断裂的多解性难题,可以更准确地确定走滑断裂的发育特征和断裂位置。
为了得到更准确的重力导数异常信息,根据布格重力异常网格化数据生成重力导数异常信息,可以按照如下步骤执行:
根据布格重力异常网格化数据生成重力垂直二次导数异常信息或重力垂直一次导数异常信息;将重力垂直二次导数异常信息或重力垂直一次导数异常信息作为重力导数异常信息。
在本发明实施例中,若布格重力异常存在背景趋势场时,求取重力垂直二次导数异常,若布格重力异常背景趋势场弱,则求取重力垂直一次导数异常,求取得到的重力垂直一次导数异常或重力垂直二次导数异常均科作为重力导数异常信息。
需要说明的是,同一网格中所有网格点应均为通过求取重力垂直一次导数异常或求取重力垂直二次导数异常得到重力导数异常信息。
为了提升滤波效果,根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据,可以按照如下步骤执行:
利用第一发育方向数据和重力导数异常信息生成目标测点在第一目标范围内的重力导数异常信息平均值;根据重力导数异常信息和目标测点的重力导数异常信息平均值生成第一滤波处理结果;利用第二发育方向数据、第一滤波结果和重力导数异常信息生成滤波处理结果。
利用第二发育方向数据、第一滤波结果和重力导数异常信息生成滤波处理结果,可以按照如下步骤执行:
获取第二发育方向数据对应的权重数据;利用第二发育方向数据、权重数据和重力导数异常信息生成补偿数据;根据第一滤波结果和补偿数据计算目标测点在第二目标范围内的第二滤波处理结果;将第二滤波处理结果作为滤波处理数据。
在本发明实施例中,第一目标范围和第二目标范围可以根据实际需求进行设置,本发明实施例对此不作具体限定。
在一个实施例中,权重数据沿第二发育方向,权重值先增大后减小,并且目标测点的权重值大于除目标观测点外的网格点的权重值。
在一个实施例中,第二目标范围以目标测点为中心,并且沿第二发育方向包括5至9个网格点;权重值沿第二发育方向依次为:1/16、2/16、3/16、4/16、3/16、2/16和1/16。
在一个实施例中,第一目标范围以目标测点为中心,并且沿第一发育方向包括5至9个网格点。
参见图3所示的某区布格重力异常图和图4所示的走滑断裂确定方法处理效果图,在某地区,按照本发明实施了确定走滑断裂的重力处理,具体方法步骤如下:
1)由布格重力异常网格化数据,若布格重力异常存在背景趋势场时,求取重力垂直二次导数异常,若布格重力异常背景趋势场弱,则求取重力垂直一次导数异常,求取得到的重力垂直一次导数异常或重力垂直二次导数异常均称为重力导数异常,本实施中求取得到了重力垂直一次导数异常;
2)确定走滑断裂可能发育的方向,根据实施例重力导数异常在平面上呈现的走向为南北向(记为Ads),按其走向的垂直方向确定为走滑断裂可能发育的方向为东西向(Azh);
3)走滑断裂信息滤波处理,计算点的走滑断裂信息为测点的重力导数异常-走向Ads上的7点平均值-方向Azh上的补偿异常,走向Ads上的(实施例为南北方向)7点平均值使用测点为中心的7点平均法,计算了“测点的重力导数异常-走向Ads上的7点平均值”之后,对其剩余值求取方向Azh上的(实施例为东西方向)补偿异常,“方向Azh上的补偿异常”使用测点为中心的7点加权平均法,加权系数分别取1/16、2/16、3/16、4/16、3/16、2/16、1/16;
4)信息与异常叠加,将步骤3)得到的走滑断裂信息叠加到步骤1)得到的重力导数异常图上,可以采用半透明彩色图上下图层叠加、不同颜色线条等值线叠加,实施例是采用,位于下层的重力导数异常图采用充填颜色的等值线显示,位于上层的走滑断裂信息采用半透明等值线彩色图;
5)确定走滑断裂,在重力导数异常沿走向追踪的扭断或错位处,信息叠加图上若有走滑断裂信息的条带状发育,追踪该条带状走滑断裂信息,以正负异常分界(0值线)为走滑断裂位置进行连线确定东西方向及近东西方向的走滑断裂的存在和发育位置。
本发明实施例提供了一种走滑断裂确定方法及装置,参见图2所示的走滑断裂确定方法实施步骤示意图,该方法是一种确定走滑断裂的重力资料处理方法,在面积性高精度重力勘探中具有重要实用价值,是现有重力资料解释研究断裂***的重要方法补充,为确定走滑断裂的重力(磁力)资料解释提供了一种新方法。
本发明实施例中还提供了一种走滑断裂确定装置,如下面的实施例所述。由于该装置解决问题的原理与走滑断裂确定方法相似,因此该装置的实施可以参见走滑断裂确定方法的实施,重复之处不再赘述。参见图5所示的一种走滑断裂确定装置结构框图,该装置包括:
第一获取模块51,用于获取布格重力异常网格化数据;计算模块52,用于根据布格重力异常网格化数据生成重力导数异常信息;第二获取模块53,用于获取第一发育方向数据和第二发育方向数据;第一发育方向数据和第二发育方向数据根据重力导数异常信息得到;生成模块54,用于根据第一发育方向数据、第二发育方向数据和重力导数异常信息生成滤波处理数据;确定模块55,用于根据重力导数异常信息和滤波处理数据生成叠加图,以利用叠加图确定走滑断裂数据。
在一个实施例中,计算模块,具体用于:根据布格重力异常网格化数据生成重力垂直二次导数异常信息或重力垂直一次导数异常信息;将重力垂直二次导数异常信息或重力垂直一次导数异常信息作为重力导数异常信息。
在一个实施例中,参见图6所示的生成模块结构框图,生成模块,包括:计算单元61,用于利用第一发育方向数据和重力导数异常信息生成目标测点在第一目标范围内的重力导数异常信息平均值;生成单元62,用于根据重力导数异常信息和目标测点的重力导数异常信息平均值生成第一滤波处理结果;滤波单元63,用于利用第二发育方向数据、第一滤波结果和重力导数异常信息生成滤波处理结果。
在一个实施例中,滤波单元,具体用于:获取所述第二发育方向数据对应的权重数据;利用所述第二发育方向数据、所述权重数据和所述重力导数异常信息生成补偿数据;根据所述第一滤波结果和补偿数据计算目标测点在第二目标范围内的第二滤波处理结果;将所述第二滤波处理结果作为滤波处理数据。
在一个实施例中,权重数据沿第二发育方向,权重值先增大后减小,并且目标测点的权重值大于除目标观测点外的网格点的权重值。
在一个实施例中,第二目标范围以目标测点为中心,并且沿第二发育方向包括5至9个网格点;权重值沿第二发育方向依次为:1/16、2/16、3/16、4/16、3/16、2/16和1/16。
在一个实施例中,第一目标范围以目标测点为中心,并且沿第一发育方向包括5至9个网格点。
本发明实施例还提供一种计算机设备,参见图7所示的计算机设备结构示意框图,该计算机设备包括存储器71、处理器72及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述任一种走滑断裂确定方法的步骤。
所属领域的技术人员可以清楚地了解到,为描述的方便和简洁,上述描述的计算机设备的具体工作过程,可以参考前述方法实施例中的对应过程,在此不再赘述。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述任一种走滑断裂确定方法的计算机程序。
本领域内的技术人员应明白,本发明的实施例可提供为方法、***、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(***)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应说明的是:以上所述实施例,仅为本发明的具体实施方式,用以说明本发明的技术方案,而非对其限制,本发明的保护范围并不局限于此,尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,其依然可以对前述实施例所记载的技术方案进行修改或可轻易想到变化,或者对其中部分技术特征进行等同替换;而这些修改、变化或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应所述以权利要求的保护范围为准。
Claims (6)
1.一种走滑断裂确定方法,其特征在于,包括:
获取布格重力异常网格化数据;
根据所述布格重力异常网格化数据生成重力导数异常信息;
获取第一发育方向数据和第二发育方向数据;所述第一发育方向数据和所述第二发育方向数据根据所述重力导数异常信息得到;根据重力导数异常信息将在平面上呈现的走向作为第一发育方向,按第一发育方向走向的垂直方向确定第二发育方向;
根据所述第一发育方向数据、所述第二发育方向数据和所述重力导数异常信息生成滤波处理数据;
根据所述重力导数异常信息和所述滤波处理数据生成叠加图,以利用所述叠加图确定走滑断裂数据;
根据所述第一发育方向数据、所述第二发育方向数据和所述重力导数异常信息生成滤波处理数据,包括:
利用所述第一发育方向数据和所述重力导数异常信息生成目标测点在第一目标范围内的重力导数异常信息平均值;
根据所述重力导数异常信息和所述目标测点的重力导数异常信息平均值生成第一滤波处理结果;
利用所述第二发育方向数据、所述第一滤波处理结果和所述重力导数异常信息生成滤波处理结果;
利用所述第二发育方向数据、所述第一滤波结果和所述重力导数异常信息生成滤波处理结果,包括:
获取所述第二发育方向数据对应的权重数据;
利用所述第二发育方向数据、所述权重数据和所述重力导数异常信息生成补偿数据;
根据所述第一滤波结果和补偿数据计算目标测点在第二目标范围内的第二滤波处理结果;
将所述第二滤波处理结果作为滤波处理数据;
所述权重数据沿第二发育方向,权重值先增大后减小,并且目标测点的权重值大于除目标观测点外的网格点的权重值;
所述第二目标范围以所述目标测点为中心,并且沿第二发育方向包括5至9个网格点;
所述权重值沿第二发育方向依次为:1/16、2/16、3/16、4/16、3/16、2/16和1/16;
所述第一目标范围以所述目标测点为中心,并且沿第一发育方向包括5至9个网格点。
2.根据权利要求1所述的方法,其特征在于,根据所述布格重力异常网格化数据生成重力导数异常信息,包括:
根据所述布格重力异常网格化数据生成重力垂直二次导数异常信息或重力垂直一次导数异常信息;
将所述重力垂直二次导数异常信息或重力垂直一次导数异常信息作为重力导数异常信息。
3.一种走滑断裂确定装置,其特征在于,包括:
第一获取模块,用于获取布格重力异常网格化数据;
计算模块,用于根据所述布格重力异常网格化数据生成重力导数异常信息;
第二获取模块,用于获取第一发育方向数据和第二发育方向数据;所述第一发育方向数据和所述第二发育方向数据根据所述重力导数异常信息得到;根据重力导数异常信息将在平面上呈现的走向作为第一发育方向,按第一发育方向走向的垂直方向确定第二发育方向;
生成模块,用于根据所述第一发育方向数据、所述第二发育方向数据和所述重力导数异常信息生成滤波处理数据;
确定模块,用于根据所述重力导数异常信息和所述滤波处理数据生成叠加图,以利用所述叠加图确定走滑断裂数据;
所述生成模块,包括:
计算单元,用于利用所述第一发育方向数据和所述重力导数异常信息生成目标测点在第一目标范围内的重力导数异常信息平均值;
生成单元,用于根据所述重力导数异常信息和所述目标测点的重力导数异常信息平均值生成第一滤波处理结果;
滤波单元,用于利用所述第二发育方向数据、所述第一滤波处理结果和所述重力导数异常信息生成滤波处理结果;
所述滤波单元,具体用于:
获取所述第二发育方向数据对应的权重数据;
利用所述第二发育方向数据、所述权重数据和所述重力导数异常信息生成补偿数据;
根据所述第一滤波结果和补偿数据计算目标测点在第二目标范围内的第二滤波处理结果;
将所述第二滤波处理结果作为滤波处理数据;
所述权重数据沿第二发育方向,权重值先增大后减小,并且目标测点的权重值大于除目标观测点外的网格点的权重值;
所述第二目标范围以所述目标测点为中心,并且沿第二发育方向包括5至9个网格点;
所述权重值沿第二发育方向依次为:1/16、2/16、3/16、4/16、3/16、2/16和1/16;
所述第一目标范围以所述目标测点为中心,并且沿第一发育方向包括5至9个网格点。
4.根据权利要求3所述的装置,其特征在于,所述计算模块,具体用于:
根据所述布格重力异常网格化数据生成重力垂直二次导数异常信息或重力垂直一次导数异常信息;
将所述重力垂直二次导数异常信息或重力垂直一次导数异常信息作为重力导数异常信息。
5.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至2任一所述走滑断裂确定方法。
6.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有计算机程序,所述计算机程序被处理器执行时实现权利要求1至2任一所述方法。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011172082.XA CN112464521B (zh) | 2020-10-28 | 2020-10-28 | 走滑断裂确定方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011172082.XA CN112464521B (zh) | 2020-10-28 | 2020-10-28 | 走滑断裂确定方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112464521A CN112464521A (zh) | 2021-03-09 |
CN112464521B true CN112464521B (zh) | 2024-05-28 |
Family
ID=74835572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011172082.XA Active CN112464521B (zh) | 2020-10-28 | 2020-10-28 | 走滑断裂确定方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112464521B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1877366A (zh) * | 2006-07-12 | 2006-12-13 | 杨辉 | 重磁延拓回返垂直导数目标优化处理技术 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN102466817A (zh) * | 2010-11-11 | 2012-05-23 | 中国石油天然气集团公司 | 归一化导数模法重力异常边界拾取方法 |
JP2016085114A (ja) * | 2014-10-26 | 2016-05-19 | 滋樹 水谷 | 地熱資源評価に有用な重力及び磁力データ総合的解析方法 |
KR101720327B1 (ko) * | 2016-10-28 | 2017-03-28 | 한국지질자원연구원 | 수중 이상체의 위치 측정 장치 및 방법 |
CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
CN110703347A (zh) * | 2019-10-24 | 2020-01-17 | 中国石油化工股份有限公司 | 基于构造背景的重力断裂影像识别方法 |
CN110989032A (zh) * | 2019-11-23 | 2020-04-10 | 北京地大惠信科技发展有限公司 | 一种基于倾斜角的重力水平总梯度断裂识别方法 |
CN111045089A (zh) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | 一种增强重磁场边界识别方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9778360B2 (en) * | 2013-12-17 | 2017-10-03 | Fugro N.V. | Method and system for generating a geoid via three computation spaces and airborne-acquired gravity data |
-
2020
- 2020-10-28 CN CN202011172082.XA patent/CN112464521B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1877366A (zh) * | 2006-07-12 | 2006-12-13 | 杨辉 | 重磁延拓回返垂直导数目标优化处理技术 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
CN102466817A (zh) * | 2010-11-11 | 2012-05-23 | 中国石油天然气集团公司 | 归一化导数模法重力异常边界拾取方法 |
JP2016085114A (ja) * | 2014-10-26 | 2016-05-19 | 滋樹 水谷 | 地熱資源評価に有用な重力及び磁力データ総合的解析方法 |
KR101720327B1 (ko) * | 2016-10-28 | 2017-03-28 | 한국지질자원연구원 | 수중 이상체의 위치 측정 장치 및 방법 |
CN108680952A (zh) * | 2018-04-03 | 2018-10-19 | 刘敬寿 | 一种走滑断层构造演化解析方法 |
CN110703347A (zh) * | 2019-10-24 | 2020-01-17 | 中国石油化工股份有限公司 | 基于构造背景的重力断裂影像识别方法 |
CN110989032A (zh) * | 2019-11-23 | 2020-04-10 | 北京地大惠信科技发展有限公司 | 一种基于倾斜角的重力水平总梯度断裂识别方法 |
CN111045089A (zh) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | 一种增强重磁场边界识别方法 |
Non-Patent Citations (1)
Title |
---|
走滑断裂和裂缝发育带的地震地质综合识别――以准噶尔盆地西北缘玛南地区为例;周卿;阎建国;黄立良;郑愈晟;;物探化探计算技术(02);第249-257页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112464521A (zh) | 2021-03-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Zhang et al. | Tectonic geomorphology of the Qilian Shan in the northeastern Tibetan Plateau: Insights into the plateau formation processes | |
Dannemiller et al. | A new method for determination of magnetization direction | |
Zheng et al. | A detailed receiver function image of the sedimentary structure in the Bohai Bay Basin | |
Camacho et al. | The 3-D gravity inversion package GROWTH2. 0 and its application to Tenerife Island, Spain | |
CN104636980B (zh) | 针对河道砂油藏类型油气汇集条件的地球物理表征方法 | |
AU2013293585B2 (en) | System and method for migration velocity modeling | |
CN103975341A (zh) | 基于4d饱和度模型和仿真模型的储层建模 | |
CN111624651B (zh) | 基于古地貌约束的储层预测方法及装置 | |
CN112464520B (zh) | 局部重力异常深度反演方法及装置 | |
Chen et al. | Interpretation of fault system in the Tana Sag, Kenya, using edge recognition techniques and Euler deconvolution | |
CN112464521B (zh) | 走滑断裂确定方法及装置 | |
CN117518233A (zh) | 一种剩余时差校正的微地震定位方法及相关设备 | |
Majdański et al. | 3D structure of the Earth's crust beneath the northern part of the Bohemian Massif | |
EP2885664A2 (en) | Methods and corresponding software module for quantifying risks or likelihoods of hydrocarbons being present in a geological basin or region | |
JP5078890B2 (ja) | 局地的かつリアルタイムの更新操作によって磁気測定データを取得し処理するための方法 | |
CN104950330B (zh) | 气顶油藏深度域成像的速度建模方法及*** | |
CN113204855B (zh) | 地层缩短速率确定方法、***、设备及可读介质 | |
CN109633760B (zh) | 一种基于自然电位等效融合成像的地下流体监测方法 | |
CN116029073A (zh) | 反演速度建模方法、装置、电子设备及介质 | |
Stewart | Scale dependence of strike and dip in sedimentary basins: Implications for field measurements and integrating subsurface datasets | |
Cox et al. | Straightening out lateral shifts in time-lapse seismic | |
CN103576188B (zh) | 一种消除速度误差影响的震源定位方法 | |
RU2563862C2 (ru) | Атрибут напряжения в горных породах | |
CN112130198B (zh) | 基于时差校正的多期次三维地震资料融合方法及装置 | |
CN113267810B (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 |