CN108362621B - 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 - Google Patents
基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 Download PDFInfo
- Publication number
- CN108362621B CN108362621B CN201810116777.2A CN201810116777A CN108362621B CN 108362621 B CN108362621 B CN 108362621B CN 201810116777 A CN201810116777 A CN 201810116777A CN 108362621 B CN108362621 B CN 108362621B
- Authority
- CN
- China
- Prior art keywords
- cal
- weak
- strong
- corrosion
- compaction
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000011435 rock Substances 0.000 title claims abstract description 27
- 238000005516 engineering process Methods 0.000 title claims abstract description 14
- 208000035126 Facies Diseases 0.000 title claims description 30
- 238000005260 corrosion Methods 0.000 claims description 306
- 230000007797 corrosion Effects 0.000 claims description 300
- 238000005056 compaction Methods 0.000 claims description 289
- 239000010453 quartz Substances 0.000 claims description 126
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 126
- 230000003628 erosive effect Effects 0.000 claims description 109
- 239000010433 feldspar Substances 0.000 claims description 68
- 239000004927 clay Substances 0.000 claims description 57
- 239000011148 porous material Substances 0.000 claims description 48
- 230000009471 action Effects 0.000 claims description 37
- 238000011160 research Methods 0.000 claims description 33
- 230000015572 biosynthetic process Effects 0.000 claims description 22
- 238000013178 mathematical model Methods 0.000 claims description 22
- 239000004568 cement Substances 0.000 claims description 20
- 238000009933 burial Methods 0.000 claims description 16
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 15
- 230000000927 lithogenic effect Effects 0.000 claims description 15
- 239000011707 mineral Substances 0.000 claims description 15
- 229910021532 Calcite Inorganic materials 0.000 claims description 14
- 230000008859 change Effects 0.000 claims description 11
- 239000010459 dolomite Substances 0.000 claims description 10
- 229910000514 dolomite Inorganic materials 0.000 claims description 10
- 238000010521 absorption reaction Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 9
- 238000004088 simulation Methods 0.000 claims description 9
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000001228 spectrum Methods 0.000 claims description 8
- 239000002245 particle Substances 0.000 claims description 7
- 238000004364 calculation method Methods 0.000 claims description 6
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 6
- 230000000694 effects Effects 0.000 claims description 5
- 239000012530 fluid Substances 0.000 claims description 5
- 238000005266 casting Methods 0.000 claims description 4
- 238000012360 testing method Methods 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 239000003292 glue Substances 0.000 claims 15
- 230000005251 gamma ray Effects 0.000 description 205
- 238000000638 solvent extraction Methods 0.000 description 8
- 238000011161 development Methods 0.000 description 6
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000035699 permeability Effects 0.000 description 2
- 239000013049 sediment Substances 0.000 description 2
- 241000239290 Araneae Species 0.000 description 1
- 229910021536 Zeolite Inorganic materials 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- HNPSIPDUKPIQMN-UHFFFAOYSA-N dioxosilane;oxo(oxoalumanyloxy)alumane Chemical compound O=[Si]=O.O=[Al]O[Al]=O HNPSIPDUKPIQMN-UHFFFAOYSA-N 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000001730 gamma-ray spectroscopy Methods 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 239000010457 zeolite Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Chemical & Material Sciences (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Dispersion Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
- Consolidation Of Soil By Introduction Of Solidifying Substances Into Soil (AREA)
Abstract
本发明公开了一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，提供了一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，该方法通过成岩相测井识别技术，对研究区成岩相进行精细研究，从而建立研究区现今成岩相剖面(Present DiageneticFacies Profile abbr.PDFP)，最后结合研究区埋藏史、热史、流体史、压力史，模拟计算井段碎屑岩孔隙度。
Description
技术领域
本发明涉及油气勘探开发技术领域，具体地指一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法。
背景技术
测井是利用岩层的电化学特征、导电性特征、放射性等地球物理特性，测量地层地球物理参数的方法。测井广泛应用于石油地质和油田勘探开发的全过程。利用测井资料，可以划分地层剖面，确定岩层厚度和埋藏深度，确定储层并识别油气水层，进行区域地层对比，而且还可以探测和研究地层主要矿物成分、孔隙度、渗透率、油气饱和度、裂缝、断层、构造特征和沉积环境与砂体的分布等，对于评价地层的储集能力、检测油气藏的开采情况，细致地分析研究油层地质特征等具有重要意义。获取地层孔隙度的测井方法通常有中子测井、孔隙度测井、声波时差测井。通过三孔隙度测井资料能连续、定量、准确获取高孔隙度或中-低孔隙度储层的孔隙值，对于深层、超深层低孔、低渗储层，三孔隙度测井计算孔隙度的准确性变差，影响了处于高成岩阶段低孔、低渗储层油气勘探开发，必须寻找出一种不同于常规方法的评价方法，来更好地、更合理地预测有效致密储层分布，最终为油气勘探开发提供有力依据。
成岩相测井是利用岩石的电性特征，在储层岩性识别的基础上对成岩作用类型和强度连续、定性、定量判定以及成岩矿物含量的定量计算。成岩相测井的目的是建立测井信息与成岩相之间关系，并对压实程度、胶结程度、溶蚀程度和孔隙类型进行判断，达到对未取心井段成岩作用、成岩相准确判别的目的。主要的成岩相测井方法有：
①利用蜘蛛网图测井识别成岩相；
②利用井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)测井数据，连续定量识别储层成岩相。
目前，通过多种测井资料数据研究了储层成岩相的判别模型，并未进一步精细研究成岩相测井以及成岩相测井数据与储层孔隙、渗透率参数之间关系。
发明内容
本发明针对现有技术存在的缺陷，提供了一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，该方法通过成岩相测井识别技术，对研究区成岩相进行精细研究，建立研究区现今成岩相剖面(Present DiageneticFacies Profile abbr.PDFP)，再结合研究区埋藏史、热史、流体史、压力史，恢复成岩相演化序列(DiageneticFaciesEvolutionary Sequence abbr.DFES)随时间发育情况，最后结合模拟算法对井段碎屑岩孔隙度进行模拟计算。
为实现上述目的，本发明提供的一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，包括以下步骤：
1)收集研究区地质资料
2)建立成岩作用数据库
在埋藏过程中，根据不同期次、不同类型的成岩作用建立成岩作用数据库；该成岩作用数据库根据不同类型成岩作用由3个小数据库组成，分为压实作用数据库、胶结作用数据库和溶蚀作用数据库：
(1)压实作用数据库D-Com
压实作用数据库根据温度和压实率建立；
(2)胶结作用数据D-Cem
胶结作用数据库根据温度和不同矿物Mi的胶结率建立；
(3)溶蚀作用数据库D-Dis
溶蚀作用数据库根据温度和不同矿物Mi溶蚀率建立；
3)建立成岩相数据库
成岩相是指在成岩环境控制下，由各种成岩作用形成并具有一定几何形态和特定成岩组构、特定成岩矿物组合的地质体。根据成岩作用，成岩相数据库可以表示为：
F＝∑fi
其中：F为成岩相；fi为成岩作用，成岩作用包括压实作用D-Com、溶蚀作用D-Dis和溶蚀作用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弱溶蚀；
4)建立成岩相孔隙演化数学模型
建立成岩作用(压实作用、石英胶结作用、长石胶结作用、碳酸盐胶结作用、黏土胶结作用、长石溶蚀作用、石英溶蚀作用、沸石溶蚀作用、碳酸盐溶蚀作用)与孔隙演化数学模型，其中各种成岩作用模型如下：
(1)压实数学模型D1
沉积物沉积后，由于上覆沉积物不断加厚，导致储层颗粒排列紧密，从而减少储层孔隙。压实作用与储层原始组构、流体排出情况、埋藏过程等紧密相关；根据压实减孔量与深度、原始结构之间的关系，可以确定不同埋深下不同碎屑岩减孔量。具体计算公式如下：
其中：颗粒密度ρs＝2650kg/m3；水密度ρw＝1.0×103kg/m3；k1为常数，k1＝0.67；c为常数，c＝3.68×10-8Pa-1；重力g＝9.8N/kg；VCom为压实作用改变孔隙体积百分量，％；
(2)溶蚀数学模型D2
a.石英溶蚀模型D2-1
其中：VQue_Dis为石英溶蚀作用改变孔隙体积百分含量，％；T为热力学温度，T＝273+T0+c*(H-H0)，T0为古地温，c为地温梯度，℃/m，H为地层埋深，H0为恒温带温度；ρH20为常温条件下水密度；
b.长石溶蚀模型D2-2
其中：t为距今时间，Ma；
为平均孔隙度，％；t1为地层温度首次达到70℃的时间，Ma；t2为地层温度首次达到90℃对应的时间，Ma；Δt＝t1-t2；VFel_Dis为长石溶蚀作用改变孔隙体积百分含量，％；
c.方解石溶蚀模型D2-3
VCal_Dis＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.442；
其中：x埋藏深度(Km)；VCal_Dis为方解石溶蚀作用改变孔隙体积百分含量，％；
d.白云石溶蚀模型D2-4
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；
其中：x埋藏深度(Km)；VDol_Dis为白云石溶蚀作用改变孔隙体积百分含量，％；
(3)胶结数学模型D3
建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
a.石英胶结模型D3-1
其中：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；
b.长石胶结模型D3-2
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；
其中：x埋藏深度(Km)；VFel_Cem为长石胶结作用改变孔隙体积百分量，％；
c.方解石胶结模型D3-3
VCal_Cem＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.44；
其中：x埋藏深度(Km)；VCal_Cem为方解石胶结作用改变孔隙体积百分量，％；
d.白云石胶结模型D3-4
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；
其中：x埋藏深度(Km)；VDol_Cem为白云石胶结作用改变孔隙体积百分量，％；
e.粘土胶结模型D3-5
VCla_Cem＝1.22932e-013*(X)^4-1.01647e-009*(X)^3+3.57019e-006*(X)^2-0.00194129*(X)+2.2677；
其中：x埋藏深度(Km)；VCla_Cem为黏土胶结作用改变孔隙体积百分量，％；
5)确定模拟地层成岩序列
(1)确定成岩相
基于测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井曲线组合，判识不同深度井段成岩相类型。
F＝∑kili
其中：F为通过测井分析识别的成岩相，成岩相为各种成岩作用之和，即：D-Com+D-Cem+D-Dis；
(2)确定地层模型参数
a.根据研究区埋藏史恢复地层L在不同演化阶段埋深：
L(Dep)＝F(X，Y，Z，t)；
b.根据研究区温度史恢复地层L在不同演化阶段温度：L(Tem)＝G(X，Y，Z，t)；6)
模拟计算
(1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
其中：Φ0为总孔隙度；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
当成岩作用为Ⅰ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；
当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
当成岩作用为Ⅰ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；
(2)将时间、深度、温度等参数代入到(1)中压实模型D1、溶蚀模型D2以及胶结模型D3即可计算目标层的孔隙度模拟值。
进一步地，所述步骤1)中，收集的地质资料包括：
1)研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料、包裹体分析测试资料等；
2)研究区测井数据，包括井径CAL、补偿中子测井CNL、密度测井(DEN)、深/浅侧向电阻率RLLD/RLLS、自然伽马GR、自然伽马能谱测井U、Th、K、声波时差AC和光电界面吸收指数PE数据；
3)研究区埋藏史、热演化史、流体史。
再进一步地，所述步骤2)中，压实作用数据库D-Com包括ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；其中，压实作用数据库划分标准如下表所示：
表1压实作用数据库划分标准
胶结作用数据库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表示矿物类型，胶结作用数据库划分标准如表下所示；
表2胶结作用数据库划分标准
溶蚀作用数据库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表示矿物类型；溶蚀作用数据库划分标准如下表所示：
表3溶蚀作用类型及划分标准
再进一步地，所述步骤5)第(1)小步中，压实作用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；
判识不同深度井段成岩相类型具体标注如下：
ⅠA强压实＝CAL∈(250，310)+CNL∈(30，40)+DEN∈(2.5，3.5)+GR∈(20，80)+AC∈(150，500)；
ⅠA中压实＝CAL∈(120，300)+CNL∈(30，35)+DEN∈(1.8，2.9)+GR∈(20，100)+AC∈(160，550)；
ⅠA弱压实＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，150)+AC∈(150，550)；
ⅠB强压实＝CAL∈(10，100)+CNL∈(30，45)+DEN∈(1.9，2.5)+GR∈(20，110)+AC∈(170，300)；
ⅠB中压实＝CAL∈(20，110)+CNL∈(30，50)+DEN∈(1.9，2.6)+GR∈(20，120)+AC∈(180，310)；
ⅠB弱压实＝CAL∈(30，120)+CNL∈(35，40)+DEN∈(1.9，2.7)+GR∈(20，130)+AC∈(190，320)；
ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
ⅡA弱压实＝CAL∈(60，150)+CNL∈(40，45)+DEN∈(1.9，3.0)+GR∈(30，150)+AC∈(220，350)；
ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
ⅡB弱压实＝CAL∈(90，180)+CNL∈(40，50)+DEN∈(1.5，2.3)+GR∈(40，70)+AC∈(250，380)；
Ⅲ强压实＝CAL∈(100，190)+CNL∈(30，45)+DEN∈(1.5，2.4)+GR∈(50，80)+AC∈(260，390)；
Ⅲ中压实＝CAL∈(110，200)+CNL∈(30，50)+DEN∈(1.5，2.5)+GR∈(60，150)+AC∈(270，400)；
Ⅲ弱压实＝CAL∈(120，210)+CNL∈(45，50)+DEN∈(1.5，2.6)+GR∈(70，150)+AC∈(280，410)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅲ-石英强胶结＝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)；
Ⅲ-石英中胶结＝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)；
Ⅲ-石英弱胶结＝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)；
ⅠA-碳酸盐强胶结＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，30)+U∈(0，100)+K∈(30，80)+PE∈(10，100)；
ⅠA-碳酸盐中胶结＝CAL∈(10，50)+CNL∈(30，35)+DEN∈(1.5，2.0)+GR∈(20，40)+U∈(10，20)+K∈(40，50)+PE∈(20，30)；
ⅠA-碳酸盐弱胶结＝CAL∈(20，60)+CNL∈(30，40)+DEN∈(1.5，2.1)+GR∈(40，50)+U∈(20，30)+K∈(50，60)+PE∈(30，40)；
ⅠB-碳酸盐强胶结＝CAL∈(30，70)+CNL∈(30，45)+DEN∈(1.5，2.2)+GR∈(30，60)+U∈(30，40)+K∈(60，70)+PE∈(40，50)；
ⅠB-碳酸盐中胶结＝CAL∈(40，80)+CNL∈(35，40)+DEN∈(1.5，2.3)+GR∈(60，70)+U∈(40，50)+K∈(70，80)+PE∈(50，100)；
ⅠB-碳酸盐弱胶结＝CAL∈(50，90)+CNL∈(35，45)+DEN∈(1.5，2.4)+GR∈(50，80)+U∈(50，60)+K∈(30，50)+PE∈(60，70)；
ⅡA-碳酸盐强胶结＝CAL∈(60，100)+CNL∈(35，50)+DEN∈(1.5，2.5)+GR∈(40，90)+U∈(60，70)+K∈(30，60)+PE∈(70，80)；
ⅡA-碳酸盐中胶结＝CAL∈(70，110)+CNL∈(40，45)+DEN∈(1.5，2.6)+GR∈(30，100)+U∈(70，80)+K∈(30，70)+PE∈(80，90)；
ⅡA-碳酸盐弱胶结＝CAL∈(80，120)+CNL∈(40，50)+DEN∈(1.5，2.7)+GR∈(90，150)+U∈(0，100)+K∈(30，80)+PE∈(90，100)；
ⅡB-碳酸盐强胶结＝CAL∈(90，130)+CNL∈(30，40)+DEN∈(1.5，2.8)+GR∈(80，110)+U∈(30，50)+K∈(30，60)+PE∈(20，50)；
ⅡB-碳酸盐中胶结＝CAL∈(100，140)+CNL∈(40，50)+DEN∈(1.5，2.9)+GR∈(70，120)+U∈(40，60)+K∈(30，70)+PE∈(30，40)；
ⅡB-碳酸盐弱胶结＝CAL∈(110，150)+CNL∈(35，50)+DEN∈(1.5，3.0)+GR∈(60，130)+U∈(50，70)+K∈(40，80)+PE∈(40，60)；
Ⅲ-碳酸盐强胶结＝CAL∈(120，160)+CNL∈(30，35)+DEN∈(1.5，3.1)+GR∈(50，140)+U∈(60，80)+K∈(50，70)+PE∈(50，70)；
Ⅲ-碳酸盐中胶结＝CAL∈(130，170)+CNL∈(40，50)+DEN∈(1.5，3.2)+GR∈40，150)+U∈(70，90)+K∈(50，80)+PE∈(60，80)；
Ⅲ-碳酸盐弱胶结＝CAL∈(140，180)+CNL∈(45，50)+DEN∈(1.5，3.3)+GR∈(30，150)+U∈(80，100)+K∈(60，80)+PE∈(70，100)；
ⅠA-黏土强胶结＝GR∈(80，150)+U∈(60，80)+Th∈(30，80)+K∈(60，90)；
ⅠA-黏土中胶结＝GR∈(80，120)+U∈(0，50)+Th∈(50，70)+K∈(30，80)；
ⅠA-黏土弱胶结＝GR∈(50，110)+U∈(70，90)+Th∈(30，80)+K∈(40，85)；
ⅠB-黏土强胶结＝GR∈(20，150)+U∈(80，100)+Th∈(30，70)+K∈(30，80)；
ⅠB-黏土中胶结＝GR∈(70，120)+U∈(0，20)+Th∈(30，80)+K∈(30，70)；
ⅠB-黏土弱胶结＝GR∈(80，110)+U∈(20，55)+Th∈(30，60)+K∈(20，70)；
ⅡA-黏土强胶结＝GR∈(30，50)+U∈(35，100)+Th∈(30，55)+K∈(30，80)；
ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
ⅡB-黏土强胶结＝GR∈(60，50)+U∈(40，60)+Th∈(30，70)+K∈(30，65)；
ⅡB-黏土中胶结＝GR∈(50，150)+U∈(0，80)+Th∈(20，80)+K∈(45，75)；
ⅡB-黏土弱胶结＝GR∈(40，50)+U∈(40，70)+Th∈(40，60)+K∈(30，50)；
Ⅲ-黏土强胶结＝GR∈(110，150)+U∈(70，100)+Th∈(30，70)+K∈(50，90)；
Ⅲ-黏土中胶结＝GR∈(60，110)+U∈(40，70)+Th∈(30，40)+K∈(30，65)；
Ⅲ-黏土弱胶结＝GR∈(20，50)+U∈(80，100)+Th∈(30，50)+K∈(30，50)；
ⅠA-石英强溶蚀＝CAL∈(280，310)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(80，150)+U∈(0，50)+Th∈(30，55)+K∈(30，70)；
ⅠA-石英中溶蚀＝CAL∈(90，150)+CNL∈(40，70)+RLLD/RLLS∈(15，20)+GR∈(20，110)+U∈(70，90)+Th∈(30，80)+K∈(30，60)；
ⅠA-石英弱溶蚀＝CAL∈(250，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(40，150)+U∈(45，90)+Th∈(35，50)+K∈(50，80)；
ⅠB-石英强溶蚀＝CAL∈(30，150)+CNL∈(30，50)+RLLD/RLLS∈(8，15)+GR∈(20，150)+U∈(60，100)+Th∈(30，80)+K∈(30，50)；
ⅠB-石英中溶蚀＝CAL∈(100，240)+CNL∈(30，50)+RLLD/RLLS∈(0.2，5)+GR∈(20，110)+U∈(50，90)+Th∈(30，55)+K∈(60，70)；
ⅠB-石英弱溶蚀＝CAL∈(90，350)+CNL∈(30，45)+RLLD/RLLS∈(11，18)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)；
ⅡA-石英强溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(14，20)+GR∈(100，150)+U∈(60，90)+Th∈(30，80)+K∈(30，50)；
ⅡA-石英中溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，8)+GR∈(100，160)+U∈(0，50)+Th∈(30，80)+K∈(40，100)；
ⅡA-石英弱溶蚀＝CAL∈(250，315)+CNL∈(30，50)+RLLD/RLLS∈(15，20)+GR∈(100，150)+U∈(55，100)+Th∈(60，80)+K∈(30，80)；
ⅡB-石英强溶蚀＝CAL∈(100，220)+CNL∈(15，45)+RLLD/RLLS∈(10，20)+GR∈(20，100)+U∈(80，100)+Th∈(50，80)+K∈(30，65)；
ⅡB-石英中溶蚀＝CAL∈(0，350)+CNL∈(30，80)+RLLD/RLLS∈(10，20)+GR∈(80，120)+U∈(50，100)+Th∈(30，80)+K∈(30，70)；
ⅡB-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，50)+RLLD/RLLS∈(0.2，8)+GR∈(50，150)+U∈(80，100)+Th∈(20，70)+K∈(30，80)；
Ⅲ-石英强溶蚀＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(11，20)+GR∈(90，140)+U∈(70，100)+Th∈(30，60)+K∈(20，65)；
Ⅲ-石英中溶蚀＝CAL∈(0，350)+CNL∈(50，90)+RLLD/RLLS∈(15，20)+GR∈(80，120)+U∈(80，100)+Th∈(30，70)+K∈(25，90)；
Ⅲ-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，80)+RLLD/RLLS∈(3，8)+GR∈(110，145)+U∈(50，100)+Th∈(20，80)+K∈(10，55)；
ⅠA-长石强溶蚀＝CAL∈(50，150)+CNL∈(30，50)+GR∈(20，80)+U∈(0，80)+AC∈(100，450)+PE∈(20，50)；
ⅠA-长石中溶蚀＝CAL∈(0，80)+CNL∈(30，40)+GR∈(20，120)+U∈(60，100)+AC∈(150，250)+PE∈(40，90)；
ⅠA-长石弱溶蚀＝CAL∈(240，325)+CNL∈(60，80)+GR∈(20，125)+U∈(40，50)+AC∈(250，550)+PE∈(20，85)；
ⅠB-长石强溶蚀＝CAL∈(90，255)+CNL∈(25，50)+GR∈(20，140)+U∈(40，80)+AC∈(450，550)+PE∈(60，80)；
ⅠB-长石中溶蚀＝CAL∈(200，235)+CNL∈(30，50)+GR∈(20，130)+U∈(45，85)+AC∈(100，550)+PE∈(45，90)；
ⅠB-长石弱溶蚀＝CAL∈(200，350)+CNL∈(10，45)+GR∈(20，70)+U∈(20，70)+AC∈(900，210)+PE∈(30，70)；
ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
ⅡA-长石中溶蚀＝CAL∈(300，350)+CNL∈(15，40)+GR∈(20，115)+U∈(35，60)+AC∈(270，450)+PE∈(70，95)；
ⅡA-长石弱溶蚀＝CAL∈(180，350)+CNL∈(30，50)+GR∈(20，85)+U∈(15，50)+AC∈(200，400)+PE∈(40，80)
ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
Ⅲ-长石强溶蚀＝CAL∈(70，140)+CNL∈(30，55)+GR∈(20，70)+U∈(50，60)+AC∈(60，450)+PE∈(70，90)；
Ⅲ-长石中溶蚀＝CAL∈(280，350)+CNL∈(0，20)+GR∈(20，125)+U∈(70，100)+AC∈(50，180)+PE∈(50，80)；
Ⅲ-长石弱溶蚀＝CAL∈(255，350)+CNL∈(30，50)+GR∈(20，150)+U∈(40，75)+AC∈(150，550)+PE∈(55，85)；
ⅠA-碳酸盐强溶蚀＝CAL∈(100，150)+U∈(80，100)+Th∈(30，55)+K∈(35，80)+AC∈(450，550)+PE∈(0，50)；
ⅠA-碳酸盐中溶蚀＝CAL∈(250，300)+U∈(70，100)+Th∈(0，38)+K∈(30，60)+AC∈(150，250)+PE∈(80，100)；
ⅠA-碳酸盐弱溶蚀＝CAL∈(280，320)+U∈(60，80)+Th∈(10，80)+K∈(40，80)+AC∈(100，300)+PE∈(45，100)；
ⅠB-碳酸盐强溶蚀＝CAL∈(250，350)+U∈(0，50)+Th∈(15，40)+K∈(65，90)+AC∈(10，100)+PE∈(0，60)；
ⅠB-碳酸盐中溶蚀＝CAL∈(240，300)+U∈(5，70)+Th∈(30，80)+K∈(0，15)+AC∈(150，240)+PE∈(20，100)；
ⅠB-碳酸盐弱溶蚀＝CAL∈(180，210)+U∈(50，90)+Th∈(30，80)+K∈(30，55)+AC∈(100，150)+PE∈(60，80)；
ⅡA-碳酸盐强溶蚀＝CAL∈(140，200)+U∈(40，100)+Th∈(60，95)+K∈(30，60)+AC∈(80，320)+PE∈(0，100)；
ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
ⅡA-碳酸盐弱溶蚀＝CAL∈(80，110)+U∈(80，100)+Th∈(40，80)+K∈(10，80)+AC∈(70，250)+PE∈(0，100)；
ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
ⅡB-碳酸盐中溶蚀＝CAL∈(200，315)+U∈(45，100)+Th∈(35，50)+K∈(30，80)+AC∈(10，80)+PE∈(10，60)；
ⅡB-碳酸盐弱溶蚀＝CAL∈(240，310)+U∈(70，90)+Th∈(30，80)+K∈(5，70)+AC∈(220，340)+PE∈(15，30)
Ⅲ-碳酸盐强溶蚀＝CAL∈(270，300)+U∈(35，80)+Th∈(20，85)+K∈(10，80)+AC∈(150，260)+PE∈(50，70)；
Ⅲ-碳酸盐中溶蚀＝CAL∈(240，350)+U∈(40，80)+Th∈(37，80)+K∈(5，50)+AC∈(310，400)+PE∈(80，95)；
Ⅲ-碳酸盐弱溶蚀＝CAL∈(110，140)+U∈(60，100)+Th∈(30，60)+K∈(30，60)+AC∈(200，375)+PE∈(55，90)。
本发明的有益效果在于：
1、预测性强，本研究基于测井成岩相的研究，对预测结果进行进一步处理，能直接对各类碎屑岩储层孔隙进行计算，尤其是低孔、低渗碎屑岩储层。
2、方法创新，本研究通过成岩相研究把测井解释和储层孔隙度预测联系在一起，有利于该方法在单井孔隙度预测上得到很好应用。
3、经济效益，该方法操作简单，避免了传统方法需要进行大量取芯和分析测试，节约了油气勘探开发成本。
附图说明
图1为井98成岩作用、成岩相的示意图；
图2为研究区埋藏史以及热史图。
具体实施方式
为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。
实施例1
基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，包括以下步骤：
1)收集研究区地质资料
(1)研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料、包裹体分析测试资料等；
(2)研究区测井数据，其中包括井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)数据；
(3)研究区埋藏史、热演化史、流体史。
2)建立成岩作用数据库
在埋藏过程中，根据不同期次、不同类型的成岩作用建立成岩作用数据库；该成岩作用数据库根据不同类型成岩作用由3个小数据库组成，分为压实作用数据库、胶结作用数据库和溶蚀作用数据库：
(1)压实作用数据库D-Com
压实作用数据库根据压实率建立，包括ⅠA强压实、ⅠA中压实、ⅠA弱压实、ⅠB强压实、ⅠB中压实、ⅠB弱压实、ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡB强压实、ⅡB中压实、ⅡB弱压实、Ⅲ强压实、Ⅲ中压实、Ⅲ弱压实；其中，压实作用数据库划分标准如表1所示。
表1压实作用数据库划分标准
(2)胶结作用数据D-Cem
胶结作用数据库根据不同矿物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所示。
表2胶结作用数据库划分标准
(3)溶蚀作用数据库D-Dis
溶蚀作用数据库根据不同矿物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所示。
表3溶蚀作用类型及划分标准
3)建立成岩相数据库
成岩相是指在成岩环境控制下，由各种成岩作用形成并具有一定几何形态和特定成岩组构、特定成岩矿物组合的地质体。根据成岩作用，成岩相数据库可以表示为：
F＝∑fi
其中：F为成岩相；fi为成岩作用，成岩作用为压实作用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弱溶蚀)。
4)建立研究区成岩相孔隙演化数学模型
建立成岩作用(压实作用、石英胶结作用、黏土胶结作用、长石溶蚀作用、方解石溶蚀作用)与孔隙演化数学模型，其中各种成岩作用模型如下：
(1)压实数学模型D1
根据压实减孔量与深度、原始结构之间的关系，确定不同埋深下不同碎屑岩减孔量，具体计算公式如下：
其中：颗粒密度ρs＝2650kg/m3；水密度ρw＝1.0×103kg/m3；k1为常数，k1＝0.67；c为常数，c＝3.68×10-8Pa-1；重力g＝9.8N/kg；VCom为压实作用改变孔隙体积百分量，％；
(2)溶蚀数学模型D2
a.长石溶蚀模型D2-2
其中：t为距今时间，Ma；
为平均孔隙度，％；t1为地层温度首次达到70℃的时间，Ma；t2为地层温度首次达到90℃对应的时间，Ma；Δt＝t1-t2；VFel_Dis为长石溶蚀作用改变孔隙体积百分含量，％；
b.方解石溶蚀模型D2-3
VCal_Dis＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.442；
其中：x埋藏深度(Km)；VCal_Dis为方解石溶蚀作用改变孔隙体积百分含量，％；
(3)胶结数学模型D3
建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
a.石英胶结模型D3-1
其中：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；
b.粘土胶结模型D3-5
VCla_Cem＝1.22932e-013*(X)^4-1.01647e-009*(X)^3+3.57019e-006*(X)^2-0.00194129*(X)+2.2677；
其中：x埋藏深度(Km)；VCla_Cem为黏土胶结作用改变孔隙体积百分量，％；
(4)主要参数选取
a.压实作用模型主要参数选取
ρs＝2650kg/m3；
ρw＝1.0×103kg/m3；
k1＝0.67；
c＝3.68×10-8Pa-1；
g＝9.8N/kg；
b.石英胶结模型主要参数选取
a＝1.98×10-22；
b＝0.022℃(-1)；
M＝60g/mol；
ρ＝0.25g/cm3；
c＝0.00000000000002℃/y；
Cn＝0.0031℃/m；
Dn＝18℃；
D＝0.02mm；
V＝1mm3；
f＝0.65；
T1＝90℃；
T2＝140℃；
c.长石溶蚀模型主要参数选取
t1＝16.2
t2＝8.2
t＝28.1
5)确定模拟地层成岩序列
(1)确定成岩相
基于测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井曲线组合，判识不同深度井段成岩相类型。
F＝∑kili
其中：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。
各成岩作用具体判识如下：
ⅠA强压实＝CAL∈(250，310)+CNL∈(30，40)+DEN∈(2.5，3.5)+GR∈(20，80)+AC∈(150，500)；
ⅠA中压实＝CAL∈(120，300)+CNL∈(30，35)+DEN∈(1.8，2.9)+GR∈(20，100)+AC∈(160，550)；
ⅠA弱压实＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，150)+AC∈(150，550)；
ⅠB强压实＝CAL∈(10，100)+CNL∈(30，45)+DEN∈(1.9，2.5)+GR∈(20，110)+AC∈(170，300)；
ⅠB中压实＝CAL∈(20，110)+CNL∈(30，50)+DEN∈(1.9，2.6)+GR∈(20，120)+AC∈(180，310)；
ⅠB弱压实＝CAL∈(30，120)+CNL∈(35，40)+DEN∈(1.9，2.7)+GR∈(20，130)+AC∈(190，320)；
ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
ⅡA弱压实＝CAL∈(60，150)+CNL∈(40，45)+DEN∈(1.9，3.0)+GR∈(30，150)+AC∈(220，350)；
ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
ⅡB弱压实＝CAL∈(90，180)+CNL∈(40，50)+DEN∈(1.5，2.3)+GR∈(40，70)+AC∈(250，380)；
Ⅲ强压实＝CAL∈(100，190)+CNL∈(30，45)+DEN∈(1.5，2.4)+GR∈(50，80)+AC∈(260，390)；
Ⅲ中压实＝CAL∈(110，200)+CNL∈(30，50)+DEN∈(1.5，2.5)+GR∈(60，150)+AC∈(270，400)；
Ⅲ弱压实＝CAL∈(120，210)+CNL∈(45，50)+DEN∈(1.5，2.6)+GR∈(70，150)+AC∈(280，410)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅲ-石英强胶结＝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)；
Ⅲ-石英中胶结＝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)；
Ⅲ-石英弱胶结＝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)；
ⅠA-碳酸盐强胶结＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，30)+U∈(0，100)+K∈(30，80)+PE∈(10，100)；
ⅠA-碳酸盐中胶结＝CAL∈(10，50)+CNL∈(30，35)+DEN∈(1.5，2.0)+GR∈(20，40)+U∈(10，20)+K∈(40，50)+PE∈(20，30)；
ⅠA-碳酸盐弱胶结＝CAL∈(20，60)+CNL∈(30，40)+DEN∈(1.5，2.1)+GR∈(40，50)+U∈(20，30)+K∈(50，60)+PE∈(30，40)；
ⅠB-碳酸盐强胶结＝CAL∈(30，70)+CNL∈(30，45)+DEN∈(1.5，2.2)+GR∈(30，60)+U∈(30，40)+K∈(60，70)+PE∈(40，50)；
ⅠB-碳酸盐中胶结＝CAL∈(40，80)+CNL∈(35，40)+DEN∈(1.5，2.3)+GR∈(60，70)+U∈(40，50)+K∈(70，80)+PE∈(50，100)；
ⅠB-碳酸盐弱胶结＝CAL∈(50，90)+CNL∈(35，45)+DEN∈(1.5，2.4)+GR∈(50，80)+U∈(50，60)+K∈(30，50)+PE∈(60，70)；
ⅡA-碳酸盐强胶结＝CAL∈(60，100)+CNL∈(35，50)+DEN∈(1.5，2.5)+GR∈(40，90)+U∈(60，70)+K∈(30，60)+PE∈(70，80)；
ⅡA-碳酸盐中胶结＝CAL∈(70，110)+CNL∈(40，45)+DEN∈(1.5，2.6)+GR∈(30，100)+U∈(70，80)+K∈(30，70)+PE∈(80，90)；
ⅡA-碳酸盐弱胶结＝CAL∈(80，120)+CNL∈(40，50)+DEN∈(1.5，2.7)+GR∈(90，150)+U∈(0，100)+K∈(30，80)+PE∈(90，100)；
ⅡB-碳酸盐强胶结＝CAL∈(90，130)+CNL∈(30，40)+DEN∈(1.5，2.8)+GR∈(80，110)+U∈(30，50)+K∈(30，60)+PE∈(20，50)；
ⅡB-碳酸盐中胶结＝CAL∈(100，140)+CNL∈(40，50)+DEN∈(1.5，2.9)+GR∈(70，120)+U∈(40，60)+K∈(30，70)+PE∈(30，40)；
ⅡB-碳酸盐弱胶结＝CAL∈(110，150)+CNL∈(35，50)+DEN∈(1.5，3.0)+GR∈(60，130)+U∈(50，70)+K∈(40，80)+PE∈(40，60)；
Ⅲ-碳酸盐强胶结＝CAL∈(120，160)+CNL∈(30，35)+DEN∈(1.5，3.1)+GR∈(50，140)+U∈(60，80)+K∈(50，70)+PE∈(50，70)；
Ⅲ-碳酸盐中胶结＝CAL∈(130，170)+CNL∈(40，50)+DEN∈(1.5，3.2)+GR∈40，150)+U∈(70，90)+K∈(50，80)+PE∈(60，80)；
Ⅲ-碳酸盐弱胶结＝CAL∈(140，180)+CNL∈(45，50)+DEN∈(1.5，3.3)+GR∈(30，150)+U∈(80，100)+K∈(60，80)+PE∈(70，100)；
ⅠA-黏土强胶结＝GR∈(80，150)+U∈(60，80)+Th∈(30，80)+K∈(60，90)；
ⅠA-黏土中胶结＝GR∈(80，120)+U∈(0，50)+Th∈(50，70)+K∈(30，80)；
ⅠA-黏土弱胶结＝GR∈(50，110)+U∈(70，90)+Th∈(30，80)+K∈(40，85)；
ⅠB-黏土强胶结＝GR∈(20，150)+U∈(80，100)+Th∈(30，70)+K∈(30，80)；
ⅠB-黏土中胶结＝GR∈(70，120)+U∈(0，20)+Th∈(30，80)+K∈(30，70)；
ⅠB-黏土弱胶结＝GR∈(80，110)+U∈(20，55)+Th∈(30，60)+K∈(20，70)；
ⅡA-黏土强胶结＝GR∈(30，50)+U∈(35，100)+Th∈(30，55)+K∈(30，80)；
ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
ⅡB-黏土强胶结＝GR∈(60，50)+U∈(40，60)+Th∈(30，70)+K∈(30，65)；
ⅡB-黏土中胶结＝GR∈(50，150)+U∈(0，80)+Th∈(20，80)+K∈(45，75)；
ⅡB-黏土弱胶结＝GR∈(40，50)+U∈(40，70)+Th∈(40，60)+K∈(30，50)；
Ⅲ-黏土强胶结＝GR∈(110，150)+U∈(70，100)+Th∈(30，70)+K∈(50，90)；
Ⅲ-黏土中胶结＝GR∈(60，110)+U∈(40，70)+Th∈(30，40)+K∈(30，65)；
Ⅲ-黏土弱胶结＝GR∈(20，50)+U∈(80，100)+Th∈(30，50)+K∈(30，50)；
ⅠA-石英强溶蚀＝CAL∈(280，310)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(80，150)+U∈(0，50)+Th∈(30，55)+K∈(30，70)；
ⅠA-石英中溶蚀＝CAL∈(90，150)+CNL∈(40，70)+RLLD/RLLS∈(15，20)+GR∈(20，110)+U∈(70，90)+Th∈(30，80)+K∈(30，60)；
ⅠA-石英弱溶蚀＝CAL∈(250，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(40，150)+U∈(45，90)+Th∈(35，50)+K∈(50，80)；
ⅠB-石英强溶蚀＝CAL∈(30，150)+CNL∈(30，50)+RLLD/RLLS∈(8，15)+GR∈(20，150)+U∈(60，100)+Th∈(30，80)+K∈(30，50)；
ⅠB-石英中溶蚀＝CAL∈(100，240)+CNL∈(30，50)+RLLD/RLLS∈(0.2，5)+GR∈(20，110)+U∈(50，90)+Th∈(30，55)+K∈(60，70)；
ⅠB-石英弱溶蚀＝CAL∈(90，350)+CNL∈(30，45)+RLLD/RLLS∈(11，18)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)；
ⅡA-石英强溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(14，20)+GR∈(100，150)+U∈(60，90)+Th∈(30，80)+K∈(30，50)；
ⅡA-石英中溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，8)+GR∈(100，160)+U∈(0，50)+Th∈(30，80)+K∈(40，100)；
ⅡA-石英弱溶蚀＝CAL∈(250，315)+CNL∈(30，50)+RLLD/RLLS∈(15，20)+GR∈(100，150)+U∈(55，100)+Th∈(60，80)+K∈(30，80)；
ⅡB-石英强溶蚀＝CAL∈(100，220)+CNL∈(15，45)+RLLD/RLLS∈(10，20)+GR∈(20，100)+U∈(80，100)+Th∈(50，80)+K∈(30，65)；
ⅡB-石英中溶蚀＝CAL∈(0，350)+CNL∈(30，80)+RLLD/RLLS∈(10，20)+GR∈(80，120)+U∈(50，100)+Th∈(30，80)+K∈(30，70)；
ⅡB-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，50)+RLLD/RLLS∈(0.2，8)+GR∈(50，150)+U∈(80，100)+Th∈(20，70)+K∈(30，80)；
Ⅲ-石英强溶蚀＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(11，20)+GR∈(90，140)+U∈(70，100)+Th∈(30，60)+K∈(20，65)；
Ⅲ-石英中溶蚀＝CAL∈(0，350)+CNL∈(50，90)+RLLD/RLLS∈(15，20)+GR∈(80，120)+U∈(80，100)+Th∈(30，70)+K∈(25，90)；
Ⅲ-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，80)+RLLD/RLLS∈(3，8)+GR∈(110，145)+U∈(50，100)+Th∈(20，80)+K∈(10，55)；
ⅠA-长石强溶蚀＝CAL∈(50，150)+CNL∈(30，50)+GR∈(20，80)+U∈(0，80)+AC∈(100，450)+PE∈(20，50)；
ⅠA-长石中溶蚀＝CAL∈(0，80)+CNL∈(30，40)+GR∈(20，120)+U∈(60，100)+AC∈(150，250)+PE∈(40，90)；
ⅠA-长石弱溶蚀＝CAL∈(240，325)+CNL∈(60，80)+GR∈(20，125)+U∈(40，50)+AC∈(250，550)+PE∈(20，85)；
ⅠB-长石强溶蚀＝CAL∈(90，255)+CNL∈(25，50)+GR∈(20，140)+U∈(40，80)+AC∈(450，550)+PE∈(60，80)；
ⅠB-长石中溶蚀＝CAL∈(200，235)+CNL∈(30，50)+GR∈(20，130)+U∈(45，85)+AC∈(100，550)+PE∈(45，90)；
ⅠB-长石弱溶蚀＝CAL∈(200，350)+CNL∈(10，45)+GR∈(20，70)+U∈(20，70)+AC∈(900，210)+PE∈(30，70)；
ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
ⅡA-长石中溶蚀＝CAL∈(300，350)+CNL∈(15，40)+GR∈(20，115)+U∈(35，60)+AC∈(270，450)+PE∈(70，95)；
ⅡA-长石弱溶蚀＝CAL∈(180，350)+CNL∈(30，50)+GR∈(20，85)+U∈(15，50)+AC∈(200，400)+PE∈(40，80)
ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
Ⅲ-长石强溶蚀＝CAL∈(70，140)+CNL∈(30，55)+GR∈(20，70)+U∈(50，60)+AC∈(60，450)+PE∈(70，90)；
Ⅲ-长石中溶蚀＝CAL∈(280，350)+CNL∈(0，20)+GR∈(20，125)+U∈(70，100)+AC∈(50，180)+PE∈(50，80)；
Ⅲ-长石弱溶蚀＝CAL∈(255，350)+CNL∈(30，50)+GR∈(20，150)+U∈(40，75)+AC∈(150，550)+PE∈(55，85)
ⅠA-碳酸盐强溶蚀＝CAL∈(100，150)+U∈(80，100)+Th∈(30，55)+K∈(35，80)+AC∈(450，550)+PE∈(0，50)；
ⅠA-碳酸盐中溶蚀＝CAL∈(250，300)+U∈(70，100)+Th∈(0，38)+K∈(30，60)+AC∈(150，250)+PE∈(80，100)；
ⅠA-碳酸盐弱溶蚀＝CAL∈(280，320)+U∈(60，80)+Th∈(10，80)+K∈(40，80)+AC∈(100，300)+PE∈(45，100)；
ⅠB-碳酸盐强溶蚀＝CAL∈(250，350)+U∈(0，50)+Th∈(15，40)+K∈(65，90)+AC∈(10，100)+PE∈(0，60)；
ⅠB-碳酸盐中溶蚀＝CAL∈(240，300)+U∈(5，70)+Th∈(30，80)+K∈(0，15)+AC∈(150，240)+PE∈(20，100)；
ⅠB-碳酸盐弱溶蚀＝CAL∈(180，210)+U∈(50，90)+Th∈(30，80)+K∈(30，55)+AC∈(100，150)+PE∈(60，80)；
ⅡA-碳酸盐强溶蚀＝CAL∈(140，200)+U∈(40，100)+Th∈(60，95)+K∈(30，60)+AC∈(80，320)+PE∈(0，100)；
ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
ⅡA-碳酸盐弱溶蚀＝CAL∈(80，110)+U∈(80，100)+Th∈(40，80)+K∈(10，80)+AC∈(70，250)+PE∈(0，100)；
ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
ⅡB-碳酸盐中溶蚀＝CAL∈(200，315)+U∈(45，100)+Th∈(35，50)+K∈(30，80)+AC∈(10，80)+PE∈(10，60)；
ⅡB-碳酸盐弱溶蚀＝CAL∈(240，310)+U∈(70，90)+Th∈(30，80)+K∈(5，70)+AC∈(220，340)+PE∈(15，30)；
Ⅲ-碳酸盐强溶蚀＝CAL∈(270，300)+U∈(35，80)+Th∈(20，85)+K∈(10，80)+AC∈(150，260)+PE∈(50，70)；
Ⅲ-碳酸盐中溶蚀＝CAL∈(240，350)+U∈(40，80)+Th∈(37，80)+K∈(5，50)+AC∈(310，400)+PE∈(80，95)；
Ⅲ-碳酸盐弱溶蚀＝CAL∈(110，140)+U∈(60，100)+Th∈(30，60)+K∈(30，60)+AC∈(200，375)+PE∈(55，90)；
(2)确定地层模型参数
a.根据研究区埋藏史恢复地层L在不同演化阶段埋深：
L(Dep)＝F(X，Y，Z，t)；
b.根据研究区温度史恢复地层L在不同演化阶段温度：L(Tem)＝G(X，Y，Z，t)；
6)模拟计算
(1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
其中：Φ0为总孔隙度；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
当成岩作用为Ⅰ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；
当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
当成岩作用为Ⅰ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；
(2)将时间、深度、温度等参数代入到(1)中压实模型D1、溶蚀模型D2以及胶结模型D3即可计算目标层的孔隙度模拟值。
实施例2
基于成岩相测井识别技术模拟计算渤中凹陷第三系湖相碎屑岩储层孔隙度的方法，
利用研究区测井数据获取成岩相数据，再根据研究区成岩作用类型、强度、埋深史、地温史数据确定各成岩作用模型参数，最后对研究区孔隙度进行模拟计算，具体步骤如下：
1)收集研究区地质资料
收集的地质资料包括以下内容：
①研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料等；
②研究区测井数据，其中包括井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)数据；
③研究区埋藏史、热演化史；
2)确定研究区目标层的成岩相
(1)确定成岩相
基于研究区测井数据(井径(CAL)、补偿中子测井(CNL)、密度测井(DEN)、深/浅侧向电阻率(RLLD/RLLS)、自然伽马(GR)、自然伽马能谱测井(U、Th、K)、声波时差(AC)和光电界面吸收指数(PE)，采用不同测井组合，判识渤中凹陷第三系不同深度井段成岩相。
F＝∑kili
其中：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。
研究区主要的成岩作用包括：ⅡA强压实、ⅡA中压实、ⅡA弱压实、ⅡA石英中胶结、ⅡA长石中胶结、ⅡA碳酸盐中溶蚀、ⅡA黏土弱胶结、ⅡB中压实、ⅡA长石强溶蚀、ⅡB长石中溶蚀。可采取下列模型对研究区成岩作用进行判识：
ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
Ⅱ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)；
ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
渤中凹陷第三系井98在2000～2400段测井数据如下：
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石英中胶结。
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)。
(2)确定地层模型参数
根据研究区地质资料确定地层在埋藏过程中埋藏史、热史。
3)模拟计算
(1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
其中：Φ0为总孔隙度，研究区为42％；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
当成岩作用为Ⅰ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；
当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀、Ⅲ-Mi中溶蚀时，ki＝1；
当成岩作用为Ⅰ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石英中胶结，其对应的孔隙预测模型为：
即井98的2008.63米深度点处的孔隙度为8.6％；
井98的2301.13米深度点处发生了Ⅱ中压实_ⅡA碳酸盐(方解石)中溶蚀，其对应的孔隙预测模型为：
其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。
Claims (4)
1.一种基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，其特征在于：包括以下步骤：
1)收集研究区地质资料
2)建立成岩作用数据库
在埋藏过程中，根据不同期次、不同类型的成岩作用建立成岩作用数据库；该成岩作用数据库根据不同类型成岩作用由3个小数据库组成，分为压实作用数据库D-Com、胶结作用数据库D-Cem和溶蚀作用数据库D-Dis；
(1)压实作用数据库D-Com
压实作用数据库根据温度和压实率建立；
(2)胶结作用数据D-Cem
胶结作用数据库根据温度和不同矿物Mi的胶结率建立；
(3)溶蚀作用数据库D-Dis
溶蚀作用数据库根据温度和不同矿物Mi溶蚀率建立；
3)建立成岩相数据库
根据成岩作用，成岩相数据库表示为：
F＝∑fi
其中：F为成岩相；fi为成岩作用，成岩作用包括压实作用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弱溶蚀；
4)建立成岩相孔隙演化数学模型
建立成岩作用与孔隙演化数学模型，其中，成岩作用包括压实作用、石英胶结作用、长石胶结作用、方解石胶结作用、白云石胶结作用、黏土胶结作用、长石溶蚀作用、石英溶蚀作用、方解石溶蚀作用和白云石溶蚀作用；各种成岩作用模型如下：
(1)压实数学模型D1
根据压实减孔量与深度、原始结构之间的关系，确定不同埋深下不同碎屑岩减孔量，具体计算公式如下：
其中：颗粒密度ρs＝2650kg/m3；水密度ρw＝1.0×103kg/m3；k1为常数，k1＝0.67；c为常数，c＝3.68×10-8Pa-1；重力g＝9.8N/kg；VCom为压实作用改变孔隙体积百分量，％；
(2)溶蚀数学模型D2
a.石英溶蚀模型D2-1
其中：VQue_Dis为石英溶蚀作用改变孔隙体积百分含量，％；T为热力学温度，T＝273+T0+c*(H-H0)，T0为古地温，c为地温梯度，℃/m，H为地层埋深，H0为恒温带温度；ρH20为常温条件下水密度；
b.长石溶蚀模型D2-2
其中：t为距今时间，Ma；
为平均孔隙度，％；t1为地层温度首次达到70℃的时间，Ma；t2为地层温度首次达到90℃对应的时间，Ma；Δt＝t1-t2；VFel_Dis为长石溶蚀作用改变孔隙体积百分含量，％；
c.方解石溶蚀模型D2-3
VCal_Dis＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.442；
其中：x埋藏深度(Km)；VCal_Dis为方解石溶蚀作用改变孔隙体积百分含量，％；
d.白云石溶蚀模型D2-4
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；
其中：x埋藏深度，Km；VDol_Dis为白云石溶蚀作用改变孔隙体积百分含量，％；
(3)胶结数学模型D3
建立起不同类型胶结物的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
a.石英胶结模型D3-1
其中：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；
b.长石胶结模型D3-2
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；
其中：x埋藏深度，Km；VFel_Cem为长石胶结作用改变孔隙体积百分量，％；
c.方解石胶结模型D3-3
VCal_Cem＝-6.169e-14*(X)^4+6.266e-10*(X)^3-2.198e-06*(X)^2+0.005019*(X)+1.44；
其中：x埋藏深度(Km)；VCal_Cem为方解石胶结作用改变孔隙体积百分量，％；
d.白云石胶结模型D3-4
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；
其中：x埋藏深度，Km；VDol_Cem为白云石胶结作用改变孔隙体积百分量，％；
e.粘土胶结模型D3-5
VCla_Cem＝1.22932e-013*(X)^4-1.01647e-009*(X)^3+3.57019e-006*(X)^2-0.00194129*(X)+2.2677；
其中：x埋藏深度，Km；VCla_Cem为黏土胶结作用改变孔隙体积百分量，％；
5)确定模拟地层成岩序列
(1)确定成岩相
基于测井数据，它包括井径CAL、补偿中子测井CNL、密度测井DEN、深/浅侧向电阻率RLLD/RLLS、自然伽马GR、自然伽马能谱测井U、Th、K、声波时差AC和光电界面吸收指数PE，采用不同测井曲线组合，判识不同深度井段成岩相类型：
F＝∑kili
其中：F为通过测井分析识别的成岩相，成岩相为各种成岩作用之和，即：D-Com+D-Cem+D-Dis；li为测井曲线数据，ki为系数；
(2)确定地层模型参数
a.根据研究区埋藏史恢复地层L在不同演化阶段埋深：
L(Dep)＝F(X，Y，Z，t)；
b.根据研究区温度史恢复地层L在不同演化阶段温度：
L(Tem)＝G(X，Y，Z，t)；
6)模拟计算
(1)根据目标层的成岩相数据，确定模拟点孔隙演化数学模型：
其中：Φ0为总孔隙度；Φ为现今孔隙值；Φi为与成岩相对应的孔隙度影响因素，0≤i≤2；i＝0，为压实影响因素；i＝1，为胶结影响因素；i＝2，为溶蚀影响因素；fi为孔隙度影响因子，fi＝0或fi＝1；ki为孔隙度影响系数，
当成岩作用为Ⅰ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；
当成岩作用为ⅠA中压实、ⅠB中压实、ⅡA中压实、ⅡB中压实、Ⅲ中压实、ⅠA-Mi中胶结、ⅠB-Mi中胶结、ⅡA-Mi中胶结、ⅡB-Mi中胶结、Ⅲ-Mi中胶结、ⅠA-Mi中溶蚀、ⅠB-Mi中溶蚀、ⅡA-Mi中溶蚀、ⅡB-Mi中溶蚀和Ⅲ-Mi中溶蚀时，ki＝1；
当成岩作用为Ⅰ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；
(2)将时间、深度和温度参数代入到(1)中压实模型D1、溶蚀模型D2以及胶结模型D3中，计算目标层的孔隙度模拟值。
2.根据权利要求1所述基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，其特征在于：所述步骤1)中，收集的地质资料包括：
1)研究区地质资料，其中包括岩心、录井信息、薄片鉴定报告、铸体薄片资料、包裹体分析测试资料；
2)研究区测井数据，包括井径CAL、补偿中子测井CNL、密度测井(DEN)、深/浅侧向电阻率RLLD/RLLS、自然伽马GR、自然伽马能谱测井U、Th、K、声波时差AC和光电界面吸收指数PE数据；
3)研究区埋藏史、热演化史、流体史。
3.根据权利要求1所述基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，其特征在于：所述步骤2)中，压实作用数据库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弱溶蚀，其中Mi表示矿物类型。
4.根据权利要求1所述基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法，其特征在于：所述步骤5)第(1)小步中，
压实作用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；
判识不同深度井段成岩相类型具体标注如下：
ⅠA强压实＝CAL∈(250，310)+CNL∈(30，40)+DEN∈(2.5，3.5)+GR∈(20，80)+AC∈(150，500)；
ⅠA中压实＝CAL∈(120，300)+CNL∈(30，35)+DEN∈(1.8，2.9)+GR∈(20，100)+AC∈(160，550)；
ⅠA弱压实＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，150)+AC∈(150，550)；
ⅠB强压实＝CAL∈(10，100)+CNL∈(30，45)+DEN∈(1.9，2.5)+GR∈(20，110)+AC∈(170，300)；
ⅠB中压实＝CAL∈(20，110)+CNL∈(30，50)+DEN∈(1.9，2.6)+GR∈(20，120)+AC∈(180，310)；
ⅠB弱压实＝CAL∈(30，120)+CNL∈(35，40)+DEN∈(1.9，2.7)+GR∈(20，130)+AC∈(190，320)；
ⅡA强压实＝CAL∈(50，130)+CNL∈(35，45)+DEN∈(1.9，2.8)+GR∈(20，140)+AC∈(200，330)；
ⅡA中压实＝CAL∈(140，350)+CNL∈(35，50)+DEN∈(1.9，2.9)+GR∈(20，90)+AC∈(210，340)；
ⅡA弱压实＝CAL∈(60，150)+CNL∈(40，45)+DEN∈(1.9，3.0)+GR∈(30，150)+AC∈(220，350)；
ⅡB强压实＝CAL∈(70，160)+CNL∈(40，50)+DEN∈(1.9，3.1)+GR∈(30，50)+AC∈(230，360)；
ⅡB中压实＝CAL∈(80，170)+CNL∈(45，50)+DEN∈(2.0，2.2)+GR∈(30，60)+AC∈(240，370)；
ⅡB弱压实＝CAL∈(90，180)+CNL∈(40，50)+DEN∈(1.5，2.3)+GR∈(40，70)+AC∈(250，380)；
Ⅲ强压实＝CAL∈(100，190)+CNL∈(30，45)+DEN∈(1.5，2.4)+GR∈(50，80)+AC∈(260，390)；
Ⅲ中压实＝CAL∈(110，200)+CNL∈(30，50)+DEN∈(1.5，2.5)+GR∈(60，150)+AC∈(270，400)；
Ⅲ弱压实＝CAL∈(120，210)+CNL∈(45，50)+DEN∈(1.5，2.6)+GR∈(70，150)+AC∈(280，410)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅰ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅱ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)；
Ⅲ-石英强胶结＝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)；
Ⅲ-石英中胶结＝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)；
Ⅲ-石英弱胶结＝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)；
ⅠA-碳酸盐强胶结＝CAL∈(0，350)+CNL∈(30，50)+DEN∈(1.5，3.5)+GR∈(20，30)+U∈(0，100)+K∈(30，80)+PE∈(10，100)；
ⅠA-碳酸盐中胶结＝CAL∈(10，50)+CNL∈(30，35)+DEN∈(1.5，2.0)+GR∈(20，40)+U∈(10，20)+K∈(40，50)+PE∈(20，30)；
ⅠA-碳酸盐弱胶结＝CAL∈(20，60)+CNL∈(30，40)+DEN∈(1.5，2.1)+GR∈(40，50)+U∈(20，30)+K∈(50，60)+PE∈(30，40)；
ⅠB-碳酸盐强胶结＝CAL∈(30，70)+CNL∈(30，45)+DEN∈(1.5，2.2)+GR∈(30，60)+U∈(30，40)+K∈(60，70)+PE∈(40，50)；
ⅠB-碳酸盐中胶结＝CAL∈(40，80)+CNL∈(35，40)+DEN∈(1.5，2.3)+GR∈(60，70)+U∈(40，50)+K∈(70，80)+PE∈(50，100)；
ⅠB-碳酸盐弱胶结＝CAL∈(50，90)+CNL∈(35，45)+DEN∈(1.5，2.4)+GR∈(50，80)+U∈(50，60)+K∈(30，50)+PE∈(60，70)；
ⅡA-碳酸盐强胶结＝CAL∈(60，100)+CNL∈(35，50)+DEN∈(1.5，2.5)+GR∈(40，90)+U∈(60，70)+K∈(30，60)+PE∈(70，80)；
ⅡA-碳酸盐中胶结＝CAL∈(70，110)+CNL∈(40，45)+DEN∈(1.5，2.6)+GR∈(30，100)+U∈(70，80)+K∈(30，70)+PE∈(80，90)；
ⅡA-碳酸盐弱胶结＝CAL∈(80，120)+CNL∈(40，50)+DEN∈(1.5，2.7)+GR∈(90，150)+U∈(0，100)+K∈(30，80)+PE∈(90，100)；
ⅡB-碳酸盐强胶结＝CAL∈(90，130)+CNL∈(30，40)+DEN∈(1.5，2.8)+GR∈(80，110)+U∈(30，50)+K∈(30，60)+PE∈(20，50)；
ⅡB-碳酸盐中胶结＝CAL∈(100，140)+CNL∈(40，50)+DEN∈(1.5，2.9)+GR∈(70，120)+U∈(40，60)+K∈(30，70)+PE∈(30，40)；
ⅡB-碳酸盐弱胶结＝CAL∈(110，150)+CNL∈(35，50)+DEN∈(1.5，3.0)+GR∈(60，130)+U∈(50，70)+K∈(40，80)+PE∈(40，60)；
Ⅲ-碳酸盐强胶结＝CAL∈(120，160)+CNL∈(30，35)+DEN∈(1.5，3.1)+GR∈(50，140)+U∈(60，80)+K∈(50，70)+PE∈(50，70)；
Ⅲ-碳酸盐中胶结＝CAL∈(130，170)+CNL∈(40，50)+DEN∈(1.5，3.2)+GR∈40，150)+U∈(70，90)+K∈(50，80)+PE∈(60，80)；
Ⅲ-碳酸盐弱胶结＝CAL∈(140，180)+CNL∈(45，50)+DEN∈(1.5，3.3)+GR∈(30，150)+U∈(80，100)+K∈(60，80)+PE∈(70，100)；
ⅠA-黏土强胶结＝GR∈(80，150)+U∈(60，80)+Th∈(30，80)+K∈(60，90)；
ⅠA-黏土中胶结＝GR∈(80，120)+U∈(0，50)+Th∈(50，70)+K∈(30，80)；
ⅠA-黏土弱胶结＝GR∈(50，110)+U∈(70，90)+Th∈(30，80)+K∈(40，85)；
ⅠB-黏土强胶结＝GR∈(20，150)+U∈(80，100)+Th∈(30，70)+K∈(30，80)；
ⅠB-黏土中胶结＝GR∈(70，120)+U∈(0，20)+Th∈(30，80)+K∈(30，70)；
ⅠB-黏土弱胶结＝GR∈(80，110)+U∈(20，55)+Th∈(30，60)+K∈(20，70)；
ⅡA-黏土强胶结＝GR∈(30，50)+U∈(35，100)+Th∈(30，55)+K∈(30，80)；
ⅡA-黏土中胶结＝GR∈(20，150)+U∈(0，100)+Th∈(20，45)+K∈(50，60)；
ⅡA-黏土弱胶结＝GR∈(70，110)+U∈(55，90)+Th∈(30，80)+K∈(30，80)；
ⅡB-黏土强胶结＝GR∈(60，50)+U∈(40，60)+Th∈(30，70)+K∈(30，65)；
ⅡB-黏土中胶结＝GR∈(50，150)+U∈(0，80)+Th∈(20，80)+K∈(45，75)；
ⅡB-黏土弱胶结＝GR∈(40，50)+U∈(40，70)+Th∈(40，60)+K∈(30，50)；
Ⅲ-黏土强胶结＝GR∈(110，150)+U∈(70，100)+Th∈(30，70)+K∈(50，90)；
Ⅲ-黏土中胶结＝GR∈(60，110)+U∈(40，70)+Th∈(30，40)+K∈(30，65)；
Ⅲ-黏土弱胶结＝GR∈(20，50)+U∈(80，100)+Th∈(30，50)+K∈(30，50)；
ⅠA-石英强溶蚀＝CAL∈(280，310)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(80，150)+U∈(0，50)+Th∈(30，55)+K∈(30，70)；
ⅠA-石英中溶蚀＝CAL∈(90，150)+CNL∈(40，70)+RLLD/RLLS∈(15，20)+GR∈(20，110)+U∈(70，90)+Th∈(30，80)+K∈(30，60)；
ⅠA-石英弱溶蚀＝CAL∈(250，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，10)+GR∈(40，150)+U∈(45，90)+Th∈(35，50)+K∈(50，80)；
ⅠB-石英强溶蚀＝CAL∈(30，150)+CNL∈(30，50)+RLLD/RLLS∈(8，15)+GR∈(20，150)+U∈(60，100)+Th∈(30，80)+K∈(30，50)；
ⅠB-石英中溶蚀＝CAL∈(100，240)+CNL∈(30，50)+RLLD/RLLS∈(0.2，5)+GR∈(20，110)+U∈(50，90)+Th∈(30，55)+K∈(60，70)；
ⅠB-石英弱溶蚀＝CAL∈(90，350)+CNL∈(30，45)+RLLD/RLLS∈(11，18)+GR∈(20，150)+U∈(0，100)+Th∈(30，80)+K∈(30，80)；
ⅡA-石英强溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(14，20)+GR∈(100，150)+U∈(60，90)+Th∈(30，80)+K∈(30，50)；
ⅡA-石英中溶蚀＝CAL∈(300，350)+CNL∈(30，50)+RLLD/RLLS∈(0.2，8)+GR∈(100，160)+U∈(0，50)+Th∈(30，80)+K∈(40，100)；
ⅡA-石英弱溶蚀＝CAL∈(250，315)+CNL∈(30，50)+RLLD/RLLS∈(15，20)+GR∈(100，150)+U∈(55，100)+Th∈(60，80)+K∈(30，80)；
ⅡB-石英强溶蚀＝CAL∈(100，220)+CNL∈(15，45)+RLLD/RLLS∈(10，20)+GR∈(20，100)+U∈(80，100)+Th∈(50，80)+K∈(30，65)；
ⅡB-石英中溶蚀＝CAL∈(0，350)+CNL∈(30，80)+RLLD/RLLS∈(10，20)+GR∈(80，120)+U∈(50，100)+Th∈(30，80)+K∈(30，70)；
ⅡB-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，50)+RLLD/RLLS∈(0.2，8)+GR∈(50，150)+U∈(80，100)+Th∈(20，70)+K∈(30，80)；
Ⅲ-石英强溶蚀＝CAL∈(0，350)+CNL∈(30，50)+RLLD/RLLS∈(11，20)+GR∈(90，140)+U∈(70，100)+Th∈(30，60)+K∈(20，65)；
Ⅲ-石英中溶蚀＝CAL∈(0，350)+CNL∈(50，90)+RLLD/RLLS∈(15，20)+GR∈(80，120)+U∈(80，100)+Th∈(30，70)+K∈(25，90)；
Ⅲ-石英弱溶蚀＝CAL∈(0，350)+CNL∈(20，80)+RLLD/RLLS∈(3，8)+GR∈(110，145)+U∈(50，100)+Th∈(20，80)+K∈(10，55)；
ⅠA-长石强溶蚀＝CAL∈(50，150)+CNL∈(30，50)+GR∈(20，80)+U∈(0，80)+AC∈(100，450)+PE∈(20，50)；
ⅠA-长石中溶蚀＝CAL∈(0，80)+CNL∈(30，40)+GR∈(20，120)+U∈(60，100)+AC∈(150，250)+PE∈(40，90)；
ⅠA-长石弱溶蚀＝CAL∈(240，325)+CNL∈(60，80)+GR∈(20，125)+U∈(40，50)+AC∈(250，550)+PE∈(20，85)；
ⅠB-长石强溶蚀＝CAL∈(90，255)+CNL∈(25，50)+GR∈(20，140)+U∈(40，80)+AC∈(450，550)+PE∈(60，80)；
ⅠB-长石中溶蚀＝CAL∈(200，235)+CNL∈(30，50)+GR∈(20，130)+U∈(45，85)+AC∈(100，550)+PE∈(45，90)；
ⅠB-长石弱溶蚀＝CAL∈(200，350)+CNL∈(10，45)+GR∈(20，70)+U∈(20，70)+AC∈(900，210)+PE∈(30，70)；
ⅡA-长石强溶蚀＝CAL∈(250，350)+CNL∈(20，55)+GR∈(20，90)+U∈(0，100)+AC∈(150，380)+PE∈(80，100)；
ⅡA-长石中溶蚀＝CAL∈(300，350)+CNL∈(15，40)+GR∈(20，115)+U∈(35，60)+AC∈(270，450)+PE∈(70，95)
ⅡA-长石弱溶蚀＝CAL∈(180，350)+CNL∈(30，50)+GR∈(20，85)+U∈(15，50)+AC∈(200，400)+PE∈(40，80)；
ⅡB-长石强溶蚀＝CAL∈(100，150)+CNL∈(15，55)+GR∈(20，130)+U∈(45，95)+AC∈(90，180)+PE∈(0，50)；
ⅡB-长石中溶蚀＝CAL∈(80，250)+CNL∈(20，75)+GR∈(20，90)+U∈(5，50)+AC∈(200，250)+PE∈(0，40)；
ⅡB-长石弱溶蚀＝CAL∈(90，150)+CNL∈(30，50)+GR∈(20，85)+U∈(0，20)+AC∈(240，380)+PE∈(30，60)；
Ⅲ-长石强溶蚀＝CAL∈(70，140)+CNL∈(30，55)+GR∈(20，70)+U∈(50，60)+AC∈(60，450)+PE∈(70，90)；
Ⅲ-长石中溶蚀＝CAL∈(280，350)+CNL∈(0，20)+GR∈(20，125)+U∈(70，100)+AC∈(50，180)+PE∈(50，80)；
Ⅲ-长石弱溶蚀＝CAL∈(255，350)+CNL∈(30，50)+GR∈(20，150)+U∈(40，75)+AC∈(150，550)+PE∈(55，85)；
ⅠA-碳酸盐强溶蚀＝CAL∈(100，150)+U∈(80，100)+Th∈(30，55)+K∈(35，80)+AC∈(450，550)+PE∈(0，50)；
ⅠA-碳酸盐中溶蚀＝CAL∈(250，300)+U∈(70，100)+Th∈(0，38)+K∈(30，60)+AC∈(150，250)+PE∈(80，100)；
ⅠA-碳酸盐弱溶蚀＝CAL∈(280，320)+U∈(60，80)+Th∈(10，80)+K∈(40，80)+AC∈(100，300)+PE∈(45，100)；
ⅠB-碳酸盐强溶蚀＝CAL∈(250，350)+U∈(0，50)+Th∈(15，40)+K∈(65，90)+AC∈(10，100)+PE∈(0，60)；
ⅠB-碳酸盐中溶蚀＝CAL∈(240，300)+U∈(5，70)+Th∈(30，80)+K∈(0，15)+AC∈(150，240)+PE∈(20，100)；
ⅠB-碳酸盐弱溶蚀＝CAL∈(180，210)+U∈(50，90)+Th∈(30，80)+K∈(30，55)+AC∈(100，150)+PE∈(60，80)；
ⅡA-碳酸盐强溶蚀＝CAL∈(140，200)+U∈(40，100)+Th∈(60，95)+K∈(30，60)+AC∈(80，320)+PE∈(0，100)；
ⅡA-碳酸盐中溶蚀＝CAL∈(275，320)+U∈(70，100)+Th∈(20，75)+K∈(30，50)+AC∈(150，550)+PE∈(0，100)；
ⅡA-碳酸盐弱溶蚀＝CAL∈(80，110)+U∈(80，100)+Th∈(40，80)+K∈(10，80)+AC∈(70，250)+PE∈(0，100)；
ⅡB-碳酸盐强溶蚀＝CAL∈(70，150)+U∈(30，60)+Th∈(0，60)+K∈(25，60)+AC∈(150，250)+PE∈(5，15)；
ⅡB-碳酸盐中溶蚀＝CAL∈(200，315)+U∈(45，100)+Th∈(35，50)+K∈(30，80)+AC∈(10，80)+PE∈(10，60)；
ⅡB-碳酸盐弱溶蚀＝CAL∈(240，310)+U∈(70，90)+Th∈(30，80)+K∈(5，70)+AC∈(220，340)+PE∈(15，30)；
Ⅲ-碳酸盐强溶蚀＝CAL∈(270，300)+U∈(35，80)+Th∈(20，85)+K∈(10，80)+AC∈(150，260)+PE∈(50，70)；
Ⅲ-碳酸盐中溶蚀＝CAL∈(240，350)+U∈(40，80)+Th∈(37，80)+K∈(5，50)+AC∈(310，400)+PE∈(80，95)；
Ⅲ-碳酸盐弱溶蚀＝CAL∈(110，140)+U∈(60，100)+Th∈(30，60)+K∈(30，60)+AC∈(200，375)+PE∈(55，90)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810116777.2A CN108362621B (zh) | 2018-02-06 | 2018-02-06 | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810116777.2A CN108362621B (zh) | 2018-02-06 | 2018-02-06 | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108362621A CN108362621A (zh) | 2018-08-03 |
| CN108362621B true CN108362621B (zh) | 2020-09-04 |
Family
ID=63004583
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810116777.2A Active CN108362621B (zh) | 2018-02-06 | 2018-02-06 | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108362621B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108828666B (zh) * | 2018-08-07 | 2019-12-27 | 中国石油大学(华东) | 一种正断层附近地层埋藏史的恢复方法 |
| US11566518B2 (en) | 2020-03-02 | 2023-01-31 | Saudi Arabian Oil Company | Method and system for diagenesis-based rock classification |
| CN115577814B (zh) * | 2022-12-08 | 2023-05-12 | 珠江水利委员会珠江水利科学研究院 | 一种溶蚀检测及修复方法及*** |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN102748016A (zh) * | 2012-07-10 | 2012-10-24 | 中国石油大学(华东) | 地质历史时期砂岩储层孔隙度演化恢复方法 |
| CN102778421A (zh) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | 地质历史时期砂岩储层渗透率演化恢复方法 |
| CN104181603A (zh) * | 2014-07-24 | 2014-12-03 | 中国石油大学(华东) | 碎屑岩沉积成岩综合相识别方法 |
| CN104331537A (zh) * | 2014-09-28 | 2015-02-04 | 长江大学 | 基于储层静态因子的井位优化设计方法 |
| CN106597548A (zh) * | 2016-12-02 | 2017-04-26 | 中国石油大学(华东) | 一种多因素定量评价地质时期3d孔隙度的方法 |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN202305523U (zh) * | 2011-09-14 | 2012-07-04 | 中国石油天然气股份有限公司 | 一种成岩作用模拟实验装置 |
| CN103278436B (zh) * | 2013-02-01 | 2015-02-04 | 西安石油大学 | 特低渗透双重介质砂岩油藏微观孔隙结构的定量表征方法 |
| US10197489B2 (en) * | 2015-07-28 | 2019-02-05 | Chevron U.S.A. Inc. | Processes and systems for characterizing and optimizing fracturing fluids |
| CN105370272A (zh) * | 2015-12-07 | 2016-03-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层的测井评价方法 |
| CN106202749B (zh) * | 2016-07-15 | 2019-03-26 | 长江大学 | 基于基准面对沉积过程控制的沉积模拟方法 |
| CN107329186B (zh) * | 2017-05-17 | 2019-08-06 | 中国石油天然气股份有限公司 | 一种确定火成岩储层次生溶蚀孔隙发育度的方法及装置 |
| CN107389528A (zh) * | 2017-07-26 | 2017-11-24 | 中国海洋石油总公司 | 一种针对超压对储层孔隙演化影响定量评价的实验方法 |
-
2018
- 2018-02-06 CN CN201810116777.2A patent/CN108362621B/zh active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562048A (zh) * | 2010-12-30 | 2012-07-11 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN102748016A (zh) * | 2012-07-10 | 2012-10-24 | 中国石油大学(华东) | 地质历史时期砂岩储层孔隙度演化恢复方法 |
| CN102778421A (zh) * | 2012-07-10 | 2012-11-14 | 中国石油大学(华东) | 地质历史时期砂岩储层渗透率演化恢复方法 |
| CN104181603A (zh) * | 2014-07-24 | 2014-12-03 | 中国石油大学(华东) | 碎屑岩沉积成岩综合相识别方法 |
| CN104331537A (zh) * | 2014-09-28 | 2015-02-04 | 长江大学 | 基于储层静态因子的井位优化设计方法 |
| CN106597548A (zh) * | 2016-12-02 | 2017-04-26 | 中国石油大学(华东) | 一种多因素定量评价地质时期3d孔隙度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108362621A (zh) | 2018-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109653725B (zh) | 基于沉积微相和岩石相的混积储层水淹程度测井解释方法 | |
| CN110318745B (zh) | 一种沉积微相约束下粒径岩性测井评价方法 | |
| CN106803021B (zh) | 一种常规与非常规储层的油气资源量的评价方法 | |
| CN106154351A (zh) | 一种低孔渗储层渗透率的估算方法 | |
| CN108362621B (zh) | 基于成岩相测井识别技术模拟计算碎屑岩储层孔隙度的方法 | |
| CN104747183A (zh) | 一种碳酸盐岩储层综合分类方法 | |
| US20180355702A1 (en) | Method of exploiting a sedimentary basin comprising hydrocarbons, by modelling the accumulation of terrestrial organic matter | |
| CN110442951B (zh) | 一种考虑密度因素的烃源岩总有机碳含量预测方法 | |
| CN109324345A (zh) | 叠合盆地油气成藏期岩石孔隙度恢复方法 | |
| CN110728402A (zh) | 一种基于地质标本分析的区域矿产资源预测*** | |
| CN112177605B (zh) | 一种海域天然气水合物有利聚集区主控因素的确定方法 | |
| CN110529106B (zh) | 一种利用测井资料确定煤层显微组分含量的方法 | |
| CN113970796A (zh) | 精确恢复沉积盆地古水深的方法 | |
| CN108345046B (zh) | 碎屑岩储层埋藏过程中成岩相划分方法 | |
| CN112185469B (zh) | 一种预测海域天然气水合物有利聚集区的方法 | |
| CN112832738B (zh) | 一种碎屑岩累计生烃强度确定方法及甜点层的识别与评价方法 | |
| CN110688781B (zh) | 一种低渗非均质气藏储层测井解释方法 | |
| Xiong et al. | Petrophysical evaluation of gas hydrate in Shenhu area, China | |
| Liu et al. | Overpressure and gas charging in tight sandstone: Xujiahe Formation, northeastern Sichuan Basin | |
| CN112746835A (zh) | 一种优化的深层页岩气地质甜点测井综合评价方法 | |
| CN112784404B (zh) | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 | |
| CN115857047A (zh) | 一种地震储层综合预测方法 | |
| CN112228050B (zh) | 一种定量评价致密油储层宏观非均质性的方法及其应用 | |
| CN108345963B (zh) | 碎屑岩储层胶结减孔量的定量预测方法 | |
| CN108446788A (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 |