黏土膨胀损害油气层的建模方法、损害程度时空演化4D定量与智能诊断方法及其系统

黏土膨胀损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
技术领域
1.本发明涉及油田勘探技术领域，具体地涉及一种黏土膨胀损害储层的建模方法与系统以及确定储层损害程度的方法与系统。


背景技术：

2.油田勘探开发的各个时期，由于受到多种内外因素影响，会导致储层原有的物理、化学、热力学和水动力学平衡状态变化，不可避免的使储层近井壁区乃至远井壁区的储层内部渗透率降低，堵塞流体流动，造成储层损害和油井产量下降，甚至“枪毙”储层。造成储层损害的原因是多样的、复杂的，特别是在生产过程中，储层岩石储渗空间、表面润湿性、水动力学场、温度场、岩石种类等不断发生变化，使损害机理随时间而变，且损害周期长、范围宽，损害更具复杂性和叠加性。储层损害一旦发生，必须根据储层损害情况采取相应的解堵措施恢复流体流动通道，以便提高油井产量和水井注入能力。因此，厘清待解堵井储层损害究竟由哪些因素造成、各损害因素所占比例如何，以及储层损害的空间分布规律和随时间变化规律对解堵措施优化设计至关重要，并直接影响解堵和增产效果好坏。
3.目前，诊断储层损害的方法可分为矿场诊断法和室内评价法。其中，所述矿场诊断法包括试井法。虽然所述试井法可以定量给出表征待诊断井的预设区域内的储层的损害程度的表皮系数、堵塞比、附加压降等重要参数，但由其表征的表皮系数与其它参数相互联系。也就是说，通过所述试井法得出的表皮系数并不仅仅反映真实储层损害特征，还是各个环节、多因素的综合表现(即其是真实损害表皮系数和由井斜表皮系数、储层形状表皮系数、打开储层不完善表皮系数、非达西流表皮系数、射孔表皮系数等组成的拟表皮系数之和)，必须进行表皮系数分解才能得到真实损害表皮系数。其中，所述室内评价法包括岩心流动实验法。所述岩心流动实验法是通过岩心驱替前后的渗透率变化来了解损害程度大小，虽然比较适合研究单因素储层损害，但难以反映较大尺度上储层损害规律。并且，因室内岩心实验条件比较理想化、评价用岩心都是原始状态岩心、无法考虑储层特性动态变化，使实验结果与井下储层真实损害出入较大。


技术实现要素：

