CN110244357A - 一种构造破坏型古油藏的判识方法 - Google Patents

一种构造破坏型古油藏的判识方法 Download PDF

Info

Publication number
CN110244357A
CN110244357A CN201910521719.2A CN201910521719A CN110244357A CN 110244357 A CN110244357 A CN 110244357A CN 201910521719 A CN201910521719 A CN 201910521719A CN 110244357 A CN110244357 A CN 110244357A
Authority
CN
China
Prior art keywords
measuring point
unit
formula
pool
mud stone
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
CN201910521719.2A
Other languages
English (en)
Other versions
CN110244357B (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 CN201910521719.2A priority Critical patent/CN110244357B/zh
Publication of CN110244357A publication Critical patent/CN110244357A/zh
Application granted granted Critical
Publication of CN110244357B publication Critical patent/CN110244357B/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. analysis, for interpretation, for correction
    • 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)
  • General Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (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-pln Hi (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号测点所在断裂未受压实作用的累积古泥岩厚度的计算参数及结果
统计单元编号 顶深/m 底深/m H<sub>i</sub>/m C<sub>i</sub> h<sub>i</sub>/m
1 2800 2820 2810 0.666 12
2 2820 2840 2830 0.665 7
3 2840 2860 2850 0.664 6
4 2860 2880 2870 0.664 5.5
5 2880 2900 2890 0.663 5.5
6 2900 2920 2910 0.662 7
7 2920 2940 2930 0.662 2
8 2940 2960 2950 0.661 9.5
9 2960 2969.6 2964.8 0.661 2.5
(四)计算现今断层泥岩比率：根据DL1地区的录井数据，利用公式(4)计算每一测点对应的现今断层泥岩比率，计算结果如表4所示。
表4 DL1地区营三段31个测点对应的SGR现今的计算参数及结果
(五)判识是否存在古油藏：当测点对应的SGR古>SGRlim且SGR现今<SGRlim时，则该测点存在古油藏，DL1地区31个测点的判识结果详见表5。
表5 DL1地区营三段31个测点古油藏的判识结果
由上述判识结果可见，在DL1地区营三段的6号测点、10号测点、15号测点、16号测点、31号测点可能存在古油藏，靠近这些测点的已知井包括DL09井、DL82井、S8井、DL80井、S1井，上述这些井位为存在古油藏的位置。
(六)采用多种方法验证上述判识结果：
(1)岩心观察法
岩心观察显示，DL09井营城组和登娄库组存在大量沥青，多数黑色的固态沥青充填于砂岩孔隙或构造裂缝中，其形态主要取决于孔、缝的形态，孤立分散，边缘有氧化环，新鲜切面散发浓烈的油气味。沥青的存在表明油气途经此处，可能曾经在该处成藏；裂缝中的油迹表明油气沿断裂发生运移。这两类现象为“古油藏遭受断裂破坏，并沿断裂运移至浅部位形成次生油气藏”的观点提供了系列证据。
(2)镜下沥青及沥青包裹体观察法
镜下观察表明，透射光下，沥青呈黑棕色或浅棕色，形态呈分散粒状、弥散状分布于颗粒表面、粒间孔隙或裂缝中，形态不规则，多数无明显的边缘，常与石英、长石等颗粒伴生。荧光下，多数沥青发绿褐色与绿色荧光，极少数不发荧光，表明其演化程度相对较低。上述特征与生物降解沥青的特征较为一致，说明研究区营城组存在古油藏，且是由于构造破坏油气的保存条件所形成的。另外，研究区营城组含沥青的烃类包裹体大量发育，主要分布于裂缝中，成条带状展布，部分沥青成膜状附着于包裹体壁上，透射光下呈现出明显的黑色边缘，荧光下呈现蓝绿色荧光，说明沥青演化程度相对较低，非高温裂解产物，证明曾经早期发生油气充注，后因地层抬升剥蚀，导致古油藏遭受破坏。
(3)GOI法
选取钻遇断裂的井位，选取断裂钻遇点上、下位置取样，进行包裹体观察，并计算GOI值，同时与测井、试油等其他类型资料综合判断古油藏是否存在。结果表明，DL09井GOI值大于5％，录井资料有油气显示，测井数据显示有差气层，但通过试油表明其仅为干层，并无油气产量，表明此处应为古油藏发育位置。据此亦可判断出S8、DL80井营城组存在古油藏。综上，该方法认为DL09井、S8井、DL80井的营城组是存在古油藏的，与本发明提供的构造破坏型古油藏的判识方法得到的古油藏位置具有一致性。
(4)地层水及原油物性对比法
对于地层水矿化度而言，其数值受地层埋深影响较大，随着地层埋深的增加而逐渐增大。对比结果表明，S1井的泉二段地层水矿化度明显比同层系邻井的矿化度高，存在异常高点，证明可能是在营城组由断层运移上来的，进而形成了次生油气藏。另外，原油物性的统计结果表明，S1和DL82井的泉头组原油密度、粘度均呈现出异常低值，指示其营城组可能存在古油藏。
利用上述多种方法对判识结果进行验证可见，本发明提供的构造破坏型古油藏的判识方法的判识结果与实钻结果具有很好的一致性，从而验证了本发明提供的构造破坏型古油藏的判识方法的可行性与准确性。

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为统计单元的总数。
CN201910521719.2A 2019-06-17 2019-06-17 一种构造破坏型古油藏的判识方法 Active CN110244357B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910521719.2A CN110244357B (zh) 2019-06-17 2019-06-17 一种构造破坏型古油藏的判识方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910521719.2A CN110244357B (zh) 2019-06-17 2019-06-17 一种构造破坏型古油藏的判识方法

Publications (2)

Publication Number Publication Date
CN110244357A true CN110244357A (zh) 2019-09-17
CN110244357B CN110244357B (zh) 2021-04-09

Family

ID=67887717

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910521719.2A Active CN110244357B (zh) 2019-06-17 2019-06-17 一种构造破坏型古油藏的判识方法

Country Status (1)

Country Link
CN (1) CN110244357B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113887078A (zh) * 2021-10-26 2022-01-04 中国科学院地质与地球物理研究所 深层-超深层海相碳酸盐岩油藏赋存深度下限预测方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113887078A (zh) * 2021-10-26 2022-01-04 中国科学院地质与地球物理研究所 深层-超深层海相碳酸盐岩油藏赋存深度下限预测方法
CN113887078B (zh) * 2021-10-26 2022-03-01 中国科学院地质与地球物理研究所 深层-超深层海相碳酸盐岩油藏赋存深度下限预测方法
US11435494B1 (en) 2021-10-26 2022-09-06 Institute Of Geology And Geophysics, Chinese Academy Of Sciences Method for predicting oil accumulation depth limit of deep and ultra-deep marine carbonate reservoirs

Also Published As

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

Similar Documents

Publication Publication Date Title
Tian et al. Multi-layered Ordovician paleokarst reservoir detection and spatial delineation: A case study in the Tahe Oilfield, Tarim Basin, Western China
CN101936165B (zh) 一种岩溶型碳酸盐岩储层布井方法
Longde et al. Petroleum exploration and development practices of sedimentary basins in China and research progress of sedimentology
Wu et al. Permeability distribution and scaling in multi-stages carbonate damage zones: Insight from strike-slip fault zones in the Tarim Basin, NW China
CN104047597A (zh) 含油气泥页岩地层测井曲线标准化方法
CN109388817A (zh) 一种储层裂缝三维建模方法
CN106503834A (zh) 一种湖相超低孔渗储层致密油甜点区的预测方法
CN105425315B (zh) 一种弱构造变形区小尺度沉积期古地貌反演方法
CN104632201A (zh) 一种页岩油气“甜点区”有机碳含量下限值的确定方法
Tinker et al. Multiple karst events related to stratigraphic cyclicity: San Andres Formation, Yates field, west Texas
Wimmers et al. Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir
Barr et al. Pre-development fracture modelling in the Clair field, west of Shetland
CN109458176A (zh) 碳酸盐岩储层压力的预测方法及其应用
Han et al. Fine-grained rock fabric facies classification and its control on shale oil accumulation: a case study from the Paleogene Kong 2 Member, Bohai Bay Basin
Zhang et al. Silurian hydrocarbon exploration breakthrough and its implications in the Shajingzi structural belt of Tarim Basin, NW China
Liu et al. The Control Theory and Application for Well Pattern Optimization of Heterogeneous Sandstone Reservoirs
Gillespie et al. Introduction: geology of fractured reservoirs
CN110244358A (zh) 一种构造破坏导致的油气逸散区的判识方法
CN110244357A (zh) 一种构造破坏型古油藏的判识方法
Jasim et al. Specifying quality of a tight oil reservoir through 3-d reservoir modeling
CN110244356A (zh) 一种构造破坏形成的后期油气藏的判识方法
Berch et al. Predicting potential unconventional production in the Tuscaloosa marine shale play using thermal modeling and log overlay analysis
Olson et al. White River Dome field: Gas production from deep coals and sandstones of the Cretaceous Williams Fork Formation
Rahman et al. Present gas reserve estimation using wireline logging data of Habiganj Gas Field, Bangladesh
Harvey et al. Shale prospectivity onshore Britain

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

Inventor after: Song Mingshui

Inventor after: Zhou Tao

Inventor after: Fan Jie

Inventor after: Lin Huixi

Inventor after: Li Guodong

Inventor after: Gao Jianxiong

Inventor after: Gong Yajun

Inventor after: Niu Jingjing

Inventor after: Zhou Huilian

Inventor after: Qin Feng

Inventor before: Song Mingshui

Inventor before: Zhou Huilian

Inventor before: Fan Jie

Inventor before: Lin Huixi

Inventor before: Li Guodong

Inventor before: Qin Feng

Inventor before: Gao Jianxiong

Inventor before: Zhou Tao

Inventor before: Gong Yajun

Inventor before: Niu Jingjing

GR01 Patent grant
GR01 Patent grant