技术领域
[0001] 本
发明涉及油气勘探开发技术领域，具体地指一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法。
背景技术
[0002] 测井是利用岩层的电化学特征、
导电性特征、
放射性等地球物理特性，测量
地层地球物理参数的方法。测井广泛应用于石油地质和油田勘探开发的全过程。利用测井资料，可以划分地层剖面，确定岩层厚度和埋藏深度，确定储层并识别油气
水层，进行区域地层对比，而且还可以探测和研究地层主要矿物成分、孔隙度、渗透率、油气
饱和度、裂缝、
断层、构造特征和沉积环境与砂体的分布等，对于评价地层的储集能
力、检测油气藏的开采情况，细致地分析研究油层地质特征等具有重要意义。获取地层孔隙度的测井方法通常有
中子测井、孔隙度测井、
声波时差测井。通过三孔隙度测井资料能连续、定量、准确获取高孔隙度或中-低孔隙度储层的孔隙值，对于深层、超深层低孔、低渗储层，三孔隙度测井计算孔隙度的准确性变差，影响了处于高成岩阶段低孔、低渗储层油气勘探开发，必须寻找出一种不同于常规方法的评价方法，来更好地、更合理地预测有效致密储层分布，最终为油气勘探开发提供有力依据。
[0003] 成岩相测井是利用
岩石的电性特征，在储层岩性识别的
基础上对成岩作用类型和强度连续、定性、定量判定以及成岩矿物含量的定量计算。成岩相测井的目的是建立测井信息与成岩相之间关系，并对
压实程度、胶结程度、溶蚀程度和孔隙类型进行判断，达到对未取心井段成岩作用、成岩相准确判别的目的。主要的成岩相测井方法有：
[0004] ①利用蜘蛛网图测井识别成岩相；
[0005] ②利用井径(CAL)、补偿中子测井(CNL)、
密度测井(DEN)、深/浅侧向
电阻率(RLLD/RLLS)、自然伽
马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)测井数据，连续定量识别储层成岩相。
[0006] 目前，通过多种测井资料数据研究了储层成岩相的判别模型，并未进一步精细研究成岩相测井以及成岩相测井数据与储层孔隙、渗透率参数之间关系。
发明内容
[0007] 本发明针对
现有技术存在的
缺陷，提供了一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，该方法通过成岩相测井识别技术，对研究区成岩相进行精细研究，建立研究区现今成岩相剖面(Present DiageneticFacies Profile abbr.PDFP)，再结合 研 究 区 埋 藏 史 、热 史 、流 体 史 、压 力 史 ，恢 复 成 岩 相 演 化 序 列(DiageneticFaciesEvolutionary Sequence abbr.DFES)随时间发育情况，最后结合模拟
算法对井段碎屑岩孔隙度进行模拟计算。
[0008] 为实现上述目的，本发明提供的一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，包括以下步骤：
[0009] 1)收集研究区地质资料
[0011] 在埋藏过程中，根据不同期次、不同类型的成岩作用建立成岩作用数据库；该成岩作用数据库根据不同类型成岩作用由3个小数据库组成，分为压实作用数据库、胶结作用数据库和溶蚀作用数据库：
[0012] (1)压实作用数据库D-Com
[0013] 压实作用数据库根据
温度和压实率建立；
[0014] (2)胶结作用数据D-Cem
[0015] 胶结作用数据库根据温度和不同矿物Mi的胶结率建立；
[0016] (3)溶蚀作用数据库D-Dis
[0017] 溶蚀作用数据库根据温度和不同矿物Mi溶蚀率建立；
[0018] 3)建立成岩相数据库
[0019] 成岩相是指在成岩
环境控制下，由各种成岩作用形成并具有一定几何形态和特定成岩组构、特定成岩矿物组合的地质体。根据成岩作用，成岩相数据库可以表示为：
[0020] F＝∑fi
[0021] 其中：F为成岩相；fi为成岩作用，成岩作用包括压实作用D-Com、溶蚀作用D-Dis和溶蚀作用D-Dis；
[0022] 压实作用D-Com包括ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实和Ⅲ弱压实；
[0023] 胶结作用D-Cem包括ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结和Ⅲ-Mi弱胶结；
[0024] 溶蚀作用D-Dis包括ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀和Ⅲ-Mi弱溶蚀；
[0025] 4)建立成岩相孔隙演化数学模型
[0026] 建立成岩作用(压实作用、
石英胶结作用、
长石胶结作用、
碳酸盐胶结作用、黏土胶结作用、长石溶蚀作用、石英溶蚀作用、沸石溶蚀作用、碳酸盐溶蚀作用)与孔隙演化数学模型，其中各种成岩作用模型如下：
[0027] (1)压实数学模型D1
[0028]
沉积物沉积后，由于上覆沉积物不断加厚，导致储层颗粒排列紧密，从而减少储层孔隙。压实作用与储层原始组构、
流体排出情况、埋藏过程等紧密相关；根据压实减孔量与深度、原始结构之间的关系，可以确定不同埋深下不同碎屑岩减孔量。具体计算公式如下：
[0029]
[0030] 其中：颗粒密度ρs＝2650kg/m3；水密度ρw＝1.0×103kg/m3；k1为常数，k1＝0.67；c为常数，c＝3.68×10-8Pa-1；重力g＝9.8N/kg；VCom为压实作用改变孔隙体积百分量，％；
[0031] (2)溶蚀数学模型D2
[0032] a.石英溶蚀模型D2-1
[0033]
[0034] 其中：VQue_Dis为石英溶蚀作用改变孔隙体积百分含量，％；T为
热力学温度，T＝273+T0+c*(H-H0)，T0为古地温，c为
地温梯度，℃/m，H为地层埋深，H0为恒温带温度；ρH20为常温条件下水密度；
[0035] b.长石溶蚀模型D2-2
[0036]
[0037] 其中：t为距今时间，Ma； 为平均孔隙度，％；t1为地层温度首次达到70℃的时间，Ma；t2为地层温度首次达到90℃对应的时间，Ma；Δt＝t1-t2；VFel_Dis为长石溶蚀作用改变孔隙体积百分含量，％；
[0038] c.方解石溶蚀模型D2-3
[0039] VCal_Dis＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.442；
[0040] 其中：x埋藏深度(Km)；VCal_Dis为方解石溶蚀作用改变孔隙体积百分含量，％；
[0042] VDol_Dis＝-0.0619(X)^6+1.1466(X)^5-7.7836(X)^4+23.231(X)^3-28.939(X)^2+14.147(X)+0.0261；
[0043] 其中：x埋藏深度(Km)；VDol_Dis为白云石溶蚀作用改变孔隙体积百分含量，％；
[0044] (3)胶结数学模型D3
[0045] 建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
[0046] a.石英胶结模型D3-1
[0047]
[0048] 其中：a，b为常数，a＝1.98×10-22，b＝0.022℃(-1)；VQue_Cem为石英胶结作用改变孔隙体积百分量，％；M为石英的摩尔
质量，g/mol；ρ为石英的密度，g/cm3；T为反应温度，℃；T＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；A0为自生石英开始形成时石英颗粒原始表面积；A0＝6fV/D；D为石英颗粒的粒径，mm；V为单位
砂岩体积，mm3；f为单位体积砂岩中石英的体积百分含量；c为平均每年地层温度的平均改变量，℃/y；
[0049] b.长石胶结模型D3-2
[0050] VFel_Cem＝(3.37798E-16*(X)^5-0.00000000000453665*(X)^4+0.0000000218747*(X)^3-0.0000457498*(X)^2+0.049135*(X)+8.67813)/5；
[0051] 其中：x埋藏深度(Km)；VFel_Cem为长石胶结作用改变孔隙体积百分量，％；
[0052] c.方解石胶结模型D3-3
[0053] VCal_Cem＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.44；
[0054] 其中：x埋藏深度(Km)；VCal_Cem为方解石胶结作用改变孔隙体积百分量，％；
[0055] d.白云石胶结模型D3-4
[0056] VDol_Cem＝-0.0619(X)^6+1.1466(X)^5-7.7836(X)4^4+23.231(X)^3-28.939(X)^2+14.147(X)+0.0261；
[0057] 其中：x埋藏深度(Km)；VDol_Cem为白云石胶结作用改变孔隙体积百分量，％；
[0058] e.粘土胶结模型D3-5
[0059] VCla_Cem＝1.22932e-013*(X)^4-1.01647e-009*(X)^3+3.57019e-006*(X)^2-0.00194129*(X)+2.2677；
[0060] 其中：x埋藏深度(Km)；VCla_Cem为黏土胶结作用改变孔隙体积百分量，％；
[0061] 5)确定模拟地层成岩序列
[0062] (1)确定成岩相
[0063] 基于测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井曲线组合，判识不同深度井段成岩相类型。
[0064] F＝∑kili
[0065] 其中：F为通过测井分析识别的成岩相，成岩相为各种成岩作用之和，即：D-Com+D-Cem+D-Dis；
[0066] (2)确定地层模型参数
[0067] a.根据研究区埋藏史恢复地层L在不同演化阶段埋深：
[0068] L(Dep)＝F(X，Y，Z，t)；
[0069] b.根据研究区温度史恢复地层L在不同演化阶段温度：L(Tem)＝G(X，Y，Z，t)；6)[0070] 模拟计算
[0071] (1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
[0072]
[0073] 其中：Φ0为总孔隙度；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
[0074] 当成岩作用为ⅠA强压实、ⅠB强压实、ⅡA强压实、ⅡB强压实、ⅠA-Mi强胶结、ⅠB-Mi强胶结、ⅡA-Mi强胶结、ⅡB-Mi强胶结、Ⅲ-Mi强胶结、ⅠA-Mi强溶蚀、ⅠB-Mi强溶蚀、ⅡA-Mi强溶蚀、ⅡB-Mi强溶蚀、Ⅲ-Mi强溶蚀时，ki＝1.6；
[0075] 当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
[0076] 当成岩作用为ⅠA弱压实、ⅠB弱压实、ⅡA弱压实、ⅡB弱压实、Ⅲ强压实、Ⅲ弱压实、ⅠA-Mi弱胶结、ⅠB-Mi弱胶结、ⅡA-Mi弱胶结、ⅡB-Mi弱胶结、Ⅲ-Mi弱胶结、ⅠA-Mi弱溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi弱溶蚀时，ki＝0.85；
[0077] (2)将时间、深度、温度等参数代入到(1)中压实模型D1、溶蚀模型D2以及胶结模型D3即可计算目标层的孔隙度模拟值。
[0078] 进一步地，所述步骤1)中，收集的地质资料包括：
[0079] 1)研究区地质资料，其中包括
岩心、录井信息、薄片鉴定报告、铸体薄片资料、包裹体分析测试资料等；
[0080] 2)研究区测井数据，包括井径CAL、补偿中子测井CNL、密度测井(DEN)、深/浅侧向电阻率RLLD/RLLS、自然伽马GR、自然伽马能谱测井U、Th、K、声波时差AC和光电界面吸收指数PE数据；
[0081] 3)研究区埋藏史、热演化史、流体史。
[0082] 再进一步地，所述步骤2)中，压实作用数据库D-Com包括ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；其中，压实作用数据库划分标准如下表所示：
[0083] 表1压实作用数据库划分标准
[0084]
[0085] 胶结作用数据库D-Cem包括ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结和Ⅲ-Mi弱胶结，其中，Mi表示矿物类型，胶结作用数据库划分标准如表下所示；
[0086] 表2胶结作用数据库划分标准
[0087]
[0088] 溶蚀作用数据库D-Dis包括ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀，其中Mi表示矿物类型；溶蚀作用数据库划分标准如下表所示：
[0089] 表3溶蚀作用类型及划分标准
[0090]
[0091]
[0092] 再进一步地，所述步骤5)第(1)小步中，压实作用D-Com包括：ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；
[0093] 胶结作用D-Cem包括：ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结、Ⅲ-Mi弱胶结；
[0094] 溶蚀作用D-Dis包括：ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀；
[0095] li为测井曲线数据，包括CAL、CNL、DEN、RLLD/RLLS、GR、U、Th、K、AC和PE测井曲线；ki为系数，ki＝0，或者ki＝1；
[0096] 判识不同深度井段成岩相类型具体标注如下：
[0097] ⅠA强压实＝CAL∈(250，310)+CNL∈(30，40)+DEN∈(2.5，3.5)+GR∈(20，80)+AC∈(150，500)；
[0098] ⅠA中压实＝CAL∈(120，300)+CNL∈(30，35)+DEN∈(1.8，2.9)+GR∈(20，100)+AC∈(160，550)；
[0099] ⅠA弱压实＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，150)+AC∈(150，550)；
[0100] ⅠB强压实＝CAL∈(10，100)+CNL∈(30，45)+DEN∈(1.9，2.5)+GR∈(20，110)+AC∈(170，300)；
[0101] ⅠB中压实＝CAL∈(20，110)+CNL∈(30，50)+DEN∈(1.9，2.6)+GR∈(20，120)+AC∈(180，310)；
[0102] ⅠB弱压实＝CAL∈(30，120)+CNL∈(35，40)+DEN∈(1.9，2.7)+GR∈(20，130)+AC∈(190，320)；
[0103] ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
[0104] ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
[0105] ⅡA弱压实＝CAL∈(60，150)+CNL∈(40，45)+DEN∈(1.9，3.0)+GR∈(30，150)+AC∈(220，350)；
[0106] ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
[0107] ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
[0108] ⅡB弱压实＝CAL∈(90，180)+CNL∈(40，50)+DEN∈(1.5，2.3)+GR∈(40，70)+AC∈(250，380)；
[0109] Ⅲ强压实＝CAL∈(100，190)+CNL∈(30，45)+DEN∈(1.5，2.4)+GR∈(50，80)+AC∈(260，390)；
[0110] Ⅲ中压实＝CAL∈(110，200)+CNL∈(30，50)+DEN∈(1.5，2.5)+GR∈(60，150)+AC∈(270，400)；
[0111] Ⅲ弱压实＝CAL∈(120，210)+CNL∈(45，50)+DEN∈(1.5，2.6)+GR∈(70，150)+AC∈(280，410)；
[0112] ⅠA-石英强胶结＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，20)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)+AC∈(150，550)；
[0113] ⅠA-石英中胶结＝CAL∈(100，130)+CNL∈(30，35)+RLLD/RLLS∈(0.2，1)+GR∈(30，40)+U∈(10，20)+Th∈(35，40)+K∈(40，45)+AC∈(200，400)；
[0114] ⅠA-石英弱胶结＝CAL∈(110，140)+CNL∈(30，40)+RLLD/RLLS∈(0.3，2)+GR∈(40，50)+U∈(20，30)+Th∈(40，45)+K∈(45，50)+AC∈(210，410)；
[0115] ⅠB-石英强胶结＝CAL∈(120，150)+CNL∈(30，45)+RLLD/RLLS∈(0.4，3)+GR∈(50，60)+U∈(30，40)+Th∈(45，50)+K∈(50，55)+AC∈(220，420)；
[0116] ⅠB-石英中胶结＝CAL∈(130，160)+CNL∈(35，40)+RLLD/RLLS∈(0.5，4)+GR∈(60，70)+U∈(40，50)+Th∈(50，55)+K∈(55，60)+AC∈(230，430)；
[0117] ⅠB-石英弱胶结＝CAL∈(140，170)+CNL∈(35，45)+RLLD/RLLS∈(0.6，5)+GR∈(70，80)+U∈(50，60)+Th∈(55，60)+K∈(60，65)+AC∈(240，440)；
[0118] ⅡA-石英强胶结＝CAL∈(150，180)+CNL∈(35，50)+RLLD/RLLS∈(0.7，6)+GR∈(80，90)+U∈(60，70)+Th∈(60，65)+K∈(65，70)+AC∈(250，450)；
[0119] ⅡA-石英中胶结＝CAL∈(160，190)+CNL∈(40，50)+RLLD/RLLS∈(0.8，7)+GR∈(90，100)+U∈(10，70)+Th∈(65，70)+K∈(40，75)+AC∈(260，460)；
[0120] ⅡA-石英弱胶结＝CAL∈(170，200)+CNL∈(40，45)+RLLD/RLLS∈(0.9，8)+GR∈(100，110)+U∈(80，90)+Th∈(70，75)+K∈(75，80)+AC∈(270，470)；
[0121] ⅡB-石英强胶结＝CAL∈(180，210)+CNL∈(45，50)+RLLD/RLLS∈(1.0，9)+GR∈(110，120)+U∈(90，100)+Th∈(75，80)+K∈(30，40)+AC∈(280，480)；
[0122] ⅡB-石英中胶结＝CAL∈(190，220)+CNL∈(30，40)+RLLD/RLLS∈(1.1，10)+GR∈(120，130)+U∈(20，40)+Th∈(40，50)+K∈(40，50)+AC∈(290，490)；
[0123] ⅡB-石英弱胶结＝CAL∈(200，230)+CNL∈(30，50)+RLLD/RLLS∈(1.2，11)+GR∈(130，140)+U∈(30，50)+Th∈(50，60)+K∈(50，60)+AC∈(300，500)；
[0124] Ⅲ-石英强胶结＝CAL∈(210，240)+CNL∈(40，50)+RLLD/RLLS∈(1.3，12)+GR∈(50，100)+U∈(40，60)+Th∈(60，70)+K∈(60，70)+AC∈(310，510)；
[0125] Ⅲ-石英中胶结＝CAL∈(220，250)+CNL∈(40，45)+RLLD/RLLS∈(1.4，13)+GR∈(60，110)+U∈(50，70)+Th∈(70，80)+K∈(70，80)+AC∈(320，520)；
[0126] Ⅲ-石英弱胶结＝CAL∈(230，260)+CNL∈(30，40)+RLLD/RLLS∈(1.5，14)+GR∈(70，120)+U∈(60，80)+Th∈(50，80)+K∈(40，80)+AC∈(330，530)；
[0127] ⅠA-碳酸盐强胶结＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，30)+U∈(0，100)+K∈(30，80)+PE∈(10，100)；
[0128] ⅠA-碳酸盐中胶结＝CAL∈(10，50)+CNL∈(30，35)+DEN∈(1.5，2.0)+GR∈(20，40)+U∈(10，20)+K∈(40，50)+PE∈(20，30)；
[0129] ⅠA-碳酸盐弱胶结＝CAL∈(20，60)+CNL∈(30，40)+DEN∈(1.5，2.1)+GR∈(40，50)+U∈(20，30)+K∈(50，60)+PE∈(30，40)；
[0130] ⅠB-碳酸盐强胶结＝CAL∈(30，70)+CNL∈(30，45)+DEN∈(1.5，2.2)+GR∈(30，60)+U∈(30，40)+K∈(60，70)+PE∈(40，50)；
[0131] ⅠB-碳酸盐中胶结＝CAL∈(40，80)+CNL∈(35，40)+DEN∈(1.5，2.3)+GR∈(60，70)+U∈(40，50)+K∈(70，80)+PE∈(50，100)；
[0132] ⅠB-碳酸盐弱胶结＝CAL∈(50，90)+CNL∈(35，45)+DEN∈(1.5，2.4)+GR∈(50，80)+U∈(50，60)+K∈(30，50)+PE∈(60，70)；
[0133] ⅡA-碳酸盐强胶结＝CAL∈(60，100)+CNL∈(35，50)+DEN∈(1.5，2.5)+GR∈(40，90)+U∈(60，70)+K∈(30，60)+PE∈(70，80)；
[0134] ⅡA-碳酸盐中胶结＝CAL∈(70，110)+CNL∈(40，45)+DEN∈(1.5，2.6)+GR∈(30，100)+U∈(70，80)+K∈(30，70)+PE∈(80，90)；
[0135] ⅡA-碳酸盐弱胶结＝CAL∈(80，120)+CNL∈(40，50)+DEN∈(1.5，2.7)+GR∈(90，150)+U∈(0，100)+K∈(30，80)+PE∈(90，100)；
[0136] ⅡB-碳酸盐强胶结＝CAL∈(90，130)+CNL∈(30，40)+DEN∈(1.5，2.8)+GR∈(80，110)+U∈(30，50)+K∈(30，60)+PE∈(20，50)；
[0137] ⅡB-碳酸盐中胶结＝CAL∈(100，140)+CNL∈(40，50)+DEN∈(1.5，2.9)+GR∈(70，120)+U∈(40，60)+K∈(30，70)+PE∈(30，40)；
[0138] ⅡB-碳酸盐弱胶结＝CAL∈(110，150)+CNL∈(35，50)+DEN∈(1.5，3.0)+GR∈(60，130)+U∈(50，70)+K∈(40，80)+PE∈(40，60)；
[0139] Ⅲ-碳酸盐强胶结＝CAL∈(120，160)+CNL∈(30，35)+DEN∈(1.5，3.1)+GR∈(50，140)+U∈(60，80)+K∈(50，70)+PE∈(50，70)；
[0140] Ⅲ-碳酸盐中胶结＝CAL∈(130，170)+CNL∈(40，50)+DEN∈(1.5，3.2)+GR∈40，150)+U∈(70，90)+K∈(50，80)+PE∈(60，80)；
[0141] Ⅲ-碳酸盐弱胶结＝CAL∈(140，180)+CNL∈(45，50)+DEN∈(1.5，3.3)+GR∈(30，150)+U∈(80，100)+K∈(60，80)+PE∈(70，100)；
[0142] ⅠA-黏土强胶结＝GR∈(80，150)+U∈(60，80)+Th∈(30，80)+K∈(60，90)；
[0143] ⅠA-黏土中胶结＝GR∈(80，120)+U∈(0，50)+Th∈(50，70)+K∈(30，80)；
[0144] ⅠA-黏土弱胶结＝GR∈(50，110)+U∈(70，90)+Th∈(30，80)+K∈(40，85)；
[0145] ⅠB-黏土强胶结＝GR∈(20，150)+U∈(80，100)+Th∈(30，70)+K∈(30，80)；
[0146] ⅠB-黏土中胶结＝GR∈(70，120)+U∈(0，20)+Th∈(30，80)+K∈(30，70)；
[0147] ⅠB-黏土弱胶结＝GR∈(80，110)+U∈(20，55)+Th∈(30，60)+K∈(20，70)；
[0148] ⅡA-黏土强胶结＝GR∈(30，50)+U∈(35，100)+Th∈(30，55)+K∈(30，80)；
[0149] ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
[0150] ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
[0151] ⅡB-黏土强胶结＝GR∈(60，50)+U∈(40，60)+Th∈(30，70)+K∈(30，65)；
[0152] ⅡB-黏土中胶结＝GR∈(50，150)+U∈(0，80)+Th∈(20，80)+K∈(45，75)；
[0153] ⅡB-黏土弱胶结＝GR∈(40，50)+U∈(40，70)+Th∈(40，60)+K∈(30，50)；
[0154] Ⅲ-黏土强胶结＝GR∈(110，150)+U∈(70，100)+Th∈(30，70)+K∈(50，90)；
[0155] Ⅲ-黏土中胶结＝GR∈(60，110)+U∈(40，70)+Th∈(30，40)+K∈(30，65)；
[0156] Ⅲ-黏土弱胶结＝GR∈(20，50)+U∈(80，100)+Th∈(30，50)+K∈(30，50)；
[0157] ⅠA-石英强溶蚀＝CAL∈(280，310)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(80，150)+U∈(0，50)+Th∈(30，55)+K∈(30，70)；
[0158] ⅠA-石英中溶蚀＝CAL∈(90，150)+CNL∈(40，70)+RLLD/RLLS∈(15，20)+GR∈(20，110)+U∈(70，90)+Th∈(30，80)+K∈(30，60)；
[0159] ⅠA-石英弱溶蚀＝CAL∈(250，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(40，150)+U∈(45，90)+Th∈(35，50)+K∈(50，80)；
[0160] ⅠB-石英强溶蚀＝CAL∈(30，150)+CNL∈(30，50)+RLLD/RLLS∈(8，15)+GR∈(20，150)+U∈(60，100)+Th∈(30，80)+K∈(30，50)；
[0161] ⅠB-石英中溶蚀＝CAL∈(100，240)+CNL∈(30，50)+RLLD/RLLS∈(0.2，5)+GR∈(20，110)+U∈(50，90)+Th∈(30，55)+K∈(60，70)；
[0162] ⅠB-石英弱溶蚀＝CAL∈(90，350)+CNL∈(30，45)+RLLD/RLLS∈(11，18)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)；
[0163] ⅡA-石英强溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(14，20)+GR∈(100，150)+U∈(60，90)+Th∈(30，80)+K∈(30，50)；
[0164] ⅡA-石英中溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，8)+GR∈(100，160)+U∈(0，50)+Th∈(30，80)+K∈(40，100)；
[0165] ⅡA-石英弱溶蚀＝CAL∈(250，315)+CNL∈(30，50)+RLLD/RLLS∈(15，20)+GR∈(100，150)+U∈(55，100)+Th∈(60，80)+K∈(30，80)；
[0166] ⅡB-石英强溶蚀＝CAL∈(100，220)+CNL∈(15，45)+RLLD/RLLS∈(10，20)+GR∈(20，100)+U∈(80，100)+Th∈(50，80)+K∈(30，65)；
[0167] ⅡB-石英中溶蚀＝CAL∈(0，350)+CNL∈(30，80)+RLLD/RLLS∈(10，20)+GR∈(80，120)+U∈(50，100)+Th∈(30，80)+K∈(30，70)；
[0168] ⅡB-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，50)+RLLD/RLLS∈(0.2，8)+GR∈(50，150)+U∈(80，100)+Th∈(20，70)+K∈(30，80)；
[0169] Ⅲ-石英强溶蚀＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(11，20)+GR∈(90，140)+U∈(70，100)+Th∈(30，60)+K∈(20，65)；
[0170] Ⅲ-石英中溶蚀＝CAL∈(0，350)+CNL∈(50，90)+RLLD/RLLS∈(15，20)+GR∈(80，120)+U∈(80，100)+Th∈(30，70)+K∈(25，90)；
[0171] Ⅲ-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，80)+RLLD/RLLS∈(3，8)+GR∈(110，145)+U∈(50，100)+Th∈(20，80)+K∈(10，55)；
[0172] ⅠA-长石强溶蚀＝CAL∈(50，150)+CNL∈(30，50)+GR∈(20，80)+U∈(0，80)+AC∈(100，450)+PE∈(20，50)；
[0173] ⅠA-长石中溶蚀＝CAL∈(0，80)+CNL∈(30，40)+GR∈(20，120)+U∈(60，100)+AC∈(150，250)+PE∈(40，90)；
[0174] ⅠA-长石弱溶蚀＝CAL∈(240，325)+CNL∈(60，80)+GR∈(20，125)+U∈(40，50)+AC∈(250，550)+PE∈(20，85)；
[0175] ⅠB-长石强溶蚀＝CAL∈(90，255)+CNL∈(25，50)+GR∈(20，140)+U∈(40，80)+AC∈(450，550)+PE∈(60，80)；
[0176] ⅠB-长石中溶蚀＝CAL∈(200，235)+CNL∈(30，50)+GR∈(20，130)+U∈(45，85)+AC∈(100，550)+PE∈(45，90)；
[0177] ⅠB-长石弱溶蚀＝CAL∈(200，350)+CNL∈(10，45)+GR∈(20，70)+U∈(20，70)+AC∈(900，210)+PE∈(30，70)；
[0178] ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
[0179] ⅡA-长石中溶蚀＝CAL∈(300，350)+CNL∈(15，40)+GR∈(20，115)+U∈(35，60)+AC∈(270，450)+PE∈(70，95)；
[0180] ⅡA-长石弱溶蚀＝CAL∈(180，350)+CNL∈(30，50)+GR∈(20，85)+U∈(15，50)+AC∈(200，400)+PE∈(40，80)
[0181] ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
[0182] ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
[0183] ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
[0184] Ⅲ-长石强溶蚀＝CAL∈(70，140)+CNL∈(30，55)+GR∈(20，70)+U∈(50，60)+AC∈(60，450)+PE∈(70，90)；
[0185] Ⅲ-长石中溶蚀＝CAL∈(280，350)+CNL∈(0，20)+GR∈(20，125)+U∈(70，100)+AC∈(50，180)+PE∈(50，80)；
[0186] Ⅲ-长石弱溶蚀＝CAL∈(255，350)+CNL∈(30，50)+GR∈(20，150)+U∈(40，75)+AC∈(150，550)+PE∈(55，85)；
[0187] ⅠA-碳酸盐强溶蚀＝CAL∈(100，150)+U∈(80，100)+Th∈(30，55)+K∈(35，80)+AC∈(450，550)+PE∈(0，50)；
[0188] ⅠA-碳酸盐中溶蚀＝CAL∈(250，300)+U∈(70，100)+Th∈(0，38)+K∈(30，60)+AC∈(150，250)+PE∈(80，100)；
[0189] ⅠA-碳酸盐弱溶蚀＝CAL∈(280，320)+U∈(60，80)+Th∈(10，80)+K∈(40，80)+AC∈(100，300)+PE∈(45，100)；
[0190] ⅠB-碳酸盐强溶蚀＝CAL∈(250，350)+U∈(0，50)+Th∈(15，40)+K∈(65，90)+AC∈(10，100)+PE∈(0，60)；
[0191] ⅠB-碳酸盐中溶蚀＝CAL∈(240，300)+U∈(5，70)+Th∈(30，80)+K∈(0，15)+AC∈(150，240)+PE∈(20，100)；
[0192] ⅠB-碳酸盐弱溶蚀＝CAL∈(180，210)+U∈(50，90)+Th∈(30，80)+K∈(30，55)+AC∈(100，150)+PE∈(60，80)；
[0193] ⅡA-碳酸盐强溶蚀＝CAL∈(140，200)+U∈(40，100)+Th∈(60，95)+K∈(30，60)+AC∈(80，320)+PE∈(0，100)；
[0194] ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
[0195] ⅡA-碳酸盐弱溶蚀＝CAL∈(80，110)+U∈(80，100)+Th∈(40，80)+K∈(10，80)+AC∈(70，250)+PE∈(0，100)；
[0196] ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
[0197] ⅡB-碳酸盐中溶蚀＝CAL∈(200，315)+U∈(45，100)+Th∈(35，50)+K∈(30，80)+AC∈(10，80)+PE∈(10，60)；
[0198] ⅡB-碳酸盐弱溶蚀＝CAL∈(240，310)+U∈(70，90)+Th∈(30，80)+K∈(5，70)+AC∈(220，340)+PE∈(15，30)
[0199] Ⅲ-碳酸盐强溶蚀＝CAL∈(270，300)+U∈(35，80)+Th∈(20，85)+K∈(10，80)+AC∈(150，260)+PE∈(50，70)；
[0200] Ⅲ-碳酸盐中溶蚀＝CAL∈(240，350)+U∈(40，80)+Th∈(37，80)+K∈(5，50)+AC∈(310，400)+PE∈(80，95)；
[0201] Ⅲ-碳酸盐弱溶蚀＝CAL∈(110，140)+U∈(60，100)+Th∈(30，60)+K∈(30，60)+AC∈(200，375)+PE∈(55，90)。
[0202] 本发明的有益效果在于：
[0203] 1、预测性强，本研究基于测井成岩相的研究，对预测结果进行进一步处理，能直接对各类碎屑岩储层孔隙进行计算，尤其是低孔、低渗碎屑岩储层。
[0204] 2、方法创新，本研究通过成岩相研究把测井解释和储层孔隙度预测联系在一起，有利于该方法在单井孔隙度预测上得到很好应用。
[0205] 3、经济效益，该方法操作简单，避免了传统方法需要进行大量取芯和分析测试，节约了油气勘探开发成本。
附图说明
[0206] 图1为井98成岩作用、成岩相的示意图；
[0207] 图2为研究区埋藏史以及热史图。
具体实施方式
[0208] 为了更好地解释本发明，以下结合具体
实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。
[0209] 实施例1
[0210] 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，包括以下步骤：
[0211] 1)收集研究区地质资料
[0212] (1)研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料、包裹体分析测试资料等；
[0213] (2)研究区测井数据，其中包括井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)数据；
[0214] (3)研究区埋藏史、热演化史、流体史。
[0215] 2)建立成岩作用数据库
[0216] 在埋藏过程中，根据不同期次、不同类型的成岩作用建立成岩作用数据库；该成岩作用数据库根据不同类型成岩作用由3个小数据库组成，分为压实作用数据库、胶结作用数据库和溶蚀作用数据库：
[0217] (1)压实作用数据库D-Com
[0218] 压实作用数据库根据压实率建立，包括ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；其中，压实作用数据库划分标准如表1所示。
[0219] 表1压实作用数据库划分标准
[0220]
[0221]
[0222] (2)胶结作用数据D-Cem
[0223] 胶结作用数据库根据不同矿物Mi的胶结率建立，包括：ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结、Ⅲ-Mi弱胶结，其中，Mi表示矿物类型，胶结作用数据库划分标准如表2所示。
[0224] 表2胶结作用数据库划分标准
[0225]
[0226] (3)溶蚀作用数据库D-Dis
[0227] 溶蚀作用数据库根据不同矿物Mi溶蚀率建立，类型包括：ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀，其中Mi表示矿物类型；溶蚀作用数据库划分标准如表3所示。
[0228] 表3溶蚀作用类型及划分标准
[0229]
[0230]
[0231] 3)建立成岩相数据库
[0232] 成岩相是指在成岩环境控制下，由各种成岩作用形成并具有一定几何形态和特定成岩组构、特定成岩矿物组合的地质体。根据成岩作用，成岩相数据库可以表示为：
[0233] F＝∑fi
[0234] 其中：F为成岩相；fi为成岩作用，成岩作用为压实作用D-Com、胶结作用D-Cem和溶蚀作用D-Dis；
[0235] 压实作用D-Com(ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实)、胶结作用D-Cem(ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结、Ⅲ-Mi弱胶结)、溶蚀作用D-Dis(ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀)。
[0236] 4)建立研究区成岩相孔隙演化数学模型
[0237] 建立成岩作用(压实作用、石英胶结作用、黏土胶结作用、长石溶蚀作用、方解石溶蚀作用)与孔隙演化数学模型，其中各种成岩作用模型如下：
[0238] (1)压实数学模型D1
[0239] 根据压实减孔量与深度、原始结构之间的关系，确定不同埋深下不同碎屑岩减孔量，具体计算公式如下：
[0240]
[0241] 其中：颗粒密度ρs＝2650kg/m3；水密度ρw＝1.0×103kg/m3；k1为常数，k1＝0.67；c为常数，c＝3.68×10-8Pa-1；重力g＝9.8N/kg；VCom为压实作用改变孔隙体积百分量，％；
[0242] (2)溶蚀数学模型D2
[0243] a.长石溶蚀模型D2-2
[0244]
[0245] 其中：t为距今时间，Ma； 为平均孔隙度，％；t1为地层温度首次达到70℃的时间，Ma；t2为地层温度首次达到90℃对应的时间，Ma；Δt＝t1-t2；VFel_Dis为长石溶蚀作用改变孔隙体积百分含量，％；
[0246] b.方解石溶蚀模型D2-3
[0247] VCal_Dis＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.442；
[0248] 其中：x埋藏深度(Km)；VCal_Dis为方解石溶蚀作用改变孔隙体积百分含量，％；
[0249] (3)胶结数学模型D3
[0250] 建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
[0251] a.石英胶结模型D3-1
[0252]
[0253] 其中：a，b为常数，a＝1.98×10-22，b＝0.022℃(-1)；VQue_Cem为石英胶结作用改变孔隙体积百分量，％；M为石英的摩尔质量，g/mol；ρ为石英的密度，g/cm3；T为反应温度，℃；T＝cn*z+dn，cn为地温梯度，℃/m；dn为地表温度或恒温带温度，℃；z为地层埋藏深度，m；A0为自生石英开始形成时石英颗粒原始表面积；A0＝6fV/D；D为石英颗粒的粒径，mm；V为单位3
砂岩体积，mm ；f为单位体积砂岩中石英的体积百分含量；c为平均每年地层温度的平均改变量，℃/y；
[0254] b.粘土胶结模型D3-5
[0255] VCla_Cem＝1.22932e-013*(X)^4-1.01647e-009*(X)^3+3.57019e-006*(X)^2-0.00194129*(X)+2.2677；
[0256] 其中：x埋藏深度(Km)；VCla_Cem为黏土胶结作用改变孔隙体积百分量，％；
[0257] (4)主要参数选取
[0258] a.压实作用模型主要参数选取
[0259] ρs＝2650kg/m3；
[0260] ρw＝1.0×103kg/m3；
[0261] k1＝0.67；
[0262] c＝3.68×10-8Pa-1；
[0263] g＝9.8N/kg；
[0264] b.石英胶结模型主要参数选取
[0265] a＝1.98×10-22；
[0266] b＝0.022℃(-1)；
[0267] M＝60g/mol；
[0268] ρ＝0.25g/cm3；
[0269] c＝0.00000000000002℃/y；
[0270] Cn＝0.0031℃/m；
[0271] Dn＝18℃；
[0272]
[0273] D＝0.02mm；
[0274] V＝1mm3；
[0275] f＝0.65；
[0276] T1＝90℃；
[0277] T2＝140℃；
[0278] c.长石溶蚀模型主要参数选取
[0279]
[0280] t1＝16.2
[0281] t2＝8.2
[0282] t＝28.1
[0283] 5)确定模拟地层成岩序列
[0284] (1)确定成岩相
[0285] 基于测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井曲线组合，判识不同深度井段成岩相类型。
[0286] F＝∑kili
[0287] 其中：F为通过测井分析识别的成岩相，成岩相为各种成岩作用之和，即：D-Com+D-Cem+D-Dis；
[0288] 压实作用D-Com包括：ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；胶结作用D-Cem包括：ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结、Ⅲ-Mi弱胶结；溶蚀作用D-Dis包括：ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀。
[0289] li为测井曲线数据，包括CAL、CNL、DEN、RLLD/RLLS、GR、U、Th、K、AC和PE测井曲线；ki为系数，ki＝0，或者ki＝1。
[0290] 各成岩作用具体判识如下：
[0291] ⅠA强压实＝CAL∈(250，310)+CNL∈(30，40)+DEN∈(2.5，3.5)+GR∈(20，80)+AC∈(150，500)；
[0292] ⅠA中压实＝CAL∈(120，300)+CNL∈(30，35)+DEN∈(1.8，2.9)+GR∈(20，100)+AC∈(160，550)；
[0293] ⅠA弱压实＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，150)+AC∈(150，550)；
[0294] ⅠB强压实＝CAL∈(10，100)+CNL∈(30，45)+DEN∈(1.9，2.5)+GR∈(20，110)+AC∈(170，300)；
[0295] ⅠB中压实＝CAL∈(20，110)+CNL∈(30，50)+DEN∈(1.9，2.6)+GR∈(20，120)+AC∈(180，310)；
[0296] ⅠB弱压实＝CAL∈(30，120)+CNL∈(35，40)+DEN∈(1.9，2.7)+GR∈(20，130)+AC∈(190，320)；
[0297] ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
[0298] ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
[0299] ⅡA弱压实＝CAL∈(60，150)+CNL∈(40，45)+DEN∈(1.9，3.0)+GR∈(30，150)+AC∈(220，350)；
[0300] ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
[0301] ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
[0302] ⅡB弱压实＝CAL∈(90，180)+CNL∈(40，50)+DEN∈(1.5，2.3)+GR∈(40，70)+AC∈(250，380)；
[0303] Ⅲ强压实＝CAL∈(100，190)+CNL∈(30，45)+DEN∈(1.5，2.4)+GR∈(50，80)+AC∈(260，390)；
[0304] Ⅲ中压实＝CAL∈(110，200)+CNL∈(30，50)+DEN∈(1.5，2.5)+GR∈(60，150)+AC∈(270，400)；
[0305] Ⅲ弱压实＝CAL∈(120，210)+CNL∈(45，50)+DEN∈(1.5，2.6)+GR∈(70，150)+AC∈(280，410)；
[0306] ⅠA-石英强胶结＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，20)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)+AC∈(150，550)；
[0307] ⅠA-石英中胶结＝CAL∈(100，130)+CNL∈(30，35)+RLLD/RLLS∈(0.2，1)+GR∈(30，40)+U∈(10，20)+Th∈(35，40)+K∈(40，45)+AC∈(200，400)；
[0308] ⅠA-石英弱胶结＝CAL∈(110，140)+CNL∈(30，40)+RLLD/RLLS∈(0.3，2)+GR∈(40，50)+U∈(20，30)+Th∈(40，45)+K∈(45，50)+AC∈(210，410)；
[0309] ⅠB-石英强胶结＝CAL∈(120，150)+CNL∈(30，45)+RLLD/RLLS∈(0.4，3)+GR∈(50，60)+U∈(30，40)+Th∈(45，50)+K∈(50，55)+AC∈(220，420)；
[0310] ⅠB-石英中胶结＝CAL∈(130，160)+CNL∈(35，40)+RLLD/RLLS∈(0.5，4)+GR∈(60，70)+U∈(40，50)+Th∈(50，55)+K∈(55，60)+AC∈(230，430)；
[0311] ⅠB-石英弱胶结＝CAL∈(140，170)+CNL∈(35，45)+RLLD/RLLS∈(0.6，5)+GR∈(70，80)+U∈(50，60)+Th∈(55，60)+K∈(60，65)+AC∈(240，440)；
[0312] ⅡA-石英强胶结＝CAL∈(150，180)+CNL∈(35，50)+RLLD/RLLS∈(0.7，6)+GR∈(80，90)+U∈(60，70)+Th∈(60，65)+K∈(65，70)+AC∈(250，450)；
[0313] ⅡA-石英中胶结＝CAL∈(160，190)+CNL∈(40，50)+RLLD/RLLS∈(0.8，7)+GR∈(90，100)+U∈(10，70)+Th∈(65，70)+K∈(40，75)+AC∈(260，460)；
[0314] ⅡA-石英弱胶结＝CAL∈(170，200)+CNL∈(40，45)+RLLD/RLLS∈(0.9，8)+GR∈(100，110)+U∈(80，90)+Th∈(70，75)+K∈(75，80)+AC∈(270，470)；
[0315] ⅡB-石英强胶结＝CAL∈(180，210)+CNL∈(45，50)+RLLD/RLLS∈(1.0，9)+GR∈(110，120)+U∈(90，100)+Th∈(75，80)+K∈(30，40)+AC∈(280，480)；
[0316] ⅡB-石英中胶结＝CAL∈(190，220)+CNL∈(30，40)+RLLD/RLLS∈(1.1，10)+GR∈(120，130)+U∈(20，40)+Th∈(40，50)+K∈(40，50)+AC∈(290，490)；
[0317] ⅡB-石英弱胶结＝CAL∈(200，230)+CNL∈(30，50)+RLLD/RLLS∈(1.2，11)+GR∈(130，140)+U∈(30，50)+Th∈(50，60)+K∈(50，60)+AC∈(300，500)；
[0318] Ⅲ-石英强胶结＝CAL∈(210，240)+CNL∈(40，50)+RLLD/RLLS∈(1.3，12)+GR∈(50，100)+U∈(40，60)+Th∈(60，70)+K∈(60，70)+AC∈(310，510)；
[0319] Ⅲ-石英中胶结＝CAL∈(220，250)+CNL∈(40，45)+RLLD/RLLS∈(1.4，13)+GR∈(60，110)+U∈(50，70)+Th∈(70，80)+K∈(70，80)+AC∈(320，520)；
[0320] Ⅲ-石英弱胶结＝CAL∈(230，260)+CNL∈(30，40)+RLLD/RLLS∈(1.5，14)+GR∈(70，120)+U∈(60，80)+Th∈(50，80)+K∈(40，80)+AC∈(330，530)；
[0321] ⅠA-碳酸盐强胶结＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，30)+U∈(0，100)+K∈(30，80)+PE∈(10，100)；
[0322] ⅠA-碳酸盐中胶结＝CAL∈(10，50)+CNL∈(30，35)+DEN∈(1.5，2.0)+GR∈(20，40)+U∈(10，20)+K∈(40，50)+PE∈(20，30)；
[0323] ⅠA-碳酸盐弱胶结＝CAL∈(20，60)+CNL∈(30，40)+DEN∈(1.5，2.1)+GR∈(40，50)+U∈(20，30)+K∈(50，60)+PE∈(30，40)；
[0324] ⅠB-碳酸盐强胶结＝CAL∈(30，70)+CNL∈(30，45)+DEN∈(1.5，2.2)+GR∈(30，60)+U∈(30，40)+K∈(60，70)+PE∈(40，50)；
[0325] ⅠB-碳酸盐中胶结＝CAL∈(40，80)+CNL∈(35，40)+DEN∈(1.5，2.3)+GR∈(60，70)+U∈(40，50)+K∈(70，80)+PE∈(50，100)；
[0326] ⅠB-碳酸盐弱胶结＝CAL∈(50，90)+CNL∈(35，45)+DEN∈(1.5，2.4)+GR∈(50，80)+U∈(50，60)+K∈(30，50)+PE∈(60，70)；
[0327] ⅡA-碳酸盐强胶结＝CAL∈(60，100)+CNL∈(35，50)+DEN∈(1.5，2.5)+GR∈(40，90)+U∈(60，70)+K∈(30，60)+PE∈(70，80)；
[0328] ⅡA-碳酸盐中胶结＝CAL∈(70，110)+CNL∈(40，45)+DEN∈(1.5，2.6)+GR∈(30，100)+U∈(70，80)+K∈(30，70)+PE∈(80，90)；
[0329] ⅡA-碳酸盐弱胶结＝CAL∈(80，120)+CNL∈(40，50)+DEN∈(1.5，2.7)+GR∈(90，150)+U∈(0，100)+K∈(30，80)+PE∈(90，100)；
[0330] ⅡB-碳酸盐强胶结＝CAL∈(90，130)+CNL∈(30，40)+DEN∈(1.5，2.8)+GR∈(80，110)+U∈(30，50)+K∈(30，60)+PE∈(20，50)；
[0331] ⅡB-碳酸盐中胶结＝CAL∈(100，140)+CNL∈(40，50)+DEN∈(1.5，2.9)+GR∈(70，120)+U∈(40，60)+K∈(30，70)+PE∈(30，40)；
[0332] ⅡB-碳酸盐弱胶结＝CAL∈(110，150)+CNL∈(35，50)+DEN∈(1.5，3.0)+GR∈(60，130)+U∈(50，70)+K∈(40，80)+PE∈(40，60)；
[0333] Ⅲ-碳酸盐强胶结＝CAL∈(120，160)+CNL∈(30，35)+DEN∈(1.5，3.1)+GR∈(50，140)+U∈(60，80)+K∈(50，70)+PE∈(50，70)；
[0334] Ⅲ-碳酸盐中胶结＝CAL∈(130，170)+CNL∈(40，50)+DEN∈(1.5，3.2)+GR∈40，150)+U∈(70，90)+K∈(50，80)+PE∈(60，80)；
[0335] Ⅲ-碳酸盐弱胶结＝CAL∈(140，180)+CNL∈(45，50)+DEN∈(1.5，3.3)+GR∈(30，150)+U∈(80，100)+K∈(60，80)+PE∈(70，100)；
[0336] ⅠA-黏土强胶结＝GR∈(80，150)+U∈(60，80)+Th∈(30，80)+K∈(60，90)；
[0337] ⅠA-黏土中胶结＝GR∈(80，120)+U∈(0，50)+Th∈(50，70)+K∈(30，80)；
[0338] ⅠA-黏土弱胶结＝GR∈(50，110)+U∈(70，90)+Th∈(30，80)+K∈(40，85)；
[0339] ⅠB-黏土强胶结＝GR∈(20，150)+U∈(80，100)+Th∈(30，70)+K∈(30，80)；
[0340] ⅠB-黏土中胶结＝GR∈(70，120)+U∈(0，20)+Th∈(30，80)+K∈(30，70)；
[0341] ⅠB-黏土弱胶结＝GR∈(80，110)+U∈(20，55)+Th∈(30，60)+K∈(20，70)；
[0342] ⅡA-黏土强胶结＝GR∈(30，50)+U∈(35，100)+Th∈(30，55)+K∈(30，80)；
[0343] ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
[0344] ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
[0345] ⅡB-黏土强胶结＝GR∈(60，50)+U∈(40，60)+Th∈(30，70)+K∈(30，65)；
[0346] ⅡB-黏土中胶结＝GR∈(50，150)+U∈(0，80)+Th∈(20，80)+K∈(45，75)；
[0347] ⅡB-黏土弱胶结＝GR∈(40，50)+U∈(40，70)+Th∈(40，60)+K∈(30，50)；
[0348] Ⅲ-黏土强胶结＝GR∈(110，150)+U∈(70，100)+Th∈(30，70)+K∈(50，90)；
[0349] Ⅲ-黏土中胶结＝GR∈(60，110)+U∈(40，70)+Th∈(30，40)+K∈(30，65)；
[0350] Ⅲ-黏土弱胶结＝GR∈(20，50)+U∈(80，100)+Th∈(30，50)+K∈(30，50)；
[0351] ⅠA-石英强溶蚀＝CAL∈(280，310)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(80，150)+U∈(0，50)+Th∈(30，55)+K∈(30，70)；
[0352] ⅠA-石英中溶蚀＝CAL∈(90，150)+CNL∈(40，70)+RLLD/RLLS∈(15，20)+GR∈(20，110)+U∈(70，90)+Th∈(30，80)+K∈(30，60)；
[0353] ⅠA-石英弱溶蚀＝CAL∈(250，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(40，150)+U∈(45，90)+Th∈(35，50)+K∈(50，80)；
[0354] ⅠB-石英强溶蚀＝CAL∈(30，150)+CNL∈(30，50)+RLLD/RLLS∈(8，15)+GR∈(20，150)+U∈(60，100)+Th∈(30，80)+K∈(30，50)；
[0355] ⅠB-石英中溶蚀＝CAL∈(100，240)+CNL∈(30，50)+RLLD/RLLS∈(0.2，5)+GR∈(20，110)+U∈(50，90)+Th∈(30，55)+K∈(60，70)；
[0356] ⅠB-石英弱溶蚀＝CAL∈(90，350)+CNL∈(30，45)+RLLD/RLLS∈(11，18)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)；
[0357] ⅡA-石英强溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(14，20)+GR∈(100，150)+U∈(60，90)+Th∈(30，80)+K∈(30，50)；
[0358] ⅡA-石英中溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，8)+GR∈(100，160)+U∈(0，50)+Th∈(30，80)+K∈(40，100)；
[0359] ⅡA-石英弱溶蚀＝CAL∈(250，315)+CNL∈(30，50)+RLLD/RLLS∈(15，20)+GR∈(100，150)+U∈(55，100)+Th∈(60，80)+K∈(30，80)；
[0360] ⅡB-石英强溶蚀＝CAL∈(100，220)+CNL∈(15，45)+RLLD/RLLS∈(10，20)+GR∈(20，100)+U∈(80，100)+Th∈(50，80)+K∈(30，65)；
[0361] ⅡB-石英中溶蚀＝CAL∈(0，350)+CNL∈(30，80)+RLLD/RLLS∈(10，20)+GR∈(80，120)+U∈(50，100)+Th∈(30，80)+K∈(30，70)；
[0362] ⅡB-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，50)+RLLD/RLLS∈(0.2，8)+GR∈(50，150)+U∈(80，100)+Th∈(20，70)+K∈(30，80)；
[0363] Ⅲ-石英强溶蚀＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(11，20)+GR∈(90，140)+U∈(70，100)+Th∈(30，60)+K∈(20，65)；
[0364] Ⅲ-石英中溶蚀＝CAL∈(0，350)+CNL∈(50，90)+RLLD/RLLS∈(15，20)+GR∈(80，120)+U∈(80，100)+Th∈(30，70)+K∈(25，90)；
[0365] Ⅲ-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，80)+RLLD/RLLS∈(3，8)+GR∈(110，145)+U∈(50，100)+Th∈(20，80)+K∈(10，55)；
[0366] ⅠA-长石强溶蚀＝CAL∈(50，150)+CNL∈(30，50)+GR∈(20，80)+U∈(0，80)+AC∈(100，450)+PE∈(20，50)；
[0367] ⅠA-长石中溶蚀＝CAL∈(0，80)+CNL∈(30，40)+GR∈(20，120)+U∈(60，100)+AC∈(150，250)+PE∈(40，90)；
[0368] ⅠA-长石弱溶蚀＝CAL∈(240，325)+CNL∈(60，80)+GR∈(20，125)+U∈(40，50)+AC∈(250，550)+PE∈(20，85)；
[0369] ⅠB-长石强溶蚀＝CAL∈(90，255)+CNL∈(25，50)+GR∈(20，140)+U∈(40，80)+AC∈(450，550)+PE∈(60，80)；
[0370] ⅠB-长石中溶蚀＝CAL∈(200，235)+CNL∈(30，50)+GR∈(20，130)+U∈(45，85)+AC∈(100，550)+PE∈(45，90)；
[0371] ⅠB-长石弱溶蚀＝CAL∈(200，350)+CNL∈(10，45)+GR∈(20，70)+U∈(20，70)+AC∈(900，210)+PE∈(30，70)；
[0372] ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
[0373] ⅡA-长石中溶蚀＝CAL∈(300，350)+CNL∈(15，40)+GR∈(20，115)+U∈(35，60)+AC∈(270，450)+PE∈(70，95)；
[0374] ⅡA-长石弱溶蚀＝CAL∈(180，350)+CNL∈(30，50)+GR∈(20，85)+U∈(15，50)+AC∈(200，400)+PE∈(40，80)
[0375] ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
[0376] ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
[0377] ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
[0378] Ⅲ-长石强溶蚀＝CAL∈(70，140)+CNL∈(30，55)+GR∈(20，70)+U∈(50，60)+AC∈(60，450)+PE∈(70，90)；
[0379] Ⅲ-长石中溶蚀＝CAL∈(280，350)+CNL∈(0，20)+GR∈(20，125)+U∈(70，100)+AC∈(50，180)+PE∈(50，80)；
[0380] Ⅲ-长石弱溶蚀＝CAL∈(255，350)+CNL∈(30，50)+GR∈(20，150)+U∈(40，75)+AC∈(150，550)+PE∈(55，85)
[0381] ⅠA-碳酸盐强溶蚀＝CAL∈(100，150)+U∈(80，100)+Th∈(30，55)+K∈(35，80)+AC∈(450，550)+PE∈(0，50)；
[0382] ⅠA-碳酸盐中溶蚀＝CAL∈(250，300)+U∈(70，100)+Th∈(0，38)+K∈(30，60)+AC∈(150，250)+PE∈(80，100)；
[0383] ⅠA-碳酸盐弱溶蚀＝CAL∈(280，320)+U∈(60，80)+Th∈(10，80)+K∈(40，80)+AC∈(100，300)+PE∈(45，100)；
[0384] ⅠB-碳酸盐强溶蚀＝CAL∈(250，350)+U∈(0，50)+Th∈(15，40)+K∈(65，90)+AC∈(10，100)+PE∈(0，60)；
[0385] ⅠB-碳酸盐中溶蚀＝CAL∈(240，300)+U∈(5，70)+Th∈(30，80)+K∈(0，15)+AC∈(150，240)+PE∈(20，100)；
[0386] ⅠB-碳酸盐弱溶蚀＝CAL∈(180，210)+U∈(50，90)+Th∈(30，80)+K∈(30，55)+AC∈(100，150)+PE∈(60，80)；
[0387] ⅡA-碳酸盐强溶蚀＝CAL∈(140，200)+U∈(40，100)+Th∈(60，95)+K∈(30，60)+AC∈(80，320)+PE∈(0，100)；
[0388] ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
[0389] ⅡA-碳酸盐弱溶蚀＝CAL∈(80，110)+U∈(80，100)+Th∈(40，80)+K∈(10，80)+AC∈(70，250)+PE∈(0，100)；
[0390] ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
[0391] ⅡB-碳酸盐中溶蚀＝CAL∈(200，315)+U∈(45，100)+Th∈(35，50)+K∈(30，80)+AC∈(10，80)+PE∈(10，60)；
[0392] ⅡB-碳酸盐弱溶蚀＝CAL∈(240，310)+U∈(70，90)+Th∈(30，80)+K∈(5，70)+AC∈(220，340)+PE∈(15，30)；
[0393] Ⅲ-碳酸盐强溶蚀＝CAL∈(270，300)+U∈(35，80)+Th∈(20，85)+K∈(10，80)+AC∈(150，260)+PE∈(50，70)；
[0394] Ⅲ-碳酸盐中溶蚀＝CAL∈(240，350)+U∈(40，80)+Th∈(37，80)+K∈(5，50)+AC∈(310，400)+PE∈(80，95)；
[0395] Ⅲ-碳酸盐弱溶蚀＝CAL∈(110，140)+U∈(60，100)+Th∈(30，60)+K∈(30，60)+AC∈(200，375)+PE∈(55，90)；
[0396] (2)确定地层模型参数
[0397] a.根据研究区埋藏史恢复地层L在不同演化阶段埋深：
[0398] L(Dep)＝F(X，Y，Z，t)；
[0399] b.根据研究区温度史恢复地层L在不同演化阶段温度：L(Tem)＝G(X，Y，Z，t)；
[0400] 6)模拟计算
[0401] (1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
[0402]
[0403] 其中：Φ0为总孔隙度；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
[0404] 当成岩作用为ⅠA强压实、ⅠB强压实、ⅡA强压实、ⅡB强压实、ⅠA-Mi强胶结、ⅠB-Mi强胶结、ⅡA-Mi强胶结、ⅡB-Mi强胶结、Ⅲ-Mi强胶结、ⅠA-Mi强溶蚀、ⅠB-Mi强溶蚀、ⅡA-Mi强溶蚀、ⅡB-Mi强溶蚀、Ⅲ-Mi强溶蚀时，ki＝1.6；
[0405] 当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
[0406] 当成岩作用为ⅠA弱压实、ⅠB弱压实、ⅡA弱压实、ⅡB弱压实、Ⅲ强压实、Ⅲ弱压实、ⅠA-Mi弱胶结、ⅠB-Mi弱胶结、ⅡA-Mi弱胶结、ⅡB-Mi弱胶结、Ⅲ-Mi弱胶结、ⅠA-Mi弱溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi弱溶蚀时，ki＝0.85；
[0407] (2)将时间、深度、温度等参数代入到(1)中压实模型D1、溶蚀模型D2以及胶结模型D3即可计算目标层的孔隙度模拟值。
[0408] 实施例2
[0409] 基于成岩相测井识别技术模拟计算渤中凹陷第三系湖相碎屑岩储层孔隙度的方法，
[0410] 利用研究区测井数据获取成岩相数据，再根据研究区成岩作用类型、强度、埋深史、地温史数据确定各成岩作用模型参数，最后对研究区孔隙度进行模拟计算，具体步骤如下：
[0411] 1)收集研究区地质资料
[0412] 收集的地质资料包括以下内容：
[0413] ①研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料等；
[0414] ②研究区测井数据，其中包括井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)数据；
[0415] ③研究区埋藏史、热演化史；
[0416] 2)确定研究区目标层的成岩相
[0417] (1)确定成岩相
[0418] 基于研究区测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井组合，判识渤中凹陷第三系不同深度井段成岩相。
[0419] F＝∑kili
[0420] 其中：F为通过测井分析识别的成岩相，成岩相为各种成岩作用之和，即：D-Com+D-Cem+D-Dis；压实作用D-Com包括：ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；胶结作用D-Cem包括：ⅠA-Mi强胶结、ⅠA-Mi中胶结、ⅠA-Mi弱胶结、ⅠB-Mi强胶结、ⅠB-Mi中胶结、ⅠB-Mi弱胶结、ⅡA-Mi强胶结、ⅡA-Mi中胶结、ⅡA-Mi弱胶结、ⅡB-Mi强胶结、ⅡB-Mi中胶结、ⅡB-Mi弱胶结、Ⅲ-Mi强胶结、Ⅲ-Mi中胶结、Ⅲ-Mi弱胶结；溶蚀作用D-Dis包括：ⅠA-Mi强溶蚀、ⅠA-Mi中溶蚀、ⅠA-Mi弱溶蚀、ⅠB-Mi强溶蚀、ⅠB-Mi中溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi强溶蚀、ⅡA-Mi中溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi强溶蚀、ⅡB-Mi中溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi强溶蚀、Ⅲ-Mi中溶蚀、Ⅲ-Mi弱溶蚀。li为测井曲线数据，包括CAL、CNL、DEN、RLLD/RLLS、GR、U、Th、K、AC和PE测井曲线；ki为系数，ki＝0，或者ki＝1。
[0421] 研究区主要的成岩作用包括：ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡA石英中胶结、ⅡA长石中胶结、ⅡA碳酸盐中溶蚀、ⅡA黏土弱胶结、ⅡB中压实、ⅡA长石强溶蚀、ⅡB长石中溶蚀。可采取下列模型对研究区成岩作用进行判识：
[0422] ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
[0423] ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
[0424] ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
[0425] ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
[0426] ⅡA-石英中胶结＝CAL∈(160，190)+CNL∈(40，50)+RLLD/RLLS∈(0.8，7)+GR∈(90，100)+U∈(10，70)+Th∈(65，70)+K∈(40，75)+AC∈(260，460)；
[0427] ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
[0428] ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
[0429] ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
[0430] ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
[0431] ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
[0432] ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
[0433] ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
[0434] ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
[0435] 渤中凹陷第三系井98在2000～2400段测井数据如下：
[0436]
[0437]
[0438]
[0439] a.选择井98的2008.63米深度7点0处测井数据CA7L、CNL3、DEN、3RLLD/R1LLS、G95R、U、Th、K、AC、PE，如果CAL∈(140，350)、CNL∈(35，50)、DEN∈(1.9，2.9)、GR∈(20，90)、AC∈10，340)，则该2008.63米深度点处成岩作用为ⅡA中压实；如果CAL∈(160，190)、CNL∈(40，
50)、RLLD/RLLS∈(0.8，7)、GR∈(90，100)、U∈(10，70)、Th∈(65，70)、K∈(40，75)、AC∈(260，460)，则2008.63米深度点处成岩作用为ⅡA-石英中胶结；则井98的2008.63米深度点处成岩相为ⅡA中压实_ⅡA石英中胶结。
[0440] b.选择井98的2301.13米深度点处测井数据CAL、CNL、DEN、RLLD/RLLS、GR、U、Th、K、AC、PE，如果CAL∈(140，350)、CNL∈(35，50)、DEN∈(1.9，2.9)、GR∈(20，90)、AC∈10，340)，则该2301.13米深度点处成岩作用为ⅡA中压实；如果CAL∈(275，320)、U∈(70，100)、Th∈(20，75)、K∈(30，50)、AC∈(150，550)以及PE∈(0，100)，则该2301.13米深度点处成岩作用为ⅡA-碳酸盐中溶蚀；则井98的2301.13米深度点处成岩相为ⅡA中压实_ⅡA碳酸盐中溶蚀(图1)。
[0441] (2)确定地层模型参数
[0442] 根据研究区地质资料确定地层在埋藏过程中埋藏史、热史。
[0443] 3)模拟计算
[0444] (1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
[0445]
[0446] 其中：Φ0为总孔隙度，研究区为42％；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
[0447] 当成岩作用为ⅠA强压实、ⅠB强压实、ⅡA强压实、ⅡB强压实、ⅠA-Mi强胶结、ⅠB-Mi强胶结、ⅡA-Mi强胶结、ⅡB-Mi强胶结、Ⅲ-Mi强胶结、ⅠA-Mi强溶蚀、ⅠB-Mi强溶蚀、ⅡA-Mi强溶蚀、ⅡB-Mi强溶蚀、Ⅲ-Mi强溶蚀时，ki＝1.6；
[0448] 当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
[0449] 当成岩作用为ⅠA弱压实、ⅠB弱压实、ⅡA弱压实、ⅡB弱压实、Ⅲ强压实、Ⅲ弱压实、ⅠA-Mi弱胶结、ⅠB-Mi弱胶结、ⅡA-Mi弱胶结、ⅡB-Mi弱胶结、Ⅲ-Mi弱胶结、ⅠA-Mi弱溶蚀、ⅠB-Mi弱溶蚀、ⅡA-Mi弱溶蚀、ⅡB-Mi弱溶蚀、Ⅲ-Mi弱溶蚀时，ki＝0.85；井98的2008.63米深度点处发生了Ⅱ中压实_ⅡA石英中胶结，其对应的孔隙
预测模型为：
[0450]
[0451] 即井98的2008.63米深度点处的孔隙度为8.6％；
[0452] 井98的2301.13米深度点处发生了Ⅱ中压实_ⅡA碳酸盐(方解石)中溶蚀，其对应的孔隙预测模型为：
[0453]
[0454] 其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。
