CN117434595A - 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法 - Google Patents

用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法 Download PDF

Info

Publication number
CN117434595A
CN117434595A CN202210824600.4A CN202210824600A CN117434595A CN 117434595 A CN117434595 A CN 117434595A CN 202210824600 A CN202210824600 A CN 202210824600A CN 117434595 A CN117434595 A CN 117434595A
Authority
CN
China
Prior art keywords
karst
distribution
leakage
reflecting
parameters
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.)
Pending
Application number
CN202210824600.4A
Other languages
English (en)
Inventor
吴兆徽
刘显太
杜玉山
吴光焕
李洪毅
刘西雷
魏超平
陈桂华
尹小梅
梁金萍
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Shengli Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN202210824600.4A priority Critical patent/CN117434595A/zh
Publication of CN117434595A publication Critical patent/CN117434595A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30Analysis
    • G01V1/301Analysis for determining seismic cross-sections or geostructures
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/612Previously recorded data, e.g. time-lapse or 4D
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/64Geostructures, e.g. in 3D data cubes
    • G01V2210/642Faults

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供一种用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,包括:步骤1,收集潜山井的井喷、泥浆漏失和放空数据;步骤2,统计潜山井泥浆漏失量及对应坐标;步骤3,编制钻遇岩溶后的漏失量平面分布图,展现当前开启未充填的岩溶分布;步骤4,用曲率参数优化后的蚂蚁地震数据体,反映潜山的裂缝及伴生岩溶;步骤5,调整曲率‑蚂蚁体参数,使之包含并反映目前岩溶漏失量的平面分布图。该用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法首次采用了潜山井泥浆漏失量平面分布图,来校正潜山缝洞地震体,从而尽最大可能如实反映了裂缝性复杂潜山内的溶蚀储集体,为准确预测尚未钻遇的潜山开启型缝洞提供了新的方法。

Description

