专利名称：地质历史时期砂岩储层渗透率演化恢复方法
技术领域：
本发明涉及石油天然气勘探与开发领域，特别是一种地质历史时期砂岩储层渗透率演化恢复方法。
背景技术：
随着中浅层油气勘探程度的不断提高和对油气储量增长的需求，油气勘探目标已逐渐转向深层，我国深层油气资源潜力巨大，并且油气勘探已取得了可喜的成果。但是，深层砂岩储层在埋藏过程中经历了复杂的成岩作用改造，地质历史时期储层孔隙度和孔隙结构经历了复杂的演化过程，导致渗透率的演化极其复杂，直接影响着中深层油气富集规律的复杂性、储层表现出“忽油忽水”、“忽储忽干”等复杂特点，使得中深层油气探勘难度大、风险高、成功率低。油气成藏时期储层渗透性的好坏，直接决定着油气储层的有效性，开展砂岩储层地质历史时期渗透率演化恢复，对于油气成藏时期储层有效性的评价，预测有利储层，指导油气勘探具有重要意义。目前国内外学者仅初步开展了地质历史时期砂岩储层孔隙度演化恢复，尚未开展地质历史时期孔隙结构与渗透率演化恢复的研究。随着油气勘探程度的不断提高，对有利储层预测精度要求越来越高，进行砂岩储层地质历史时期渗透率演化恢复显得十分必要。
发明内容
本发明的目的在于以地质历史时期砂岩储层渗透率演化恢复方法为目标，在孔隙度演化恢复的基础上，采用图像分析的方法，恢复地质历史时期砂岩储层孔隙结构，然后进行孔隙结构约束下的砂岩储层渗透率演化恢复。本发明的技术方案为一种地质历史时期砂岩储层渗透率演化恢复方法，具体步骤如下(I)地质历史时期砂岩储层孔隙度演化恢复第一步建立砂岩储层成岩作用演化序列及确定成岩作用发生时储层相应的古埋深对现今铸体薄片进行观察、扫描电镜分析和流体包裹体分析，结合成岩环境演化分析，建立砂岩储层成岩作用演化序列，根据建立的成岩演化序列，确定各成岩作用开始及结束的时间，将这些时间投影到所要恢复样品点的单井埋藏史上，获得成岩作用发生时储层相应的古埋深；第二步建立正常压实图版第一，分析沉积物粒径、分选及刚性颗�：慷匝故底饔玫挠跋斐潭龋范ㄓ跋煅故底饔玫闹骺匾蛩兀坏诙友芯壳钟惺挡馕镄缘难分校粞》系夭阊沽ξＱ�、胶结物含量小于5%、次生孔隙含量小于1%的样品，作为正常压实条件下的样品；第三，将所选的正常压实条件下的样品，分不同主控因素进行深度与孔隙度拟合，建立孔隙度演化曲线，即可作为研究区储层正常压实图版；第四，结合本地区地温梯度演化史，确定目的层演化到现今经历的平均地温梯度，建立的正常压实图版可认为是该平均地温梯度控制下的正常压实图版；第三步建立孔隙度与面孔率之间的函数关系根据物理学原理，人眼可分辨明视距处的最小直线距离为约0. Imm,因此，在显微镜下0. 5 y m的直线距离放大200倍时即可被人眼分辨明视，即200倍镜下人眼能分辨明视的孔隙半径为0. 25 u m,基于上述情况，将200倍镜下孔隙半径小于0. 25 y m的孔隙视为微孔隙，在统计面孔率与孔隙度的关系时不予考虑；借助于压汞资料确定出实测孔隙度中半径大于0. 25 y m的孔隙含量，便可求取200倍镜下人眼能分辨明视孔隙度，称为显孔隙度，然后借助现今铸体薄片，利用偏光显微镜及计算机图像分析技术统计其对应的面孔率，将面孔率与对应的显孔隙度进行拟合，就可以建立面孔率与显孔隙度之间的函数关系；第四步成岩序列约束下的砂岩储层孔隙度反演回剥以现今铸体薄片孔隙特征为基�。猿裳已莼蛄形际幼詈笠黄诔裳易饔每蓟匕指锤髌诔裳易饔每记暗目紫短卣鳎捎眉扑慊枷穹治黾际酰考扑愀髌诔裳易饔茫故底饔贸猓源⒉忝婵茁视跋毂浠浚⒏菝婵茁视胂钥紫抖戎涞墓叵担纸裰灞∑械拿婵茁首杂Φ南钥紫抖龋指锤鞒裳易饔每�、结束时的孔隙度，进而获得各古埋深下的孔隙度；第五步机械压实作用和热压实作用校正根据成岩演化序列，在主要胶结作用发生之前，样品为正常压实阶段，其在各埋深下的压实损孔量可根据与样品具相同特征的正常压实图版读取；而在主要胶结时期开始后，胶结作用抑制了压实作用的正常进行，各时期的压实减孔量可采用主要胶结作用开始后的压实总减孔量按正常压实图版上各时期压实减孔量的比例对其进行分配而确定各成岩作用时期的压实减孔校正量，对步骤四所得的反演回剥孔隙度进行机械压实作用校正；然后，根据研究区埋藏史及地温梯度演化史，确定各成岩作用时期经历的平均地温梯度及古埋深，利用地温场与砂岩孔隙度关系图，对步骤四所得的反演回剥孔隙度进行热压实作用校正；第六步在上述步骤的基础上，确定各主要成岩阶段真实孔隙度，建立地质历史时期砂岩储层孔隙度演化曲线。(2)地质历史时期砂岩储层孔隙结构演化恢复第一步利用反演回剥法，在成岩作用演化序列的约束下，将胶结作用及溶蚀作用成岩事件进行反演回剥；第二步对压实损失的孔隙面积进行等比例放大(如图I所示)，具体恢复方法如下①参量设置待恢复状态设视域长为L,宽为D岩石体积V2z孔隙体积V2P视域面积S2z (未知)孔隙面积S2P (未知)、
面孔率A2 = S2P/S2z*100% (可以由①于A的关系求取，A2 = F(O2))孔隙度02= V2P/V2z*100% (已知)已知状态设视域长为L+AL,宽为D+AD岩石体积Vlz孔隙体积:V1P 视域面积Slz (已知)孔隙面积Sip (已知)面孔率=A1= S1P/Slz*100% (已知)孔隙度=Ol= V1P/Vlz*100% (已知)设由待恢复状态到已知状态压实减小的面积为Sy，取正直设由待恢复状态到已知状态胶结减小的面积为Sj，取正直设由待恢复状态到已知状态溶解增加的面积为S 取正直②推导过程S2p-Sy-Sj+Sr = Slp 公式 IA2 = S2P/S2z公式 2S2z-Sy = Slz公式 3把公式3带入公式2可得，A2 = S2P/S2z = S2P/ (Slz+Sy)S2p = A2* (Slz+Sy)公式 4把公式4带入公式I :A2* (Slz+Sy) -Sy-Sj+Sr = SlpA2*Slz+A2*Sy-Sy-Sj+Sr = SlpSy (I-A2) = A2*Slz+Sr-Sj-SlpSy = (A2*Slz-Slp+Sr-Sj) / (I-A2) 公式 5Sy= (AdSlz-AdS1JSr-Sj)/(1-A2)Sy = (Slz (A2-A1)+Sr-Sj) / (1~A2) 公式 5’Sy = (A2*Slz-Slp_Sj+Sr) / (I-A2)设由待恢复状态到已知状态压实作用使总视域长宽成等比例缩小，即L/D = (L+ AL)/(D+ AD) 公式 6Slz+Sy = (L+ AL)*(D+ AD) 公式 7将公式6带入公式7可得，Slz+Sy = (L/D) * (D+ A D)2D (Slz+Sy) /L = (D+ A D)2AD = [D(Slz+Sy)/I] 1/2-D 公式 8A L = L* (D+ A D) /D-LA L = L* A D/D 公式 9求得AD、AL后，选择中心颗粒不动，从两端开始按比例移动颗粒，假设某一颗�：嵯虻阶菹蛑行南叩木嗬胛皆蛐柘虮咭贫木嗬階d = 2 AL d/ ( A L+L)；
第三步建立地质历史时期砂岩储层孔隙结构演化图版，并进行各时期孔隙结构分类；(3)地质历史时期砂岩储层渗透率演化恢复第一步利用压汞资料，根据排驱压力(Pd)、毛管压力中值(P50)等参数，建立每一类孔隙结构的K与K/0关系，其中K为渗透率，O为孔隙度；第二步根据每一类孔隙结构K与K/O的关系，结合孔隙结构恢复与分类结果，求取地质历史时期砂岩储层渗透率；第四步；建立地质历史时期砂岩储层渗透率演化曲线。本发明的有益效果为本发明解决了地质历史时期砂岩储层渗透率恢复的问题，为研究油气成藏时期砂岩储层渗透性的好坏，评价成藏时期储层的有效性，预测有利储层， 指导油气勘探具有重要意义。
图I为孔隙结构恢复示意图。图2为恢复样品点的单井埋藏史。图3为胜利油田某地区古近系砂岩不同分选下储层正常压实图版。图4为胜利油田某地区面孔率与显孔隙度之间的函数关系。图5为利用反演回剥法恢复的各成岩作用时期的物性状态示意图。图6为胜利油田某地区机械压实校正示意图。图7为地温场与砂岩的孔隙度关系图。图8为某井3431. 25m地质历史时期储层孔隙度演化曲线。图9为某井3431. 25m深度孔隙结构演化恢复图版。图10为某地区孔隙结构分类及物性关系拟合图。图11为某井3431. 25m地质历史时期储层渗透率演化曲线。
具体实施例方式以胜利油田某地区地质历史时期砂岩储层渗透率演化恢复方法来说明该发明的具体技术方案(I)砂岩储层地质历史时期孔隙度演化恢复第一步建立砂岩储层成岩作用演化序列及成岩作用发生时储层相应的古埋深利用薄片观察、扫描电镜分析及流体包裹体分析等技术，结合成岩环境演化分析，认为胜利油田某地区沙四段成岩演化序列为距今44. 5Ma-距今42. 5Ma，成岩作用以压实作用为主；距今42. 5Ma-距今32Ma，压实作用/长石溶解/石英加大；距今32Ma_距今24. 6Ma，压实作用/石英溶解/碳酸盐胶结；距今24. 6Ma_距今2Ma，压实作用/碳酸盐胶结物溶解；距今2Ma-现今，压实作用/晚期黄铁矿。根据上述建立的成岩演化序列，确定各成岩作用开始及结束的时间，将这些时间投影到所要恢复样品点的单井埋藏史(图2)上，可获得各期成岩作用开始和结束的古埋深。以某井3431. 25m为例，距今约45Ma前，开始沉积，古埋深为Om ;距今42. 5Ma，埋深约1350m ;距今32Ma，埋深约2520m ;距今24. 6Ma,埋深约2800m ;距今2Ma，埋深约3360m ;现今埋深 3431. 25m。第二步建立正常压实图版 选取胜利油田某地区符合正常压实的样品点进行正常压实图版的建立，建立了分选系数为I. 5-1. 75、I. 75-2,2-2. 5及大于2. 5的各种分选下的正常压实曲线(图3)。第三步建立孔隙度与面孔率之间的函数关系利用大量样品实测孔隙度，结合压汞资料分析实测孔隙度中人眼能分辨明视的显孔隙度含量，求取显孔隙度，然后将显孔隙度与对应铸体薄片实测面孔率进行函数拟合，建立面孔率与显孔隙度之间的函数关系(图4)。第四步成岩序列约束下的砂岩储层孔隙度反演回剥以某井3431. 25m为例(分选系数I. 59)，选取典型视域，进行3行X 3列共9个视域的照片进行拼接，现今实测孔隙度为11. 7% ;碳酸盐胶结物溶解产生I. 26%的面孔率(转化为孔隙度为3. 06% );碳酸盐胶结损失6. 5%的面孔率(转化为孔隙度为12. 27%),石英溶解产生0. 15%的面孔率(转化为孔隙度为0. 5 % )，石英加大损失面孔率0. 06%(0. 23% );长石和岩屑溶解产生2. 34%的面孔率(转化为孔隙度为5. 17% )，各成岩作用时期面孔率特征见图5。第五步机械压实作用和热压实作用校正I.机械压实作用校正以胜利油田某地区为例进行说明，假设反演回剥长石溶解/石英加大前(距今42. 5Ma)孔隙度为Ct11，此时正常压实情况下的孔隙度为，此时真实孔隙度为(K，这一阶段压实损失孔隙度为第一期酸性溶解结束时(距今32Ma)反演回剥孔隙度为小21，对应正常压实孔隙度为A2，这一阶段压实损失孔隙度为此时真实孔隙度为cK ;碳酸盐胶结后及第二期酸性溶解前(距今24. 6Ma)通过反演回剥恢复的孔隙度为031，对应相 同深度正常压实情况下孔隙度为小3，这一阶段压实损失孔隙度为此时真实孔隙度为小。，第二期酸性溶解结束(距今2Ma)通过反演回剥恢复的孔隙度为A41，对应相同深度正常压实情况下孔隙度为A4，这一阶段压实损失孔隙度为此时真实孔隙度为ctd，距今2Ma至今压实损失孔隙度为小5压损，对应真实孔隙度t，即叭=小a今(图6)。由上述可知，由于反演回剥法得出的孔隙度是没有考虑压实损失的情况下得出的，此时真实的孔隙度要大于这个值，即叭> A11，由于这一时期没有经历胶结，基本为正常压实阶段。因此，可认为0a= Ct1，而经历了第一期酸性溶解后(长石溶解、石英加大)的真实孔隙度(^应该为该时期正常压实孔隙度A2加上该时期长石溶解孔隙度，减去石英加大损失孔隙度，这一阶段基本也可以看作是正常压实。经历了碱性环境后，发生了碳酸盐的胶结，少量石英的溶解，该时期由于碳酸盐胶结物一般含量较高，大多会抑制正常压实的进行，不能利用正常压实曲线进行求取该时期的真实孔隙度。同时后期酸性溶解这一阶段也不能用正常压实曲线求取该时期的真实孔隙度，从第一期酸性溶解结束、碱性环境开始到现今(距今32Ma_现今)的压实损失的总孔隙度(小2-5_)可以计算出来，即 4) 2_5压损=b_ ^ _ ^+ ^石難》+ ^細盐胶结溶解而第一期酸性溶解结束-石英溶解/碳酸盐胶结结束(距今32Ma_距今24. 6Ma)这一阶段压实损失的孔隙度与碳酸盐胶结结束至距今2Ma (距今24. 6Ma-2Ma)压实损失的孔隙度距今2Ma至现今压实损失的孔隙度，无法精确求出，根据正常压实曲线上这两个阶段孔隙压实损失的比例对进行分配，求的
^ 2-3压损、办3—4压损及
由4-5压损进而求取七、:= 4>b~4>細滅结+ ^石英溶解"^ 2-3压损 (J) d = (J) b- (J)碳酸盐胶结+ 石英溶解_小2-3压损_小3-4压损 (J) e = (J) b_ (J)碳酸盐胶结+小石英溶解_小2-3压损_ ^ 3-4压损_ ^ 4-5压损=^现今某井3431. 25m沉积初期孔隙度为49. I %，长石溶解石英加大之前，主要为正常压实，此时孔隙度为对应正常压实曲线上的孔隙度32. 6% (埋深1350m);长石溶解/石英加大之后碳酸盐胶结/石英加大之时(埋深为2520m)对应正常压实曲线上孔隙度(19. 3%) +长石溶解孔隙度(5. 17%)-石英加大孔隙度(0. 23%)，即为24. 24%;碳酸盐胶结/石英溶解结束后，由于碳酸盐胶结抑制了压实作用的进行，各成岩阶段压实减孔量按正常压实图版上各时期压实减孔量的比例对其进行分配，而确定各成岩作用时期的压实减孔校正量，碳酸盐胶结之后各时期压实总减孔量为3. 56%，按2. 8 5 0.5的比例进行分配；碳酸盐胶结物溶解之前(埋深2800m)的孔隙度为长石溶解/石英加大之后碳酸盐胶结/石英加大之前对应的孔隙度(24. 24%)-碳酸盐胶结孔隙度(12. 27%)-压实损失孔隙度(1.2%)，即为10. 77%;碳酸盐胶结物溶解之后(埋深3360m)的孔隙度为碳酸盐胶结物溶解之前孔隙度(10. 77%) +碳酸盐胶结物溶解孔隙度(3. 06%)-压实损失孔隙度(2. 14%)，即为
11.69%;现今(埋深3431. 25m)的孔隙度为碳酸盐胶结物溶解之后孔隙度(11.69%)-压实损失孔隙度(0. 22%)，即为11. 47%。由此恢复出各成岩作用开始(结束)时的孔隙度。2.热压实作用校正某井3431. 25m在沉积后到长石溶解/石英加大之前(距今42. 5Ma，埋深1350m)经历了沙四上亚段-沙三下亚段沉积时期，经历的平均地温梯度为4. 3°C /100m，在埋深为1350m时，与地温梯度为3. 87V /IOOm的情况下，孔隙度将多压实3. 5%左右(图6)，SP此时真实的孔隙度为，地温梯度为3. 870C /IOOm的情况下的32. 6 % -3. 5 % = 29. I % ;同理，至碳酸盐胶结/石英溶解之前(距今32Ma，埋深2520m左右)，经历的平均地温梯度为
4.18°C /100m(沙四上亚段-沙一段沉积时期的平均地温梯度)，与3. 87V /IOOm下相比，2520m埋深下需要较正的孔隙度值为3% (图6)，即此时真实孔隙度应为24.24% -3% =21. 24%，其余各段平均地温梯度与建立正常压实图版的平均地温梯度3. 87V /IOOm相差不大，因此，可以忽略热压实作用的影响，不进行热压实作用校正。第六步最终恢复各主要成岩作用时期真实孔隙度(表I)，并建立了地质历史时期砂岩储层孔隙度演化曲线(图7)。表I某井3239. Im处各成岩作用时期储层孔隙度值恢复结果
权利要求
1.一种地质历史时期砂岩储层渗透率演化恢复方法，具体步骤如下 (1)地质历史时期砂岩储层孔隙度演化恢复； (2)地质历史时期砂岩储层孔隙结构演化恢复 第一步利用反演回剥法，在成岩作用演化序列的约束下，将胶结作用及溶蚀作用成岩事件进行反演回剥； 第二步对压实损失的孔隙面积进行等比例放大，具体恢复方法如下 ①参量设置 待恢复状态 设视域长为L,宽为D 岩石体积V2z 孔隙体积V2P 视域面积S2z (未知) 孔隙面积S2P(未知) 面孔率A2 = S2p/S2z*100% (可以由①于A的关系求�。珹2 = F(O2)) 孔隙度 A2 = V2P/V2z*100% (已知) 已知状态 设视域长为L+ A L，宽为D+ A D 岩石体积Vlz 孔隙体积V1P 视域面积Slz(已知) 孔隙面积S1P(已知) 面孔率=A1 = S1P/Slz*100% (已知) 孔隙度 Al = V1P/Vlz*100% (已知) 设由待恢复状态到已知状态压实减小的面积为Sy，取正直 设由待恢复状态到已知状态胶结减小的面积为=Sj,取正直 设由待恢复状态到已知状态溶解增加的面积为=Sp取正直 ②推导过程 S2p-Sy-Sj+Sr = Slp 公式 I A2 = S2P/S2z公式 2 S2z-Sy = Slz公式 3 把公式3带入公式2可得， A2 = S2P/S2z = S2P/ (Slz+Sy) S2p = A2* (Slz+Sy)公式 4 把公式4带入公式I :A2*(Slz+Sy)-Sy-Sj+Sr = SlpA2*Slz+A2*Sy-Sy-Sj+Sr = SlpSy(I-A2) = A2*Slz+Sr-S厂SlpSy = (A2*Slz-Slp+Sr-Sj) / (I-A2) 公式 5Sy = (A2^S1-A1^S1+S-Sj) / (I-A2)Sy = (Slz (A2-A1) +Sr-Sj) / (1-A2) 公式 5’Sy = (A2*Slz-Slp-Sj+Sr)/(I-A2) 设由待恢复状态到已知状态压实作用使总视域长宽成等比例缩�。� L/D = (L+ AL)/(D+ AD) 公式 6 Slz+Sy = (L+ AL)*(D+ AD) 公式 I 将公式6带入公式7可得，Slz+Sy = (L/D)*(D+AD)2D(Slz+Sy)/L = (D+AD)2AD= [D(Slz+Sy)/LJ172-D 公式 8A L = L* (D+ A D) /D-LA L = L* A D/D 公式 9 求得AD、AL后，选择中心颗粒不动，从两端开始按比例移动颗粒，假设某一颗�：嵯虻阶菹蛑行南叩木嗬胛猟，则需向边移动的距离Ad = 2 AL d/(AL+L)； 第三步建立地质历史时期砂岩储层孔隙结构演化图版，并进行各时期孔隙结构分类； (3)地质历史时期砂岩储层渗透率演化恢复 第一步利用压汞资料，根据排驱压力(Pd)、毛管压力中值(P50)等参数，建立每一类孔隙结构的K与K/0关系，其中K为渗透率，O为孔隙度； 第二步根据每一类孔隙结构K与K/O的关系，结合孔隙结构恢复与分类结果，求取地质历史时期砂岩储层渗透率； 第四步；建立地质历史时期砂岩储层渗透率演化曲线。
2.如权利要求I所述的地质历史时期砂岩储层渗透率演化恢复方法，其特征在于地质历史时期砂岩储层孔隙度演化恢复方法，具体步骤如下 第一步建立砂岩储层成岩作用演化序列及确定成岩作用发生时储层相应的古埋深 第二步建立正常压实图版 第一，分析沉积物粒径、分选及刚性颗�：慷匝故底饔玫挠跋斐潭龋范ㄓ跋煅故底饔玫闹骺匾蛩兀坏诙友芯壳钟惺挡馕镄缘难分校粞》系夭阊沽ξＱ埂⒔航嵛锖啃∮�5%、次生孔隙含量小于1%的样品，作为正常压实条件下的样品；第三，将所选的正常压实条件下的样品，分不同主控因素进行深度与孔隙度拟合，建立孔隙度演化曲线，作为研究区储层正常压实图版；第四，结合本地区地温梯度演化史，确定目的层演化到现今经历的平均地温梯度，建立的正常压实图版可认为是该平均地温梯度控制下的正常压实图版； 第三步建立孔隙度与面孔率之间的函数关系 借助于压汞资料确定出实测孔隙度中半径大于0. 25 的孔隙含量，求取200倍镜下人眼能分辨明视孔隙度，称为显孔隙度，然后借助现今铸体薄片，利用偏光显微镜及计算机图像分析技术统计其对应的面孔率，将面孔率与对应的显孔隙度进行拟合，就可以建立面孔率与显孔隙度之间的函数关系； 第四步成岩序列约束下的砂岩储层孔隙度反演回剥 以现今铸体薄片孔隙特征为基�。猿裳已莼蛄形际幼詈笠黄诔裳易饔每蓟匕指锤髌诔裳易饔每记暗目紫短卣鳎捎眉扑慊枷穹治黾际酰考扑愀髌诔裳易饔茫故底饔贸猓源⒉忝婵茁视跋毂浠浚⒏菝婵茁视胂钥紫抖戎涞墓叵担纸裰灞∑械拿婵茁首杂Φ南钥紫抖龋指锤鞒裳易饔每�、结束时的孔隙度，进而获得各古埋深下的孔隙度； 第五步机械压实作用和热压实作用校正 根据成岩演化序列，在主要胶结作用发生之前，样品为正常压实阶段，其在各埋深下的压实损孔量可根据与样品具相同特征的正常压实图版读�。欢谥饕航崾逼诳己螅航嶙饔靡种屏搜故底饔玫恼＝校魇逼诘难故导蹩琢靠刹捎弥饕航嶙饔每己蟮难故底芗蹩琢堪凑Ｑ故低及嫔细魇逼谘故导蹩琢康谋壤云浣蟹峙涠范ǜ鞒裳易饔檬逼诘难故导蹩仔Ｕ浚圆街杷乃玫姆囱莼匕紫抖冉谢笛故底饔眯Ｕ蝗缓螅菅芯壳癫厥芳暗匚绿荻妊莼罚范ǜ鞒裳易饔檬逼诰钠骄匚绿荻燃肮怕裆睿玫匚鲁∮肷把铱紫抖裙叵低迹圆街杷乃玫姆囱莼匕紫抖冉腥妊故底饔眯Ｕ坏诹皆谏鲜霾街璧幕∩希范ǜ髦饕裳医锥握媸悼紫抖龋⒌刂世肥逼谏把掖⒉憧紫抖妊莼�。
3.如权利要求2所述的地质历史时期砂岩储层渗透率演化恢复方法，其特征在于 第一步中对现今铸体薄片进行观察、扫描电镜分析和流体包裹体分析，结合成岩环境演化分析，建立砂岩储层成岩作用演化序列，根据建立的成岩演化序列，确定各成岩作用开始及结束的时间，将这些时间投影到所要恢复样品点的单井埋藏史上，获得成岩作用发生时储层相应的古埋深。
全文摘要
本发明涉及一种地质历史时期砂岩储层渗透率演化恢复方法，其步骤为(1)地质历史时期砂岩储层孔隙度演化恢复；(2)地质历史时期砂岩储层孔隙结构演化恢复；(3)地质历史时期砂岩储层渗透率演化恢复。本发明解决了地质历史时期砂岩储层渗透率恢复的问题，为研究油气成藏时期砂岩储层渗透性的好坏，评价成藏时期储层的有效性，预测有利储层，指导油气勘探具有重要意义。
文档编号G01N15/08GK102778421SQ20121023534
公开日2012年11月14日 申请日期2012年7月10日 优先权日2012年7月10日
发明者刘惠民, 张善文, 操应长, 王永诗, 王艳忠, 葸克来, 陈林, 马奔奔 申请人:中国石油大学(华东)
专业数控铣床产品供应商
20年生产经验，实力工厂，精工品质，质量有保障