4.本发明的目的是提供一种黏土膨胀损害储层的建模方法与系统以及确定储层损害程度的方法与系统，其可定量模拟由黏土膨胀引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
5.为了实现上述目的，本发明第一方面提供一种黏土膨胀损害储层的建模方法，所述建模方法包括：确定待诊断井的预设区域内的储层中的流体的达西表观速度；根据所述
流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程；根据菲克扩散定律，建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程；以及根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程，其中所述黏土为所述岩石的组成成分。
6.优选地，所述确定待诊断井的预设区域内的储层中的流体的达西表观速度包括：建立所述流体进入所述储层的压力传导方程；以及根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
7.优选地，所述建立所述流体中的水分子的质量平衡方程包括：根据所述流体的达西表观速度u与所述水分子的扩散系数d
w
，建立下式表示的所述质量平衡方程：其中，φ0为所述储层的孔隙度的初始值；为所述储层内的孔隙的含水体积分数；以及为所述储层内的任意点的空间位置。
8.优选地，所述流体中的水分子的质量平衡方程的初始条件为以及所述流体中的水分子的质量平衡方程的边界条件为其中，φ0为所述储层的孔隙度的初始值；r
w
为所述待诊断井的井筒半径；以及s
wc
为所述储层中的束缚水饱和度。
9.优选地，所述确定黏土膨胀损害储层的时空演化模拟方程包括：根据所述质量平衡方程及该质量平衡方程的边界条件与初始条件，确定所述储层内的孔隙的含水体积分数；根据所述扩散方程、所述扩散方程的边界条件、所述扩散方程的初始条件及所述储层内的孔隙的含水体积分数，确定所述储层中的岩石的吸水速率；以及根据所述储层中的岩石的吸水速率，确定所述黏土膨胀损害储层的时空演化模拟方程。
10.优选地，所述建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程包括：根据所述菲克扩散定律，建立下式表示的所述流体中的水分子向储层内部扩散的扩散方程：其中，n为以储层中的位置为原点且以所述流体与所述岩石的界面在位置处的法线方向为坐标轴的方向建立的一维坐标系中的坐标；t为时间；d
w
为所述水分子的扩散系数；以及c(n,t)为所述储层中的岩石的含水体积分数。
11.优选地，所述确定所述储层中的岩石的吸水速率包括：根据所述扩散方程、所述扩散方程的边界条件、所述扩散方程的初始条件及所述储层内的孔隙的含水体积分数，确定下式表示的所述储层的吸水速率
12.其中，所述扩散方
程的初始条件为c(z,t＝0)＝c0；所述扩散方程的边界条件为；所述扩散方程的边界条件为与d
w
为所述水分子的扩散系数；c(z,t)为所述储层中的岩石的含水体积分数；以及k
f
为膜交换系数。
13.优选地，所述确定所述黏土膨胀损害储层的时空演化模拟方程包括：根据所述储层的吸水速率确定由下式表示的所述黏土膨胀损害储层的时空演化模拟方程：其中，为所述储层的孔隙度；以及λ为黏土膨胀系数。
14.通过上述技术方案，本发明创造性地确定待诊断井的预设区域内的储层中的流体的达西表观速度；根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程；建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程；以及根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程。由此，通过所确定的时空演化模拟方程可定量模拟由黏土膨胀引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
15.本发明第二方面提供一种确定储层损害程度的方法，所述方法包括：基于根据所述的黏土膨胀损害储层的建模方法建立的时空演化模拟方程，确定所述储层的孔隙度；以及基于所确定的所述储层的孔隙度，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
16.优选地，所述特征参数为所述储层的渗透率和/或所述储层的滤失系数，相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：基于所述储层的孔隙度及公式确定所述储层的渗透率和/或基于所述储层的孔隙度及公式确定所述储层的滤失系数其中，φ0为所述储层的孔隙度的初始值；为所述储层的孔隙度；φ
dmax
为所述储层的最大孔隙度；m
k
与m
k
分别为第一经验值与第二经验值；为所述储层的渗透率的初始值；以及所述储层的滤失系数的初始值。
17.优选地，所述特征参数为所述储层的表皮系数，相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：基于所述储层的孔隙度及公式确定所述储层的渗透率以及基于所述储层的渗透率
及公式确定所述储层的表皮系数其中，为所述储层的渗透率的初始值，r
w
为所述待诊断井的井筒半径，以及r
sw
为所述储层的损害半径。
18.通过上述技术方案，本发明创造性地通过所确定的时空演化模拟方程可确定所述储层的孔隙度，再根据所述储层的孔隙度可确定表征待诊断井的预设区域内的储层的损害程度的特征参数(例如所述储层的渗透率和/或表皮系数)，由此可定量模拟由黏土膨胀引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
19.相应地，本发明第三方面还提供一种黏土膨胀损害储层的建模系统，所述建模系统包括：速度确定装置，用于确定待诊断井的预设区域内的储层中的流体的达西表观速度；第一建立装置，用于根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程，第二建立装置，用于根据菲克扩散定律，建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程；以及模拟方程确定装置，用于根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程，其中所述黏土为所述岩石的组成成分。
20.所述黏土膨胀损害储层的建模系统与上述黏土膨胀损害储层的建模方法相对于现有技术所具有的优势相同，在此不再赘述。
21.相应地，本发明第四方面还提供一种确定储层损害程度的系统，所述系统包括：孔隙度确定装置，用于基于所述的黏土膨胀损害储层的建模系统建立的时空演化模拟方程，确定所述储层的孔隙度；以及特征参数确定装置，用于基于所确定的所述储层的孔隙度，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
22.所述确定储层损害程度的系统与上述确定储层损害程度的方法相对于现有技术所具有的优势相同，在此不再赘述。
23.相应地，本发明第五方面还提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行所述的黏土膨胀损害储层的建模方法和/或所述的确定储层损害程度的方法。
24.本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
25.附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
26.图1是本发明一实施例提供的黏土膨胀损害储层的建模方法的流程图；
27.图2是本发明一实施例提供的储层中的孔隙内的水分子向岩石内部扩散的示意
图；
28.图3是本发明一实施例提供的确定黏土膨胀损害储层的时空演化模拟方程的流程图；
29.图4是本发明一实施例提供的确定储层损害程度的方法的流程图；
30.图5是本发明一实施例提供的吸水速率随时间演化的示意图；
31.图6是本发明一实施例提供的表皮系数随时间演化的示意图；
32.图7是本发明一实施例提供的由储层渗透率损害率表征的在第500天黏土膨胀损害储层的半径的示意图；
33.图8是本发明一实施例提供的黏土膨胀损害储层的建模系统的结构图；以及
34.图9是本发明一实施例提供的确定储层损害程度的系统的结构图。
具体实施方式
35.以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
36.外来流体(例如注入水)的水分子通过固
‑
液(储层中的岩石与流体之间的)界面向固相介质(储层中的岩石)内部的扩散过程可认为是固
‑
液界面的任意局部区域的情况，对于该界面上的任意足够小尺度的局域区域而言，可认为水分子的扩散方向垂直于该区域的某点(例如o点)的切线方向(即扩散方向与该区域所在的平面相垂直)，如图2(其中阴影部分表示岩石，其他空白部分表示储层中的孔隙)所示。所述储层中的岩石包括黏土，在水分子向岩石扩散的过程中，所述黏土会发生膨胀，进而可能会导致储层的渗透率降低(甚至堵塞)。因此，本发明各实施例的核心是要建立水分子向岩石内部扩散和储层中的孔隙内的水含量变化的动力学模型(即水分子通过固
‑
液界面由孔隙中的液相向固相内部扩散的扩散方程与孔隙中的流体的对流扩散方程)。具体地，基于菲克扩散定律以及储层中的孔隙内的流体的对流扩散关系等建立黏土膨胀影响待诊断井周围的储层中的孔隙度分布的时空演化控制唯象模型(该模型包含储层内的孔隙的含水体积分数c1与储层中的岩石的含水体积分数的初始值c0)，再结合储层的孔隙度和渗透率等储层损害特征参数间的关系，就可诊断渗透率等储层损害特征参数的时空场分布。
37.需要说明的是，为了简单描述起见，在本发明的各个实施例中的随时空演化的物理量、化学量可省略变量例如可简写为c1；以及可简写为
38.图1是本发明一实施例提供的黏土膨胀损害储层的建模方法的流程图。所述建模方法可包括步骤s101
‑
s104。
39.步骤s101，确定待诊断井的预设区域内的储层中的流体的达西表观速度。
40.其中，所述待诊断井可例如为注水井或采油井。
41.对于步骤s101，所述确定储层中的流体的速度可包括：建立所述流体进入所述储层的压力传导方程；以及根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
42.具体地，压力是驱动固
‑
液混合液从注水井的井筒持续侵入周围的储层的动力，由此可建立如公式(1)的所述流体进入储层的压力传导方程：
[0043][0044]
再根据公式(1)及达西公式(如下式(2))可确定所述流体的达西表观速度，
[0045][0046]
其中，为所述流体的压力；μ为流体粘度；c
t
为流体
‑
岩石综合压缩系数以及为所述储层的渗透率。
[0047]
步骤s102，根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程。
[0048]
在储层条件下，所述储层中的孔隙内的不同位置的含水量满足质量守恒方程。其中，所述储层内的含水量变化主要由对流与扩散两个过程决定。具体地，对于步骤s102，所述建立所述流体中的水分子的质量平衡方程可包括：根据所述流体的达西表观速度u与所述水分子的扩散系数d
w
，建立下式(3)表示的所述质量平衡方程：
[0049][0050]
其中，φ0为所述储层的孔隙度的初始值；为所述储层内的孔隙的含水体积分数；以及为所述储层内的任意点的空间位置。
[0051]
所述流体中的水分子的质量平衡方程的初始条件为以及所述流体中的水分子的质量平衡方程的边界条件为(也就是说，在注水井的井壁处的储层孔隙完全被水充满，即孔隙中的水饱和度为1)。其中，φ0为所述储层的孔隙度的初始值；r
w
为所述待诊断井的井筒半径；以及s
wc
为所述储层中的束缚水饱和度。
[0052]
步骤s103，根据菲克扩散定律，建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程。
[0053]
需要说明的是，c(n，t)是储层中的岩石在时刻t、在以位置(例如图2中的o点)为原点建立的一维坐标系(坐标轴的方向为固
‑
液界面在该位置处的法线方向)中的坐标n处的含水体积分数；相应地，储层中的孔隙在位置处的含水体积分数为
[0054]
可建立垂直固
‑
液界面指向固相内部的一维坐标系n，界面处有n＝0，固相内部n>0，如图2所示。对于步骤s103，所述建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程可包括：根据所述菲克扩散定律，建立下式(4)表示的所述流体中的水分子向储层内部扩散的扩散方程：
[0055][0056]
其中，n为以储层中的位置为原点且以所述流体与所述岩石的界面在位置处的法线方向为坐标轴的方向建立的一维坐标系中的坐标；t为时间；d
w
为所述水分子的扩散系
数；以及c(n,t)为所述储层中的岩石的含水体积分数。
[0057]
其中，所述扩散方程的初始条件为c(n,t＝0)＝c0；以及所述扩散方程的边界条件为与d
w
为所述水分子的扩散系数；c(n,t)为所述储层中的岩石的含水体积分数；以及k
f
为膜交换系数。
[0058]
步骤s104，根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程。
[0059]
对于步骤s104，如图3所示，所述确定黏土膨胀损害储层的时空演化模拟方程可包括步骤s301
‑
s303。
[0060]
步骤s301，根据所述质量平衡方程、该质量平衡方程的边界条件与初始条件，确定所述储层内的孔隙的含水体积分数。
[0061]
根据上述公式(3)及所述流体中的水分子的质量平衡方程的初始条件与边界条件，可求解得到所述储层内的孔隙的含水体积分数
[0062]
步骤s302，根据所述扩散方程、所述扩散方程的边界条件、所述扩散方程的初始条件及所述储层内的孔隙的含水体积分数，确定所述储层中的岩石的吸水速率。
[0063]
首先，根据上述公式(4)及所述扩散方程的边界条件、所述扩散方程的初始条件，可求解得到所述储层中的岩石的含水体积分数c(n，t)，
[0064][0065]
其中，即余误差函数，l
‑1{
·
}表示逆拉普拉斯变换。
[0066]
然后，根据c(n，t)及所述储层中的岩石的吸水速率的定义可求解得到
[0067][0068]
具体地，c1、c0的相对大小决定了的正负。若c1＞c0，则说明孔隙中水分含量大于固相中水分含量，将向固相内扩散；反之，若c1＜c0，两式将给出这意味着岩石固相失水。最后，在储层条件下，孔隙中的含水量总是要大于或等于固相中的含水量，需对加限制条件，使其在c1＜c0时等于0。因此，表达式为：
[0069]
[0070]
步骤s303，根据所述储层中的岩石的吸水速率，确定所述黏土膨胀损害储层的时空演化模拟方程。
[0071]
对于步骤s103，所述确定所述黏土膨胀损害储层的时空演化模拟方程可包括：根据所述储层的吸水速率确定由下式(6)表示的所述黏土膨胀损害储层的时空演化模拟方程：
[0072][0073]
其中，为所述储层的孔隙度；以及λ为黏土膨胀系数。
[0074]
具体地，黏土膨胀系数其中cc为岩石中的黏土质量百分数；pi为岩石的塑性系数(无量纲)，若pi＜1～2，该岩石为脆性岩石，若2＜pi＜6，该岩石为塑脆性岩石，若pi＞6，该岩石为塑性岩石；k
′
为经验参数。
[0075]
根据上述公式(6)可求解得到若(即吸水速率为正)，黏土就会膨胀，使得即孔隙度下降。
[0076]
综上所述，本发明创造性地确定待诊断井的预设区域内的储层中的流体的达西表观速度；根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程；建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程；以及根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程。由此，通过所确定的时空演化模拟方程可定量模拟由黏土膨胀引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
[0077]
图4是本发明一实施例提供的确定储层损害程度的方法的流程图。如图4所示，所述确定储层损害程度的方法可包括步骤s401
‑
s402。
[0078]
步骤s401，基于所述的黏土膨胀损害储层的建模方法建立的时空演化模拟方程，确定所述储层的孔隙度。
[0079]
对于上述公式(6)所示的黏土膨胀损害储层的时空演化模拟方程的求取，需要根据公式(3)计算得到对于公式(3)而言，在一维情形下，该类方程可以整理为以下的一般形式：
[0080][0081]
其中，a
a
,b
b
,c
c
可为常数(如扩散系数)，也可为函数(如所述流体的速度)；f可为压力、物质浓度(例如体积分数)、应力等。对时间采用向后差分，空间采用中心差分。则上述方程可以有如下差分格式：
[0082][0083]
其中，i＝1,2,3...n
i
，t＝nδt，n
i
为离散空间点个数。
[0084]
求解区间为x∈(0,x
max
)，δx、δt为空间、时间步长。同时，考虑初始条件f
in
|
n＝0
＝f
i0
，i＝1，2，3...，n
i
和边界条件(f
in
|
i＝1
＝f0,n＝1,2,3...(井壁处)以及)(构造了一个虚拟网格i+1，预设范围的边界处或距井壁数米处)。
[0085]
首先，对于i＝2,3,...,n
i
‑
1整理上述差分格式有：
[0086][0087]
其中，a1
i
，a2
i
，a3
i
分别为，
[0088][0089]
同时，据公式(3)可确定a
i
、b
i
与c
i
。并将所确定的a
i
、b
i
与c
i
代入公式(10)可得到迭代关系式(9)的具体表现形式，由于该迭代关系式(9)的具体表现形式复杂，故在此不对其进行列出。然后，利用初始条件和边界条件进行迭代计算就可得到场f的值。
[0090]
接着，对说明边界条件的差分求解过程进行说明。
[0091]
上述迭代关系式(9)适用于非边界网格。而对于i＝1(井壁处)而言，因为采用的是点中心网格，且其为狄利克雷(dirichlet)边界条件，故直接可得到以下关系式：
[0092]
f
1n
＝f0(常数)，i＝1
ꢀꢀ
(11)
[0093]
对于i＝n(预设范围的边界处距井壁数米处)而言，其为诺伊曼或第二类(neumann)边界条件，增加一个虚拟网格i＝n
i
+1，由知将其代入式(9)可知：
[0094]
[0095]
根据上述过程可求解出场函数f的时空变化情况。由于上述数值模型是针对待诊断井(注水井)的井筒附近储层建立的，在求解某物理量f在井周的分布时，需要采用柱坐标系。由此，式需要变换为这种形式不利于等距差分，可以引入坐标变换：r＝r
w
e
x
′
，其中，r
w
为井筒半径，x
′
为一个无量纲的空间坐标。将这个变换代入一般方程中，可以得到关于x
′
的方程：
[0096][0097]
如果将和作为新的方程系数，则上式和相比，本质上是一样的。因此，便可以在x
′
坐标进行等距差分并沿用前述的迭代格式。计算完f的值后，再将空间坐标从x
′
映射回r即可得到f(r,t)。
[0098]
在通过上述方法计算得到所述储层内的孔隙的含水体积分数c1(r，t)之后，再根据上述公式(5)可计算得到所述储层的吸水速率(图5示出了储层中的特定位置r处的吸水速率随时间变化情况)，由此通过上述黏土膨胀损害储层的建模方法建立的时空演化模拟方程综合考虑了黏土膨胀时多种物理化学因素对储层损害的影响，由此通过该步骤s401求解得到的储层的孔隙度非常精确。
[0099]
步骤s402，基于所确定的所述储层的孔隙度，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
[0100]
其中，所述特征参数可为所述储层的渗透率。
[0101]
对于步骤s402，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数可包括：基于所述储层的孔隙度及公式(14)，确定所述储层的渗透率
[0102][0103]
其中，φ0为孔隙度的初始值；m
k
为第二经验值；以及为所述储层的渗透率的初始值。
[0104]
其中，所述特征参数可为所述储层的滤失系数。
[0105]
对于步骤s402，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数可包括：基于所述储层的孔隙度及公式(15)，确定所述储层的滤失系数
[0106][0107]
其中，φ0为孔隙度的初始值；φ
dmax
为所述储层的最大孔隙度；m
k
为第一经验值；以及所述储层的滤失系数的初始值。
[0108]
在所述特征参数为所述储层的渗透率与所述储层的滤失系数的情况下，可通过公式(14)确定所述储层的渗透率，并通过公式(15)所述储层的渗透率所述储层的渗透率。
[0109]
其中，所述特征参数可为所述储层的表皮系数。
[0110]
对于步骤s402，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数可包括：基于所述储层的孔隙度φ(r,t)及公式确定所述储层的渗透率以及基于所述储层的渗透率及公式(16)，确定所述储层的表皮系数
[0111][0112]
其中，为所述储层的渗透率的初始值；以及r
w
为所述待诊断井的井筒半径，以及r
sw
为所述储层的损害半径。
[0113]
通过该步骤s402得到的特征参数(例如所述储层的渗透率与表皮系数)是时空演化4d定量模拟的结果(如图6所示)。更具体地，图7示出了由储层渗透率损害率(基于所述储层的渗透率及公式确定所述储层的渗透率损害率i(r
i
,t)，其中为的最大值)表征的在第40天黏土膨胀损害储层的半径的示意图(如箭头所指示的半径)，相关工作人员可通过该图7直观地确认储层被损害的程度。因此，可根据渗透率或表皮系数的演化特点进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义。
[0114]
综上所述，本发明创造性地通过所确定的时空演化模拟方程可确定所述储层的孔隙度，再根据所述储层的孔隙度可确定表征待诊断井的预设区域内的储层的损害程度的特征参数(例如所述储层的渗透率和/或表皮系数)，由此可定量模拟由黏土膨胀引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
[0115]
图8是本发明一实施例提供的黏土膨胀损害储层的建模系统的结构图。如图8所示，所述建模系统包括：速度确定装置10，用于确定待诊断井的预设区域内的储层中的流体的达西表观速度；第一建立装置20，用于根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述流体中的水分子的质量平衡方程，第二建立装置30，用于根据菲克扩散定律，建立所述流体中的水分子向储层中的岩石内部扩散的扩散方程；以及模拟方程确定装置40，用于根据所述扩散方程及所述质量平衡方程，确定黏土膨胀损害储层的时空演化模拟方程，其中所述黏土为所述岩石的组成成分。
[0116]
可选的，所述速度确定装置10包括：压力传导方程建立模块(未示出)，用于所述流
体进入所述储层的压力传导方程；以及速度确定模块(未示出)，用于根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
[0117]
可选的，所述模拟方程确定装置40包括：含水量确定模块(未示出)，用于根据所述质量平衡方程及该质量平衡方程的边界条件与初始条件，确定所述储层内的孔隙的含水体积分数；吸水速率确定模块(未示出)，用于根据所述扩散方程、所述扩散方程的边界条件、所述扩散方程的初始条件及所述储层内的孔隙的含水体积分数，确定所述储层中的岩石的吸水速率；以及模拟方程确定模块(未示出)，用于根据所述储层中的岩石的吸水速率，确定所述黏土膨胀损害储层的时空演化模拟方程。
[0118]
所述黏土膨胀损害储层的建模系统与上述黏土膨胀损害储层的建模方法相对于现有技术所具有的优势相同，在此不再赘述。
[0119]
图9是本发明一实施例提供的确定储层损害程度的系统的结构图。如图9所示，所述系统可包括：孔隙度确定装置50，用于基于所述的黏土膨胀损害储层的建模系统建立的时空演化模拟方程，确定所述储层的孔隙度；以及特征参数确定装置60，用于基于所确定的所述储层的孔隙度，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
[0120]
可选的，所述特征参数为所述储层的渗透率和/或所述储层的滤失系数，相应地，所述特征参数确定装置60包括：渗透率计算模块(未示出)，用于基于所述储层的孔隙度φ(r,t)及公式确定所述储层的渗透率和/或滤失系数计算模块(未示出)，用于基于所述储层的孔隙度及公式确定所述储层的滤失系数其中，φ0为所述储层的孔隙度的初始值；为所述储层的孔隙度；φ
dmax
为所述储层的最大孔隙度；m
k
与m
k
分别为第一经验值与第二经验值；为所述储层的渗透率的初始值；以及所述储层的滤失系数的初始值。
[0121]
可选的，所述特征参数为所述储层的表皮系数，相应地，所述特征参数确定装置60包括：渗透率计算模块(未示出)，用于基于所述储层的孔隙度φ(r,t)及公式确定所述储层的渗透率以及表皮系数计算模块(未示出)，用于基于所述储层的渗透率及公式确定所述储层的表皮系数其中，为所述储层的渗透率的初始值，r
w
为所述待诊断井的井筒半径，以及r
sw
为所述储层的损害半径。
[0122]
所述确定储层损害程度的系统与上述确定储层损害程度的方法相对于现有技术所具有的优势相同，在此不再赘述。
[0123]
相应地，本发明一实施例还提供一种机器可读存储介质，该机器可读存储介质上
存储有指令，该指令用于使得机器执行所述的黏土膨胀损害储层的建模方法和/或所述的确定储层损害程度的方法。
[0124]
所述机器可读存储介质包括但不限于相变内存(相变随机存取存储器的简称，phase change random access memory，pram，亦称为rcm/pcram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体(flash memory)或其他内存技术、只读光盘只读存储器(cd
‑
rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备等各种可以存储程序代码的介质。
[0125]
上述步骤s101
‑
s104、步骤s301
‑
s303及步骤s401
‑
s402均可通过计算机来执行。并且，步骤s101
‑
s104所涉及的各种物理化学量的处理过程实现了对黏土膨胀损害储层的时空演化场的模拟，步骤s301
‑
s303所涉及的各种物理化学量的处理过程实现了对黏土膨胀损害储层的时空演化场的具体模拟，以及步骤s401
‑
s402所涉及的各种物理化学量的处理过程实现了对黏土膨胀损害储层的时空演化的预测。
[0126]
以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
[0127]
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。
[0128]
此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。