用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法
技术领域
本发明涉及钻井-地球物理技术的交叉技术领域,特别是涉及到一种用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法。
背景技术
如果钻井过程中存在泥浆漏失现象,必然钻遇了无充填物的开启式岩溶。由于漏失量能反映存在岩溶,漏失量分布图可以作为验证其它方法是否正确的新方法,比如目前新兴的古岩溶图像学习法。
目前还未出现过利用“泥浆漏失量平面分布图”之类的钻井漏失量、放空量、井喷量等平面分布图,来反映裂缝性复杂潜山缝洞分布的方法。由于这类图件的原创性,“泥浆漏失量平面分布图”无论在百度、Bing、CNKI、万方等任何搜索引擎里,都是无法搜素到的。对于这种能反映开启型岩溶缝洞的平面图,以前之所以从未有笔者所在课题组之外的人编制过,原因是多方面的。
泥浆漏失量资料不属于开发路常用资料库能下载到的地质资料,而是属于跨学科的钻井资料,可通过联系泥浆公司或文献中找到该类资料。常规Direct等制图建模软件中,不能像加载测井数据那样,加载这类漏失量资料,也就是这类平面图是不易编制的。这些难题也促使地质和开发人员觉得“漏失量”这类难度大且名字晦气的图件没有必要编制。现代岩溶分布看上去随机性也很大,规律不易摸清。但是,岩溶泥浆漏失量平面分布图却是最可靠的反映地下存在岩溶的资料,即存在漏失量必然存在溶洞。
在申请号:CN202011276684.X的中国专利申请中,涉及到泥浆漏失风险评估方法及装置。所述方法包括以下步骤:1)提取叠前-叠后地震数据体频率信息,计算所使用数据体针对断层—裂缝体系延伸长度分辨率,分别提取叠后相干等几何属性、沿层高精度曲率和叠前波动力学属性;2)分别网格化叠后相干等几何属性、沿层高精度曲率属性和叠前波动力学属性,并将属性数据进行归一化处理,利用条件聚类方法逐步刻画出断裂-裂缝分布轮廓;3)叠合断裂纵向延伸长度>λ/4的分布图、纵向延伸长度介于λ/8~λ/4之间的断裂分布图、纵向延伸长度<λ/8断裂分布图。该发明目的在于提供泥浆漏失风险评估方法及装置,为油气勘探开发提供切实可行的技术支撑。
在申请号:CN202210032430.6的中国专利申请中,涉及到基于主控因素类比的钻前泥浆漏失预测方法及系统,对漏失区域数据资料进行分析,确定影响漏失的主控因素集以及漏失位置的泥浆漏失量;通过相关性分析得到主控因素集中各因素与泥浆漏失量之间的相关关系;获取待预测区域的各主控因素值,并根据得到的各因素与泥浆漏失量之间的相关关系,得到待预测区域的泥浆漏失量预测结果。实现复杂井况钻井前的泥浆漏失预测,通过在钻井前对地层进行井漏预测,对钻井提供指导,避开钻井漏失层段,提高钻探成功率。
在申请号:CN202210013790.1的中国专利申请中,涉及到一种基于BP神经网络的复杂井况钻前泥浆漏失预测方法,所述方法为:S1、收集目标区域历史地震属性数据、漏失案例相应的工程及地质数据资料;S2、对收集数据进行标准化处理,并且为了方便方法的实施,对所有数据进行转换处理;S3、以预处理后的历史地震数据为输入,以漏失情况为输出,以真实漏失状态为标准值,监督训练并优化得到漏失预测神经网络模型;S4、输入目标区域地震即时地震数据,由模型自动评判出此区域各深度对应的漏失情况。该发明解决了现有复杂井况钻前难以准确预测井内状况、容易出现泥浆漏失的问题。
以上现有技术均与本发明有较大区别,未能解决我们想要解决的技术问题,为此我们发明了一种新的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法。
发明内容
本发明的目的是提供一种反映了裂缝性潜山缝洞分布的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法。
本发明的目的可通过如下技术措施来实现:用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,该用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法包括:
步骤1,收集潜山井的井喷、泥浆漏失和放空数据;
步骤2,统计潜山井泥浆漏失量及对应坐标;
步骤3,编制钻遇岩溶后的漏失量平面分布图,展现当前开启未充填的岩溶分布;
步骤4,用曲率参数优化后的蚂蚁地震数据体,反映潜山的裂缝及伴生岩溶;
步骤5,调整曲率-蚂蚁体参数,使断裂分布图与漏失量平面图反映的位置信息一致。
本发明的目的还可通过如下技术措施来实现:
在步骤2中,对于井喷资料多的情况,编制井喷等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与井喷资料吻合。
在步骤2中,当放空资料大于5口时,编制放空资料等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与放空资料吻合。
在步骤3中,通过统计潜山井泥浆漏失量及对应坐标,编制出钻遇岩溶后的漏失量平面分布图。
在步骤5中,利用钻遇岩溶后的漏失量平面分布图,调整曲率参数优化蚂蚁追踪数据体。
在步骤5中,钻井资料中的井喷、放空数据,收集并在分析缝洞分布时是予以考虑的,但是,由于井喷和放空数据较少,这两类数据往往不足以编制平面图件,因此并不制平面图;仅进行泥浆漏失量制图。
在步骤5中,调整蚂蚁体平滑度SigmaXYZ参数越大越平滑。
本发明中的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,首次采用钻遇岩溶的“泥浆漏失量平面分布图”来调整优化地震参数,利用钻井资料调整曲率-蚂蚁体反映裂缝性潜山开启型缝洞分布,首次采用了潜山井泥浆漏失量平面分布图,来校正潜山缝洞地震体,从而尽最大可能如实反映了裂缝性复杂潜山内的溶蚀储集体,为准确预测尚未钻遇的潜山开启型缝洞提供了新的方法。
附图说明
图1为本发明的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法的一具体实施例的流程图;
图2为本发明的一具体实施例中钻遇岩溶井的泥浆漏失量等值线平面图;
图3为本发明的一具体实施例中草古1潜山及周边地震工区的蚂蚁追踪体平面图;
图4为本发明的一具体实施例中遇岩溶后的漏失量分布图、缝洞分布图的叠加吻合示意图;
图5为本发明的一具体实施例中最终反映溶蚀缝洞空间分布的三维地震缝洞体的示意图。
具体实施方式
应该指出,以下详细说明都是示例性的,旨在对本发明提供进一步的说明。除非另有指明,本文使用的所有技术和科学术语具有与本发明所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本发明的示例性实施方式。如在这里所使用的,除非上下文另外明确指出,否则单数形式也意图包括复数形式,此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作和/或它们的组合。
本发明提供了用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,包括:步骤1,收集潜山井的井喷、泥浆漏失和放空的钻井数据;步骤2,专门统计潜山井泥浆漏失量及对应坐标数据;步骤3,编制钻遇岩溶后的漏失量平面分布图,从而展现当前开启未充填的岩溶分布;步骤4,用曲率参数优化后的某蚂蚁地震数据体,反映潜山的裂缝及伴生岩溶;步骤5,调整曲率-蚂蚁体参数,使之尽可能反映并包含岩溶漏失量平面分布图。该用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,首次采用了潜山井泥浆漏失量平面分布图,来校正潜山缝洞地震体,从而尽最大可能如实反映了裂缝性复杂潜山内的溶蚀储集体,为准确预测尚未钻遇的潜山开启型缝洞提供了新的方法。
以下为应用本发明的几个具体实施例
实施例1
在应用本发明的一具体实施例1中,该用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法包括了以下步骤:
步骤1,收集潜山井的井喷、泥浆漏失和放空数据;
步骤2,统计潜山井泥浆漏失量及对应坐标;井喷资料虽少,遇到井喷资料多的情况,可编制井喷等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与井喷资料吻合。
放空资料一般很少,遇到放空资料大于5口的情况,可编制放空资料等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与放空资料吻合。
步骤3,编制钻遇岩溶后的漏失量平面分布图,从而展现当前开启未充填的岩溶分布;首次通过统计潜山井泥浆漏失量及对应坐标,编制出钻遇岩溶后的漏失量平面分布图。
步骤4,用曲率参数优化后的某蚂蚁地震数据体,反映潜山的裂缝及伴生岩溶;
步骤5,调整曲率-蚂蚁体参数,使之尽可能包含并反映目前岩溶漏失量的平面分布图。
本发明首次利用钻遇岩溶后的漏失量平面分布图,调整曲率参数优化蚂蚁追踪数据体。
对井喷、放空数据收集考虑,不制平面图,仅对漏失量收集并制平面图。这是由于井喷和放空数据常较少,这两类数据往往不足以编制平面图件。
由于可调参数仅6种且可选数据不多,通过对各类参数穷举尝试,使裂缝分布图、漏失量分布图两者尽可能一致。由于是断裂分布图、漏失量分布图两种复杂图像的对比,只能采用人工比对的方式。调整参数的依据为:溶洞主要呈现沿断裂分布的一般规律,断裂密集部位的泥浆漏失量更大,断裂稀疏区域的泥浆漏失量小。
实施例2
在应用本发明的一具体实施例2中,如图1所示,图1为本发明的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法的流程图。在草桥地区,用漏失量分布图来调整地震参数、反映潜山缝洞分布,具体包括以下步骤:
步骤1,收集潜山井的井喷、泥浆漏失和放空数据;
步骤2,统计潜山井泥浆漏失量及对应坐标;
步骤3,利用表1编制钻遇岩溶后的漏失量平面分布图,从而展现当前开启未充填的岩溶分布;
表1草桥局部井漏失量统计数据表
步骤4,用曲率参数优化后的某蚂蚁体,反映潜山的裂缝及伴生岩溶;
步骤5,调整曲率-蚂蚁体参数,使之调整后的平面图可能包含并反映目前岩溶漏失量的平面分布图。
此处调整蚂蚁体平滑度SigmaXYZ参数(1~3之间)越大越平滑,此处设置为2;
方差体Variance(1~11之间)滤波去噪参数设置为5;
垂向平滑去噪参数(0~200ms)设置为50ms;
构造平滑参数Structural smoothing中的Inline,Crossline,Vertical scale均设置1.5;
Vertical radius设置为25;Inline/Xline radius设置为1。
由于各参数可选值不多,通过尝试调整各类参数,从而在显示上,使断裂分布图的断裂密集区域、漏失量平面图的漏失量较大区域两者应尽可能一致,并使断裂稀疏区域、漏失量较小区域两者也应尽可能一致。从而针对代表性工区采用了以上参数。
从而使草桥地区的蚂蚁断裂-缝洞追踪体、遇岩溶后的泥浆漏失平面图尽可能存在包含和吻合关系。
实施例3
在应用本发明的一具体实施例3中,在草桥潜山边部地区应用本发明的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,具体包括:
步骤1,收集潜山井的泥浆漏失量数据(表2);
表2草桥潜山边部井漏失量统计数据表
步骤2,专门统计潜山井泥浆漏失量及对应坐标数据;
步骤3,编制钻遇岩溶后的漏失量平面分布图(图2),从而展现当前开启未充填的岩溶分布;
步骤4,用曲率参数优化后的某蚂蚁地震数据体(图3),反映潜山的裂缝及伴生岩溶;
步骤5,人工对比漏失量平面分布图和蚂蚁体断裂分布图(图4),调整曲率-蚂蚁体参数,使蚂蚁体(图5)尽可能反映并包含岩溶漏失量平面分布图。
最后应说明的是:以上所述仅为本发明的优选实施例而已,并不用于限制本发明,尽管参照前述实施例对本发明进行了详细的说明,对于本领域技术人员来说,其依然可以对前述实施例记载的技术方案进行修改,或者对其中部分技术特征进行等同替换。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
除说明书所述的技术特征外,均为本专业技术人员的已知技术。

