CN108364096B - 基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 - Google Patents
基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 Download PDFInfo
- Publication number
- CN108364096B CN108364096B CN201810116730.6A CN201810116730A CN108364096B CN 108364096 B CN108364096 B CN 108364096B CN 201810116730 A CN201810116730 A CN 201810116730A CN 108364096 B CN108364096 B CN 108364096B
- Authority
- CN
- China
- Prior art keywords
- delta
- facies
- river
- lake
- erosion
- 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
- 230000003628 erosive effect Effects 0.000 title claims abstract description 160
- 239000011435 rock Substances 0.000 title claims abstract description 60
- 238000000034 method Methods 0.000 title claims abstract description 59
- 230000008569 process Effects 0.000 title claims abstract description 34
- 239000011148 porous material Substances 0.000 title claims abstract description 31
- 230000004044 response Effects 0.000 title claims abstract description 13
- 208000035126 Facies Diseases 0.000 claims description 186
- 238000005260 corrosion Methods 0.000 claims description 83
- 230000007797 corrosion Effects 0.000 claims description 83
- 230000009471 action Effects 0.000 claims description 53
- 239000012530 fluid Substances 0.000 claims description 39
- 230000000694 effects Effects 0.000 claims description 38
- 230000015572 biosynthetic process Effects 0.000 claims description 26
- 239000010433 feldspar Substances 0.000 claims description 25
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 20
- 239000010453 quartz Substances 0.000 claims description 19
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 19
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 18
- 238000009933 burial Methods 0.000 claims description 15
- 239000004576 sand Substances 0.000 claims description 15
- 229910052500 inorganic mineral Inorganic materials 0.000 claims description 13
- 239000011707 mineral Substances 0.000 claims description 13
- 229910021532 Calcite Inorganic materials 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 11
- 238000013178 mathematical model Methods 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 10
- 229910000514 dolomite Inorganic materials 0.000 claims description 9
- 239000010459 dolomite Substances 0.000 claims description 9
- 239000002131 composite material Substances 0.000 claims description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 7
- 238000004090 dissolution Methods 0.000 claims description 6
- 230000002378 acidificating effect Effects 0.000 claims description 4
- 238000005253 cladding Methods 0.000 claims description 4
- 238000012886 linear function Methods 0.000 claims description 4
- 230000009467 reduction Effects 0.000 claims description 4
- 239000004927 clay Substances 0.000 claims description 3
- 239000013505 freshwater Substances 0.000 claims description 3
- 239000004615 ingredient Substances 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 239000011034 rock crystal Substances 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 19
- 230000008021 deposition Effects 0.000 description 16
- 230000008859 change Effects 0.000 description 5
- 230000000704 physical effect Effects 0.000 description 4
- 238000004062 sedimentation Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 3
- 239000000919 ceramic Substances 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000002734 clay mineral Substances 0.000 description 1
- 238000005056 compaction Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 231100000069 corrosive reaction Toxicity 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000002195 soluble material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Strategic Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- General Business, Economics & Management (AREA)
- General Physics & Mathematics (AREA)
- Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Health & Medical Sciences (AREA)
- Development Economics (AREA)
- Agronomy & Crop Science (AREA)
- General Health & Medical Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Primary Health Care (AREA)
- Life Sciences & Earth Sciences (AREA)
- Animal Husbandry (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，该方法提供了一种基于成岩过程的碎屑岩储层溶蚀增孔量的定量预测方法，该方法通过查清不同地质过程的关键地质因素及其差异以及不同地质过程中各地质因素对储层孔隙度的影响，构建孔隙响应定量方程，依据定量响应方程模拟储层溶蚀作用下孔隙演化过程，再现沉积物埋藏成岩过程中溶蚀作用对储层物性变化规律影响，为油藏评价提供依据。
Description
技术领域
本发明涉及油气勘探开发技术领域，具体地指一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法。
背景技术
碎屑岩骨架颗粒在埋藏成岩过程中的溶解作用及其与碎屑岩储层中次生孔隙形成机制之间关系的研究，一直是沉积学界所关注的热点问题之一。在储层演化过程中，由于孔隙水成分的改变，导致岩石骨架以及颗粒之间溶蚀物大量溶解，形成次生溶蚀孔隙，使储层孔隙度增大。这种次生溶蚀孔隙对改善储层物性的重要性近年来受到愈来愈多的重视。影响溶解作用的因素很多，如碎屑岩储层中岩石颗粒粒度的大小、孔隙渗透的好坏、可溶性物质的多少、流体的酸碱性、流体流速的快慢等都影响次生孔隙形成。通常情况下，较粗的颗粒、较好的渗透物性、多的可溶性物质、酸性流体介质、流体流速较快有利于次生孔隙形成。
前人在研究碎屑岩储层溶蚀作用的时候，多是在实验室条件下关注单个矿物(石英、长石、方解石)溶蚀作用研究，由于受实验条件限制，模拟参数的选择多为温度、压力、酸碱度，而不能综合考虑多种地质因素作用下溶蚀作用的变化规律，因此其地质模型的合理性值得商榷。通过岩石薄片观测，能简单、有效、客观地观察各种溶蚀现象，由于受样品数量的限制，其结论往往都是以偏概全，不能对整个研究区溶蚀作用作一个客观全面的反映。由于成岩过程中，溶蚀类型有多有少，溶蚀时间有先有后，溶蚀的强度有强有弱，从而使得定量预测碎屑岩储层溶蚀量存在很大难度。如果通过研究碎屑岩储层演化过程中影响和控制孔隙值变化地质因素及其差异性，构建影响孔隙演化定量响应关系模型，再现沉积物埋藏成岩过程中溶蚀作用对储层物性变化规律影响，就能对碎屑岩溶蚀增孔量进行很好定量预测，从而为盆地储层物性研究提供新的思路和依据。
成岩作用是一切储层形成必经过程，有什么样的作用过程，就有什么样的响应特征。储层中物性变化是地质历史中压实、溶蚀以及溶蚀作用等的共同结果。考虑储层演化过程中不同时期储层发生的各种溶蚀作用，特别是不同时期溶蚀作用发生时的条件及其差异，以及这些条件及差异作用下的溶蚀结果，可以为碎屑岩储层孔隙预测提供很好的手段。
发明内容
本发明针对上述存在的问题，提供了一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，该方法基于碎屑岩成岩过程中不同埋深、不同温度、不同演化阶段、不同成岩条件下矿物溶蚀作用研究，建立起多种地质参数下矿物溶蚀响应机制以及溶蚀数学模型，确定储层演化过程中原始孔隙减少情况，为油藏评价提供依据。
为了实现上述目的，本发明提供了一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测
方法，包括以下步骤：
1)收集研究区地质资料
包括：研究区构造等值线图、研究区沉积相图、研究区岩性分布图、研究区不同时期流体分布图、研究区溶蚀作用数据。
2)建立原始资料数据库
(1)沉积相类型数据库F_m
对碎屑岩储层沉积相类型进行概括和分类建立沉积相类型数据库F_m；
(2)岩性数据库R_n
根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库；
(3)流体数据库P_o
根据储层孔隙中流体性质，建立流体数据库；
(4)成岩阶段数据库S_q
(5)建立溶蚀作用数据库D_c
根据岩石中不同矿物Mi的溶蚀率建立溶蚀作用数据库D_c；
(6)建立溶蚀级别数学模型B_k
建立起不同类型溶蚀物Mi的溶蚀级别数学模型B_k，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
3)建立成岩阶段预测模型
(1)确定演化期次
根据目的层上覆地层发育情况划分目标层在地质历史时期所经历的演化阶段，目的层L上覆地层有i层，即从上到下依次标记为：L1、L2、L3……Li-1、Li，则目的层在地质历史时期演化阶段共有i个，按照时间演化，目的层L演化期次先后顺序依次计为Li、Li-1、Li-2……L2、L1；
(2)确定不同期次埋深
在(1)的基础下，计算目的层L不同演化期次的埋深，其计算公式如下：
Dep(Li)＝H(L)-H(Li)；
Dep(Li-1)＝H(L)-H(Li-1)；
Dep(Li-2)＝H(L)-H(Li-2)；
……
Dep(L2)＝H(L)-H2；
Dep(L1)＝H(L)-H1；
其中：H(L)为目的层L顶界面，H(Li)为上覆层Li顶界面，Dep(Li)为目的层Li阶段埋深；
(3)确定不同期次地层温度
地质体在埋藏的过程中，温度的大小可表示为与深度的线性函数关系，通过该模型，可计算目的层L在不同时期、不同深度、不同位置地层温度；
温度计算模型公式：
T＝T0+c*(D(Ti)-H0)
其中T0为常温带温度，c为常数，Dep(Ti)为目的层Ti阶段埋深，H0为恒温带埋深，为常数，T为目标层古地温；
(4)确定不同期次地层成岩阶段
当T∈[20～65)，DS为早成岩阶段早期，即为ⅠA；
当T∈[65～85)，DS为早成岩阶段晚期，即为ⅠB；
当T∈[85～140)，DS为中成岩阶段早期，即ⅡA；
当T∈[140～175)，DS为中成岩阶段晚期，即为ⅡB；
当T∈[175～200)，DS为晚成岩阶段，即为Ⅲ；
4)建立不同成岩阶段溶蚀作用预测模型
(1)研究区储层网格化
对研究区储层网格化，即为将空间上不均匀分布的数据按一定方法(如滑动平均法、克里格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程；研究区储层的每个网格可用Wi(X，Y)表示；
(2)确定网格属性
a.根据步骤3)确定网格Wi(X，Y)不同演化期次成岩阶段S_q；
b.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；
c.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；
d.根据研究区流体数据确定网格Wi(X，Y)的流体属性P_o；
e.根据上述步骤a、b、c和d确定每个网格不同演化期次Wi(X，Y)的综合属性G(S_q，F_m，R_n，P_o)，即为：
Wi(X，Y)＝S_q+F_m+R_n+P_o；
(3)基于研究区地质资料，确定研究区已知井Hj所在网格Wj(X，Y)不同演化期次的成岩阶段S_q_j以及溶蚀作用D_s_j；
(4)确定已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段的综合属性Gj(F_m_j，R_n_j，P_o_j)，建立已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段溶蚀作用D_s与综合属性G(F_m_j，R_n_j，P_o_j)的对应关系，即为：
Wj(S_q_j，D_s_j)＝F_m_j+R_n_j+P_o_j；
(5)任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)成岩阶段S_q时的综合属性G(F_m，R_n，P_o)，即为：
G i＝(F_m_i，R_n_i，P_o_i)；
(6)将网格Wi(X，Y)成岩阶段S_q时的溶蚀作用综合属性G i与井Hj所在网格Wj(X，Y)成岩阶段S_q的综合属性G j相比较，即为：
Gi-Gj＝(F_m_i，R_n_i，P_o_i)-(F_m_j+R_n_j+P_o_j)；
如果满足G i-G j＝0，则未知网格Wi(X，Y)的溶蚀作用D_s_i与井j具有相同的溶蚀作用D_s_j；
如果G i-G j≠0，则按不同属性优先级顺序，即：沉积相F_m一级、岩性R_n二级、流体性质P_o三级对未知网格Wi(X，Y)的溶蚀作用进行判识，即：
ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，P_o_i-P_o_j≠0；
ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，P_o_i-P_o_j≠0；
未知网格Wi(X，Y)的溶蚀作用D_s_i满足条件i时，未知网格Wi(X，Y)与具有相同属性F_j的网格Wj(X，Y)的溶蚀作用相同；
D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j的网格Wj(X，Y)的溶蚀作用相同；
(7)根据每个网格Wi(X，Y)不同演化期次的溶蚀作用确定研究区溶蚀作用演化；
5)计算不同演化期次溶蚀量
(1)根据成岩序列数据，确定已知井所在网格Wj(X，Y)的成岩阶段S_q所对应的溶蚀作用下的溶蚀级别B_k_j，即Wj(S_q_j，D_s_j)＝B_k_j；
(2)将Wi(X，Y)的Wi(S_q_j，D_s_j)与Wj(S_q_j，D_s_j)比对，满足Wi(S_q_j，D_s_j)＝Wj(S_q_j，D_s_j)，则Wi(X，Y)与Wj(S_q_j，D_s_j)具有相同的溶蚀级别B_k，即B_k_i＝B_k_j；如果出现同一未知网格Wi(X，Y)与多个已知网格Wj(S_q_j，D_s_j)相同，则B_k_i＝1/n∑B_k_j，n为相同网格Wj(S_q_j，D_s_j)的网格数；
(3)根据网格Wi(X，Y)不同演化期次的溶蚀级别B_k_i计算每个网格Wi(X，Y)不同演化期次的溶蚀率PCem；
(4)根据每个网格Wi(X，Y)的溶蚀率RCem计算溶蚀量VCem，
其中
为原始孔隙；
(5)根据每个网格Wi(X，Y)不同演化期次的溶蚀量确定目标储层溶蚀量，即溶蚀增孔量。
进一步地，所述步骤1)中，1)收集研究区地质资料包括：研究区构造等值线图、研究区沉积相图、研究区岩性分布图、研究区不同时期流体分布图、研究区溶蚀作用数据。
再进一步地，所述步骤2)中，
(1)沉积相类型数据库F_m包括：冲积扇-旱扇-扇根F_1、冲积扇-旱扇-扇中F_2、冲积扇-旱扇-扇缘F_3、冲积扇-湿扇-扇根F_4、冲积扇-湿扇-扇中F_5、冲积扇-湿扇-扇缘F_6、河流相-曲流河-河床亚相-河床滞留F_7、河流相-曲流河-河床亚相-边滩F_8、河流相-曲流河-堤岸亚相-天然堤F_9、河流相-曲流河-堤岸亚相-决口扇F_10、河流相-曲流河-河漫亚相-河漫滩F_11、河流相-曲流河-泛滥盆地-河漫湖泊F_12、河流相-曲流河-河漫沼泽F_13、河流相-辫状河-牛轭湖F_14、河流相-辫状河-河床滞留F_15、河流相-辫状河-心滩F_16、河道F_17、河流相-辫状河-泛滥平原F_18、湖泊相-断陷型-湖成三角洲F_19、湖泊相-断陷型-滨湖F_20、湖泊相-断陷型-浅湖F_21、湖泊相-断陷型-半深湖F_22、湖泊相-断陷型-深湖F_23、湖泊相-断陷型-湖湾F_24、湖泊相-坳陷型-湖成三角洲F_25、湖泊相-坳陷型-滨湖F_26、湖泊相-坳陷型-浅湖F_27、湖泊相-坳陷型-半深湖F_28、湖泊相-坳陷型-深湖F_29、湖泊相-坳陷型-湖湾F_30、湖泊相-前陆型-湖成三角洲F_31、湖泊相-前陆型-滨湖F_32、湖泊相-前陆型-浅湖F_33、湖泊相-前陆型-半深湖F_34、湖泊相-前陆型-深湖F_35、湖泊相-前陆型-湖湾F_36、三角洲相-辫状河三角洲-三角洲平原-分支F_37、三角洲相-辫状河三角洲-三角洲平原-河道F_38、三角洲相-辫状河三角洲-三角洲平原-陆上天然堤F_39、三角洲相-辫状河三角洲-三角洲平原-决口扇F_40、三角洲相-辫状河三角洲-三角洲平原-沼泽F_41、三角洲相-辫状河三角洲-三角洲平原-淡水湖泊F_42、三角洲相-辫状河三角洲-三角洲前缘-水下分支河道F_43、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤F_44、三角洲相-辫状河三角洲-三角洲前缘-支流间湾F_45、三角洲相-辫状河三角洲-三角洲前缘-分支河口砂坝F_46、三角洲相-辫状河三角洲-三角洲前缘-远砂坝F_47、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂F_48、三角洲相-辫状河三角洲-前三角洲F_49、三角洲相-扇三角洲-扇三角洲平原-泥石流F_50、三角洲相-扇三角洲-扇三角洲平原-河道充填F_51、三角洲相-扇三角洲-扇三角洲平原-漫滩F_52、三角洲相-扇三角洲-扇三角洲前缘-碎屑流F_53、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道F_54、三角洲相-扇三角洲-扇三角洲前缘-支流间湾F_55、三角洲相-扇三角洲-扇三角洲前缘-河口砂坝F_56、三角洲相-扇三角洲-扇三角洲前缘-远砂坝F_57、三角洲相-扇三角洲-前扇三角洲F_58；
(2)岩性数据库R_n包括细砾岩R_1、中砾岩R_2、粗砾岩R_3、巨砾岩R_4、石英砂岩R_5、长石质石英砂岩R_6、岩屑质石英砂岩R_7、长石岩屑质石英砂岩R_8、长石砂岩R_9、岩屑质长石砂岩R_10、岩屑长石砂岩R_11、岩屑砂岩R_12、长石质岩屑砂岩R_13、长石岩屑砂岩R_14、黏土R_15、泥岩R_16、页岩R_17；
(3)流体数据库P_o包括淡水环境P_1、酸性水环境P_2、碱性水环境P_3；
(4)成岩阶段数据库S_q包括：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；
(5)溶蚀作用数据库D_c中溶蚀作用类型包括：
石英强溶蚀作用D1-1：其溶蚀率RCem＞10％；
石英中溶蚀作用D1-2：其溶蚀率5％＜RCem≤10％；
石英弱溶蚀作用D1-3：其溶蚀率RCem≤5％；
长石强溶蚀作用D2-1：其溶蚀率RCem＞10％；
长石中溶蚀作用D2-2：其溶蚀率5％＜RCem≤10％；
长石弱溶蚀作用D2-3：其溶蚀率RCem≤5％；
方解石强溶蚀作用D3-1：其溶蚀率RCem＞10％；
方解石中溶蚀作用D3-2：其溶蚀率5％＜RCem≤10％；
方解石弱溶蚀作用D3-3：其溶蚀率RCem≤5％；
白云石强溶蚀作用D4-1：其溶蚀率RCem＞10％；
白云石中溶蚀作用D4-2：其溶蚀率5％＜RCem≤10％；
白云石弱溶蚀作用D4-3：其溶蚀率RCem≤5％；
(6)溶蚀级别数学模型B_k包括以下几类：
B_1＝Ⅰ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/10；
B_2＝Ⅱ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/8；
B_3＝Ⅲ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/6；
B_4＝Ⅳ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/5；
B_5＝Ⅴ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/3；
B_6＝Ⅵ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.6；
B_7＝Ⅶ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.5；
B_8＝Ⅷ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.4；
B_9＝Ⅸ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.2；
B_10＝Ⅹ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1；
其中，X埋藏时间，Ma，Y溶蚀率，％。
本发明有益效果
提供了一种基于过程响应的碎屑岩储层溶蚀增孔量的定量预测方法，该方法基于碎屑岩成岩过程中不同演化阶段、不同成岩条件下矿物溶蚀作用研究，建立起不同条件下矿物溶蚀响应机制以及溶蚀数学模型，再现沉积物埋藏成岩过程中溶蚀作用对储层物性变化规律影响，确定储层演化过程中原始孔隙增加情况，能简单、有效、客观地为油藏评价提供依据，从而降低油气勘探风险。
附图说明
图1为研究区平原组构造等值线图
图2为研究区明化镇组构造等值线图
图3为研究区馆陶组构造等值线图
图4为研究区东营组二段构造等值线图
图5为研究区东营组二段沉积相图
图6为研究区东营组二段岩性分布图
图7为研究区馆陶组沉积时东营组二段流体分布图
图8为研究区明化镇组沉积时东营组二段流体分布图
图9为研究区平原组沉积时东营组二段流体分布图
图10为研究区馆陶组沉积时东营组二段地温
图11为研究区明化镇组沉积时东营组二段地温
图12为研究区平原组沉积时东营组二段地温
图13为研究区馆陶组沉积时东营组二段成岩阶段图
图14为研究区明化镇组沉积时东营组二段成岩阶段图
图15为研究区平原组沉积时东营组二段成岩阶段图
图16为研究区馆陶组沉积时东营组二段溶蚀作用图
图17为研究区明化镇组沉积时东营组二段溶蚀作用图
图18为研究区平原组沉积时东营组二段溶蚀作用图
图19为研究区馆陶组沉积时东营组二段溶蚀量分布图
图20为研究区明化镇组沉积时东营组二段溶蚀量分布图
图21为研究区平原组沉积时东营组二段溶蚀量分布图
具体实施方式
为了更好地解释本发明，以下结合具体实施例进一步阐明本发明的主要内容，但本发明的内容不仅仅局限于以下实施例。
实施例1
一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，包括以下步骤：
1)收集研究区地质资料，它包括：
①研究区东营组二段溶蚀作用表；
②研究区平原组构造等值线图；
③研究区明化镇组构造等值线图；
④研究区馆陶组构造等值线图；
⑤研究区东营组二段构造等值线图；
⑥研究区东营组二段沉积相图；
⑦研究区东营组二段岩性分布图；
⑧研究区馆陶组沉积时东营组二段流体分布图；
⑨研究区明化镇组沉积时东营组二段流体分布图；
⑩研究区平原组沉积时东营组二段流体分布图；
2)建立原始资料数据库
(1)沉积相类型数据库F_m
对碎屑岩储层沉积相类型进行概括和分类，其中包括：冲积扇-旱扇-扇根F_1、冲积扇-旱扇-扇中F_2、冲积扇-旱扇-扇缘F_3、冲积扇-湿扇-扇根F_4、冲积扇-湿扇-扇中F_5、冲积扇-湿扇-扇缘F_6、河流相-曲流河-河床亚相-河床滞留F_7、河流相-曲流河-河床亚相-边滩F_8、河流相-曲流河-堤岸亚相-天然堤F_9、河流相-曲流河-堤岸亚相-决口扇F_10、河流相-曲流河-河漫亚相-河漫滩F_11、河流相-曲流河-泛滥盆地-河漫湖泊F_12、河流相-曲流河-河漫沼泽F_13、河流相-辫状河-牛轭湖F_14、河流相-辫状河-河床滞留F_15、河流相-辫状河-心滩F_16、河道F_17、河流相-辫状河-泛滥平原F_18、湖泊相-断陷型-湖成三角洲F_19、湖泊相-断陷型-滨湖F_20、湖泊相-断陷型-浅湖F_21、湖泊相-断陷型-半深湖F_22、湖泊相-断陷型-深湖F_23、湖泊相-断陷型-湖湾F_24、湖泊相-坳陷型-湖成三角洲F_25、湖泊相-坳陷型-滨湖F_26、湖泊相-坳陷型-浅湖F_27、湖泊相-坳陷型-半深湖F_28、湖泊相-坳陷型-深湖F_29、湖泊相-坳陷型-湖湾F_30、湖泊相-前陆型-湖成三角洲F_31、湖泊相-前陆型-滨湖F_32、湖泊相-前陆型-浅湖F_33、湖泊相-前陆型-半深湖F_34、湖泊相-前陆型-深湖F_35、湖泊相-前陆型-湖湾F_36、三角洲相-辫状河三角洲-三角洲平原-分支F_37、三角洲相-辫状河三角洲-三角洲平原-河道F_38、三角洲相-辫状河三角洲-三角洲平原-陆上天然堤F_39、三角洲相-辫状河三角洲-三角洲平原-决口扇F_40、三角洲相-辫状河三角洲-三角洲平原-沼泽F_41、三角洲相-辫状河三角洲-三角洲平原-淡水湖泊F_42、三角洲相-辫状河三角洲-三角洲前缘-水下分支河道F_43、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤F_44、三角洲相-辫状河三角洲-三角洲前缘-支流间湾F_45、三角洲相-辫状河三角洲-三角洲前缘-分支河口砂坝F_46、三角洲相-辫状河三角洲-三角洲前缘-远砂坝F_47、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂F_48、三角洲相-辫状河三角洲-前三角洲F_49、三角洲相-扇三角洲-扇三角洲平原-泥石流F_50、三角洲相-扇三角洲-扇三角洲平原-河道充填F_51、三角洲相-扇三角洲-扇三角洲平原-漫滩F_52、三角洲相-扇三角洲-扇三角洲前缘-碎屑流F_53、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道F_54、三角洲相-扇三角洲-扇三角洲前缘-支流间湾F_55、三角洲相-扇三角洲-扇三角洲前缘-河口砂坝F_56、三角洲相-扇三角洲-扇三角洲前缘-远砂坝F_57、三角洲相-扇三角洲-前扇三角洲F_58；
(2)岩性数据库R_n
根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库，它包括细砾岩R_1、中砾岩R_2、粗砾岩R_3、巨砾岩R_4、石英砂岩R_5、长石质石英砂岩R_6、岩屑质石英砂岩R_7、长石岩屑质石英砂岩R_8、长石砂岩R_9、岩屑质长石砂岩R_10、岩屑长石砂岩R_11、岩屑砂岩R_12、长石质岩屑砂岩R_13、长石岩屑砂岩R_14、黏土R_15、泥岩R_16、页岩R_17；
(3)流体数据库P_o
根据储层孔隙中流体性质，建立流体数据库：它包括淡水环境P_1、酸性水环境P_2、碱性水环境P_3；
(4)成岩阶段数据库S_q
成岩阶段数据库S_q包括：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；
(5)建立溶蚀作用数据库D_c
根据岩石中不同矿物Mi的溶蚀率建立溶蚀作用数据库，溶蚀作用类型包括：
石英强溶蚀作用D1-1：其溶蚀率RCem＞10％；
石英中溶蚀作用D1-2：其溶蚀率5％＜RCem≤10％；
石英弱溶蚀作用D1-3：其溶蚀率RCem≤5％；
长石强溶蚀作用D2-1：其溶蚀率RCem＞10％；
长石中溶蚀作用D2-2：其溶蚀率5％＜RCem≤10％；
长石弱溶蚀作用D2-3：其溶蚀率RCem≤5％；
方解石强溶蚀作用D3-1：其溶蚀率RCem＞10％；
方解石中溶蚀作用D3-2：其溶蚀率5％＜RCem≤10％；
方解石弱溶蚀作用D3-3：其溶蚀率RCem≤5％；
白云石强溶蚀作用D4-1：其溶蚀率RCem＞10％；
白云石中溶蚀作用D4-2：其溶蚀率5％＜RCem≤10％；
白云石弱溶蚀作用D4-3：其溶蚀率RCem≤5％；
(6)建立溶蚀级别数学模型B_k
建立起不同类型溶蚀物Mi的数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量。
B_1＝Ⅰ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/10；
B_2＝Ⅱ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/8；
B_3＝Ⅲ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/6；
B_4＝Ⅳ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/5；
B_5＝Ⅴ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/3；
B_6＝Ⅵ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.6；
B_7＝Ⅶ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.5；
B_8＝Ⅷ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.4；
B_9＝Ⅸ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1.2；
B_10＝Ⅹ：Y＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/1；
其中，X埋藏时间，Ma，Y溶蚀率，％；
3)建立成岩阶段预测模型
(1)确定演化期次
根据目的层上覆地层发育情况划分目标层在地质历史时期所经历的演化阶段，目的层L上覆地层有i层，即从上到下依次标记为：L1、L2、L3……Li-1、Li，则目的层在地质历史时期演化阶段共有i个，按照时间演化，目的层L演化期次先后顺序依次计为Li、Li-1、Li-2……L2、L1；
(2)确定不同期次埋深
在(1)的基础下，计算目的层L不同演化期次的埋深，其计算公式如下：
Dep(Li)＝H(L)-H(Li)；
Dep(Li-1)＝H(L)-H(Li-1)；
Dep(Li-2)＝H(L)-H(Li-2)；
……
Dep(L2)＝H(L)-H2；
Dep(L1)＝H(L)-H1；
其中：H(L)为目的层L顶界面，H(Li)为上覆层Li顶界面，Dep(Li)为目的层Li阶段埋深；
(3)确定不同期次地层温度
地质体在埋藏的过程中，温度的大小可表示为与深度的线性函数关系，通过该模型，可计算目的层L在不同时期、不同深度、不同位置地层温度；
温度计算模型公式：
T＝T0+c*(D(Ti)-H0)
其中T0为常温带温度，c为常数，Dep(Ti)为目的层Ti阶段埋深，H0为恒温带埋深，为常数，T为目标层古地温；
(4)确定不同期次地层成岩阶段
当T∈[20～65)，DS为早成岩阶段早期，即为ⅠA；
当T∈[65～85)，DS为早成岩阶段晚期，即为ⅠB；
当T∈[85～140)，DS为中成岩阶段早期，即ⅡA；
当T∈[140～175)，DS为中成岩阶段晚期，即为ⅡB；
当T∈[175～200)，DS为晚成岩阶段，即为Ⅲ；
4)建立不同成岩阶段溶蚀作用预测模型
(1)研究区储层网格化
对研究区储层网格化，即为将空间上不均匀分布的数据按一定方法(如滑动平均法、克里格法或其他适当的数值推算方法)归算成规则网格中的代表值(趋势值)的过程；研究区储层的每个网格可用Wi(X，Y)表示；
(2)确定网格属性
a.根据步骤3)确定网格Wi(X，Y)不同演化期次成岩阶段S_q；
b.根据研究区沉积相数据确定网格Wi(X，Y)的沉积相属性F_m；
c.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；
d.根据研究区流体数据确定网格Wi(X，Y)的流体属性P_o；
e.根据上述步骤a、b、c和d确定每个网格不同演化期次Wi(X，Y)的综合属性G(S_q，F_m，R_n，P_o)，即为：
Wi(X，Y)＝S_q+F_m+R_n+P_o；
(3)基于研究区地质资料，确定研究区已知井Hj所在网格Wj(X，Y)不同演化期次的成岩阶段S_q_j以及溶蚀作用D_s_j；
(4)确定已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段的综合属性Gj(F_m_j，R_n_j，P_o_j)，建立已知井Hj所在网格Wj(X，Y)不同演化期次不同演化阶段溶蚀作用D_s与综合属性G(F_m_j，R_n_j，P_o_j)的对应关系，即为：
Wj(S_q_j，D_s_j)＝F_m_j+R_n_j+P_o_j；
(5)任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)成岩阶段S_q时的综合属性G(F_m，R_n，P_o)，即为：
G i＝(F_m_i，R_n_i，P_o_i)；
(6)将网格Wi(X，Y)成岩阶段S_q时的溶蚀作用综合属性G i与井Hj所在网格Wj(X，Y)成岩阶段S_q的综合属性G j相比较，即为：
Gi-Gj＝(F_m_i，R_n_i，P_o_i)-(F_m_j+R_n_j+P_o_j)；
如果满足G i-G j＝0，则未知网格Wi(X，Y)的溶蚀作用D_s_i与井j具有相同的溶蚀作用D_s_j；
如果G i-G j≠0，则按不同属性优先级顺序，即：沉积相F_m一级、岩性R_n二级、流体性质P_o三级对未知网格Wi(X，Y)的溶蚀作用进行判识，即：
ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，P_o_i-P_o_j≠0；
ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，P_o_i-P_o_j≠0；
未知网格Wi(X，Y)的溶蚀作用D_s_i满足条件i时，未知网格Wi(X，Y)与具有相同属性F_j的网格Wj(X，Y)的溶蚀作用相同；
D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j的网格Wj(X，Y)的溶蚀作用相同；
(7)根据每个网格Wi(X，Y)不同演化期次的溶蚀作用确定研究区溶蚀作用演化；
5)计算不同演化期次溶蚀量
(1)根据成岩序列数据，确定已知井所在网格Wj(X，Y)的成岩阶段S_q所对应的溶蚀作用下的溶蚀级别B_k_j，即Wj(S_q_j，D_s_j)＝B_k_j；
(2)将Wi(X，Y)的Wi(S_q_j，D_s_j)与Wj(S_q_j，D_s_j)比对，满足Wi(S_q_j，D_s_j)＝Wj(S_q_j，D_s_j)，则Wi(X，Y)与Wj(S_q_j，D_s_j)具有相同的溶蚀级别B_k，即B_k_i＝B_k_j；如果出现同一未知网格Wi(X，Y)与多个已知网格Wj(S_q_j，D_s_j)相同，则B_k_i＝1/n∑B_k_j，n为相同网格Wj(S_q_j，D_s_j)的网格数；
(3)根据网格Wi(X，Y)不同演化期次的溶蚀级别B_k_i计算每个网格Wi(X，Y)不同演化期次的溶蚀率PCem；
(4)根据每个网格Wi(X，Y)的溶蚀率RCem计算溶蚀量VCem，
其中
为原始孔隙；
(5)根据每个网格Wi(X，Y)不同演化期次的溶蚀量确定目标储层溶蚀量，即溶蚀增孔量；
实施例2
基于成岩过程的碎屑岩储层溶蚀增孔量的定量预测方法预测渤中凹陷东营组二段储层，具体步骤如下：
1)收集研究区地质资料
收集渤中凹陷东营组二段储层的地质资料，并按照原始资料数据库对各种类型数据进行分类。收集的地质资料包括以下内容：研究区东营组二段溶蚀作用表(表1)、研究区平原组构造等值线图(图1)、研究区明化镇组构造等值线图(图2)、研究区馆陶组构造等值线图(图3)、研究区东营组二段构造等值线图(图4)、研究区东营组二段沉积相图(图5)、研究区东营组二段岩性分布图(图6)、研究区馆陶组沉积时东营组二段流体分布图(图7)、研究区明化镇组沉积时东营组二段流体分布图(图8)、研究区平原组沉积时东营组二段流体分布图(图9)。
表1研究区东营组二段溶蚀作用表
2)确定成岩阶段
(1)确定演化期次
东营组二段上覆地层从上到下依次是平原组、明化镇组以及馆陶组，即从上到下依次标记为L1、L2、L3，则研究区演化期次共有3个，按时间演化，东营组二段L的演化期次先后顺序依次计为L3、L2、L1；
(2)确定不同期次埋深
在(1)的基础下，计算东营组二段L不同演化期次的埋深，其计算公式如下：
Dep(Li)＝H(L)-H(Li)；
Dep(Li-1)＝H(L)-H(Li-1)；
Dep(Li-2)＝H(L)-H(Li-2)；
……
Dep(L2)＝H(L)-H2；
Dep(L1)＝H(L)-H1；
其中：H(L)为目的层L顶界面，H(Li)为上覆层Li顶界面，Dep(Li)为目的层Li阶段埋深；
(3)确定不同期次地层温度
地质体在埋藏的过程中，温度的大小可表示为与深度的线性函数关系，通过该模型，可计算目的层L在不同时期、不同深度、不同位置地层温度；
温度计算模型公式：
T＝T0+c*(D(Ti)-H0)
其中T0＝21℃，c＝0.0031℃/m，Dep(Ti)为目的层Ti阶段埋深，H0＝30m，T为目标层古地温；
(4)确定不同期次地层成岩阶段
当T∈[20～65)，DS为早成岩阶段早期，即为ⅠA；
当T∈[65～85)，DS为早成岩阶段晚期，即为ⅠB；
当T∈[85～140)，DS为中成岩阶段早期，即ⅡA；
当T∈[140～175)，DS为中成岩阶段晚期，即为ⅡB；
当T∈[175～200)，DS为晚成岩阶段，即为Ⅲ；
3)确定演化期次L3溶蚀作用
(1)将东营组二段L网格化，每个网格可用Wi(X，Y)表示；
(2)确定演化期次L3网格属性
a.根据步骤2确定网格Wi(X，Y)演化期次L3成岩阶段S_q；
b.根据研究区沉积相数据确定网格Wi(X，Y)演化期次L3的沉积相属性F_m；
c.根据研究区岩性数据确定网格Wi(X，Y)的岩性属性R_n；
d.根据研究区流体数据确定网格Wi(X，Y)演化期次L3的流体属性P_o；
e.根据上述步骤a、b、c和d确定每个网格Wi(X，Y)演化期次L3的综合属性G(S_q，F_m，R_n，P_o)，即为：
Wi(X，Y)＝S_q+F_m+R_n+P_o；
(3)基于研究区地质资料，确定研究区已知井Hj所在网格Wj(X，Y)演化期次L3的成岩阶段S_q_j以及溶蚀作用D_s_j；
(4)确定已知井Hj所在网格Wj(X，Y)演化期次L3的综合属性Gj(F_m_j，R_n_j，P_o_j)，建立已知井Hj所在网格Wj(X，Y)演化期次L3溶蚀作用D_s与综合属性G(F_m_j，R_n_j，P_o_j)的对应关系，即为：
Wj(S_q_j，D_s_j)＝F_m_j+R_n_j+P_o_j；
(5)任取一未知网格Wi(X，Y)，确定未知网格Wi(X，Y)成岩阶段S_q时的综合属性G(F_m，R_n，P_o)，即为：
G i＝(F_m_i，R_n_i，P_o_i)；
(6)将网格Wi(X，Y)成岩阶段S_q时的溶蚀作用综合属性G i与井Hj所在网格Wj(X，Y)成岩阶段S_q的综合属性G j相比较，即为：
Gi-Gj＝(F_m_i，R_n_i，P_o_i)-(F_m_j+R_n_j+P_o_j)；
如果满足G i-G j＝0，则未知网格Wi(X，Y)的溶蚀作用D_s_i与井j具有相同的溶蚀作用D_s_j；
如果G i-G j≠0，则按不同属性优先级顺序，即：沉积相F_m一级、岩性R_n二级、流体性质P_o三级对未知网格Wi(X，Y)的溶蚀作用进行判识，即：
ⅰ：F_m_i-F_m_j＝0，R_n_i-R_n_j≠0，P_o_i-P_o_j≠0；
ⅱ：F_m_i-F_m_j＝0，R_n_i-R_n_j＝0，P_o_i-P_o_j≠0；
未知网格Wi(X，Y)的溶蚀作用D_s_i满足条件i时，未知网格Wi(X，Y)与具有相同属性F_j的网格Wj(X，Y)的溶蚀作用相同；
D_s_i满足条件ⅱ时，未知网格Wi(X，Y)与具有相同属性F_m_j、R_n_j的网格Wj(X，Y)的溶蚀作用相同；
如在时间t＝5.0Ma，埋深H＝250m处，网格W30＝(S_1，F_43，R_6，P_2)-(S_1_CFD17-3-1，F_m_CFD17-3-1，R_n_CFD17-3-1+P_o_CFD17-3-1)＝0，即网格W30与井CFD17-3-1具有相同的溶蚀作用，也即：W30＝D3-2；
如在时间t＝3.6Ma，埋深H＝180m处，网格W40＝(S_1，F_43，R_10，P_2)-(S_1_CFD18-2E-1，F_CFD18-2E-1，R_n_CFD18-2E-1，P_o_CFD18-2E-1)＝0，即网格W40与井CFD18-2E-1具有相同的溶蚀作用，也即：W40＝D3-3；
如在时间t＝5.3Ma，埋深H＝265m处，网格W50＝(S_1，F_26，R_11，P_2)网格属性满足条件ⅱ，即：F_m-F_m_BZ-13-1-3＝0，R_n-R_n_BZ-13-1-3＝0，P_o-P_o_BZ-13-1-3≠0；则W50与井BZ-13-1-3具有相同溶蚀作用，也即：W50＝D5-3；
(7)根据每个网格Wi(X，Y)演化期次L3的溶蚀作用确定研究区溶蚀作用演化(图16)。
4)计算演化期次L3溶蚀量
(1)根据成岩序列数据，确定已知井所在网格Wj(X，Y)的成岩阶段S_q所对应的溶蚀作用下的溶蚀级别B_k_j，即为：
Wj(S_q_j，D_s_j)＝B_k_j；
(2)将Wi(X，Y)的Wi(S_q_j，D_s_j)与Wj(S_q_j，D_s_j)比对，满足Wi(S_q_j，D_s_j)＝Wj(S_q_j，D_s_j)，则Wi(X，Y)与Wj(S_q_j，D_s_j)具有相同的溶蚀级别B_k，即B_k_i＝B_k_j；如果出现同一未知网格Wi(X，Y)与多个已知网格Wj(S_q_j，D_s_j)相同，则B_k_i＝1/n∑B_k_j，n为相同网格Wj(S_q_j，D_s_j)的网格数；
(3)根据网格Wi(X，Y)演化期次L3的溶蚀级别B_k_i计算每个网格Wi(X，Y)不同演化期次的溶蚀率PCem；
(4)根据每个网格Wi(X，Y)的溶蚀率RCem计算溶蚀量VCem，
其中
为原始孔隙；
如网格W30＝(S_1，D3-2)＝W(S_1_CFD17-3-1，D3-3_CFD17-3-1)，则网格W30演化期次L3时方解石矿物溶蚀级别B_k＝B_4，其计算公式如下：
RCal＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/8
＝[-4.859e-06*(5)^4+0.001028*(5)^3-0.0745*(5)^2+2.862*(5)+10.6]/8
＝2.9
如网格W40＝(S_1，D3-3)＝W(S_1_CFD18-2E-1，D3-3_CFD18-2E-1)＝W(S_1_CFD24-1-1，D3-3_CFD24-1-1)，则网格W40演化期次L3方解石矿物溶蚀级别B_k＝1/2∑(B_2+B_2)＝B_2；
RCal＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/8
＝[-4.859e-06*(3.6)^4+0.001028*(3.6)^3-0.0745*(3.6)^2+2.862*(3.6)+10.6]/8
＝2.5
如网格W50＝(S_1，D5-3)＝W(S_1_CFD16-1-2D，D5-3_CFD16-1-2D)＝W(S_1_CFD23-3-1，D5-3_CFD23-3-1)＝W(S_1_CFD16-2-1，D5-3_CFD16-2-1)＝W(S_1_BZ13-1-3，D5-3_BZ13-1-3)，则网格W50演化期次L3黏土矿物溶蚀级别B_k＝1/4∑(B_2+B_2+B_2+B_2)＝B_2；
RCla＝[-4.859e-06*(X)^4+0.001028*(X)^3-0.0745*(X)^2+2.862*(X)+10.6]/8
＝[-4.859e-06*(5.3)^4+0.001028*(5.3)^3-0.0745*(5.3)^2+2.862*(5.3)+10.6]/8
＝3.0
(5)根据每个网格Wi(X，Y)不同演化期次的溶蚀量确定目标储层溶蚀量，即溶蚀增孔量(图19)；
5)确定演化期次L2、L1溶蚀作用
利用步骤3，采用相同办法，确定东营组二段L演化期次L2、L1的溶蚀作用(图17、图18)；
6)确定演化期次L2、L1溶蚀量
利用步骤4)，采用相同办法，确定东营组二段L演化期次L2、L1的溶蚀量(图20、图21)。
其它未详细说明的部分均为现有技术。尽管上述实施例对本发明做出了详尽的描述，但它仅仅是本发明一部分实施例，而不是全部实施例，人们还可以根据本实施例在不经创造性前提下获得其他实施例，这些实施例都属于本发明保护范围。
Claims (3)
1.一种基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，其特征在于：包括以下步骤：
1）收集研究区地质资料
2）建立原始资料数据库
（1）沉积相类型数据库
对碎屑岩储层沉积相类型进行概括和分类建立沉积相类型数据库；
（2）岩性数据库
根据组成岩石颗粒的颗粒大小及矿物成分，将碎屑岩岩性进行划分，建立岩性数据库；
（3）流体数据库
根据储层孔隙中流体性质，建立流体数据库；
（4）成岩阶段数据库
（5）建立溶蚀作用数据库根据岩石中不同矿物的溶蚀率建立溶蚀作用数据库；
（6）建立溶蚀级别数学模型建立起不同类型溶蚀物的溶蚀级别数学模型，从而定量刻画碎屑岩在成岩过程中孔隙减小量；
3）建立成岩阶段预测模型
（1）确定演化期次
根据目的层上覆地层发育情况划分目标层在地质历史时期所经历的演化阶段，目的层L上覆地层有i层，即从上到下依次标记为：L1、L2、L3……Li-1、Li，则目的层在地质历史时期演化阶段共有i个，按照时间演化，目的层L演化期次先后顺序依次计为Li、Li-1、Li-2……L2、L1；
（2）确定不同期次埋深
在（1）的基础下，计算目的层L不同演化期次的埋深，其计算公式如下：
Dep(Li)=H(L)-H(Li)；
Dep (Li-1)=H(L)-H(Li-1)；
Dep (Li-2) =H(L)-H(Li-2)；
……
Dep (L2)=H(L)-H2；
Dep (L1)=H(L)-H1；
其中：H(L)为目的层L顶界面，H(Li) 为上覆层Li顶界面，Dep (Li)为目的层Li 阶段埋深；
（3）确定不同期次地层温度
地质体在埋藏的过程中，温度的大小可表示为与深度的线性函数关系，通过该模型，可计算目的层L在不同时期、不同深度、不同位置地层温度；
温度计算模型公式：
T = T0 + c*（D(Ti)-H0）
其中T0为常温带温度，c为常数，Dep(Ti)为目的层Ti 阶段埋深，H0为恒温带埋深，为常数，T为目标层古地温；
（4）确定不同期次地层成岩阶段
当T∈[20～65），DS为早成岩阶段早期，即为ⅠA；
当T∈[65～85），DS为早成岩阶段晚期，即为ⅠB；
当T∈[85～140），DS为中成岩阶段早期，即ⅡA；
当T∈[140～175），DS为中成岩阶段晚期，即为ⅡB；
当T∈[175～200），DS为晚成岩阶段，即为Ⅲ；
4）建立不同成岩阶段溶蚀作用预测模型
（1）研究区储层网格化
对研究区储层网格化，研究区储层的每个网格用
（X，Y）表示；
（2）确定网格属性
a. 根据步骤3）确定网格Wi（X，Y）不同演化期次成岩阶段S_q；
b. 根据研究区沉积相数据确定网格Wi（X，Y）的沉积相属性F_m；
c. 根据研究区岩性数据确定网格Wi（X，Y）的岩性属性R_n；
d. 根据研究区流体数据确定网格Wi（X，Y）的流体属性P_o；
e. 根据上述步骤a、b、c和d确定每个网格不同演化期次Wi（X，Y）的综合属性G（S_q，F_m，R_n，P_o），即为：
Wi（X，Y）= S_q +F_m+ R_n+P_o；
（3）基于研究区地质资料，确定研究区已知井Hj所在网格Wj（X，Y）不同演化期次的成岩阶段S_q_j以及溶蚀作用D_s_j；
（4）确定已知井Hj所在网格Wj（X，Y）不同演化期次不同演化阶段的综合属性Gj（F_m_j，R_n_j，P_o_j），建立已知井Hj所在网格Wj（X，Y）不同演化期次不同演化阶段溶蚀作用D_s_j与综合属性G（F_ m_j，R_n_j，P_o_j）的对应关系，即为：
Wj（S_q_j，D_s_j）= F_ m_j+R_n_j+P_o_j；
（5）任取一未知网格Wi（X，Y），确定未知网格Wi（X，Y）成岩阶段S_q时的综合属性G（F_m，R_n，P_o），即为：
G i = （F_m_i，R_n_i，P_o_i）；
（6）将网格Wi（X，Y）成岩阶段S_q时的溶蚀作用综合属性G i与井Hj所在网格Wj（X，Y）成岩阶段S_q的综合属性G j相比较，即为：
Gi-Gj=（F_ m_i，R_n_i，P_o_i）-（F_m_j+ R_n_j+P_o_j）；
如果满足G i-G j= 0，则未知网格Wi（X，Y）的溶蚀作用D_s_i与井j具有相同的溶蚀作用D_s_j；
如果G i-G j≠0，则按不同属性优先级顺序，即：沉积相属性F_m一级、岩性岩性属性R_n二级、流体属性P_o三级对未知网格Wi（X，Y）的溶蚀作用进行判识，即：
ⅰ：F_ m_i-F_ m_j = 0，R_n_i -R_n_j≠0，P_o_i - P_o_j≠0；
ⅱ：F_ m_i-F_ m_j = 0，R_n_i -R_n_j= 0，P_o_i - P_o_j≠0；
未知网格Wi（X，Y）的溶蚀作用D_s_i满足条件i时，未知网格Wi（X，Y）与具有相同属性F_j的网格Wj（X，Y）的溶蚀作用相同；
D_s_i 满足条件ⅱ时，未知网格Wi（X，Y）与具有相同属性F_ m_j、R_n_j 的网格Wj（X，Y）的溶蚀作用相同；
（7）根据每个网格Wi（X，Y）不同演化期次的溶蚀作用确定研究区溶蚀作用演化；
5）计算不同演化期次溶蚀量
（1）根据成岩序列数据，确定已知井所在网格Wj（X，Y）的成岩阶段S_q所对应的溶蚀作用下的溶蚀级别B_k_j，即Wj（S_q_j，D_s_j）= B_k_j；
（2）将Wi（X，Y）的Wi（S_q_j，D_s_j）与Wj（S_q_j，D_s_j）比对，满足Wi（S_q_j，D_s_j）=Wj（S_q_j，D_s_j），则Wi（X，Y）与Wj（S_q_j，D_s_j）具有相同的溶蚀级别B_k，即B_k_i=B_k_j；如果出现同一未知网格Wi（X，Y）与多个已知网格Wj（S_q_j，D_s_j）相同，则B_k_i=1/n∑B_k_j，n为相同网格Wj（S_q_j，D_s_j）的网格数；
（3）根据网格Wi（X，Y）不同演化期次的溶蚀级别B_k_i计算每个网格Wi（X，Y）不同演化期次的溶蚀率PCem；
（4）根据每个网格Wi（X，Y）的溶蚀率RCem计算溶蚀量VCem，VCem=∅0*RCem，其中∅0为原始孔隙；
（5）根据每个网格Wi（X，Y）不同演化期次的溶蚀量确定目标储层溶蚀量，即溶蚀增孔量。
2.根据权利要求1所述基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，其特征在于：所述步骤1）中，收集研究区地质资料由以下组成：研究区构造等值线图、研究区沉积相图、研究区岩性分布图、研究区不同时期流体分布图、研究区溶蚀作用数据。
3.根据权利要求1所述基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法，其特征在于：所述步骤2）中，
（1）沉积相类型数据库由以下组成：冲积扇-旱扇-扇根F_1、冲积扇-旱扇-扇中F_2、冲积扇-旱扇-扇缘F_3、冲积扇-湿扇-扇根F_4、冲积扇-湿扇-扇中F_5、冲积扇-湿扇-扇缘F_6、河流相-曲流河-河床亚相-河床滞留F_7、河流相-曲流河-河床亚相-边滩F_8、河流相-曲流河-堤岸亚相-天然堤F_9、河流相-曲流河-堤岸亚相-决口扇F_10、河流相-曲流河-河漫亚相-河漫滩F_11、河流相-曲流河-泛滥盆地-河漫湖泊F_12、河流相-曲流河-河漫沼泽F_13、河流相-辫状河-牛轭湖F_14、河流相-辫状河-河床滞留F_15、河流相-辫状河-心滩F_16、河道F_17、河流相-辫状河-泛滥平原F_18、湖泊相-断陷型-湖成三角洲F_19、湖泊相-断陷型-滨湖F_20、湖泊相-断陷型-浅湖F_21、湖泊相-断陷型-半深湖F_22、湖泊相-断陷型-深湖F_23、湖泊相-断陷型-湖湾F_24、湖泊相-坳陷型-湖成三角洲F_25、湖泊相-坳陷型-滨湖F_26、湖泊相-坳陷型-浅湖F_27、湖泊相-坳陷型-半深湖F_28、湖泊相-坳陷型-深湖F_29、湖泊相-坳陷型-湖湾F_30、湖泊相-前陆型-湖成三角洲F_31、湖泊相-前陆型-滨湖F_32、湖泊相-前陆型-浅湖F_33、湖泊相-前陆型-半深湖F_34、湖泊相-前陆型-深湖F_35、湖泊相-前陆型-湖湾F_36、三角洲相-辫状河三角洲-三角洲平原-分支F_37、三角洲相-辫状河三角洲-三角洲平原-河道F_38、三角洲相-辫状河三角洲-三角洲平原-陆上天然堤F_39、三角洲相-辫状河三角洲-三角洲平原-决口扇F_40、三角洲相-辫状河三角洲-三角洲平原-沼泽F_41、三角洲相-辫状河三角洲-三角洲平原-淡水湖泊F_42、三角洲相-辫状河三角洲-三角洲前缘-水下分支河道F_43、三角洲相-辫状河三角洲-三角洲前缘-水下天然堤F_44、三角洲相-辫状河三角洲-三角洲前缘-支流间湾F_45、三角洲相-辫状河三角洲-三角洲前缘-分支河口砂坝F_46、三角洲相-辫状河三角洲-三角洲前缘-远砂坝F_47、三角洲相-辫状河三角洲-三角洲前缘-前缘席状砂F_48、三角洲相-辫状河三角洲-前三角洲F_49、三角洲相-扇三角洲-扇三角洲平原-泥石流F_50、三角洲相-扇三角洲-扇三角洲平原-河道充填F_51、三角洲相-扇三角洲-扇三角洲平原-漫滩F_52、三角洲相-扇三角洲-扇三角洲前缘-碎屑流F_53、三角洲相-扇三角洲-扇三角洲前缘-水下分流河道F_54、三角洲相-扇三角洲-扇三角洲前缘-支流间湾F_55、三角洲相-扇三角洲-扇三角洲前缘-河口砂坝F_56、三角洲相-扇三角洲-扇三角洲前缘-远砂坝F_57、三角洲相-扇三角洲-前扇三角洲F_58；
（2）岩性数据库由以下组成：细砾岩R_1、中砾岩R_2、粗砾岩R_3、巨砾岩R_4、石英砂岩R_5、长石质石英砂岩R_6、岩屑质石英砂岩R_7、长石岩屑质石英砂岩R_8、长石砂岩R_9、岩屑质长石砂岩R_10、岩屑长石砂岩R_11、岩屑砂岩R_12、长石质岩屑砂岩R_13、长石岩屑砂岩R_14、黏土R_15、泥岩R_16、页岩R_17；
（3）流体数据库由以下组成：淡水环境P_1、酸性水环境P_2、碱性水环境P_3；
（4）成岩阶段数据库由以下组成：早成岩阶段A期ⅠA，即为：S_1；早成岩阶段B期ⅠB，即为：S_2；中成岩阶段A期ⅡA，即为：S_3；中成岩阶段B期ⅡB，即为S_4；晚成岩阶段Ⅲ，即为：S_5；
（5）溶蚀作用数据库中溶蚀作用类型由以下组成：
石英强溶蚀作用D1-1：其溶蚀率RCem＞10%；
石英中溶蚀作用D1-2：其溶蚀率5%＜RCem≤10%；
石英弱溶蚀作用D1-3：其溶蚀率RCem≤5%；
长石强溶蚀作用D2-1：其溶蚀率RCem＞10%；
长石中溶蚀作用D2-2：其溶蚀率5%＜RCem≤10%；
长石弱溶蚀作用D2-3：其溶蚀率RCem≤5%；
方解石强溶蚀作用D3-1：其溶蚀率RCem＞10%；
方解石中溶蚀作用D3-2：其溶蚀率5%＜RCem≤10%；
方解石弱溶蚀作用D3-3：其溶蚀率RCem≤5%；
白云石强溶蚀作用D4-1：其溶蚀率RCem＞10%；
白云石中溶蚀作用D4-2：其溶蚀率5%＜RCem≤10%；
白云石弱溶蚀作用D4-3：其溶蚀率RCem≤5%；
（6）溶蚀级别数学模型由以下几类组成：
B_1=Ⅰ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/10；
B_2=Ⅱ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/8；
B_3=Ⅲ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6] /6；
B_4=Ⅳ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6] /5；
B_5=Ⅴ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6] /3；
B_6=Ⅵ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/1.6；
B_7=Ⅶ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/1.5；
B_8=Ⅷ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/1.4；
B_9=Ⅸ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/1.2；
B_10=Ⅹ：Y =[-4.859e-06* (X)^4+0.001028* (X)^3-0.0745* (X)^2 +2.862* (X)+10.6]/1；
其中，X为埋藏时间，单位：Ma，Y为溶蚀率，单位：%。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810116730.6A CN108364096B (zh) | 2018-02-06 | 2018-02-06 | 基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810116730.6A CN108364096B (zh) | 2018-02-06 | 2018-02-06 | 基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108364096A CN108364096A (zh) | 2018-08-03 |
| CN108364096B true CN108364096B (zh) | 2021-10-15 |
Family
ID=63004594
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810116730.6A Active CN108364096B (zh) | 2018-02-06 | 2018-02-06 | 基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108364096B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111751903A (zh) * | 2019-03-29 | 2020-10-09 | 中国石油化工股份有限公司 | 一种碳酸盐岩气藏储层溶孔成因类型的薄片识别方法 |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102912000A (zh) * | 2012-11-07 | 2013-02-06 | 中国科学院地球化学研究所 | 一种定量微藻对灰岩的生物溶蚀作用方法 |
| CN102562048B (zh) * | 2010-12-30 | 2014-08-13 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN108345963A (zh) * | 2018-02-06 | 2018-07-31 | 长江大学 | 碎屑岩储层胶结减孔量的定量预测方法 |
-
2018
- 2018-02-06 CN CN201810116730.6A patent/CN108364096B/zh active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102562048B (zh) * | 2010-12-30 | 2014-08-13 | 长江大学 | 一种预测高成岩阶段低孔低渗碎屑岩有效储层的方法 |
| CN102912000A (zh) * | 2012-11-07 | 2013-02-06 | 中国科学院地球化学研究所 | 一种定量微藻对灰岩的生物溶蚀作用方法 |
| CN108345963A (zh) * | 2018-02-06 | 2018-07-31 | 长江大学 | 碎屑岩储层胶结减孔量的定量预测方法 |
Non-Patent Citations (1)
| Title |
|---|
| 海相碳酸盐岩优质储层形成机理与分布预测;何治亮等;《石油与天然气地址》;20110831;第32卷(第4期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108364096A (zh) | 2018-08-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bhattacharya et al. | Estimation of source area, river paleo-discharge, paleoslope, and sediment budgets of linked deep-time depositional systems and implications for hydrocarbon potential | |
| Friend et al. | Sedimentation of the Wood Bay Formation (Devonian) of Spitsbergen: regional analysis of a late orogenic basin | |
| Roy | Heavy mineral beach placers in southeastern Australia; their nature and genesis | |
| Prosser et al. | Holocene valley aggradation and gully erosion in headwater catchments, south‐eastern highlands of Australia | |
| Callot et al. | Paleo-fluids characterisation and fluid flow modelling along a regional transect in Northern United Arab Emirates (UAE) | |
| CN108344854B (zh) | 基于成岩过程的碎屑岩储层胶结减孔量的定量预测方法 | |
| Dou et al. | Petroleum geology of the Melut basin and the great Palogue field, Sudan | |
| Droste | High-resolution seismic stratigraphy of the Shu'aiba and Natih formations in the Sultanate of Oman: implications for Cretaceous epeiric carbonate platform systems | |
| Grauch et al. | New perspectives on the geometry of the Albuquerque Basin, Rio Grande rift, New Mexico: Insights from geophysical models of rift-fill thickness | |
| CN108376189B (zh) | 碎屑岩储层埋藏过程中成岩相演化的恢复方法 | |
| CN108364096B (zh) | 基于过程响应的碎屑岩储层溶蚀增孔量定量预测方法 | |
| Liu et al. | Geomorphological evolution and sediment dispersal processes in strike-slip and extensional composite basins: A case study in the Liaodong Bay Depression, Bohai Bay Basin, China | |
| Kadlec et al. | Cenozoic history of the Moravian Karst (northern segment): cave sediments and karst morphology | |
| CN108446788B (zh) | 基于成岩过程的碎屑岩储层孔隙度预测方法 | |
| CN108388708B (zh) | 碎屑岩储层溶蚀增孔量的定量预测方法 | |
| CN108345963B (zh) | 碎屑岩储层胶结减孔量的定量预测方法 | |
| Alqahtani et al. | Depositional environments of the Miocene-Pliocene siliciclastic succession, Al Rehaili area, north of Jeddah, Saudi Arabia | |
| CN113341480B (zh) | 分频rgb切片与多属性融合的砂质水合物储层预测方法 | |
| Ajakaiye et al. | Three-dimensional gravity interpretation of the dutsen-wai complex, nigerian younger granite province | |
| CN108388709A (zh) | 基于成岩相预测的储层孔隙度定量模拟方法 | |
| CN108469634B (zh) | 基于井的碎屑岩储层成岩相横向预测方法 | |
| Reman | Assessing Stratigraphic Controls on the Secondary Detrital Footprint from Buried Mineralization and Alteration at the Highland Valley Copper Mine, British Columbia | |
| CN108388960B (zh) | 压实减孔量的定量预测方法 | |
| Williams | Depositional history of the Black Mountain conglomerate, Mohave County, Arizona: Sedimentary response to Miocene extension | |
| Meunier | Reservoir characterization of the realgrunnen subgroup in Wisting Central III (7324/8-3), SW Barents Sea |
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 |