本公开提供了一种重力异常的提取方法、装置、设备和存储介质，属于重力勘探技术领域。所述方法包括：对目标区域的布格重力异常进行地形校正和曲化平处理，得到第一重力异常；对第一重力异常进行常密度重力剥层处理，得到第二重力异常；对第二重力异常进行变密度重力剥层处理，得到第三重力异常，变密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度在沿层横向上是变化的；确定目标区域的第一重力区域场；根据第三重力异常和第一重力区域场提取目标地质体的重力异常。对第二重力异常进行变密度重力剥层处理，可以消除地层沿层横向上的变化密度对目标地质体的重力异常的影响，从而可以提高提取到的目标重力异常的准确性。准确性。准确性。
【技术实现步骤摘要】
重力异常的提取方法、装置、设备和存储介质
[0001]本公开涉及重力勘探
，特别涉及一种重力异常的提取方法、装置、设备和存储介质。
技术介绍
[0002]重力勘探是油气勘探中的一种重要的方法，在区域构造、断线盆地、局部构造带、潜山、裂谷、山前带等研究中具有重要作用。在重力勘探中，由地下岩矿体密度分布不均匀所引起的实际观测重力值与理论上的正常重力值之间的偏差，称为重力异常。通过对重力异常进行定性或定量解释，可以推测出目标区域的地质体构造以及岩矿体的分布情况。
[0003]相关技术中，重力异常的提取方法包括：对目标区域的布格重力异常进行地形校正和曲化平处理，得到第一重力异常；对第一重力异常进行常密度重力剥层处理，得到第二重力异常，该常密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度均为常数；确定目标区域的第一重力区域场；将第二重力异常与第一重力区域场之间的差值，确定为目标地质体的重力异常。
[0004]在复杂区域中，采用上述方法提取到的目标地质体的重力异常的准确性较低，可能会出现错误的地质解释。
技术实现思路
[0005]本公开实施例提供了一种重力异常的提取方法、装置、设备和存储介质，能够提高提取的目标地质体的重力异常的准确性。所述技术方案如下：
[0006]第一方面，提供了一种重力异常的提取方法，所述方法包括：对目标区域的布格重力异常进行地形校正和曲化平处理，得到第一重力异常；对所述第一重力异常进行常密度重力剥层处理，得到第二重力异常，所述常密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度均为常数；根据所述目标区域的测井资料和电法勘探资料，确定所述目标区域的地层沿层横向密度分布模型，所述地层沿层横向密度分布模型中，所述目标地质体所在地层到所述地表之间的多层地层的密度在沿层横向上是变化的；根据所述地层沿层横向密度分布模型，对所述第二重力异常进行变密度重力剥层处理，得到第三重力异常；确定所述目标区域的第一重力区域场；根据所述第三重力异常和所述第一重力区域场，提取所述目标地质体的重力异常。
[0007]可选地，所述根据所述地层沿层横向密度分布模型，对所述第二重力异常进行地层变密度重力剥层处理，得到第三重力异常，包括：根据所述地层沿层横向密度分布模型，计算每个所述地层的上界面与下界面之间的密度差；根据所述密度差，计算每个所述地层的变密度重力剥层值；将所述第二重力异常与所述多层地层的变密度重力剥层值的差值，作为所述第三重力异常。
[0008]可选地，所述根据所述第三重力异常和所述第一重力区域场，提取所述目标地质体的重力异常，包括：对所述第三重力异常进行湿度重力校正，得到第四重力异常；将所述
第四重力异常与所述第一重力区域场之间的差值，确定为所述目标地质体的重力异常。
[0009]可选地，所述对所述第三重力异常进行湿度重力校正，得到第四重力异常，包括：根据所述目标区域的电阻率数据，确定潜水面；根据所述潜水面的上层地层与下层地层之间的密度差，计算湿度重力校正值；将所述第三重力异常与所述湿度重力校正值之间的差值，作为所述第四重力异常。
[0010]可选地，所述确定所述目标区域的第一重力区域场，包括：确定所述目标地质体的磁力异常；确定所述目标区域的多个插值切割半径；计算每个所述插值切割半径对应的第二重力区域场；根据所述第二重力区域场与所述磁力异常的相似性，从所述多个插值切割半径中确定目标插值切割半径，所述目标插值切割半径对应的第二重力区域场与所述目标磁力异常的相似性最大；将所述目标插值切割半径对应的第二重力区域场，确定为所述目标区域的第一重力区域场。
[0011]第二方面，提供了一种重力异常的提取装置，所述装置包括：地形校正和曲化平处理模块，用于对目标区域的布格重力异常进行地形校正和曲化平处理，得到第一重力异常；常密度重力剥层处理模块，用于对所述第一重力异常进行常密度重力剥层处理，得到第二重力异常，所述常密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度均为常数；第一确定模块，用于根据所述目标区域的测井资料和电法勘探资料，确定所述目标区域的地层沿层横向密度分布模型，所述地层沿层横向密度分布模型中，所述目标地质体所在地层到所述地表之间的多层地层的密度在沿层上是横向变化的；变密度重力剥层处理模块，用于根据所述地层沿层横向密度分布模型，对所述第二重力异常进行变密度重力剥层处理，得到第三重力异常；第二确定模块，用于确定所述目标区域的第一重力区域场；重力异常提取模块，用于根据所述第三重力异常和所述第一重力区域场，提取所述目标地质体的重力异常。
[0012]可选地，所述变密度重力剥层处理模块用于根据所述地层沿层横向密度分布模型，计算每个所述地层的上界面与下界面之间的密度差；根据所述密度差，计算每个所述地层的变密度重力剥层值；将所述第二重力异常与所述多层地层的变密度重力剥层值的差值，作为所述第三重力异常。
[0013]可选地，所述重力异常提取模块用于对所述第三重力异常进行湿度重力校正，得到第四重力异常；将所述第四重力异常与所述第一重力区域场之间的差值，确定为所述目标地质体的重力异常。
[0014]可选地，所述重力异常提取模块用于根据所述目标区域的电阻率数据，确定潜水面；根据所述潜水面的上层地层与下层地层之间的密度差，计算湿度重力校正值；将所述第三重力异常与所述湿度重力校正值之间的差值，作为所述第四重力异常。
[0015]可选地，所述第二确定模块用于确定所述目标地质体的磁力异常；确定所述目标区域的多个插值切割半径；计算每个所述插值切割半径对应的第二重力区域场；根据所述第二重力区域场与所述磁力异常的相似性，从所述多个插值切割半径中确定目标插值切割半径，所述目标插值切割半径对应的第二重力区域场与所述目标磁力异常的相似性最大；将所述目标插值切割半径对应的第二重力区域场，确定为所述目标区域的第一重力区域场。
[0016]第三方面，提供了一种计算机设备，包括：处理器；用于存储处理器可执行指令的
存储器；其中，所述处理器被配置为执行第一方面所述的方法。
[0017]第四方面，提供了一种计算机可读介质，当计算机可读介质中的指令由计算机设备的处理器执行时，使得计算机设备能够执行第一方面所述的方法。
[0018]第五方面，提供了一种计算机程序产品，包括计算机程序/指令，其特征在于，所述计算机程序/指令被处理器执行时实现第一方面所述的方法。
[0019]本公开实施例提供的技术方案带来的有益效果是：
[0020]本公开实施例中，对目标区域的布格重力异常进行地形校正和曲化平处理得到第一重力异常；对第一重力异常进行常密度重力剥层处理，得到第二重力异常，常密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度为常数；根据目标区域的测井资料和电法勘探资料，确定目标区域的地层沿层横向密度分布模型；根据地层沿层横向密度分布模型，对第二重力异常进行变密度重力剥层处理，得到第三重力异常；根据第三本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种重力异常的提取方法，其特征在于，所述方法包括：对目标区域的布格重力异常进行地形校正和曲化平处理，得到第一重力异常；对所述第一重力异常进行常密度重力剥层处理，得到第二重力异常，所述常密度重力剥层处理中，目标地质体所在地层到地表之间的多层地层的密度均为常数；根据所述目标区域的测井资料和电法勘探资料，确定所述目标区域的地层沿层横向密度分布模型，所述地层沿层横向密度分布模型中，所述目标地质体所在地层到所述地表之间的多层地层的密度在沿层横向上是变化的；根据所述地层沿层横向密度分布模型，对所述第二重力异常进行变密度重力剥层处理，得到第三重力异常；确定所述目标区域的第一重力区域场；根据所述第三重力异常和所述第一重力区域场，提取所述目标地质体的重力异常。2.根据权利要求1所述的方法，其特征在于，所述根据所述地层沿层横向密度分布模型，对所述第二重力异常进行地层变密度重力剥层处理，得到第三重力异常，包括：根据所述地层沿层横向密度分布模型，计算每个所述地层的上界面与下界面之间的密度差；根据所述密度差，计算每个所述地层的变密度重力剥层值；将所述第二重力异常与所述多层地层的变密度重力剥层值的差值，作为所述第三重力异常。3.根据权利要求1所述的方法，其特征在于，所述根据所述第三重力异常和所述第一重力区域场，提取所述目标地质体的重力异常，包括：对所述第三重力异常进行湿度重力校正，得到第四重力异常；将所述第四重力异常与所述第一重力区域场之间的差值，确定为所述目标地质体的重力异常。4.根据权利要求3所述的方法，其特征在于，所述对所述第三重力异常进行湿度重力校正，得到第四重力异常，包括：根据所述目标区域的电阻率数据，确定潜水面；根据所述潜水面的上层地层与下层地层之间的密度差，计算湿度重力校正值；将所述第三重力异常与所述湿度重力校正值之间的差值，作为所述第四重力异常。5.根据权利要求1至4任一项所述的方法，其特征在于，所述确定所述目标区域的第一重力区域场，包括：确定所述目标地质体的磁力异常；确定所述目标区域的多个插值切割半径；计算每个所述插值切割半径对应的第二重力区域...
【专利技术属性】
技术研发人员：刘云祥，赵文举，王财富，
申请(专利权)人：中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。