Claims (7)

1.用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,该用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法包括:
步骤1,收集潜山井的井喷、泥浆漏失和放空数据;
步骤2,统计潜山井泥浆漏失量及对应坐标;
步骤3,编制钻遇岩溶后的漏失量平面分布图,展现当前开启未充填的岩溶分布;
步骤4,用曲率参数优化后的蚂蚁地震数据体,反映潜山的裂缝及伴生岩溶;
步骤5,调整曲率-蚂蚁体参数,使之包含并反映目前岩溶漏失量的平面分布图。
2.根据权利要求1所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤2中,对于井喷资料多的情况,编制井喷等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与井喷资料吻合。
3.根据权利要求1所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤2中,当放空资料大于5口时,编制放空资料等值线图,以便调整蚂蚁体SigmaXYZ参数,使缝洞体与放空资料吻合。
4.根据权利要求1所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤3中,通过统计潜山井泥浆漏失量及对应坐标,编制出钻遇岩溶后的漏失量平面分布图。
5.根据权利要求1所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤5中,利用钻遇岩溶后的漏失量平面分布图,调整曲率参数优化蚂蚁追踪数据体。
6.根据权利要求5所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤5中,钻井资料中的井喷、放空数据,收集并在分析缝洞分布时是予以考虑的,但是,由于井喷和放空数据较少,这两类数据往往不足以编制平面图件,因此并不制平面图;仅进行泥浆漏失量制图。
7.根据权利要求5所述的用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法,其特征在于,在步骤5中,调整蚂蚁体平滑度SigmaXYZ参数越大越平滑。
CN202210824600.4A 2022-07-13 2022-07-13 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法 Pending CN117434595A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210824600.4A CN117434595A (zh) 2022-07-13 2022-07-13 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210824600.4A CN117434595A (zh) 2022-07-13 2022-07-13 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法

