CN110244358A - 一种构造破坏导致的油气逸散区的判识方法 - Google Patents

一种构造破坏导致的油气逸散区的判识方法 Download PDF

Info

Publication number
CN110244358A
CN110244358A CN201910521732.8A CN201910521732A CN110244358A CN 110244358 A CN110244358 A CN 110244358A CN 201910521732 A CN201910521732 A CN 201910521732A CN 110244358 A CN110244358 A CN 110244358A
Authority
CN
China
Prior art keywords
measuring point
unit
formula
mud stone
oil
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
Application number
CN201910521732.8A
Other languages
English (en)
Other versions
CN110244358B (zh
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 Henan Oilfield Branch Co
Exploration and Development Research Institute of Sinopec Shengli Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Henan Oilfield Branch 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 Henan Oilfield Branch Co filed Critical China Petroleum and Chemical Corp
Priority to CN201910521732.8A priority Critical patent/CN110244358B/zh
Publication of CN110244358A publication Critical patent/CN110244358A/zh
Application granted granted Critical
Publication of CN110244358B publication Critical patent/CN110244358B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • 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
    • 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/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Remote Sensing (AREA)
  • General Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提出一种构造破坏导致的油气逸散区的判识方法，属于油气勘探技术领域，其通过判断不同成藏期时的断‑盖配置关系确定油气的垂向保存条件，进而识别油气逸散区。该构造破坏导致的油气逸散区的判识方法，包括如下步骤：在研究区的目的层内，选取均匀覆盖于研究区范围内断层的多个测点；根据研究区的录井、钻井及地震资料，计算每一测点对应的断层泥岩比率下限值SGRlim；根据研究区的录井、钻井及地震资料，计算每一测点对应的成藏期断层泥岩比率SGR古；根据研究区录井数据，计算每一测点对应的现今断层泥岩比率SGR现今；当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，则该测点存在油气逸散。

Description

一种构造破坏导致的油气逸散区的判识方法
技术领域
本发明属于油气勘探技术领域，尤其涉及一种构造破坏导致的油气逸散区的判识方法。
背景技术
实际勘探结果表明，根据流体包裹体法、GOI测试等方法确定某井位处有油气成藏，但是常常出现钻探效果差、产量低的现象。经研究发现，对于两期或多期成藏的油气藏而言，油气成藏后，若各成藏期时的断-盖配置条件均较差，则油气藏由于构造抬升或断裂开启遭到破坏，油气极易沿断层向上运移发生逸散，最终残留油气极少，产量较低，钻探效果差。
然而，针对油气成藏期的研究，目前主要方法包括生排烃期法、流体包裹体法和圈闭时间法等，这些方法仅能确定成藏期，并不能确定油气藏是否遭到破坏，可能会误导油气勘探预测结果，影响油气的勘探进程。
因而，识别由构造破坏导致的油气逸散区，从而在勘探定井时进行规避，对于提高高产探井的钻遇成功率具有重要意义。
发明内容
本发明提出一种构造破坏导致的油气逸散区的判识方法，通过判断不同成藏期时的断-盖(断层与盖层)配置关系确定油气的垂向保存条件，进而识别油气逸散区，应用于井位勘探部署中，能够有效增加钻遇高产探井的成功率，有利于减少钻井成本。
为了达到上述目的，本发明采用的技术方案为：
一种构造破坏导致的油气逸散区的判识方法，包括如下步骤：
选取测点：在研究区的目的层内，选取均匀覆盖于研究区范围内断层的多个测点；
计算断层泥岩比率下限值：以研究区内具有实测岩石排替压力数据的检测点为样品点，根据样品点的实测数据，利用公式(1)拟合研究区内岩石排替压力与岩石的埋深和泥质含量乘积之间的函数关系，获得常数a和b；再根据研究区的录井、钻井及地震资料，利用公式(2)计算每一测点对应的断层泥岩比率下限值，公式(1)和(2)的表达式如下：
其中，Pd样为样品点的岩石排替压力，单位为MPa；a、b均为拟合常数；Z样为样品点的岩石埋深，单位为m；V样为样品点的岩石泥质含量，单位为％；SGRlim为测点的断层泥岩比率下限值；Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；ρw为地层水密度，单位为g/cm3；g为重力加速度，单位为m/s2；θ为测点所在断裂的断裂倾角，单位为°；V储为测点对应储层的泥质含量，单位为％；
计算成藏期断层泥岩比率：根据研究区的录井、钻井及地震资料，利用公式(3)计算每一测点对应的成藏期断层泥岩比率，公式(3)的表达式如下：
其中，SGR古为测点的成藏期断层泥岩比率；d古为测点所在断裂在成藏期时的古断距，单位为m；H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m；
计算现今断层泥岩比率：根据研究区录井数据，利用公式(4)计算每一测点对应的现今断层泥岩比率，公式(4)的表达式如下：
其中，SGR现今为测点的现今断层泥岩比率；d现今为测点所在断裂的现今断距，单位为m；H为测点所在断裂的现今累积泥岩厚度，单位为m；
判识是否存在油气逸散：当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，则该测点存在油气逸散。
作为优选，计算断层泥岩比率下限值的步骤中，根据研究区的三维地震资料和钻井资料，利用公式(5)计算获得公式(2)中的参数Z断'，公式(5)的表达式如下：
其中，Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；T断为测点所在断裂的断层岩开始成岩的时间，单位为Ma；T储为测点对应储层开始成岩的时间，单位为Ma；ρr为上覆沉积岩石密度，单位为g/cm3；ρw为地层水密度，单位为g/cm3；θ为测点所在断裂的断裂倾角，单位为°。
作为优选，计算断层泥岩比率下限值的步骤中，公式(2)中参数V储的具体计算步骤如下：根据靠近测点的已知井的录井和GR测井数据，利用公式(6)和公式(7)计算获得V储，公式(6)和公式(7)的表达式如下：
其中，Vsh为测点对应储层的泥质含量指数，单位为％；GR目的层为靠近该测点的已知井在目的层的GR测井曲线值，单位为API；GR纯泥岩层为靠近该测点的已知井的纯泥岩层的GR测井曲线值，单位为API；GR纯砂岩层为靠近该测点的已知井的纯砂岩层的GR测井曲线值，单位为API；V储为测点对应储层的泥质含量，单位为％；GCUR为与地层有关的经验系数，古近系和新近系地层的GCUR取值为3.7，老地层的GCUR取值为2。
作为优选，公式(6)中的参数GR目的层、GR纯泥岩层和GR纯砂岩层的具体获取步骤如下：根据目的层的深度，以所述已知井的GR测井曲线在相同深度对应的GR测井值作为GR目的层，选取所述已知井中大套泥岩的GR测井最高值作为GR纯泥岩层，选取所述已知井中大套砂岩的GR测井最低值作为GR纯砂岩层。
作为优选，计算成藏期断层泥岩比率的步骤中，公式(3)中的参数d古的具体获取步骤如下：切取垂直于该断裂不同部位的地震剖面，读取该断裂成藏期时不同部位的断距，取其中的最大断距作为该断裂自成藏期至现今的活动量，将该断裂的现今断距与其自成藏期至现今的活动量作差，得到该断裂在成藏期时的古断距d古。
作为优选，计算成藏期断层泥岩比率的步骤中，公式(3)中参数H0的具体计算步骤如下：将目的层中在古断距d古内沉积的地层等分为多个统计单元，分别统计每个统计单元对应的现今泥岩累积厚度，根据去压实校正法利用公式(8)计算获得H0，公式(8)的表达式如下：
其中，H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m；hi为第i个统计单元对应的现今泥岩累积厚度，单位为m；Ci为第i个统计单元对应的泥岩压实系数，利用去压实校正法获得；α为地层倾角，单位为°；n为统计单元的总数。
与现有技术相比，本发明的优点和积极效果在于：
1、本发明提供的构造破坏导致的油气逸散区的判识方法，从构造对油气成藏的破坏作用入手，通过评价不同时期断-盖配置的有效性，对比现今与各成藏期时的断-盖配置关系，明确构造破坏导致的油气逸散位置，其判识结果与实钻结果具有很好的一致性，判识准确性高，为油气成藏过程提供了技术支持，可以广泛适用于油气资源地质勘探及有利区优选等多种领域。
2、本发明提供的构造破坏导致的油气逸散区的判识方法，仅通过地震、钻井、测井等基础资料即可完成判识，无需做大量的实验，操作简单且判识成本低，为油田勘探选区降低风险，对提高钻遇高产探井成功率具有重要的理论价值。
附图说明
图1为本发明实施例所提供的构造破坏导致的油气逸散区的判识方法的流程示意图；
图2为本发明实施例所提供的断-盖配置下储层中某点的力学示意图；
图3为本发明实施例1提供的松辽盆地某断陷DL1地区的岩石排替压力与岩石的埋深和泥质含量乘积之间的关系图。
具体实施方式
下面将对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
如图1所示，本发明实施例提供了一种构造破坏导致的油气逸散区的判识方法，包括如下步骤：
S1选取测点：在研究区的目的层内，选取均匀覆盖于研究区范围内断层的多个测点。在本步骤中，需要说明的是，当断层延伸长度较大时，可在同一条断层选取多个测点。
S2计算断层泥岩比率下限值：以研究区内具有实测岩石排替压力数据的检测点为样品点，根据样品点的实测数据，利用公式(1)拟合研究区内岩石排替压力与岩石的埋深和泥质含量乘积之间的函数关系，获得常数a和b；再根据研究区的录井、钻井及地震资料，利用公式(2)计算每一测点对应的断层泥岩比率下限值，公式(1)和(2)的表达式如下：
其中，Pd样为样品点的岩石排替压力，单位为MPa；a、b均为拟合常数；Z样为样品点的岩石埋深，单位为m；V样为样品点的岩石泥质含量，单位为％；SGRlim为测点的断层泥岩比率下限值；Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；ρw为地层水密度，单位为g/cm3；g为重力加速度，单位为m/s2；θ为测点所在断裂的断裂倾角，单位为°；V储为测点对应储层的泥质含量，单位为％。
本步骤依据断-盖配置的临界条件，推导出断层泥岩比率下限值的计算公式(2)，其推导原理如下：
如图2所示，根据断-盖配置下储层中某点的力学分析可知，该点所受的动力为浮力和储层排替压力，阻力为断过盖层的断层岩的排替压力，当处于临界条件下时，即动力和阻力平衡，可得如下关系式：
Pd储+F浮sinθ＝Pd断lim (9)
其中，Pd储为储层排替压力，单位为MPa；F浮为浮力，单位为MPa；Pd断lim为断层岩排替压力下限，单位为MPa；θ为断裂倾角，单位为°。
由于岩石的排替压力与岩石的埋深和泥质含量具有一定的函数关系，因而，Pd储可表示如下：
对于断层岩而言，泥质含量利用断层的SGR表征，因而，Pd断lim可表示如下：
而F浮可表示如下：
F浮＝ρwgZ (12)
将公式(10)-公式(12)代入公式(9)，即可得上述公式(2)。
需要说明的是，上述公式(2)中的参数Z断'的计算方法，优选为：根据研究区的三维地震资料和钻井资料，利用公式(5)计算获得公式(2)中的参数Z断'，公式(5)的表达式如下：
其中，Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；T断为测点所在断裂的断层岩开始成岩的时间，单位为Ma；T储为测点对应储层开始成岩的时间，单位为Ma；ρr为上覆沉积岩石密度，单位为g/cm3；ρw为地层水密度，单位为g/cm3；θ为测点所在断裂的断裂倾角，单位为°。
还需要说明的是，上述公式(2)中的参数V储的计算方法，优选为：根据靠近测点的已知井的录井和GR测井数据，利用公式(6)和公式(7)计算获得V储，公式(6)和公式(7)的表达式如下：
其中，Vsh为测点对应储层的泥质含量指数，单位为％；GR目的层为靠近该测点的已知井在目的层的GR测井曲线值，单位为API；GR纯泥岩层为靠近该测点的已知井的纯泥岩层的GR测井曲线值，单位为API；GR纯砂岩层为靠近该测点的已知井的纯砂岩层的GR测井曲线值，单位为API；V储为测点对应储层的泥质含量，单位为％；GCUR为与地层有关的经验系数，古近系和新近系地层的GCUR取值为3.7，老地层的GCUR取值为2。其中，参数GR目的层、GR纯泥岩层和GR纯砂岩层的具体获取步骤如下：根据目的层的深度，以所述已知井的GR测井曲线在相同深度对应的GR测井值作为GR目的层，选取所述已知井中大套泥岩的GR测井最高值作为GR纯泥岩层，选取所述已知井中大套砂岩的GR测井最低值作为GR纯砂岩层。需要说明的是，大套泥岩是指录井资料中显示的厚度大且岩性分布稳定的泥岩，大套砂岩是指录井资料中显示的厚度大且岩性分布稳定的砂岩。此外，还需要说明的是，当无靠近测点的已知井时，可选取测点附近的两个已知井A和B，依据于测点到两个已知井之间的距离，以及井A和井B的GR测井值，通过线性插值法计算获得该测点的GR测井值。
S3计算成藏期断层泥岩比率：根据研究区的录井、钻井及地震资料，利用公式(3)计算每一测点对应的成藏期断层泥岩比率，公式(3)的表达式如下：
其中，SGR古为测点的成藏期断层泥岩比率；d古为测点所在断裂在成藏期时的古断距，单位为m；H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m。
本步骤通过恢复古断距和恢复古断距范围内的累积古泥岩厚度，计算得到测点对应的成藏期断层泥岩比率。
需要说明的是，上述公式(3)中的参数d古的获取步骤，优选为：切取垂直于该断裂不同部位的地震剖面，读取该断裂成藏期时不同部位的断距，取其中的最大断距作为该断裂自成藏期至现今的活动量，将该断裂的现今断距与其自成藏期至现今的活动量作差，得到该断裂在成藏期时的古断距d古。
上述公式(3)中参数H0的计算步骤，优选为：将目的层中在古断距d古内沉积的地层等分为多个统计单元，分别统计每个统计单元对应的现今泥岩累积厚度，根据去压实校正法利用公式(8)计算获得H0，公式(8)的表达式如下：
其中，H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m；hi为第i个统计单元对应的现今泥岩累积厚度，单位为m；Ci为第i个统计单元对应的泥岩压实系数，利用去压实校正法获得；α为地层倾角，单位为°；n为统计单元的总数。
需要说明的是，泥岩压实系数Ci可利用经验公式法或孔隙度-岩性法计算获得。其中，利用经验公式法计算泥岩压实系数Ci的公式如下：
Ci＝q-plnHi (13)
其中，Ci为第i个统计单元对应的泥岩压实系数；Hi为第i个统计单元深度跨度的中值，单位为m；对于泥岩，取p＝0.1、q＝1.46。
利用孔隙度-岩性法计算泥岩压实系数Ci的公式如下：
其中，当Ci为第i个统计单元对应的泥岩压实系数；Hi为第i个统计单元深度跨度的中值，单位为m；φ为深度为Hi处的孔隙度，单位为％；φ0为地表初始孔隙度，单位为％。
可以理解的是，本领域技术人员也可以采用其他的方式恢复古断距和恢复古断距范围内的累积古泥岩厚度。
S4计算现今断层泥岩比率：根据研究区录井数据，利用公式(4)计算每一测点对应的现今断层泥岩比率，公式(4)的表达式如下：
其中，SGR现今为测点的现今断层泥岩比率；d现今为测点所在断裂的现今断距，单位为m；H为测点所在断裂的现今累积泥岩厚度，单位为m。
S5判识是否存在油气逸散：当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，则该测点存在油气逸散。本步骤中，需要说明的是，当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，该测点油气虽然在各成藏期发生充注，但是由于断-盖配置关系较差导致油气遭到破坏，油气大量逸散，因而，该测点的现今断-盖配置关系和成藏期时断-盖配置关系均较差，可确定该测点存在构造破坏导致的油气逸散，不利于成藏。本步骤从构造对油气成藏的破坏作用入手，根据断-盖配置有效性的实质，通过分析盖层与断层在不同成藏期时的配置关系，确定是否存在油气逸散。
本发明提供的上述构造破坏导致的油气逸散区的判识方法，从构造对油气成藏的破坏作用入手，通过评价不同时期断-盖配置的有效性，对比现今与各成藏期时的断-盖配置关系，明确构造破坏导致的油气逸散位置，其判识结果与实钻结果具有很好的一致性，判识准确性高，为油气成藏过程提供了技术支持，可以广泛适用于油气资源地质勘探及有利区优选等多种领域。同时，本发明提供的上述构造破坏导致的油气逸散区的判识方法，仅通过地震、钻井、测井等基础资料即可完成判识，无需做大量的实验，操作简单且判识成本低，为油田勘探选区降低风险，对提高钻遇高产探井成功率具有重要的理论价值。
为了更清楚详细地介绍本发明实施例所提供的构造破坏导致的油气逸散区的判识方法，下面将结合具体实施例进行描述。
实施例1
以松辽盆地某断陷DL1地区为研究区，研究区主要发育白垩纪地层，自下而上分别为：沙河子组(K1sh)、营城组(K1yc)、登娄库组(K1d)、泉头组(K1q)、青山口组(K2qn)、姚家组(K2y)和嫩江组(K2n)等。勘探成果表明，从沙河子组至泉头组均有油气分布，其中泉头组和营城组是主力含油气层段。油气成藏条件分析结果认为，DL1地区发育两套区域性泥岩盖层，分别位于营三段和泉二段。由于营三段盖层经历了两期油气成藏，泉二段盖层仅经历了后期油气成藏，因而，选定营三段为目的层。营三段中构造破坏导致的油气逸散区的判识方法，包括如下步骤：
(一)选取测点：在营三段选取了均匀覆盖于DL1地区断层的31个测点。
(二)计算断层泥岩比率下限值：
以DL1地区内具有实测岩石排替压力数据的检测点为样品点，根据样品点的实测数据，利用公式(1)拟合DL1地区岩石排替压力与岩石的埋深和泥质含量乘积之间的函数关系(关系图详见图3)，拟合得到关系式为R2＝0.828，获得常数a＝0.3772、b＝0.0014。
再根据DL1地区的录井、钻井及地震资料，利用公式(2)计算每一测点对应的断层泥岩比率下限值，计算结果如表1所示。需要说明的是，计算时，公式(2)中的ρw取1.00g/cm3；g取9.8m/s2；参数Z断'的计算步骤为：根据DL1地区的三维地震资料和钻井资料，利用公式(5)计算获得，计算时ρr取2.52g/cm3，计算结果如表1所示；参数V储的具体计算步骤为：根据靠近测点的已知井的录井和GR测井数据，利用公式(6)和公式(7)计算获得，其中，由于营城组为老地层，GCUR取值为2；根据目的层的深度，以所述已知井的GR测井曲线在相同深度对应的GR测井值作为GR目的层，选取所述已知井中大套泥岩的GR测井最高值作为GR纯泥岩层，选取所述已知井中大套砂岩的GR测井最低值作为GR纯砂岩层，计算结果如表1所示。
表1 DL1地区营三段31个测点对应的SGRlim的计算参数及结果
(三)计算成藏期断层泥岩比率：
根据DL1地区的录井、钻井及地震资料，利用公式(3)计算每一测点对应的成藏期断层泥岩比率，计算结果如表2所示。
表2 DL1地区营三段31个测点对应的SGR古的计算参数及结果
需要说明的是，公式(3)中的参数d古的具体获取步骤如下：切取垂直于该断裂不同部位的地震剖面，读取该断裂成藏期时不同部位的断距(即登娄库组沉积末期时的断距—T3断距)，取其中的最大断距作为该断裂自成藏期至现今的活动量，将该断裂的现今断距与其自成藏期至现今的活动量作差，得到该断裂在成藏期时的古断距d古，计算结果如表2所示。
公式(3)中参数H0的具体计算步骤如下：将目的层中在古断距d古内沉积的地层等分为多个统计单元，分别统计每个统计单元对应的现今泥岩累积厚度，根据去压实校正法利用公式(8)计算获得H0。以12号测点为例，12号测点的古断距为169.6m，按照每20米为一个统计单元，可划分为9个统计单元(2800～2820、2820～2840、……2960～2969.6)，分别统计每个统计单元对应的现今泥岩累积厚度，利用公式(8)计算获得H0，其中，Ci利用经验公式法计算获得，地层倾角α为4°，计算结果如表3所示。
表3 12号测点所在断裂未受压实作用的累积古泥岩厚度的计算参数及结果
(四)计算现今断层泥岩比率：根据DL1地区的录井数据，利用公式(4)计算每一测点对应的现今断层泥岩比率，计算结果如表4所示。
表4 DL1地区营三段31个测点对应的SGR现今的计算参数及结果
编号 Z/m d<sub>现今</sub>/m H/m SGR<sub>现今</sub>
1 2960 59.4 28.7 0.484
2 2890 82.5 39.9 0.484
3 2850 52.8 25.6 0.484
4 2600 33 20.2 0.613
5 2370 128.7 69.2 0.538
6 2270 66 26.9 0.407
7 1930 148.5 60.4 0.407
8 3060 69.3 33.5 0.484
9 2920 108.9 52.7 0.484
10 2500 32.9 21.9 0.665
11 2250 112.2 74.6 0.665
12 2800 237.6 120.9 0.509
13 2300 36.3 18.5 0.509
14 2230 135.4 82.1 0.606
15 2080 75.9 20.0 0.264
16 1900 26.4 7.0 0.264
17 2885 132 63.9 0.484
18 2710 155.1 78.9 0.509
19 2550 194.7 139.0 0.714
20 2400 96 40.9 0.426
21 2180 185.8 118.5 0.638
22 2080 59.4 42.4 0.714
23 2020 6.6 2.0 0.303
24 1920 29.7 14.5 0.488
25 1875 13.2 7.4 0.561
26 3420 221.1 107.0 0.484
27 3120 112.2 54.3 0.484
28 2760 39.6 20.0 0.505
29 2590 76.9 39.1 0.509
30 2500 56.1 40.1 0.714
31 2460 29.7 8.5 0.286
(五)判识是否存在油气逸散：当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，则该测点存在油气逸散，DL1地区31个测点的判识结果详见表5。
表5 DL1地区营三段31个测点油气逸散的判识结果
由上述判识结果可见，在DL1地区营三段的1-5、7-9、12-13、17-18、22-30号测点均存在油气逸散。
(6)上述判识结果的验证：
根据现今的勘探成果可知，DL地区整体上油气富集程度较低，产量较低，勘探效果并不理想。另外，在上述确定的油气逸散位置，油气成藏规模较小甚至几乎没有油气，附近的S2、SN100、DL100、DL85等井位在该层段内仅有油气显示，通过镜下观察包裹体，存在烃类包裹体，且通过测试包裹体均一温度可知，其烃类包裹体形成于两期，但是包裹体数量较少，GOI值均小于1％。另外，从勘探现状来看，这几口井产量极低，甚至只有油花，未形成规模储量，为油气逸散区。因而，本发明提供的构造破坏导致的油气逸散区的判识方法的判识结果与实际勘探结果具有较高的一致性，从而验证了本发明提供的构造破坏导致的油气逸散区的判识方法的可行性与准确性。

Claims (6)

1.一种构造破坏导致的油气逸散区的判识方法，其特征在于，包括如下步骤：
选取测点：在研究区的目的层内，选取均匀覆盖于研究区范围内断层的多个测点；
计算断层泥岩比率下限值：以研究区内具有实测岩石排替压力数据的检测点为样品点，根据样品点的实测数据，利用公式(1)拟合研究区内岩石排替压力与岩石的埋深和泥质含量乘积之间的函数关系，获得常数a和b；再根据研究区的录井、钻井及地震资料，利用公式(2)计算每一测点对应的断层泥岩比率下限值，公式(1)和(2)的表达式如下：
其中，Pd样为样品点的岩石排替压力，单位为MPa；a、b均为拟合常数；Z样为样品点的岩石埋深，单位为m；V样为样品点的岩石泥质含量，单位为％；SGRlim为测点的断层泥岩比率下限值；Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；ρw为地层水密度，单位为g/cm3；g为重力加速度，单位为m/s2；θ为测点所在断裂的断裂倾角，单位为°；V储为测点对应储层的泥质含量，单位为％；
计算成藏期断层泥岩比率：根据研究区的录井、钻井及地震资料，利用公式(3)计算每一测点对应的成藏期断层泥岩比率，公式(3)的表达式如下：
其中，SGR古为测点的成藏期断层泥岩比率；d古为测点所在断裂在成藏期时的古断距，单位为m；H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m；
计算现今断层泥岩比率：根据研究区录井数据，利用公式(4)计算每一测点对应的现今断层泥岩比率，公式(4)的表达式如下：
其中，SGR现今为测点的现今断层泥岩比率；d现今为测点所在断裂的现今断距，单位为m；H为测点所在断裂的现今累积泥岩厚度，单位为m；
判识是否存在油气逸散：当测点对应的SGR古<SGRlim且SGR现今<SGRlim时，则该测点存在油气逸散。
2.根据权利要求1所述的构造破坏导致的油气逸散区的判识方法，其特征在于，计算断层泥岩比率下限值的步骤中，根据研究区的三维地震资料和钻井资料，利用公式(5)计算获得公式(2)中的参数Z断'，公式(5)的表达式如下：
其中，Z断'为测点所在断裂的断层岩开始成岩时的埋深，单位为m；Z为测点埋深，单位为m；T断为测点所在断裂的断层岩开始成岩的时间，单位为Ma；T储为测点对应储层开始成岩的时间，单位为Ma；ρr为上覆沉积岩石密度，单位为g/cm3；ρw为地层水密度，单位为g/cm3；θ为测点所在断裂的断裂倾角，单位为°。
3.根据权利要求1所述的构造破坏导致的油气逸散区的判识方法，其特征在于，计算断层泥岩比率下限值的步骤中，公式(2)中参数V储的具体计算步骤如下：根据靠近测点的已知井的录井和GR测井数据，利用公式(6)和公式(7)计算获得V储，公式(6)和公式(7)的表达式如下：
其中，Vsh为测点对应储层的泥质含量指数，单位为％；GR目的层为靠近该测点的已知井在目的层的GR测井曲线值，单位为API；GR纯泥岩层为靠近该测点的已知井的纯泥岩层的GR测井曲线值，单位为API；GR纯砂岩层为靠近该测点的已知井的纯砂岩层的GR测井曲线值，单位为API；V储为测点对应储层的泥质含量，单位为％；GCUR为与地层有关的经验系数，古近系和新近系地层的GCUR取值为3.7，老地层的GCUR取值为2。
4.根据权利要求3所述的构造破坏导致的油气逸散区的判识方法，其特征在于，公式(6)中的参数GR目的层、GR纯泥岩层和GR纯砂岩层的具体获取步骤如下：根据目的层的深度，以所述已知井的GR测井曲线在相同深度对应的GR测井值作为GR目的层，选取所述已知井中大套泥岩的GR测井最高值作为GR纯泥岩层，选取所述已知井中大套砂岩的GR测井最低值作为GR纯砂岩层。
5.根据权利要求1所述的构造破坏导致的油气逸散区的判识方法，其特征在于，计算成藏期断层泥岩比率的步骤中，公式(3)中的参数d古的具体获取步骤如下：切取垂直于该断裂不同部位的地震剖面，读取该断裂成藏期时不同部位的断距，取其中的最大断距作为该断裂自成藏期至现今的活动量，将该断裂的现今断距与其自成藏期至现今的活动量作差，得到该断裂在成藏期时的古断距d古。
6.根据权利要求5所述的构造破坏导致的油气逸散区的判识方法，其特征在于，计算成藏期断层泥岩比率的步骤中，公式(3)中参数H0的具体计算步骤如下：将目的层中在古断距d古内沉积的地层等分为多个统计单元，分别统计每个统计单元对应的现今泥岩累积厚度，根据去压实校正法利用公式(8)计算获得H0，公式(8)的表达式如下：
其中，H0为测点所在断裂在古断距d古范围内的未受压实作用的累积古泥岩厚度，单位为m；hi为第i个统计单元对应的现今泥岩累积厚度，单位为m；Ci为第i个统计单元对应的泥岩压实系数，利用去压实校正法获得；α为地层倾角，单位为°；n为统计单元的总数。
CN201910521732.8A 2019-06-17 2019-06-17 一种构造破坏导致的油气逸散区的判识方法 Active CN110244358B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910521732.8A CN110244358B (zh) 2019-06-17 2019-06-17 一种构造破坏导致的油气逸散区的判识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910521732.8A CN110244358B (zh) 2019-06-17 2019-06-17 一种构造破坏导致的油气逸散区的判识方法

Publications (2)

Publication Number Publication Date
CN110244358A true CN110244358A (zh) 2019-09-17
CN110244358B CN110244358B (zh) 2021-04-09

Family

ID=67887720

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910521732.8A Active CN110244358B (zh) 2019-06-17 2019-06-17 一种构造破坏导致的油气逸散区的判识方法

Country Status (1)

Country Link
CN (1) CN110244358B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946774A (zh) * 2021-01-08 2021-06-11 中国石油大学(华东) 一种以单层泥岩为计算单元的断层侧向封闭性评价方法
CN112943229A (zh) * 2021-02-04 2021-06-11 大庆油田有限责任公司 一种储气库盖层突破压力连续预测方法

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012045001A2 (en) * 2010-09-30 2012-04-05 Vanderbilt University Influenza virus antibodies and immunogens and uses therefor
CN104111483A (zh) * 2014-03-26 2014-10-22 中国石油化工股份有限公司 一种断陷湖盆斜坡带油藏分布的勘探方法
CN106650042A (zh) * 2016-12-01 2017-05-10 中国石油天然气股份有限公司 断裂输导能力参数的确定方法和装置
KR101853334B1 (ko) * 2016-04-25 2018-04-30 충북대학교 산학협력단 순차층서 분석을 통한 비전통자원 개발 구간 예측 방법
CN108241181A (zh) * 2018-01-30 2018-07-03 中国石油化工股份有限公司 一种断层封闭性的评价方法
CN109239778A (zh) * 2018-09-14 2019-01-18 中国石油大学(华东) 一种断层侧向封闭性定量评价方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012045001A2 (en) * 2010-09-30 2012-04-05 Vanderbilt University Influenza virus antibodies and immunogens and uses therefor
CN104111483A (zh) * 2014-03-26 2014-10-22 中国石油化工股份有限公司 一种断陷湖盆斜坡带油藏分布的勘探方法
KR101853334B1 (ko) * 2016-04-25 2018-04-30 충북대학교 산학협력단 순차층서 분석을 통한 비전통자원 개발 구간 예측 방법
CN106650042A (zh) * 2016-12-01 2017-05-10 中国石油天然气股份有限公司 断裂输导能力参数的确定方法和装置
CN108241181A (zh) * 2018-01-30 2018-07-03 中国石油化工股份有限公司 一种断层封闭性的评价方法
CN109239778A (zh) * 2018-09-14 2019-01-18 中国石油大学(华东) 一种断层侧向封闭性定量评价方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
王超 等: "油气成藏期断层古侧向封闭性预测方法及其应用", 《地球科学》 *
范婕 等: "恢复不整合剥蚀厚度的旋回分析法", 《中国矿业大学学报》 *
范婕 等: "长岭断陷龙凤山地区断裂与油气运聚的关系", 《地球科学》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112946774A (zh) * 2021-01-08 2021-06-11 中国石油大学(华东) 一种以单层泥岩为计算单元的断层侧向封闭性评价方法
CN112946774B (zh) * 2021-01-08 2022-02-11 中国石油大学(华东) 一种以单层泥岩为计算单元的断层侧向封闭性评价方法
CN112943229A (zh) * 2021-02-04 2021-06-11 大庆油田有限责任公司 一种储气库盖层突破压力连续预测方法

Also Published As

Publication number Publication date
CN110244358B (zh) 2021-04-09

Similar Documents

Publication Publication Date Title
CN109061765A (zh) 非均质薄砂岩互层油藏的圈闭评价方法
CN109653725A (zh) 基于沉积微相和岩石相的混积储层水淹程度测井解释方法
McCulloh Mass properties of sedimentary rocks and gravimetric effects of petroleum and natural-gas reservoirs
CN110056346B (zh) 一种基于趋势变化函数的油藏三维原始含水饱和度模拟方法
CN104514552A (zh) 一种煤层气储层识别与丰度预测的方法
Abdideh et al. Cluster analysis of petrophysical and geological parameters for separating the electrofacies of a gas carbonate reservoir sequence
Wimmers et al. Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir
CN114910964B (zh) 一种断陷湖盆陡坡带砂砾岩体甜点区的预测方法
CN109458176A (zh) 碳酸盐岩储层压力的预测方法及其应用
Esmaeili et al. Developing a saturation-height function for reservoir rock types and comparing the results with the well log-derived water saturation, a case study from the Fahliyan formation, Dorood oilfield, Southwest of Iran
CN110244358A (zh) 一种构造破坏导致的油气逸散区的判识方法
CN111206921A (zh) 一种适用于火山岩溢流相有利储层的描述方法
Xiong et al. Petrophysical evaluation of gas hydrate in Shenhu area, China
Jasim et al. Specifying quality of a tight oil reservoir through 3-d reservoir modeling
CN103698494B (zh) 一种确定岩性圈闭含油气饱和度的方法及装置
Davies Permeability Modelling of a Sandstone Reservoir in Parts of the Niger Delta
CN110244356A (zh) 一种构造破坏形成的后期油气藏的判识方法
Rahman et al. Present gas reserve estimation using wireline logging data of Habiganj Gas Field, Bangladesh
CN110244357A (zh) 一种构造破坏型古油藏的判识方法
Oluwajana Petrophysical Analysis of Reservoirs in Well-Y, Coastal Swamp Depobelt, Niger Delta Basin, Nigeria.
Sokhal et al. An Example from the Algerian Sahara Illustrates the Use of the Hydraulic Flow Unit Technique to Discriminate Fluid Flow Routes in Confined Sand Reservoirs
CN110297280A (zh) 一种碳酸盐岩超压空间分布特征的预测方法
Hamdulla et al. Building geological model for tertiary reservoir of exploration Ismail oil field, North Iraq
Anyiam et al. Formation evaluation of an onshore appraisal well ‘KG-5’,“green field”, Niger Delta, Nigeria
RU2225020C1 (ru) Способ геофизической разведки для определения нефтепродуктивности трещинных глинистых коллекторов в межскважинном пространстве

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
CB03 Change of inventor or designer information

Inventor after: Fan Jie

Inventor after: Ma Ji

Inventor after: Wang Jinduo

Inventor after: Zeng Zhiping

Inventor after: Mao Minghai

Inventor after: Li Songtao

Inventor after: Zhao Leqiang

Inventor after: Zhang Yi

Inventor after: Zhou Tao

Inventor after: Gao Jianxiong

Inventor before: Fan Jie

Inventor before: Xue Zelei

Inventor before: Lin Huixi

Inventor before: Qin Feng

Inventor before: Zhou Tao

Inventor before: Mao Minghai

Inventor before: Li Songtao

Inventor before: Liu Dezhi

Inventor before: Zhang Yi

Inventor before: Li Ling

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant