[0019] ②在保证模拟退火降温过程随机性的基础上，计算交换单位体像素点26个空间占位对邻域不稳定性的贡献程度S；当S>Sd时，认为该点的不稳定程度较高，可作为系统更新的交换点；当S
[0020] ③计算交换单位体像素点后系统的相关参数，包括单点概率函数、自相关函数、线性路径函数、分形函数和能量值，计算与未交换前系统的能量差值ΔE；当ΔE<0时，更新系统；当ΔE>0时，根据Metropolis准则来判断系统是否更新，即在一定的概率条件下接受系统更新；如果判断后不满足系统更新条件，则返回步骤②；
[0021] ④判断内循环终止条件，即判断在同一温度条件下系统能量差值是否小于设定最小能量差值；同时为避免系统刚降温，系统能量上升而立刻导致内循环结束而产生的降温，通过设定系统更新的失败率ff来避免该现象的出现，其中：
[0022]
[0023] 式中，Nf为导致系统能量回升的更新失败的次数；N为系统更新的总次数；
[0024] 当ff大于一定值后，则进行降温处理，降温过程采取等比降温方案，并返回步骤②；
[0025] ⑤当模拟过程温度降低到最终设定温度时或与上次降温的系统能量差值ΔE小于设定值时，整个模拟过程终止；
[0026] 作为约束条件，模拟退火算法中使用的统计函数包括：单点概率函数P(r)、自相关函数、线性路径函数和分形函数，利用自相关函数和线性路径函数对初始系统进行退火模拟，当模型具备一定分形特征后，引入分形函数进一步约束重建模型；
[0027] 第四步，将混合算法重建后初始数字岩心模型中的类球岩石颗粒，与过程法中构建的基础数字岩心模型的原始球形岩石颗粒相比较并取二者补集，将初始数字岩心模型初步划分为岩石骨架相、孔隙相和黏土矿物相三大类；
[0028] 第五步，通过Hoshen-Kopelman算法对初始数字岩心模型中的黏土矿物基团进行统计和划分，其中被M相占据的概率为c，被T相占据的概率为1-c，对于晶格中的每一个占位i，当其被M相占据时，则给该占位赋予一个基团标记 其中α是基团标记的特征符号，t为基团标记的次数，某一离散点的标记由一系列自然数表示：
[0029]
[0030] 在这一系列自然数中只有一个自然数是基团α的准确标记，该标记为 且该值是集合(3)中所有自然数的最小值，其它各基团标记之间的关系则由以下整数集给出：
[0031]
[0032] 其中，只有 是正整数元素，该值为基团中M相的个数，当进行第t次标记时，若基团中M相个数少于上次标记过程基团α的M相个数，则将该差值表示为相应t次的基团α的T相个数，(4)中的其它元素皆为负整数，反映了 与其它基团标记 的关系， 与 的关系用式(5)表示：
[0033]
[0034] 检查被判断离散点是否有被扫描过的相邻离散点，若相邻离散点为T相，则将当前被判断离散点赋予新基团的标记；如果有一个相邻离散点已经赋予基团标记，则将当前网格与相邻离散点赋予相同的标记；如果有一个以上的相邻离散点已经赋予基团标记，且基团标记各不相同，则将基团中所有离散点赋予相同的标记，最后统计并划分模型中黏土矿物相基团的个数及尺寸；
[0035] 第六步，较大的连通基团为黏土相中基团尺寸大于相邻基质颗粒尺寸的黏土矿物基团，通过K-means算法对初始数字岩心模型中黏土矿物相基团尺寸较大的黏土矿物基团进行划分，具体步骤如下：
[0036] ①读取数据样本的集合；
[0037] ②设定样本聚类的个数k，随机的选取k个数据样本作为初始的数据样本聚类中心；
[0038] ③计算欧氏距离，计算数据样本中每个数据到各聚类中心的欧式几何距离，然后根据最小误差平方和准则函数将数据按照远近距离划分到相应的不同聚类中心所对应的聚类当中；
[0039] ④更新聚类中心，将每个聚类中所有数据的均值作为各个聚类新的中心，并以最小误差平方和准则重新计算新的聚类中心的值；
[0040] ⑤
迭代判别，将步骤④中计算得到的数值与前一次计算得到的数值相比较，如果两者差值小于或等于预先设定的临界值，则停止迭代，否则重新进行步骤③进行迭代；
[0041] ⑥输出数据样本及聚类结果，包括每个聚类的聚类中心、大小；
[0042] 第七步，当黏土矿物基团边界的离散点为单个岩石颗粒时，则将该黏土矿物基团划分为交代形式，交代形式主要分布于岩石颗粒内，呈单个离散点的形式分布；当黏土矿物基团边界的相邻离散点为单个岩石骨架颗粒及孔隙时，则将该黏土矿物相基团划分为颗粒表面充填形式；
[0043] 当黏土矿物基团边界的相邻离散点为多个岩石骨架颗粒及孔隙时，则将该黏土矿物基团划分为粒间充填形式；
[0044] 将交代形式、颗粒表面充填形式和粒间充填形式的黏土矿物基团分别标记为A、B、C；最终得到不同结构黏土矿物基团分布和不同类型的黏土矿物基团分布；
[0045] 第八步，基于Hoshen-Kopelman算法和K-means算法得到初始数字岩心模型中黏土矿物基团大小及数量分布，以及按结构划分得到的黏土矿物基团类型及数量分布，结合真实储层黏土含量及分布以及主要的黏土矿物结构特点，按黏土矿物相基团大小和结构特点将模型中的黏土矿物赋予相应的黏土性质，得到含多组分岩石矿物分布的数字岩心模型。
[0046] 所述步骤2中，室内实验包括为模拟原始储层条件及不同生产条件下，当储层温度、压力、层内流体和注入流体性质发生变化时，不同岩石矿物的转化及体积变化情况，具体包括不同类型岩石矿物在不同模拟条件下(不同温度，不同润湿环境条件，不同水类型)的
膨胀率、溶蚀率和转化率。
[0047] 所述步骤3中，对不同模拟条件(不同温度，不同润湿环境，不同水类型)下岩石矿物的膨胀过程进行模拟的具体步骤如下：
[0048] ①基于步骤2中室内实验得到的膨胀率以及步骤1中所构建含多种岩石矿物组分数字岩心模型的
分辨率，计算不同类型岩石矿物在不同模拟条件下的膨胀体积；
[0049] ②读取步骤1中得到的含多种岩石矿物组分数字岩心模型(中不同岩石矿物基团的大小；
[0050] ③判断孔隙空间占位是否能用于放置膨胀后的岩石矿物；
[0051] ④计算空间占位的稳定性；
[0052] ⑤计算空间占位与岩石矿物基团表面的最小距离；
[0053] ⑥将可用于放置膨胀岩石矿物的空间占位进行排序；
[0054] ⑦将膨胀岩石矿物放置在优先等级较高的空间占位上。
[0055] 所述步骤4中，对不同模拟条件(不同温度，不同润湿环境，不同水类型)下岩石矿物的溶蚀过程进行模拟的具体步骤如下：
[0056] ①基于室内实验得到的溶蚀率以及步骤1中含多种岩石矿物组分数字岩心模型的分辨率，计算不同类型岩石矿物在不同模拟条件下的溶蚀体积；
[0057] ②读取步骤1得到的含多种岩石矿物组分数字岩心模型中不同岩石矿物基团的大小；
[0058] ③判断岩石矿物基团边界能否用于溶蚀；
[0059] ④计算边界占位的稳定性；
[0060] ⑤将可用于发生溶蚀的边界占位进行排序；
[0061] ⑥按照各岩石矿物基团的大小选取可用于溶蚀的边界占位；
[0062] ⑦在优先等级较高的边界占位上进行溶蚀。
[0063] 所述步骤5中，对不同模拟条件(不同温度，不同润湿环境，不同水类型)下岩石矿物的转化过程进行模拟的具体步骤如下：
[0064] ①基于室内实验得到的转化率以及步骤1中含多种岩石矿物组分数字岩心模型的分辨率，计算不同类型岩石矿物在不同模拟条件下的转化体积与被转化体积；
[0065] ②读取步骤1得到的含多种岩石矿物组分数字岩心模型中不同岩石矿物基团的大小；
[0066] ③判断转化相岩石矿物的占位是否能够用于放置被转化相岩石矿物；
[0067] ④计算空间占位的稳定性；
[0068] ⑤将可用于转化的岩石矿物空间占位按所属基团的大小进行排序；
[0069] ⑥选择优先等级较高的转化岩石矿物放置被转化相岩石矿物；
[0070] ⑦将转化岩石矿物放置在优先等级较高的空间占位上。
[0071] 通过对比不同储层敏感性伤害前后数字岩心模型中岩石矿物的体积变化、数字岩心模型的孔隙度和渗透率的变化，进一步研究储层敏感性伤害对储层微观结构的影响。
[0072] 与现有技术相比，本发明的一种基于数字岩心模型的储层敏感性伤害模拟方法，通过将有限的矿场资料与室内实验结果相结合，提出了一种基于数字岩心技术的储层敏感性伤害模拟研究的新手段，该方法首先利用有限的矿场资料构建了包含多种岩石矿物类型与产状的数字岩心模型；其次，在不同模拟条件(温度、压力和流体性质等)下研究了不同类型岩石矿物的转化及体积变化规律，以此为基础，最后通过基于模型离散点稳定性判别和形态学中膨胀和腐蚀算法等方法将室内研究结果与数字岩心紧密结合，实现了基于数字岩心技术对储层敏感性伤害的模拟，最后通过对储层伤害前后数字岩心模型微观结构及孔渗变化的分析，得到不同模拟条件下，岩石矿物敏感性对储层伤害的程度，该方法的提出进一步拓展了数字岩心技术在油气田开发领域的应用，也为岩石矿物敏感性伤害的研究提供了新手段。
附图说明
[0073] 图1是含多种岩石矿物数字岩心模型重建
流程图。
[0074] 图2是岩石矿物膨胀运算流程图。
[0075] 图3是岩石矿物腐蚀运算流程图。
[0076] 图4是岩石矿物转化运算流程图。
[0077] 图5是含多种岩石矿物数字岩心模型。
[0078] 图6是含多种岩石矿物数字岩心模型中黏土矿物基团的分布。
[0079] 图7是不同类型的黏土矿物模型及其分布。
[0080] 图8是高温高pH模拟冷凝液条件下的储层敏感性伤害模型。
[0081] 图9是黏土膨胀溶蚀前后特征。
[0082] 图10是高温高pH模拟冷凝液条件下各类黏土基团大小的分布。
[0083] 图11是高温高pH模拟冷凝液条件下黏土基团的分布。
具体实施方式
[0084] 下面结合附图和
实施例对本发明作详细说明，但本发明不限于下列的实施例。
[0085] 本发明具体涉及到一种基于数字岩心模型的储层敏感性伤害模拟方法，该实施例的高温高pH模拟冷凝液条件下储层敏感性伤害的模拟过程是通过以下步骤来实现的：
[0086] 实施例中所采用的是基于储层岩石二维信息的含多种岩石矿物的数字岩心重建技术，其中包含的储层岩石二维信息主要包括储层的粒度分布、铸体资料、孔隙度、岩石矿物含量及产状分布等。
[0087] (1)含多种岩石矿物原始数字岩心模型的构建
[0088] 含多种岩石矿物的数字岩心模型的构建按图1所示的过程进行构建，其中利用过程法构建基础数字岩心模型，在满足高能环境和重力势能梯度最大的下落模拟原则的基础上模拟沉积过程，并结合真实岩心孔隙度，选择压实因子控制数字岩心的孔隙度。为构建含多组分岩石矿物的数字岩心模型，在沉积过程中，根据真实粒度随机选择球体颗粒的半径时，考虑了其他类型岩石矿物所占体积，因此沉积颗粒的尺寸不但由原始的粒度分布决定，同时额外考虑了其他类型岩石矿物与储层砂岩颗粒之间的比例。在运用混合算法构建初始数字岩心模型时，将单位体像素点的空间占位(点、线和面)按其对邻域不稳定性的贡献程度赋予权值，其中面为5，边为3，点为2；在选取交换体像素点时，计算该体像素点与邻域占位点、线和面上的不稳定性贡献程度S，并基于模拟退火算法中能量值下降的过程，引入交换单位体像素点对其邻域不稳定性的贡献程度参数Sd，对交换点的可交换性进行判断，提高交换单位体像素点的有效性，其中Sd为与模拟过程中系统能量相关的无因次值：
[0089] Sd＝N×β(E0-Ei/△Emax)(1)
[0090] 式中，N为单位体像素点影响的邻域接触点的个数，无量纲；β为单位体像素点对邻域不稳定性系数，无量纲；E0为系统的初始能量，无量纲；Ei为第i次降温后系统的能量，无量纲；ΔEmax为初始模型和参考模型系统的能量差值，无量纲。
[0091] 作为约束条件，模拟退火算法中常用的统计函数包括：单点概率函数、自相关函数、线性路径函数和分形函数等，利用自相关函数和线性路径函数对初始系统进行退火模拟，当模型具备一定分形特征后，引入分形函数进一步约束重建模型。
[0092] 将混合算法重建后初始数字岩心模型中的类球岩石颗粒，与过程法中重建基础模型的原始球形岩石颗粒相比较，并将模型初步划分为岩石骨架相(R)、孔隙相(P)和黏土矿物相(C)三大类，其中C相在空间中以大小不同的不规则离散基团的形式分布，其中模型运算之前需要将
三维重建模型划分为两相，孔隙相和岩石骨架相需合并为一相T，利用Hoshen-Kopelman算法得到模型中黏土基团的编号、尺寸和数量，由于统计划分的黏土矿物基团中存在部分尺寸较大的连通黏土基团，而在真实储层中不同种类的黏土矿物在岩石颗粒表面也存在连通、接触的情况；利用K-means聚类算法可以将对岩石颗粒(聚类中心)周围的黏土矿物按所属关系进行聚类。因此，本发明以岩石颗粒的球心和所有较大尺寸的黏土矿物基团作为K-means算法的数据样本对Hoshen-Kopelman算法划分后的黏土矿物基团进行有效划分。对于较大尺寸的连通黏土矿物基团可按照K-means算法划分为多个附着于岩石颗粒表面的有效黏土基团。
[0093] 储层中常见的黏土矿物包括蒙脱石、伊利石、伊蒙混层、绿泥石、
高岭石；常见的分布形式为粒间孔隙充填、颗粒包壳、交代和包壳衬边等，且不同黏土的分布特点也各不相同。重建模型中黏土矿物的填充形式主要为颗粒表面填充(单个黏土表面填充，多个黏土表面填充和层状黏土表面填充)、颗粒间填充(双颗粒间黏土填充、多颗粒间黏土填充)和颗粒内部填充，故在进行数字岩心模型的黏土矿物
构建时，结合实际黏土的分布形式，按照单个黏土基团与岩石骨架颗粒的相邻关系将黏土矿物基团分布的主要形式划分为：粒间充填、颗粒表面充填和交代作用。
[0094] 基于Hoshen-Kopelman算法得到模型中黏土矿物基团大小及数量分布，以及按结构划分得到的重建模型中黏土基团类型及数量分布，结合真实储层黏土含量分布以及主要的黏土矿物结构特点，按黏土矿物基团大小和结构特点将模型中的黏土矿物赋予相应的黏土性质，得到含不同类型岩石矿物分布的数字岩心模型，如图5所示。
[0095] (2)数字岩心模型岩石矿物分布特征
[0096] 参考模型储层孔隙度为26.38％，渗透率0.614μm2，泥质含量12.36％；其中黏土含量分布为：蒙脱石40.8％，高岭石19.1％，绿泥石27.4％，伊利石6.3％。其中蒙脱石产状主要以颗粒包壳为主，存在部分粒间充填形式；高岭石以粒间孔隙充填，呈分散质点状集合体分布；绿泥石以包壳衬边，粒间充填和交代状分布；伊利石的分布形式包括粒间充填、交代和
薄膜式分布。
[0097] ①重建模型中的黏土基团分布
[0098] 基于Hoshen-Kopelman算法得到的数字岩心模型中不同黏土矿物基团的分布情况如图6所示，其中最大的黏土矿物基团的尺寸为27953个
体素，最小的黏土基团大小为1个体素(基团个数为9432)。基团大小小于11个体素的黏土基团仅占所有黏土基团的1.91％；而主要的黏土基团则分布在10000个体素到25000个体素之间，占总黏土体积的97.29％。整体的黏土基团分布呈现“大基团为主，小基团分散”的特点，这与实际储层中黏土矿物的分布形式相近似。
[0099] ②结构划分后模型中黏土矿物的统计
[0100] 通过黏土基团的结构判别，所有的黏土矿物基团按产状被划分为三种主要类型：表面充填，粒间充填和交代作用。其中以粒间充填形式分布的黏土矿物基团共有4685个基团，含量占黏土总体积的67.13％；表面充填产状的黏土矿物基团共有4530个，含量占黏土总体积的32.30％；而交代作用的黏土矿物则零星的分布于岩石颗粒当中，其含量仅占到黏土总体积的0.28％；模型中较大的黏土矿物基团主要为表面充填和粒间充填的形式。
[0101] 通过对数字岩心模型中黏土矿物基团的划分和结构判别，每个黏土矿物基团都标记了不同的属性(包括基团大小、序号、产状等)。因此，结合真实储层的相关信息(包括黏土含量、黏土类型、黏土的产状等)，模型中的黏土矿物按黏土含量和产状特征被划分为不同的黏土类型，蒙脱石是模型中含量最多的黏土矿物，含量占黏土矿物总体积的40.84％，绿泥石占27.43％，高岭石占19.11％，伊利石占6.28％。且对于不同尺寸的黏土矿物基团，各种黏土矿物基团的分布相对均匀。
[0102] ③含多种岩石矿物数字岩心模型中黏土矿物的分布
[0103] 从各层中黏土矿物的分布情况来看，黏土矿物分布中包括部分基团大小小于5的黏土矿物颗粒，同时大基团黏土矿物的分布以粒间充填(双颗粒间和多粒间)，颗粒表面充填(蚀变类黏土、包壳衬边、薄膜式)为主，存在少量交代式分布的黏土矿物。
[0104] 由图7a和图7e可以看出，蒙脱石在模型中由于含量较高，主要成连片充填并附着于岩石基质表面，同时蒙脱石黏土基团主要以粒间充填和表面充填的形式分布于重建模型当中，基团数量分别为2117个和1935个；粒间充填类蒙脱石和表面充填类蒙脱石的含量占黏土矿物总体积的41.41％和58.39％；其中最大的粒间充填类蒙脱石基团大小为22716个体素，最大的表面充填类蒙脱石基团大小为21273个体素；由图7b和7f可以看出，绿泥石以环状和部分连片的基团分布于模型当中，粒间充填类绿泥石共有900个基团，占黏土矿物总体积的62,53％；表面充填类绿泥石共有975个基团，占黏土矿物总体积的37.14％，最大的粒间充填类和表面充填类绿泥石基团大小分别为22767个体素和21193个体素；由图7c和图7g可以看出，高岭石一般以粒间充填的形式分布于储层岩石当中，而由模型中黏土矿物的分布可知，粒间充填是模型中高岭石的主要分布形式，占黏土总体积的98.58％，其中最大的粒间充填类黏土基团的大小为27953个体素；由图7d和7h可以看出，伊利石在模型中的产状包括粒间充填，表面充填和交代形式，其中表面充填类和粒间充填类伊利石分别占
41.32％和58.12％。交代作用在四种黏土矿物中均有分布，且主要以零星分布的形式分布于岩石颗粒当中，蒙脱石、绿泥石、高岭石和伊利石中交代状黏土基团的个数分别为504、
619、61和244个，所构建的含黏土三维多孔介质模型与真实储层的黏土矿物分布、产状较为吻合。
[0105] (3)储层岩石矿物性质变化造成的储层伤害模型构建
[0106] 储层中的黏土矿物在不同的储层环境中，经过复杂的水岩反应，容易造成黏土矿物的转化，由室内模拟实验结果可知：高温高pH模拟冷凝液反应后，1％的蒙脱石向其它类型岩石矿物转化，6％的高岭石向蒙脱石转化，1％的高岭石向伊利石转化，7％的伊利石向蒙脱石转化，16％的绿泥石向其它类型岩石矿物转化。
[0107] 蒙脱石及一些混层黏土矿物是常见的水敏性黏土矿物，在不同的环境中，水敏性黏土矿物的膨胀体积和对储层造成的伤害程度亦不尽相同。针对不同模拟条件下黏土矿物的体积膨胀变化，开展了相应的实验研究，在高温高pH的模拟冷凝液环境中，体积膨胀较为明显的黏土矿物为蒙脱石，膨胀率为159.13％；非膨胀性的伊利石向膨胀性的蒙脱石发生了部分转化造成了体积膨胀，膨胀率为58.8％；模型中蒙脱石的基团表面发生了明显的膨胀；随着蒙脱石的膨胀，模型的孔隙度不断下降，同时模型的孔喉结构也发生了较为明显的变化，进一步加剧了水敏性黏土矿物对储层渗透性能带来的影响；绿泥石是典型的酸敏性黏土矿物，然而在高温高pH的模拟冷凝液环境中绿泥石在水岩反应作用下，也发生了较少量的溶蚀作用。
[0108] 储层岩石矿物性质变化造成的储层伤害，包括了因岩石矿物性质变化造成岩石矿物的体积膨胀、溶蚀、转化等一系列过程，不同模拟条件下黏土矿物的体积变化和各类黏土矿物的转化过程是在岩石的水岩反应过程中同时发生的，且黏土转化过程并不改变黏土矿物的体积，而黏土矿物的体积变化均反映了不同类型黏土矿物转化后的体积变化。因此由于储层岩石矿物性质造成的储层伤害模型应首先综合考虑岩石矿物的转化和体积变化，在此基础上，基于黏土矿物的转化关系，对体积变化后的储层伤害模型做黏土转化过程的模拟，具体地，高温高pH模拟冷凝液条件下岩石矿物的膨胀过程进行模拟的具体步骤按图2所示的流程进行模拟：
[0109] ①基于室内实验得到的岩石矿物膨胀率以及所构建含多种岩石矿物组分数字岩心模型的分辨率，计算不同类型岩石矿物在高温高pH模拟冷凝液条件下的膨胀体积；
[0110] ②读取含多种岩石矿物组分数字岩心模型中不同岩石矿物基团的大小；
[0111] ③判断孔隙空间占位是否能用于放置膨胀后的岩石矿物；
[0112] ④计算空间占位的稳定性；
[0113] ⑤计算空间占位与岩石矿物基团表面的最小距离；
[0114] ⑥将可用于放置膨胀岩石矿物的空间占位进行排序；
[0115] ⑦将膨胀岩石矿物放置在优先等级较高的空间占位上。
[0116] 高温高pH模拟冷凝液条件下岩石矿物的溶蚀过程进行模拟的具体步骤按图3所示的流程进行模拟：
[0117] ①基于室内实验得到的岩石矿物溶蚀率以及所构建含多种岩石矿物组分数字岩心模型的分辨率，计算不同类型岩石矿物在高温高pH模拟冷凝液条件下的溶蚀体积；
[0118] ②读取含多种岩石矿物组分数字岩心模型中不同岩石矿物基团的大小；
[0119] ③判断岩石矿物基团边界能否用于溶蚀；
[0120] ④计算边界占位的稳定性；
[0121] ⑤将可用于发生溶蚀的边界占位进行排序；
[0122] ⑥按照各岩石矿物基团的大小选取可用于溶蚀的边界占位；
[0123] ⑦在优先等级较高的边界占位上进行溶蚀。
[0124] 高温高pH模拟冷凝液条件下岩石矿物的转化过程进行模拟的具体步骤按图4所示的流程进行模拟：
[0125] ①基于室内实验得到的岩石矿物转化率以及所构建含多种岩石矿物组分数字岩心模型的分辨率，计算不同类型岩石矿物在高温高pH模拟冷凝液条件下的转化体积与被转化体积；
[0126] ②读取含多种岩石矿物组分数字岩心模型中不同岩石矿物基团的大小；
[0127] ③判断转化相岩石矿物的占位是否能够用于放置被转化相岩石矿物；
[0128] ④计算空间占位的稳定性；
[0129] ⑤将可用于转化的岩石矿物空间占位按所属基团的大小进行排序；
[0130] ⑥选择优先等级较高的转化岩石矿物放置被转化相岩石矿物；
[0131] ⑦将转化岩石矿物放置在优先等级较高的空间占位上。
[0132] (4)储层岩石矿物性质变化造成的储层伤害特征研究
[0133] 利用含多组分数字岩心模型构建了由于储层岩石矿物性质变化造成的储层伤害模型，由于在不同模拟条件下不同岩石矿物的转化及体积变化会造成黏土矿物产状、含量的变化，同时黏土矿物的变化会进一步造成储层岩石物性的变化。因此，本例研究了由于高温高pH模拟冷凝液条件下，储层岩石矿物性质变化造成的黏土矿物含量及体积变化，以及模型的孔隙度及渗透率的变化。
[0134] ①黏土矿物含量与体积的变化
[0135] 如图8所示为高温高pH模拟冷凝液条件下的储层敏感性伤害模型，在黏土转化过程中，虽然总的黏土矿物含量及体积并没有发生变化，但是不同类型的黏土矿物之间发生了相互转化，各类黏土矿物的含量和体积都发生了变化；而在黏土矿物的溶蚀和膨胀过程中，由于体积的变化是从黏土矿物的表面开始的，因此在变化过程中，黏土矿物的含量及体积，甚至黏土矿物的产状都发生了明显的变化，如图9所示，其中图9a和图9c分别为黏土矿物膨胀前后的示意图，图9b和图9d分别为黏土矿物溶蚀前后的示意图。
[0136] 表1高温高pH模拟冷凝液条件下黏土矿物变化
[0137]
[0138]
[0139] 如表1所示为高温高pH模拟冷凝液条件下黏土矿物的变化情况。模型中的黏土含量发生了明显的变化，其中模型的黏土含量由原始的12.29％上升为高温高pH模拟冷凝液条件下19.00％，黏土矿物含量的变化主要是由于模拟反应条件下黏土矿物的膨胀所导致的，其中蒙脱石由原始模型的40.84％上升为高温高pH模拟冷凝液反应后的66.06％。岩石矿物经历转化、溶蚀、膨胀后，相比原始模型，高温高pH模拟冷凝液模型中，蒙脱石和伊利石的膨胀占主导作用，而绿泥石和高岭石向蒙脱石、伊利石和其它矿物的转化作用以及二者的溶蚀作用使得模拟反应后绿泥石和高岭石的含量相比原始模型有所降低。
[0140] 由图10和图11所示，图10a、10b、10c和10d分别为高温高pH模拟冷凝液储层伤害数字岩心模型和原始含多种岩石矿物数字岩心模型中蒙脱石、高岭石、伊利石和绿泥石的基团大小；图11a和图11e为高温高pH模拟冷凝液伤害前后蒙脱石的分布，图11b和图11f为高温高pH模拟冷凝液伤害前后高岭石的分布，图11c和图11g为高温高pH模拟冷凝液伤害前后伊利石的分布，图11d和图11h为高温高pH模拟冷凝液伤害前后绿泥石的分布，不同类型黏土矿物基团在储层伤害过程前后，基团的大小及基团的数量都发生了明显的变化，从黏土的结构上来看，由于模拟过程中交代类黏土矿物均位于其它类黏土矿物中，因此在储层伤害过程中，交代类黏土矿物的含量基本不发生变化，对于膨胀性黏土矿物蒙脱石，发生敏感性伤害后，表面充填和粒间充填蒙脱石的基团大小都有所上升，且由图10可以看出，在膨胀类黏土矿物的膨胀过程中，部分表面充填的黏土矿物随着体积的不断增加，附着于岩石颗粒表面的黏土不断向孔隙空间及邻近的岩石颗粒表面膨胀，因此膨胀过程中表面类黏土的增幅低于粒间类黏土的增幅。绿泥石由于在不同的模拟反应过程中都会发生少量的溶蚀，因此不同类型的绿泥石基团大小都有所下降。
[0141] ②孔隙度和渗透率变化
[0142] 表2高温高pH模拟冷凝液条件下模型孔渗参数变化
[0143]
[0144] 由表2可以看出，储层伤害后黏土矿物含量上升而孔隙度下降，高温高pH模拟冷凝液条件下岩石矿物性质造成的储层伤害明显，黏土含量由原始模型的12.29％上升为模拟反应后的19.00％，增幅为50.08％，而孔隙度由原始的26.38％下降为20.56％；同时渗透率由原始模型的589.76×103μm2下降为279.20×103μm2，由于岩石矿物性质的差异对原始的储层造成了极大的伤害，储层孔渗下降明显。