Publications (1)

Publication Number Publication Date
CN117434595A true CN117434595A (zh) 2024-01-23

Family

ID=89544909

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210824600.4A Pending CN117434595A (zh) 2022-07-13 2022-07-13 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法

Country Status (1)

Country Link
CN (1) CN117434595A (zh)

Similar Documents

Publication Publication Date Title
AU2009260453B2 (en) Heterogeneous earth models for a reservoir field
CN103245971B (zh) 井中断点引导的地震小断层解释方法及装置
CA2823710C (en) Methods and systems regarding models of underground formations
CN104329079A (zh) 一种识别气测录井油气层的方法以及系统
CN111611673B (zh) 一种碳酸盐岩油藏古暗河型储层建模方法
CN104155701A (zh) 一种利用叠前地震资料及井信息的多尺度裂缝预测方法
CN104636980A (zh) 针对河道砂油藏类型油气汇集条件的地球物理表征方法
CN118273768B (zh) 基于gis底座的煤矿水害全息本源模态预警方法和系统
CN117474340B (zh) 用于地铁盾构施工沉降的风险评价方法及系统
CN112489208A (zh) 基于蚂蚁算法的裂缝片提取方法和三维地质模型构建方法
CN115877447A (zh) 一种直平联合井网条件下地震约束三维地质建模的储层预测方法
CN112084553B (zh) 一种用于隧道规划的勘测方法
CN111830558B (zh) 一种断裂带雕刻方法
CN114581556B (zh) 一种区域地质调查中的数字填图方法
CN117434595A (zh) 用岩溶漏失分布图调整地震参数反映潜山缝洞分布的方法
CN112285773B (zh) 一种隐伏断层二元结构定量表征方法
Lyu et al. Classification, modeling and characterization of marine carbonate paleokarst reservoirs in Tahe Oilfield, Tarim Basin, China
CN114060015A (zh) 一种致密砂岩含气性的评价方法及装置
CN105467447A (zh) 相控趋势能量匹配的地震保幅评价方法
CN115479540B (zh) 一种基于激光扫描的岩溶突水灾害量级评估方法和系统
CN108510586B (zh) 一种构建地质导向地层模型的方法
CN112147698B (zh) 裂缝发育带识别与特征确定的方法及系统
Burger et al. Geostatistics and the polygonal method: a re-examination
CN115903028A (zh) 一种井震联合描述碳酸盐岩滩相储层分布的方法、装置、设备及计算机可读存储介质
CN118506916A (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