CN109324345A - 叠合盆地油气成藏期岩石孔隙度恢复方法 - Google Patents
叠合盆地油气成藏期岩石孔隙度恢复方法 Download PDFInfo
- Publication number
- CN109324345A CN109324345A CN201811463682.4A CN201811463682A CN109324345A CN 109324345 A CN109324345 A CN 109324345A CN 201811463682 A CN201811463682 A CN 201811463682A CN 109324345 A CN109324345 A CN 109324345A
- Authority
- CN
- China
- Prior art keywords
- porosity
- formula
- oil
- model
- pool
- 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
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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
- G01V1/50—Analysing data
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (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,建立面孔率‑孔隙度模型;步骤6,建立埋藏过程中的有机酸溶蚀孔隙度模型;步骤7,根据步骤4建立的压实孔隙度模型及步骤6建立的有机酸溶蚀孔隙度模型,建立成藏期孔隙度模型。该叠合盆地油气成藏期岩石孔隙度恢复方法丰富了多期次成藏过程中溶蚀孔隙度的定量评价,可广泛应用于叠合盆地多期次成藏过程中的孔隙度评价,对油气勘探具有重要的指导意义。
Description
技术领域
本发明涉及油气地质勘探技术领域,特别是涉及到一种叠合盆地油气成藏期岩石孔隙度恢复方法。
背景技术
低渗透油气资源在我国油气资源中占有十分重要的地位,也是未来中国油气勘探的主要对象。低渗透储层的形成主要受沉积作用和后期成岩作用的影响,沉积作用主控储层物性,而成岩作用中的压实和胶结作用则对储层物性起破坏作用,溶蚀作用则能够改善储层物性。现有储层评价体系多属于静态评价,真正决定油气是否大规模进入储层的成藏关键期的储层物性则是进行储层物性评价的关键。
在低勘探程度区,钻井少,分析数据少等制约了对区内古孔隙度的恢复,仅通过少量井的测井数据、普通薄片、铸体薄片数据等有限的资料对储层古孔隙度进行分析是解决勘探程度低的地区物性分析的一种较为有效的手段。成藏关键时期储层物性很大程度上决定了油气能否大规模进入储层,先充注后致密或者先致密后充注在勘探开发潜力上存在较大差异,后期工程上需要采取的措施也千差万别。
众多学者已对储层孔隙演化情况开展过研究,他们利用储层物性及铸体薄片等资料对成藏期储层恢复进行了大量的尝试,并且在孔隙度恢复过程中,多考虑压实作用及胶结作用等对孔隙度演化的影响,而忽略了地层水中有机酸对孔隙的溶蚀作用对储层物性的影响。
目前恢复储层孔隙演化的研究方法主要有压实趋势法和反演回剥法。其中,第一种:压实趋势法,通过利用测井解释孔隙度建立储层现今孔隙度-深度关系曲线,在明确成藏期古埋深的基础上,根据现今孔-深关系,读取储层在该埋深下对应的孔隙度,近似认为是储层在主成藏期的孔隙度。压实趋势法的缺点是:没有考虑其他因素的影响,且在孔-深关系上,同一深度下孔隙度分布范围较大,给古孔隙度的恢复造成了困难。第二种:反演回剥法,步骤包括:①计算原始孔隙度;②建立面孔率与孔隙度函数关系,线性拟合方法;③定量统计胶结作用对面孔率的贡献量,采用“反演回剥”原理恢复不同阶段的孔隙度;④根据面孔率反演结果,恢复各主要胶结作用时期储集层孔隙结构。反演回剥法的缺点是:该方法在实际应用中具有很大的局限性,首先,这一方法工作量较大,所要求的资料较为多样,仅适合于勘探程度相对较高的地区;其次,该方法对于储层也有一定的要求,如储层的杂基含量相对较低,裂缝应不发育,且粒径要适中;第三,这一方法目前仍存在一些难点及问题,如在多旋回成岩作用类型的期次区分上,面孔率-孔隙度转化问题上,粘土矿物胶结物对储层物性的影响如何计算等。为此我们发明了一种新的叠合盆地油气成藏期岩石孔隙度恢复方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种可完成在多次油气充注过程中有机酸对孔隙的溶蚀作用的定量化,恢复成藏期时油气孔隙度的叠合盆地油气成藏期岩石孔隙度恢复方法。
本发明的目的可通过如下技术措施来实现:叠合盆地油气成藏期岩石孔隙度恢复方法,该叠合盆地油气成藏期岩石孔隙度恢复方法包括:步骤1,进行地层剥蚀量的恢复;步骤2,建立热演化史模型;步骤3,确定区内目的层油气成藏期次及各个成藏期的地层温度;步骤4,建立埋藏过程中的压实孔隙度模型;步骤5,建立面孔率-孔隙度模型;步骤6,建立埋藏过程中的有机酸溶蚀孔隙度模型;步骤7,根据步骤4建立的压实孔隙度模型及步骤6建立的有机酸溶蚀孔隙度模型,建立成藏期孔隙度模型。
本发明的目的还可通过如下技术措施来实现:
在步骤1,利用声波时差曲线,挑选厚度不小于1.5m的泥岩段对应的声波时差曲线进行曲线拟合,进行剥蚀量恢复,进而恢复不同地质历史时期地层的剥蚀量。
步骤1包括:
步骤a1,选取测井曲线:统计区内探井的测井曲线,剔除非正常缩径、扩径段以及曲线上非正常跳跃段数值;
步骤a2,采集合适的测井曲线数据:选取步骤a1中纯泥岩段厚度不小于1.5m的声波时差曲线,读取其中间值;
步骤a3,剥蚀量恢复:首先,对步骤a2所得数据进行回归拟合,横坐标为声波时差数据的对数坐标,纵坐标为深度:
y=a1·x+b 公式(1),
式中,a1、b为拟合系数;x为声波时差曲线值,μs/m;y为深度,m;
其次,将压实曲线外延,与对数横坐标相交于600μs/m,则该点的水平延伸线即为古地表H0,由古地表至现今残留目的层顶面即为地层剥蚀量H:
H=H1-H0 公式(2),
式中,H0为古地表,m;H1为古地表到现今残留目的层顶面的厚度,m;H为地层剥蚀量,m。
在步骤2,通过步骤1确定的不同地层的剥蚀量,恢复古埋深,并结合古地温梯度数据,建立热演化史模型。
步骤2包括:
步骤b1,埋藏史的建立:①整理单井岩性剖面的数据,包括井上的地层分层及对应的地质年代数据;②通过步骤1确定不同地层的剥蚀量以及地层遭受剥蚀的时间;③建立直角坐标系,横坐标表示地质时间,纵坐标表示叠加地层的深度,按照地质时间由新到老进行地层及剥蚀量的叠加;④将同一地层底界相连则为该地层整个地质历史时期的古埋深,获得埋藏史图,为后期的热演化史的建立奠定基础;
步骤b2,热演化史模型的建立:首先收集研究区的古地温梯度数据,地表温度默认为20℃,其次,由于地温梯度是深度每增加100m温度的增加值,在已知地表温度为20℃及步骤b1,得到的古埋深的条件下,计算地质历史时期某一深度的地层温度,将相同温度相连即可建立热演化史模型。
在步骤3,利用流体包裹体确定成藏期次及各个期次的最大峰温值,进而在步骤2建立的生烃演化史模型的基础上确定成藏时间,并确定其古埋深。
步骤3包括:
步骤c1,单井模拟:统计单井基础资料,包括井深h、分层、古地温梯度数据、地层剥蚀量H,进行热演化史的恢复;
步骤c2,成藏期的确定:根据步骤c1的单井模拟演化结果,在生烃层中与均一温度相同的古地温线的第一个交点所对应的时间,即为油气的成藏期n;
步骤c3,流体包裹体均一温度的测定:选取新鲜样品,通过制片、岩石学观察及冷热台测试,确定第n个成藏期对应的地层流体包裹体的均一温度的主峰值Tn。
在步骤4,在明确原始孔隙度和古埋深的基础上,建立压实孔隙度模型,获得不同时期的压实孔隙度。
步骤4包括:
步骤d1,原始孔隙度计算
利用原始孔隙度φ0与分选系数S0之间的经验函数关系式:
φ0=20.91+22.9/S0 公式(3),
其中,S0=(P25/P75)1/2 公式(4),
式中,S0为分选系数;φ0为原始孔隙度,%;P25,P75分别代表粒度累计曲线上颗粒含量在25%和75%处所对应的颗粒直径,μm;
步骤d2,利用机械压实孔隙度与埋深的指数关系,来恢复经历压实作用后成藏期末的孔隙度,即:
φs=φ0·e-0.0002·h 公式(5),
其中,φs为剩余孔隙度,%;φ0为恢复砂岩的原始孔隙度,%;h为埋深,m;
步骤d3,根据Tn在埋藏史图上标出目的层在成藏期的埋深,将埋深代入公式(5),得到成藏期的压实孔隙度,
在步骤5,根据薄片、铸体薄片观察开展镜下溶蚀孔识别及划分,结合物性分析结果,统计溶蚀孔的面孔率与孔隙度的关系,建立面孔率-孔隙度模型。
在步骤6,利用温度与浓度的差异确定有机酸充注系数,进而利用步骤5的现今孔隙度与现今溶蚀面孔率之间的关系确定不同成藏期次的溶蚀孔隙度模型。
步骤6包括:
步骤f1,利用面孔率和孔隙度之间的关系,建立线性关系模型如下:
φ=k·S+k1 公式(7),
式中,φ为孔隙度,%;S为溶蚀面孔率,%;k、k1为拟合系数;
利用井上现今孔隙度和现今溶蚀面孔率确定k、k1这两个拟合系数;
φf=k·Sf+k1 公式(8),
式中,φf为现今孔隙度,%;Sf为现今溶蚀面孔率,%;
步骤f2,地层水中的有机酸和CO2是储层发生溶蚀、形成次生孔隙的重要因素,而这些有机酸和CO2又是烃源岩在成熟过程中的产物,有机酸浓度随温度变化曲线近似于抛物线,因此建立有机酸浓度与温度的关系:
C=a2·T2+b·T+c1,
式中,C为有机酸充注高峰浓度,mg/L;a2、b、c1为常数,T为有机酸充注高峰时的地层温度;
通过模拟实验,建立温度与有机酸分布曲线如下:
C=-0.00385·T9+1.39+·T-9.54 公式(9);
同等条件下,有机酸浓度与其溶蚀面孔率成正比,则:
Sn=m·Sf 公式(10);
m=Cn/(C1+C2+…+Cn+…+Cz) 公式(11);
式中,第n个原油充注高峰时溶蚀孔隙的面孔率为Sn,%;m为第n期成藏时的有机酸充注系数;
将Sn所得值代入公式(7)中,得到该期充注的溶蚀孔隙度,
φn=k·m·Sf+k1 公式(12);
式中,φn为第n次成藏时有机酸溶蚀孔隙度,%。
在步骤7,单个成藏期孔隙度模型为地层压实孔隙度+有机酸溶蚀孔隙度,即:
φ成藏期=φs+φn 公式(13);
式中,φs为剩余孔隙度,%;φn为第n次成藏时有机酸溶蚀孔隙度,%。
本发明中的叠合盆地油气成藏期岩石孔隙度恢复方法,提出了成藏期中油气充注时有机酸溶蚀孔隙度确定的量化方法,丰富了多期次成藏过程中溶蚀孔隙度的定量评价方法,可广泛应用于叠合盆地多期次成藏过程中的孔隙度评价,对油气勘探具有重要的指导意义。
附图说明
图1为本发明的叠合盆地油气成藏期岩石孔隙度恢复方法的一具体实施例的流程图;
图2为本发明的一具体实施例中恢复的生烃演化史图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图1所示,图1为本发明的叠合盆地油气成藏期岩石孔隙度恢复方法的流程图。
步骤101,地层剥蚀量的恢复
由于构造运动的影响,地层在抬升过程中会发生剥蚀,因此,需要对地层进行剥蚀量的恢复。利用声波时差曲线,挑选厚度不小于1.5m的泥岩段对应的声波时差曲线进行曲线拟合,进行剥蚀量恢复,进而恢复不同地质历史时期地层的剥蚀量;
其中,步骤101的具体步骤如下:
步骤a1,选取测井曲线:统计区内探井的测井曲线,剔除非正常缩径、扩径段以及曲线上非正常跳跃段数值;
步骤a2,采集合适的测井曲线数据:选取步骤a1中纯泥岩段厚度不小于1.5m的声波时差曲线,读取其中间值;
步骤a3,剥蚀量恢复:首先,对步骤a2所得数据进行回归拟合,横坐标为声波时差数据的对数坐标,纵坐标为深度:
y=a1·x+b 公式(1),
式中,a1、b为拟合系数;x为声波时差曲线值,μs/m;y为深度,m;
其次,将压实曲线外延,与对数横坐标相交于600μs/m,则该点的水平延伸线即为古地表H0,由古地表至现今残留目的层顶面即为地层剥蚀量H:
H=H1-H0 公式(2),
式中,H0为古地表,m;H1为古地表到现今残留目的层顶面的厚度,m;H为地层剥蚀量,m。
步骤102,建立热演化史模型
通过步骤101确定的不同地层的剥蚀量,恢复古埋深,并结合古地温梯度数据,建立热演化史模型;
其中,步骤102的具体步骤如下:
建立热演化史模型前先建立埋藏史;
步骤b1,埋藏史的建立:①整理单井岩性剖面的数据,包括井上的地层分层及对应的地质年代数据;②通过步骤101确定不同地层的剥蚀量以及地层遭受剥蚀的时间;③建立直角坐标系,横坐标表示地质时间,纵坐标表示叠加地层的深度,按照地质时间由新到老进行地层及剥蚀量的叠加;④将同一地层底界相连则为该地层整个地质历史时期的古埋深,获得埋藏史图,为后期的热演化史的建立奠定基础;
步骤b2,热演化史模型的建立:首先收集研究区的古地温梯度数据,地表温度默认为20℃,其次,由于地温梯度是深度每增加100m温度的增加值,在已知地表温度为20℃及步骤b1,得到的古埋深的条件下,计算地质历史时期某一深度的地层温度,将相同温度相连即可建立热演化史模型。
步骤103,确定区内目的层油气成藏期次及各个成藏期的地层温度
利用流体包裹体确定成藏期次及各个期次的最大峰温值,进而在步骤102建立的生烃演化史模型的基础上确定成藏时间,并确定其古埋深;
其中,步骤103的具体步骤如下:
步骤c1,单井模拟:统计单井基础资料,包括井深h、分层、古地温梯度数据、地层剥蚀量H,进行热演化史的恢复;
步骤c2,成藏期的确定:根据步骤c1的单井模拟演化结果,在生烃层中与均一温度相同的古地温线的第一个交点所对应的时间,即为油气的成藏期n;
步骤c3,流体包裹体均一温度的测定:选取新鲜样品,通过制片、岩石学观察及冷热台测试,确定第n个成藏期对应的地层流体包裹体的均一温度的主峰值Tn。
步骤104,建立埋藏过程中的压实孔隙度模型
在明确原始孔隙度和古埋深的基础上,建立压实孔隙度模型,获得不同时期的压实孔隙度;
其中,步骤104的具体步骤如下:
步骤d1,原始孔隙度计算
利用原始孔隙度φ0与分选系数S0之间的经验函数关系式:
φ0=20.91+22.9/S0 公式(3),
其中,S0=(P25/P75)1/2 公式(4),
式中,S0为分选系数;φ0为原始孔隙度,%;P25,P75分别代表粒度累计曲线上颗粒含量在25%和75%处所对应的颗粒直径,μm;
步骤d2,利用机械压实孔隙度与埋深的指数关系,来恢复经历压实作用后成藏期末的孔隙度,即:
φs=φ0·e-0.0002·h 公式(5),
其中,φs为剩余孔隙度,%;φ0为恢复砂岩的原始孔隙度,%;h为埋深,m;
步骤d3,根据Tn在埋藏史图上标出目的层在成藏期的埋深,将埋深代入公式(5),得到成藏期的压实孔隙度,
步骤105,建立面孔率-孔隙度模型
根据薄片、铸体薄片观察开展镜下溶蚀孔识别及划分,结合物性分析结果,统计溶蚀孔的面孔率与孔隙度的关系,建立面孔率-孔隙度模型;
步骤106,建立埋藏过程中的有机酸溶蚀孔隙度模型
利用温度与浓度的差异确定有机酸充注系数,进而利用步骤105的现今孔隙度与现今溶蚀面孔率之间的关系确定不同成藏期次的溶蚀孔隙度模型;其中,步骤106的具体步骤如下:
步骤f1,利用面孔率和孔隙度之间的关系,建立线性关系模型如下:
φ=k·S+k1 公式(7),
式中,φ为孔隙度,%;S为溶蚀面孔率,%;k、k1为拟合系数;
利用井上现今孔隙度和现今溶蚀面孔率可以确定k、k1这两个拟合系数;
φf=k·Sf+k1 公式(8),
式中,φf为现今孔隙度,%;Sf为现今溶蚀面孔率,%;
步骤f2,地层水中的有机酸和CO2是储层发生溶蚀、形成次生孔隙的重要因素,而这些有机酸和CO2又是烃源岩在成熟过程中的产物,有机酸浓度随温度变化曲线近似于抛物线,因此可以建立有机酸浓度与温度的关系:
C=a2·T2+b·T+c1,
式中,C为有机酸充注高峰浓度,mg/L;a2、b、c1为常数,T为有机酸充注高峰时的地层温度;
通过模拟实验,建立温度与有机酸分布曲线如下:
C=-0.00385·T9+1.39·T-9.54 公式(9);
同等条件下,有机酸浓度与其溶蚀面孔率成正比,则:
Sn=m·Sf 公式(10);
m=Cn/(C1+C2+…+Cn+…+Cz) 公式(11);
式中,第n个原油充注高峰时溶蚀孔隙的面孔率为Sn,%;m为第n期成藏时的有机酸充注系数;
将Sn所得值代入公式(7)中,得到该期充注的溶蚀孔隙度,
φn=k·m·Sf+k1 公式(12);
式中,φn为第n次成藏时有机酸溶蚀孔隙度,%。
步骤107,建立成藏期孔隙度模型
根据步骤104建立的压实孔隙度模型及步骤106建立的溶蚀孔隙度模型,建立成藏期孔隙度模型。单个成藏期孔隙度模型为地层压实孔隙度+有机酸溶蚀孔隙度,即:
φ成藏期=φs+φn 公式(13);
式中,φs为剩余孔隙度,%;φn为第n次成藏时有机酸溶蚀孔隙度,%。
在应用本发明的一具体实施例中,以木垒凹陷为例,对凹陷内二叠系平地泉组的成藏期(充注期)的孔隙度进行恢复。包括以下步骤:
S01地层剥蚀量的恢复
选取木参1井,利用声波时差曲线进行剥蚀量恢复,选择该井泥岩段大于1.5m的声波时差曲线值,进行回归拟合,并根据公式(1)得出声波时差与埋深关系为:
y=-6.03·x+2208 公式(14)
将压实曲线外延至对数横坐标交于600μs/m处,则古地表埋深H0为-1410m,现今残留目的层顶面埋深H1为192m,根据公式(12)得出该地层剥蚀厚度H为1602m。
S02成藏期的确定
选取油气显示较为丰富的二叠系平地泉组砂岩段进行流体包裹体分析,确定均一化温度主峰值T1为90℃和T2为110℃。对木参1井进行单井模拟,收集包括该井井深、地层分层、古地温梯度、地层剥蚀量H等数据,进行热演化史的恢复(如图2所示)。结合均一化温度主峰值确定区内共有两次成藏,第一次成藏期为距今约170Ma左右的晚侏罗世,第二次成藏期为距今约20Ma的新近纪。
S03埋藏过程中压实孔隙度模型的建立
S03-1原始孔隙度计算
粒度分析得出P25、P75分别为9μm和11μm,利用公式(3)、(4)计算原始孔隙度φ0为46.23%;
S03-2压实孔隙度计算
A、在第一次成藏期时,二叠系平地泉组埋深h为2100m,根据公式(6)得:本次成藏期的压实孔隙度φs为30.37%;
B、在第二次成藏期时,二叠系平地泉组埋深h为2900m,根据公式(6)得:本次成藏期的压实孔隙度φs为25.88%。
S04埋藏过程中的溶蚀增孔模型的建立
统计井上现今孔隙度和现今溶蚀面孔率数据,利用公式(8)求得k、k1,建立孔隙度和溶蚀面孔率的线性关系模型:
φn=0.3197·Sn+0.5183 公式(15)。
两次充注高峰的主峰温度T1、T2分别为90℃和110℃,因此,利用公式(9)、(10)、(11)得第一次充注的有机酸充注系数为0.42,现今二叠系平地泉组所测得的溶蚀孔平均面孔率Sf为4%,第一次充注时的溶蚀孔面孔率S1为1.68%,利用公式(12)得出第一次充注时的溶蚀孔隙度为0.02%。
第二次充注的有机酸充注系数为0.58,现今二叠系平地泉组所测得的溶蚀孔平均面孔率Sf为4%,利用公式(9)、(10)、(11)得第二次充注时的溶蚀孔的面孔率S2为2.32%,利用公式(12)得出第二次充注时的溶蚀孔隙度为0.22%。
则根据公式(13)可得:
第一次成藏的孔隙度为:
φ成藏期=φs+φn=30.37%+0.02%=30.39%。
第二次成藏的孔隙度为:
φ成藏期=φs+φn=25.88%+0.22%=26.10%。
本实施例中,木垒凹陷二叠系平地泉组两次成藏赋值见表1。
表1木垒凹陷二叠系平地泉组两次成藏赋值表
以上所述是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员来说,在不脱离本发明所述原理的前提下,还可以作出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (13)
1.叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,该叠合盆地油气成藏期岩石孔隙度恢复方法包括:
步骤1,进行地层剥蚀量的恢复;
步骤2,建立热演化史模型;
步骤3,确定区内目的层油气成藏期次及各个成藏期的地层温度;
步骤4,建立埋藏过程中的压实孔隙度模型;
步骤5,建立面孔率-孔隙度模型;
步骤6,建立埋藏过程中的有机酸溶蚀孔隙度模型;
步骤7,根据步骤4建立的压实孔隙度模型及步骤6建立的有机酸溶蚀孔隙度模型,建立成藏期孔隙度模型。
2.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤1,利用声波时差曲线,挑选厚度不小于1.5m的泥岩段对应的声波时差曲线进行曲线拟合,进行剥蚀量恢复,进而恢复不同地质历史时期地层的剥蚀量。
3.根据权利要求2所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,步骤1包括:
步骤a1,选取测井曲线:统计区内探井的测井曲线,剔除非正常缩径、扩径段以及曲线上非正常跳跃段数值;
步骤a2,采集合适的测井曲线数据:选取步骤a1中纯泥岩段厚度不小于1.5m的声波时差曲线,读取其中间值;
步骤a3,剥蚀量恢复:首先,对步骤a2所得数据进行回归拟合,横坐标为声波时差数据的对数坐标,纵坐标为深度:
y=a1·x+b 公式(1),
式中,a1、b为拟合系数;x为声波时差曲线值,μs/m;y为深度,m;
其次,将压实曲线外延,与对数横坐标相交于600μs/m,则该点的水平延伸线即为古地表H0,由古地表至现今残留目的层顶面即为地层剥蚀量H:
H=H1-H0 公式(2),
式中,H0为古地表,m;H1为古地表到现今残留目的层顶面的厚度,m;H为地层剥蚀量,m。
4.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤2,通过步骤1确定的不同地层的剥蚀量,恢复古埋深,并结合古地温梯度数据,建立热演化史模型。
5.根据权利要求4所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,步骤2包括:
步骤b1,埋藏史的建立:①整理单井岩性剖面的数据,包括井上的地层分层及对应的地质年代数据;②通过步骤1确定不同地层的剥蚀量以及地层遭受剥蚀的时间;③建立直角坐标系,横坐标表示地质时间,纵坐标表示叠加地层的深度,按照地质时间由新到老进行地层及剥蚀量的叠加;④将同一地层底界相连则为该地层整个地质历史时期的古埋深,获得埋藏史图,为后期的热演化史的建立奠定基础;
步骤b2,热演化史模型的建立:首先收集研究区的古地温梯度数据,地表温度默认为20℃,其次,由于地温梯度是深度每增加100m温度的增加值,在已知地表温度为20℃及步骤b1,得到的古埋深的条件下,计算地质历史时期某一深度的地层温度,将相同温度相连即可建立热演化史模型。
6.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤3,利用流体包裹体确定成藏期次及各个期次的最大峰温值,进而在步骤2建立的生烃演化史模型的基础上确定成藏时间,并确定其古埋深。
7.根据权利要求6所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,步骤3包括:
步骤c1,单井模拟:统计单井基础资料,包括井深h、分层、古地温梯度数据、地层剥蚀量H,进行热演化史的恢复;
步骤c2,成藏期的确定:根据步骤c1的单井模拟演化结果,在生烃层中与均一温度相同的古地温线的第一个交点所对应的时间,即为油气的成藏期n;
步骤c3,流体包裹体均一温度的测定:选取新鲜样品,通过制片、岩石学观察及冷热台测试,确定第n个成藏期对应的地层流体包裹体的均一温度的主峰值Tn。
8.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤4,在明确原始孔隙度和古埋深的基础上,建立压实孔隙度模型,获得不同时期的压实孔隙度。
9.根据权利要求8所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,步骤4包括:
步骤d1,原始孔隙度计算
利用原始孔隙度φ0与分选系数S0之间的经验函数关系式:
φ0=20.91+22.9/S0 公式(3),
其中,S0=(P25/P75)1/2 公式(4),
式中,S0为分选系数;φ0为原始孔隙度,%;P25,P75分别代表粒度累计曲线上颗粒含量在25%和75%处所对应的颗粒直径,μm;
步骤d2,利用机械压实孔隙度与埋深的指数关系,来恢复经历压实作用后成藏期末的孔隙度,即:
φs=φ0·e-0.0002·h 公式(5),
其中,φs为剩余孔隙度,%;φ0为恢复砂岩的原始孔隙度,%;h为埋深,m;
步骤d3,根据Tn在埋藏史图上标出目的层在成藏期的埋深,将埋深代入公式(5),得到成藏期的压实孔隙度,
10.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤5,根据薄片、铸体薄片观察开展镜下溶蚀孔识别及划分,结合物性分析结果,统计溶蚀孔的面孔率与孔隙度的关系,建立面孔率-孔隙度模型。
11.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤6,利用温度与浓度的差异确定有机酸充注系数,进而利用步骤5的现今孔隙度与现今溶蚀面孔率之间的关系确定不同成藏期次的溶蚀孔隙度模型。
12.根据权利要求11所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,步骤6包括:
步骤f1,利用面孔率和孔隙度之间的关系,建立线性关系模型如下:
φ=k·S+k1 公式(7),
式中,φ为孔隙度,%;S为溶蚀面孔率,%;k、k1为拟合系数;
利用井上现今孔隙度和现今溶蚀面孔率确定k、k1这两个拟合系数;
φf=k·Sf+k1 公式(8),
式中,φf为现今孔隙度,%;Sf为现今溶蚀面孔率,%;
步骤f2,地层水中的有机酸和CO2是储层发生溶蚀、形成次生孔隙的重要因素,而这些有机酸和CO2又是烃源岩在成熟过程中的产物,有机酸浓度随温度变化曲线近似于抛物线,因此建立有机酸浓度与温度的关系:
C=a2·T2+b·T+c1,
式中,C为有机酸充注高峰浓度,mg/L;a2、b、c1为常数,T为有机酸充注高峰时的地层温度;
通过模拟实验,建立温度与有机酸分布曲线如下:
C=-0.00365·T2+1.39·T-9.54 公式(9);
同等条件下,有机酸浓度与其溶蚀面孔率成正比,则:
Sn=m·Sf 公式(10);
m=Cn/(C1+C2+…+Cn+…+Cz) 公式(11);
式中,第n个原油充注高峰时溶蚀孔隙的面孔率为Sn,%;m为第n期成藏时的有机酸充注系数;
将Sn所得值代入公式(7)中,得到该期充注的溶蚀孔隙度,
φn=k·m·Sf+k1 公式(12);
式中,φn为第n次成藏时有机酸溶蚀孔隙度,%。
13.根据权利要求1所述的叠合盆地油气成藏期岩石孔隙度恢复方法,其特征在于,在步骤7,单个成藏期孔隙度模型为地层压实孔隙度+有机酸溶蚀孔隙度,即:
φ成藏期=φs+φn 公式(13);
式中,φs为剩余孔隙度,%;φn为第n次成藏时有机酸溶蚀孔隙度,%。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811463682.4A CN109324345B (zh) | 2018-12-03 | 2018-12-03 | 叠合盆地油气成藏期岩石孔隙度恢复方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811463682.4A CN109324345B (zh) | 2018-12-03 | 2018-12-03 | 叠合盆地油气成藏期岩石孔隙度恢复方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109324345A true CN109324345A (zh) | 2019-02-12 |
CN109324345B CN109324345B (zh) | 2020-07-14 |
Family
ID=65256173
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811463682.4A Active CN109324345B (zh) | 2018-12-03 | 2018-12-03 | 叠合盆地油气成藏期岩石孔隙度恢复方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109324345B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111009179A (zh) * | 2019-10-23 | 2020-04-14 | 上海同继地质科技有限公司 | 剥蚀厚度的确定方法及装置 |
CN112528463A (zh) * | 2020-11-06 | 2021-03-19 | 中国石油天然气股份有限公司 | 一种古老碳酸盐岩储层成岩-孔隙演化恢复方法 |
CN112528226A (zh) * | 2020-11-06 | 2021-03-19 | 中国石油天然气股份有限公司 | 一种深层碳酸盐岩储层孔隙演化恢复方法 |
CN112965114A (zh) * | 2021-02-07 | 2021-06-15 | 中国石油大学(华东) | 一种海上深层天然气储层甜点评价方法 |
CN113138414A (zh) * | 2020-07-24 | 2021-07-20 | 中国石油天然气股份有限公司 | 深层砂岩储层孔隙度的预测方法及装置 |
CN113309511A (zh) * | 2021-07-13 | 2021-08-27 | 西安石油大学 | 一种致密油成藏动力的获取方法 |
CN114428088A (zh) * | 2020-09-08 | 2022-05-03 | 中国石油化工股份有限公司 | 烃源流体溶蚀增孔能力评价图版的生成方法 |
CN114660269A (zh) * | 2022-03-23 | 2022-06-24 | 中国海洋石油集团有限公司 | 一种古潜山天然气动态成藏过程恢复方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105334150A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 前陆盆地砂岩储层成岩过程与孔隙演化过程的确定方法 |
US20160103245A1 (en) * | 2014-10-09 | 2016-04-14 | Chevron U.S.A. Inc. | Conditioning of object or event based reservior models using local multiple-point statistics simulations |
CN106368694A (zh) * | 2016-11-16 | 2017-02-01 | 中国海洋石油总公司 | 一种针对复杂领域储层孔隙演化恢复与物性预测方法 |
CN108363115A (zh) * | 2018-02-05 | 2018-08-03 | 西安石油大学 | 一种致密石英砂岩孔隙度恢复模拟方法 |
-
2018
- 2018-12-03 CN CN201811463682.4A patent/CN109324345B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160103245A1 (en) * | 2014-10-09 | 2016-04-14 | Chevron U.S.A. Inc. | Conditioning of object or event based reservior models using local multiple-point statistics simulations |
CN105334150A (zh) * | 2015-11-24 | 2016-02-17 | 中国石油大学(北京) | 前陆盆地砂岩储层成岩过程与孔隙演化过程的确定方法 |
CN106368694A (zh) * | 2016-11-16 | 2017-02-01 | 中国海洋石油总公司 | 一种针对复杂领域储层孔隙演化恢复与物性预测方法 |
CN108363115A (zh) * | 2018-02-05 | 2018-08-03 | 西安石油大学 | 一种致密石英砂岩孔隙度恢复模拟方法 |
Non-Patent Citations (1)
Title |
---|
王艳忠 等: ""碎屑岩储层地质历史时期孔隙度演化恢复方法———以济阳坳陷东营凹陷沙河街组四段上亚段为例"", 《石油学报》 * |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111009179A (zh) * | 2019-10-23 | 2020-04-14 | 上海同继地质科技有限公司 | 剥蚀厚度的确定方法及装置 |
CN111009179B (zh) * | 2019-10-23 | 2022-01-25 | 上海同继地质科技有限公司 | 剥蚀厚度的确定方法及装置 |
CN113138414A (zh) * | 2020-07-24 | 2021-07-20 | 中国石油天然气股份有限公司 | 深层砂岩储层孔隙度的预测方法及装置 |
CN114428088A (zh) * | 2020-09-08 | 2022-05-03 | 中国石油化工股份有限公司 | 烃源流体溶蚀增孔能力评价图版的生成方法 |
CN114428088B (zh) * | 2020-09-08 | 2024-02-23 | 中国石油化工股份有限公司 | 烃源流体溶蚀增孔能力评价图版的生成方法 |
CN112528463A (zh) * | 2020-11-06 | 2021-03-19 | 中国石油天然气股份有限公司 | 一种古老碳酸盐岩储层成岩-孔隙演化恢复方法 |
CN112528226A (zh) * | 2020-11-06 | 2021-03-19 | 中国石油天然气股份有限公司 | 一种深层碳酸盐岩储层孔隙演化恢复方法 |
CN112965114A (zh) * | 2021-02-07 | 2021-06-15 | 中国石油大学(华东) | 一种海上深层天然气储层甜点评价方法 |
CN113309511A (zh) * | 2021-07-13 | 2021-08-27 | 西安石油大学 | 一种致密油成藏动力的获取方法 |
CN113309511B (zh) * | 2021-07-13 | 2023-08-25 | 西安石油大学 | 一种致密油成藏动力的获取方法 |
CN114660269A (zh) * | 2022-03-23 | 2022-06-24 | 中国海洋石油集团有限公司 | 一种古潜山天然气动态成藏过程恢复方法 |
CN114660269B (zh) * | 2022-03-23 | 2023-09-19 | 中国海洋石油集团有限公司 | 一种古潜山天然气动态成藏过程恢复方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109324345B (zh) | 2020-07-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
CN105866835B (zh) | 一种基于地应力分布的断层三维封闭性定量评价方法 | |
CN106295095B (zh) | 基于常规测井资料预测低渗透砂岩储层产能的方法 | |
Evans et al. | Microseismicity and permeability enhancement of hydrogeologic structures during massive fluid injections into granite at 3 km depth at the Soultz HDR site | |
CN103256046B (zh) | 非常规油气藏水平井全缝长压裂参数模拟的方法及装置 | |
CN104564041B (zh) | 基于开发渗透率下限的低渗透碎屑岩储层有效性评价方法 | |
CN102562048B (zh) | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 | |
CN106154351A (zh) | 一种低孔渗储层渗透率的估算方法 | |
CN106569288B (zh) | 基于储层裂缝有效性聚类分析的裂缝性储层质量评价方法 | |
CN109946147B (zh) | 一种煤系复合储层压裂物理模拟试件制备方法 | |
CN105629325B (zh) | 前陆盆地冲积扇精细刻画与预测方法 | |
CN105651962B (zh) | 成岩相识别方法 | |
Bekele et al. | Modeling secondary oil migration with core-scale data: Viking Formation, Alberta Basin | |
CN111006987B (zh) | 一种储层裂缝孔渗参数预测方法 | |
CN109509111A (zh) | 探井地层压力的预测方法及系统 | |
CN105134185A (zh) | 储层流体性质识别方法 | |
CN107301483A (zh) | 快速综合评价未动用储量经济有效动用可行性的方法 | |
CN105467466B (zh) | 一种基于多尺度信息的致密储层成岩相预测方法 | |
Wimmers et al. | Integration of sedimentology, petrophysics and rock typing as key to understanding a tight gas reservoir | |
CN112177605A (zh) | 一种海域天然气水合物有利聚集区主控因素的确定方法 | |
Sun et al. | Quantitative evaluation the physical properties evolution of sandstone reservoirs constrained by burial and thermal evolution reconstruction: A case study from the Lower Cretaceous Baxigai Formation of the western Yingmaili Area in the Tabei Uplift, Tarim Basin, NW China | |
CN108304959A (zh) | 提高地层流体压力预测精度的方法 | |
Yang et al. | The coupling of dynamics and permeability in the hydrocarbon accumulation period controls the oil-bearing potential of low permeability reservoirs: a case study of the low permeability turbidite reservoirs in the middle part of the third member of Shahejie Formation in Dongying Sag | |
CN108362621B (zh) | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 | |
CN112185469B (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 |