CN113703069B - 贾敏损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 - Google Patents

贾敏损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 Download PDF

Info

Publication number
CN113703069B
CN113703069B CN202110991049.8A CN202110991049A CN113703069B CN 113703069 B CN113703069 B CN 113703069B CN 202110991049 A CN202110991049 A CN 202110991049A CN 113703069 B CN113703069 B CN 113703069B
Authority
CN
China
Prior art keywords
reservoir
equation
pores
determining
permeability
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
Application number
CN202110991049.8A
Other languages
English (en)
Other versions
CN113703069A (zh
Inventor
蒋官澄
李奕政
贺垠博
杨丽丽
彭春耀
骆小虎
罗绪武
梁兴
谭宾
冉启华
刘小波
程荣超
王增林
陈刚
董腾飞
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Publication of CN113703069A publication Critical patent/CN113703069A/zh
Application granted granted Critical
Publication of CN113703069B publication Critical patent/CN113703069B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G01V20/00—
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/04—Investigating osmotic effects
    • 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
    • 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/082—Investigating permeability by forcing a fluid through a sample
    • G01N15/0826—Investigating permeability by forcing a fluid through a sample and measuring fluid flow rate, i.e. permeation rate or pressure change
    • 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
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00—Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N2013/003—Diffusion; diffusivity between liquids
    • 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
    • G01N2015/0023—Investigating dispersion of liquids
    • G01N2015/0034—Investigating dispersion of liquids in solids

Landscapes

  • Chemical & Material Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Dispersion Chemistry (AREA)
  • Fluid Mechanics (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及油田勘探技术领域，公开一种贾敏效应损害储层的建模方法与系统及确定储层损害程度的方法与系统。所述建模方法包括：确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；根据达西表观速度与流体中的水分子的扩散系数，建立水相运动方程；根据储层的孔隙的孔径分布特征及预设渗透率模型，建立储层的渗透率的分布方程；及根据所述分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程。本发明可定量模拟由贾敏引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演。

Description

贾敏损害油气层的建模方法、损害程度时空演化4D定量与智 能诊断方法及其系统
技术领域
本发明涉及油田勘探技术领域，具体地涉及一种贾敏效应损害储层的建模方法、损害程度确定方法及系统。
背景技术
油田勘探开发的各个时期，由于受到多种内外因素影响，会导致储层原有的物理、化学、热力学和水动力学平衡状态变化，不可避免的使储层近井壁区乃至远井壁区的储层内部渗透率降低，堵塞流体流动，造成储层损害和油井产量下降，甚至“枪毙”储层。造成储层损害的原因是多样的、复杂的，特别是在生产过程中，储层岩石储渗空间、表面润湿性、水动力学场、温度场、岩石种类等不断发生变化，使损害机理随时间而变，且损害周期长、范围宽，损害更具复杂性和叠加性。储层损害一旦发生，必须根据储层损害情况采取相应的解堵措施恢复流体流动通道，以便提高油井产量和水井注入能力。具体地，贾敏效应主要发生在低渗油藏水驱油过程中，储层岩石喉道的大小决定了贾敏效应的强弱程度。对于低渗透储层来说，其岩石喉道狭小，更容易引发贾敏效应，随着油井含水上升，储层内含油饱和度逐渐降低，大量油珠产生的贾敏效应会给低渗透油藏水驱油过程造成很大阻力。例如，不仅会使大量的油滴被捕集成为残余油，而且可能“锁死”已经形成的水驱油通道。因此定量评价各类因素对贾敏效应的影响程度，并对储层贾敏损害程度进行准确的预测尤为重要。
因此，厘清待解堵井储层损害究竟由哪些因素造成、各损害因素所占比例如何，以及储层损害的空间分布规律和随时间变化规律对解堵措施优化设计至关重要，并直接影响解堵和增产效果好坏。
目前，诊断储层损害的方法可分为矿场诊断法和室内评价法。其中，所述矿场诊断法包括试井法。虽然所述试井法可以定量给出表征待诊断井的预设区域内的储层的损害程度的表皮系数、堵塞比、附加压降等重要参数，但由其表征的表皮系数与其它参数相互联系。也就是说，通过所述试井法得出的表皮系数并不仅仅反映真实储层损害特征，还是各个环节、多因素的综合表现(即其是真实损害表皮系数和由井斜表皮系数、储层形状表皮系数、打开储层不完善表皮系数、非达西流表皮系数、射孔表皮系数等组成的拟表皮系数之和)，必须进行表皮系数分解才能得到真实损害表皮系数。其中，所述室内评价法包括岩心流动实验法。所述岩心流动实验法是通过岩心驱替前后的渗透率变化来了解损害程度大小，虽然比较适合研究单因素储层损害，但难以反映较大尺度上储层损害规律。并且，因室内岩心实验条件比较理想化、评价用岩心都是原始状态岩心、无法考虑储层特性动态变化，使实验结果与井下储层真实损害出入较大。
发明内容
本发明的目的是提供一种贾敏效应损害储层的建模方法与系统以及确定储层损害程度的方法与系统，其可定量模拟由贾敏效应引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
为了实现上述目的，本发明第一方面提供一种贾敏效应损害储层的建模方法，所述建模方法包括：确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程。
优选地，所述确定待诊断井的预设区域内的储层中的流体的达西表观速度包括：建立所述流体进入所述储层的压力传导方程；以及根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
优选地，所述建立储层的水相运动方程包括：根据所述流体的达西表观速度u与所述水分子的扩散系数Dw，建立下式表示的所述质量平衡方程：
Figure BDA0003232423690000031
其中，φ0为所述储层的孔隙度的初始值；
Figure BDA0003232423690000032
为所述储层内的孔隙被水相占据的绝对孔隙度；以及
Figure BDA0003232423690000033
为所述储层内的任意点的空间位置；以及根据所述质量平衡方程及所述储层的水相饱和度的时空分布函数
Figure BDA0003232423690000034
建立下式表示的所述水相运动方程：
Figure BDA0003232423690000035
优选地，所述建立所述储层的渗透率的分布方程包括：根据所述储层的孔隙的孔径分布特征，确定孔径为λ的孔隙的体积密度函数与所述储层的水相饱和度的孔径分布方程；以及根据所述预设渗透率模型、所述孔径为λ的孔隙的体积密度函数及所述水相饱和度的孔径分布方程，建立所述储层的渗透率的分布方程。
优选地，在所述储层的孔隙的孔径分布特征为所述储层中的孔径大于λ的孔隙的数量N(＞λ)满足
Figure BDA0003232423690000036
的情况下，所述确定孔径为λ的孔隙的体积密度函数包括：根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定所述储层中的孔隙的总体积为
Figure BDA0003232423690000041
以及根据所述储层中的孔隙的总体积Φmax及所述储层中的孔径小于λ的孔隙的体积
Figure BDA0003232423690000042
确定所述孔径为λ的孔隙的体积密度函数为
Figure BDA0003232423690000043
其中，D为所述孔隙的分形维数；以及λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
优选地，确定所述储层的水相饱和度的孔径分布方程包括：根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定所述被非水相占据的孔隙的体积为
Figure BDA0003232423690000044
其中λpc为非水相受到所述最小毛细管阻力处的孔隙的特定孔径；以及根据所述储层中的孔隙的总体积Φmax与被非水相占据的孔隙的体积Φnw(λ)，确定下式表示的所述水相饱和度的孔径分布方程，
Figure BDA0003232423690000045
其中，D为所述孔隙的分形维数；λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
优选地，在所述储层的预设渗透率模型满足
Figure BDA0003232423690000046
的情况下，所述建立所述储层的渗透率的分布方程包括：根据所述储层的预设渗透率模型
Figure BDA0003232423690000047
所述孔径为λ的孔隙的体积密度函数dξ及所述水相饱和度的孔径分布方程，建立下式表示的所述渗透率的分布方程，
Figure BDA0003232423690000048
通过上述技术方案，本发明创造性地确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程。由此，通过所确定的时空演化模拟方程可定量模拟由贾敏效应引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
本发明第二方面提供一种确定储层损害程度的方法，所述方法包括：基于所述的贾敏效应损害储层的建模方法建立的时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
优选地，所述特征参数为所述储层的渗透率，相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000051
优选地，所述特征参数为所述储层的表皮系数，相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000052
以及基于所述储层的渗透率
Figure BDA0003232423690000053
及公式
Figure BDA0003232423690000054
确定所述储层的表皮系数
Figure BDA0003232423690000055
其中，
Figure BDA0003232423690000056
为所述储层的渗透率的初始值，rw为所述待诊断井的井筒半径，以及rsw为所述储层的损害半径。
通过上述技术方案，通过所确定的时空演化模拟方程可确定表征待诊断井的预设区域内的储层的损害程度的特征参数(例如所述储层的渗透率)，由此可定量模拟由贾敏效应引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
相应地，本发明第三方面还提供一种贾敏效应损害储层的建模系统，所述建模系统包括：速度确定装置，用于确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；第一建立装置，用于根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；第二建立装置，用于根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及模拟方程确定装置，用于根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程。
所述贾敏效应损害储层的建模系统与上述贾敏效应损害储层的建模方法相对于现有技术所具有的优势相同，在此不再赘述。
相应地，本发明第四方面还提供一种确定储层损害程度的系统，所述系统包括：接收装置，用于接收基于所述的贾敏效应损害储层的建模系统建立的时空演化模拟方程；以及特征参数确定装置，用于基于所述时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
所述确定储层损害程度的系统与上述确定储层损害程度的方法相对于现有技术所具有的优势相同，在此不再赘述。
相应地，本发明第五方面还提供一种机器可读存储介质，所述机器可读存储介质上存储有指令，该指令用于使得机器执行所述的贾敏效应损害储层的建模方法和/或所述的确定储层损害程度的方法。
本发明实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本发明实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本发明实施例，但并不构成对本发明实施例的限制。在附图中：
图1是本发明一实施例提供的贾敏效应损害储层的建模方法的流程图；
图2是本发明一实施例提供的建立所述储层的渗透率的分布方程的流程图；
图3是本发明一实施例提供的确定孔径为λ的孔隙的体积密度函数的流程图；
图4是本发明一实施例提供的确定所述储层的水相饱和度的孔径分布方程的流程图；
图5是本发明一实施例提供的渗透率随含水饱和度演化的示意图；
图6是本发明一实施例提供的含水饱和度随空间演化的示意图；
图7是本发明一实施例提供的渗透率损害率随含水饱和度的示意图；
图8是本发明一实施例提供的贾敏效应损害储层的建模系统的结构图；以及
图9是本发明一实施例提供的确定储层损害程度的系统的结构图。
具体实施方式
以下结合附图对本发明的具体实施方式进行详细说明。应当理解的是，此处所描述的具体实施方式仅用于说明和解释本发明，并不用于限制本发明。
在低渗透储层内的油、水两相流动过程中，通常含有大量分散的油滴和气泡。当这些个体运动到储层中的狭小孔喉(例如孔径为2.5μm的孔隙)时，由于这些个体的直径大于孔喉直径，其流动受到阻力，迫使油滴或气泡等发生变形，故将油滴或气泡等在直径分布不均匀的储层中运动时受到的附加的毛细管阻力的现象称为贾敏效应，如图1所示。与水锁效应相比，微小孔喉中受到阻力被迫变形的是油滴或者气泡，此时润湿相是水，非润湿相是油滴或气泡。
通常情况下，产生贾敏效应需同时满足以下条件：首先，原油应该处于分散状态，不能以连续的油流形式运动；其次，地质储层的孔道直径应该比较细小。当分散的油滴运动到油层狭小孔道时，由于其直径大于孔道直径，就会受到阻力而致使其无法继续运动，或者油滴在变直径孔道中流动时会产生毛细管效应，这样就会导致贾敏损害(因贾敏效应导致储层的渗透率下降的情形)。
具体地，如图1所示，在储层左右两侧的驱替压力P的作用下，液滴从左向右运动。当油滴运动到孔隙1处时，因孔隙1的阻挡，对于油滴前端(2R1＝λ，其中R1为前端液面的曲率半径，λ为孔隙1的孔径)，其受到泡前压力P1的作用；对于油滴后端(其直径
Figure BDA0003232423690000081
其中R2为后端液面的曲率半径，
Figure BDA0003232423690000082
为孔隙的平均孔径)，其受到P3的作用；以及对于油滴内部，其泡内压力为P2，并且，P1、P2、P3满足以下两个等式：
Figure BDA0003232423690000083
Figure BDA0003232423690000084
由此，因贾敏效应油滴受到的附加阻力为：
Figure BDA0003232423690000085
考虑液滴前端
Figure BDA0003232423690000086
液滴后端
Figure BDA0003232423690000087
Figure BDA0003232423690000088
其中σ为液滴的表面张力。当驱替压力(即孔喉外压力差)至少达到Pc时，气泡或液滴才能通过孔喉，否则就会被堵住。
贾敏效应受储层的孔隙结构、岩性、物性、侵入流体等多种因素的影响。贾敏损害同岩石的孔隙介质的几何特性关系密切，不同的孔喉结构分布模式和复杂性会导致水润湿相在岩石中的展布模式显著变化，从而影响储层的渗透率。因此，本发明各实施例的核心是要建立储层中的孔隙内的水相饱和度变化的动力学模型(即储层中的水相运动方程及渗透率的分布方程)。具体地，基于储层中的孔隙内的流体的对流扩散关系(即质量平衡方程)等建立储层中的水相运动方程；基于孔隙的孔径分布特征及预设渗透率模型，建立渗透率的分布方程；然后根据所述水相运动方程及所述渗透率的分布方程，确定贾敏效应影响待诊断井周围的储层中的渗透率分布的时空演化控制唯象模型，由此就可诊断渗透率等储层损害特征参数的时空场分布。
需要说明的是，为了简单描述起见，在本发明的各个实施例中的随时空演化的物理量、化学量可省略变量
Figure BDA0003232423690000091
例如
Figure BDA0003232423690000092
可简写为K。
图1是本发明一实施例提供的贾敏效应损害储层的建模方法的流程图。所述建模方法可包括步骤S101-S104。
步骤S101，确定待诊断井的预设区域内的储层中的流体的达西表观速度。
其中，所述储层的渗透率低于预设渗透率；以及所述待诊断井可例如为采油井。
对于步骤S101，所述确定储层中的流体的速度可包括：建立所述流体进入所述储层的压力传导方程；以及根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
具体地，压力是驱动固-液混合液从注水井的井筒持续侵入待诊断井周围的储层的动力，由此可建立如公式(1)的所述流体进入储层的压力传导方程：
Figure BDA0003232423690000101
其中，所述压力传导方程的初始条件为
Figure BDA0003232423690000102
以及边界条件为
Figure BDA0003232423690000103
(也就是说，在注水井的井壁处的压力为Pw)。
根据公式(1)及达西公式(如下式(2))可确定所述流体的达西表观速度，
Figure BDA0003232423690000104
其中，
Figure BDA0003232423690000105
为所述流体的压力；φ0为所述储层的孔隙度的初始值；μ为流体粘度；ct为流体-岩石综合压缩系数以及
Figure BDA0003232423690000106
为所述储层的渗透率。
步骤S102，根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程。
在储层条件下，所述储层中的孔隙内的不同位置的含水量满足质量守恒方程。其中，外来水相在所述储层内的运动主要由对流与扩散两个过程决定。具体地，对于步骤S102，所述建立所述储层的水相运动方程可包括：根据所述流体的达西表观速度u与所述水分子的扩散系数Dw，建立下式表示的所述质量平衡方程：
Figure BDA0003232423690000107
其中，φ0为所述储层的孔隙度的初始值；
Figure BDA0003232423690000108
为所述储层内的孔隙被水相占据的绝对孔隙度；以及
Figure BDA0003232423690000109
为所述储层内的任意点的空间位置(例如，以所述待诊断井的中心为原点)。
根据所述质量平衡方程及所述储层的水相饱和度的时空分布函数
Figure BDA00032324236900001010
建立下式(3)表示的所述水相运动方程：
Figure BDA0003232423690000111
所述水相运动方程的初始条件为
Figure BDA0003232423690000112
以及所述水相运动方程的边界条件为
Figure BDA0003232423690000113
(也就是说，在注水井的井壁处的储层孔隙完全被水充满，即孔隙中的水相饱和度为1)。其中，φ0为所述储层的孔隙度的初始值；rw为所述待诊断井的井筒半径；以及Swc为所述储层中的束缚水饱和度。
步骤S103，根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程。
对于步骤S103，如图2所示，所述建立所述储层的渗透率的分布方程可包括步骤S201-S202。
步骤S201，根据所述储层的孔隙的孔径分布特征，确定孔径为λ的孔隙的体积密度函数与所述储层的水相饱和度的孔径分布方程。
储层岩石的孔隙结构具有强烈的不规则性，受到多种因素的控制，使得经典几何学为基础的孔隙结构模型与真实的孔隙相差太远。为了定量描述储层的孔隙结构，采用分形理论来研究孔隙结构中两相流的贾敏效应，并忽略水相在孔隙通道内表面的润湿薄膜影响。根据分形的几何原理，若储层的孔径分布具有分形特征，则储层中孔径大于λ的孔隙的数量N(>λ)与λ有如下幂函数关系：
Figure BDA0003232423690000114
其中，D为孔隙的分形维数(2<D<3)。
在所述储层的孔隙的孔径分布特征为所述储层中的孔径大于λ的孔隙的数量N(＞λ)满足上式
Figure BDA0003232423690000115
的情况下，如图3所示，所述步骤S201中的确定孔径为λ的孔隙的体积密度函数可包括步骤S301-S302。
步骤S301，根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定所述储层中的孔隙的总体积为
Figure BDA0003232423690000121
具体地，根据上式
Figure BDA0003232423690000122
可得到储层中的孔隙的总数目N(＞λmin)及孔径小于λ的孔隙的数目：
Figure BDA0003232423690000123
Figure BDA0003232423690000124
根据公式(4)与(5)可得到下式(6)
Figure BDA0003232423690000125
上式(6)中的f(λ)为储层的孔径分布密度函数，则孔径小于λpc的孔隙的数量N(＜λpc)与λ之间的关系为下式(7)所表示的幂函数关系：
Figure BDA0003232423690000126
上述相关公式中的λ、λmin、λmax与λpc分别为所述孔隙的孔径、最小孔径、最大孔径(λmin与λmax可以由平均孔径与孔径分布标准差得出，一般而言有
Figure BDA0003232423690000127
)与发生贾敏效应的最大孔径(即非水相受到所述最小毛细管阻力处的孔隙的特定孔径)；a为比例常数。
接着，由公式(7)可得到储层的孔径分布密度函数f(λ)满足下式(8)：
Figure BDA0003232423690000128
式中a'＝-Da为比例常数。
由上式(8)表示的所述孔径分布密度函数可以得到储层中的孔隙的总体积分形表达式：
Figure BDA0003232423690000131
其中α为与孔隙的形状有关的常数(孔隙的形状为立方体时α＝1，或者孔隙的形状为球体时α＝π/6)，通过积分可得：
Figure BDA0003232423690000132
类似地，所述储层中的孔径小于λ的孔隙的体积
Figure BDA0003232423690000133
步骤S302，根据所述储层中的孔隙的总体积Φmax及所述储层中的孔径小于λ的孔隙的体积
Figure BDA0003232423690000134
确定所述孔径为λ的孔隙的体积密度函数为
Figure BDA0003232423690000135
其中，D为所述孔隙的分形维数；以及λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)(常数)。
步骤S202，根据所述预设渗透率模型、所述孔径为λ的孔隙的体积密度函数及所述水相饱和度的孔径分布方程，建立所述储层的渗透率的分布方程。
如图4所示，步骤S202中的确定所述储层的水相饱和度的孔径分布方程可包括步骤S401-S402。
步骤S401，根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定被非水相占据的孔隙的体积为
Figure BDA0003232423690000136
其中，λpc为非水相受到所述最小毛细管阻力处的孔隙的特定孔径(即发生贾敏效应时的最大孔喉直径)。
在储层为水湿(即亲水性)的情况下，向所述储层的两侧施加一个驱替压力P(从零逐渐增加到P)，根据贾敏效应产生条件，油滴或气泡可顺利通过孔径处于λmax到λpc范围之间的孔隙(即连通孔径由最大λmax到λpc逐次打开)。当压力P未到Pc时，λ大于λpc的孔喉完全被水占据，而λ小于λpc的孔喉还被油滴或气泡堵塞。
由于根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，可确定上式(8)确定的孔径分布密度函数f(λ)，则进而可根据f(λ)得到压力P未到Pc时的非水相占据的孔隙体积Φnw(λ)，
Figure BDA0003232423690000141
步骤S402，根据所述储层中的孔隙的总体积Φmax与所述被非水相占据的孔隙的体积Φnw(λ)，确定下式表示的所述水相饱和度的孔径分布方程
Figure BDA0003232423690000142
其中，D为所述孔隙的分形维数；λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
具体地，根据公式(10)与公式(12)可确定水相占据的孔隙体积为：
Figure BDA0003232423690000143
然后再结合公式(10)可得到下式：
Figure BDA0003232423690000144
由于λmin＜＜λmax，
Figure BDA0003232423690000145
故上式可写为下式(13)，即可确定所述水相饱和度的孔径分布方程(13)，
Figure BDA0003232423690000146
根据流体力学理论，在压差ΔP作用下，粘度为μ的流体流过毛细管束总流量可以由Hagen-poiseuille方程描述，
Figure BDA0003232423690000151
上式中，λi为孔隙通道i的直径(即孔隙i的孔径)，Vi为孔隙通道i的体积，Vi＝Φ0MLmξi，Φ0为储层的孔隙度，ξi为孔隙通道i所占孔隙体积的分数，Lc为孔隙通道i的弯曲长度，Lm为孔隙通道i的直线长度，M为平均孔喉截面积
Figure BDA0003232423690000152
为孔隙通道i的迂曲度(同时可由经验公式
Figure BDA0003232423690000153
计算)，将上述各量代入上式得到，
Figure BDA0003232423690000154
结合达西定律，可得到毛管束模型的渗透率可表示为
Figure BDA0003232423690000155
在本发明实施例中，可将储层的渗透通道看做是多个毛管束累积起来的。由于孔隙的尺寸分布具有连续性，因此可以将所述毛管束模型的渗透率的表达式用积分形式写出：
Figure BDA0003232423690000156
在所述储层的预设渗透率模型满足
Figure BDA0003232423690000157
的情况下，所述建立所述储层的渗透率的分布方程可包括：根据所述储层的预设渗透率模型
Figure BDA0003232423690000158
所述孔径为λ的孔隙的体积密度函数dξ及所述水相饱和度的孔径分布方程，建立下式表示的所述渗透率的分布方程，
Figure BDA0003232423690000159
其中，所述压力传导方程的初始条件为K(t＝0)＝1。K(Sw)为渗透率关于含水饱和度的函数(如图5所示，该函数可简称为无因次渗透率)；φ0为储层的初始孔隙度，Sw为含水饱和度；D为所述孔隙的分形维数，其可由
Figure BDA0003232423690000161
计算(D越大说明小孔喉所占比例越大，常规储层λmax＞20μm，D<2.6，低渗储层λmax＜0.2μm，D>2.8)，τ为迂曲度，λmax为最大孔喉直径，λmin为最小孔喉直径)。
具体地，首先，逐渐增加驱替压力直至P，此时直径大于λ的孔喉被水水占据，被水相占据的绝对孔隙度为φw，则将公式(11)代入公式(14)可得出所述储层的渗透率的孔径分布函数：
Figure BDA0003232423690000162
然后，根据式(13)中的λ与含水饱和度Sw(λ)之间的关系式将上述渗透率的孔径分布函数K(λ)中的变量λ替换为Sw，并结合φw＝φ0Sw，可得到公式(15)(即得到无因次渗透率的表达式，该无因次渗透率为含水饱和度的函数)。
步骤S104，根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程。
具体地，根据公式(3)可求解到所述储层的水相饱和度的时空分布函数
Figure BDA0003232423690000163
(如图6所示)，并将
Figure BDA0003232423690000164
代入公式(15)可得到所述储层的渗透率在四维时空上的分布形式，即贾敏效应损害储层的时空演化模拟方程。
综上所述，本发明创造性地确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；建立所述储层的渗透率的分布方程；以及根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程。由此，通过所确定的时空演化模拟方程可定量模拟由贾敏效应引起的储层损害特征的四维时空演化过程，从而对未发生储层损害的井进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义，以及对已损害井优化设计解堵措施、提高或恢复油井产量和水井注水能力，以及提高油藏数值模拟精度都具有十分重大意义。
相应地，本发明另一实施例还提供一种确定储层损害程度的方法，所述方法可包括：基于所述的贾敏效应损害储层的建模方法建立的时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
对于上述公式(15)所示的贾敏效应损害储层的时空演化模拟方程的求取，需要根据公式(3)计算得到
Figure BDA0003232423690000171
对于公式(3)而言，在一维情形下，该类方程可以整理为以下的一般形式：
Figure BDA0003232423690000172
其中，aa,bb,cc可为常数(如扩散系数)，也可为函数(如所述流体的速度)；f可为压力、水相饱和度、物质浓度(例如体积分数)、应力等。对时间采用向后差分，空间采用中心差分。则上述方程可以有如下差分格式：
Figure BDA0003232423690000173
其中，i＝1,2,3…Ni，
Figure BDA0003232423690000174
n＝1,2,3…，t＝nΔt，Ni为离散空间点个数。
求解区间为x∈(0,xmax)，Δx、Δt为空间、时间步长。同时，考虑初始条件fi n|n＝0＝fi 0,i＝1,2,3…,Ni和边界条件(fi n|i＝1＝f0,n＝1,2,3...(井壁处)以及
Figure BDA0003232423690000175
n＝1,2,3...)(构造了一个虚拟网格i+1，预设范围的边界处或距井壁数米处)。
首先，对于i＝2,3,...,Ni-1整理上述差分格式有：
Figure BDA0003232423690000176
Figure BDA0003232423690000181
其中，A1i，A2i，A3i分别为，
Figure BDA0003232423690000182
同时，据公式(3)可确定ai、bi与ci。并将所确定的ai、bi与ci代入公式(19)可得到迭代关系式(18)的具体表现形式，由于该迭代关系式(18)的具体表现形式复杂，故在此不对其进行列出。然后，利用初始条件和边界条件进行迭代计算就可得到场f的值。
接着，对说明边界条件的差分求解过程进行说明。
上述迭代关系式(18)适用于非边界网格。而对于i＝1(井壁处)而言，因为采用的是点中心网格，且其为狄利克雷(Dirichlet)边界条件，故直接可得到以下关系式：
f1 n＝f0(常数)，i＝1 (20)
对于i＝N(预设范围的边界处距井壁数米处)而言，其为诺伊曼或第二类(Neumann)边界条件，增加一个虚拟网格i＝Ni+1，由
Figure BDA0003232423690000183
n＝1,2,3...知
Figure BDA0003232423690000184
将其代入式(18)可知：
Figure BDA0003232423690000185
根据上述过程可求解出场函数f的时空变化情况。由于上述数值模型是针对待诊断井(注水井)的井筒附近储层建立的，在求解某物理量f在井周的分布时，需要采用柱坐标系。由此，式
Figure BDA0003232423690000191
需要变换为
Figure BDA0003232423690000192
这种形式不利于等距差分，可以引入坐标变换：r＝rwex′，其中，rw为井筒半径，x′为一个无量纲的空间坐标。将这个变换代入一般方程中，可以得到关于x′的方程：
Figure BDA0003232423690000193
如果将
Figure BDA0003232423690000194
和
Figure BDA0003232423690000195
作为新的方程系数，则上式和
Figure BDA0003232423690000196
相比，本质上是一样的。因此，便可以在x′坐标进行等距差分并沿用前述的迭代格式。计算完f的值后，再将空间坐标从x′映射回r即可得到f(r,t)。
在通过上述方法计算得到所述储层的水相饱和度
Figure BDA0003232423690000197
之后，再根据上述公式(15)可计算得到所述储层的渗透率
Figure BDA0003232423690000198
由此通过上述贾敏效应损害储层的建模方法建立的时空演化模拟方程综合考虑了贾敏损害时多种物理化学因素对储层损害的影响，由此通过本实施例求解得到的储层的渗透率非常精确。
在得到所述储层的渗透率的基础上，表征待诊断井的预设区域内的储层的损害程度的特征参数。
在一实施例中，所述特征参数可为所述储层的渗透率损害率。
相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数可包括：
基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000199
以及基于所述储层的渗透率
Figure BDA00032324236900001910
及公式(23)，确定所述储层的渗透率损害率
Figure BDA00032324236900001911
Figure BDA0003232423690000201
其中，
Figure BDA0003232423690000202
为
Figure BDA0003232423690000203
的最大值。
在另一实施例中，所述特征参数可为所述储层的表皮系数。
所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数可包括：基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000204
以及基于所述储层的渗透率
Figure BDA0003232423690000205
及公式(24)，确定所述储层的表皮系数
Figure BDA0003232423690000206
Figure BDA0003232423690000207
其中，
Figure BDA0003232423690000208
为所述储层的渗透率的初始值，以及
Figure BDA0003232423690000209
rw为所述待诊断井的井筒半径，以及rsw为所述储层的损害半径。
通过上述各个实施例得到的特征参数(例如所述储层的渗透率损害率I(Sw)是时空演化4D定量模拟的结果(分别如图7所示)。因此，可根据渗透率或表皮系数的演化特点进行储层损害定量预测和损害规律时空推演，为预防或避免储层损害、制定油藏的开发方案以及之后增产措施具有科学指导意义。
图9是本发明一实施例提供的贾敏效应损害储层的建模系统的结构图。如图9所示，所述建模系统包括：速度确定装置10，用于确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；第一建立装置20，用于根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；第二建立装置30，用于根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及模拟方程确定装置40，用于根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程。
可选的，所述速度确定装置10包括：压力传导方程建立模块(未示出)，用于所述流体进入所述储层的压力传导方程；以及速度确定模块(未示出)，用于根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
可选的，所述第一建立装置20包括：质量平衡方程建立模块(未示出)，用于根据所述流体的达西表观速度u与所述水分子的扩散系数Dw，建立下式表示的所述质量平衡方程：
Figure BDA0003232423690000211
其中，φ0为所述储层的孔隙度的初始值；
Figure BDA0003232423690000212
为所述储层内的孔隙被水相占据的绝对孔隙度；以及
Figure BDA0003232423690000213
为所述储层内的任意点的空间位置；以及水相运动方程建立模块(未示出)，用于根据所述质量平衡方程及所述储层的水相饱和度的时空分布函数
Figure BDA0003232423690000214
建立下式表示的所述水相运动方程：
Figure BDA0003232423690000215
可选的，所述第二建立装置30包括：第一分布方程建立模块(未示出)，用于根据所述储层的孔隙的孔径分布特征，确定孔径为λ的孔隙的体积密度函数与所述储层的水相饱和度的孔径分布方程；以及第二分布方程建立模块(未示出)，用于根据所述预设渗透率模型、所述孔径为λ的孔隙的体积密度函数及所述水相饱和度的孔径分布方程，建立所述储层的渗透率的分布方程。
可选的，在所述储层的孔隙的孔径分布特征为所述储层中的孔径大于λ的孔隙的数量N(＞λ)满足
Figure BDA0003232423690000216
的情况下，所述第一分布方程建立模块(未示出)包括：第一体积确定单元，用于根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定所述储层中的孔隙的总体积为
Figure BDA0003232423690000221
以及体积密度函数确定单元，用于根据所述储层中的孔隙的总体积Φmax及所述储层中的孔径小于λ的孔隙的体积
Figure BDA0003232423690000222
确定所述孔径λ的孔隙的体积密度函数为
Figure BDA0003232423690000223
其中，D为所述孔隙的分形维数；以及λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
可选的，所述第一分布方程建立模块(未示出)包括：第二分布方程建立模块(未示出)，用于根据所述储层中的孔径大于λ的孔隙的数量N(＞λ)，确定被非水相占据的孔隙的体积为
Figure BDA0003232423690000224
其中λpc为非水相受到所述最小毛细管阻力处的孔隙的特定孔径；以及水相饱和度分布方程确定单元，用于根据所述储层中的孔隙的总体积Φmax与所述被非水相占据的孔隙的体积Φnw(λ)，确定下式表示的所述水相饱和度的孔径分布方程，
Figure BDA0003232423690000225
其中，D为所述孔隙的分形维数；λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
所述贾敏效应损害储层的建模系统与上述贾敏效应损害储层的建模方法相对于现有技术所具有的优势相同，在此不再赘述。
图8是本发明一实施例提供的确定储层损害程度的系统的结构图。如图8所示，所述系统可包括：接收装置50，用于接收基于所述的贾敏效应损害储层的建模系统建立的时空演化模拟方程；以及特征参数确定装置60，用于基于所述时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
可选的，所述特征参数为所述储层的渗透率。相应地，所述特征参数确定装置60包括：渗透率计算模块(未示出)，用于基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000231
可选的，所述特征参数为所述储层的表皮系数。相应地，所述特征参数确定装置60包括：渗透率计算模块(未示出)，用于基于所述时空演化模拟方程，确定所述储层的渗透率
Figure BDA0003232423690000232
以及表皮系数计算模块(未示出)，用于基于所述储层的渗透率
Figure BDA0003232423690000233
及公式
Figure BDA0003232423690000234
确定所述储层的表皮系数
Figure BDA0003232423690000235
所述确定贾敏效应程度的系统与上述确定贾敏效应程度的方法相对于现有技术所具有的优势相同，在此不再赘述。
相应地，本发明一实施例还提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令用于使得机器执行所述的贾敏效应损害储层的建模方法和/或所述的确定储层损害程度的方法。
所述机器可读存储介质包括但不限于相变内存(相变随机存取存储器的简称，Phase Change Random Access Memory，PRAM，亦称为RCM/PCRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体(Flash Memory)或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁盘存储或其他磁性存储设备等各种可以存储程序代码的介质。
上述步骤S101-S104、步骤S201-S202、步骤S301-S302及步骤S401-S402均可通过计算机来执行。并且，步骤S101-S104所涉及的各种物理化学量的处理过程实现了对贾敏效应损害储层的时空演化场的模拟，步骤S201-S202、步骤S301-S302及步骤S401-S402所涉及的各种物理化学量的处理过程实现了对贾敏效应损害储层的时空演化场的具体模拟。
以上结合附图详细描述了本发明的优选实施方式，但是，本发明并不限于上述实施方式中的具体细节，在本发明的技术构思范围内，可以对本发明的技术方案进行多种简单变型，这些简单变型均属于本发明的保护范围。
另外需要说明的是，在上述具体实施方式中所描述的各个具体技术特征，在不矛盾的情况下，可以通过任何合适的方式进行组合。为了避免不必要的重复，本发明对各种可能的组合方式不再另行说明。
此外，本发明的各种不同的实施方式之间也可以进行任意组合，只要其不违背本发明的思想，其同样应当视为本发明所公开的内容。

Claims (10)

1.一种贾敏效应损害储层的建模方法，其特征在于，所述建模方法包括：
确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；
根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；
根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及
根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程，
其中，所述建立所述储层的渗透率的分布方程包括：
根据所述储层的孔隙的孔径分布特征，确定孔径为λ的孔隙的体积密度函数与所述储层的水相饱和度的孔径分布方程；以及
根据所述预设渗透率模型、所述孔径为λ的孔隙的体积密度函数及所述水相饱和度的孔径分布方程，建立所述储层的渗透率的分布方程，
其中，在所述储层的孔隙的孔径分布特征为所述储层中的孔径大于λ的孔隙的数量N(>λ)满足
Figure FDA0003864056580000011
的情况下，所述确定孔径为λ的孔隙的体积密度函数包括：
根据所述储层中的孔径大于λ的孔隙的数量N(>λ)，确定所述储层中的孔隙的总体积为
Figure FDA0003864056580000012
以及
根据所述储层中的孔隙的总体积Φmax及所述储层中的孔径小于λ的孔隙的体积
Figure FDA0003864056580000013
确定所述孔径为λ的孔隙的体积密度函数为
Figure FDA0003864056580000021
其中，D为所述孔隙的分形维数；以及λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)，
其中，确定所述储层的水相饱和度的孔径分布方程包括：
根据所述储层中的孔径大于λ的孔隙的数量N(>λ)，确定所述被非水相占据的孔隙的体积为
Figure FDA0003864056580000022
其中λpc为非水相受到最小毛细管阻力处的孔隙的特定孔径；以及
根据所述储层中的孔隙的总体积Φmax与被非水相占据的孔隙的体积Φnw(λ)，确定下式表示的所述水相饱和度的孔径分布方程，
Figure FDA0003864056580000023
其中，D为所述孔隙的分形维数；λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
2.根据权利要求1所述的贾敏效应损害储层的建模方法，其特征在于，所述确定待诊断井的预设区域内的储层中的流体的达西表观速度包括：
建立所述流体进入所述储层的压力传导方程；以及
根据所述压力传导方程及达西公式，确定所述流体的达西表观速度。
3.根据权利要求1所述的贾敏效应损害储层的建模方法，其特征在于，所述建立储层的水相运动方程包括：
根据所述流体的达西表观速度u与所述水分子的扩散系数Dw，建立下式表示的质量平衡方程：
Figure FDA0003864056580000024
其中，φ0为所述储层的孔隙度的初始值；
Figure FDA0003864056580000031
为所述储层内的孔隙被水相占据的绝对孔隙度；以及
Figure FDA0003864056580000032
为所述储层内的任意点的空间位置；以及t为时间；以及
根据所述质量平衡方程及所述储层的水相饱和度的时空分布函数
Figure FDA0003864056580000033
建立下式表示的所述水相运动方程：
Figure FDA0003864056580000034
4.根据权利要求1所述的贾敏效应损害储层的建模方法，其特征在于，在所述储层的预设渗透率模型满足的情况下，所述建立所述储层的渗透率的分布方程包括：
根据所述储层的预设渗透率模型
Figure FDA0003864056580000035
所述孔径为λ的孔隙的体积密度函数dξ及所述水相饱和度的孔径分布方程，建立下式表示的所述渗透率的分布方程，
Figure FDA0003864056580000036
5.一种确定储层损害程度的方法，其特征在于，所述方法包括：
基于根据权利要求1-4中任一项权利要求所述的贾敏效应损害储层的建模方法建立的时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
6.根据权利要求5所述的确定储层损害程度的方法，其特征在于，所述特征参数为所述储层的渗透率，
相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：
基于所述时空演化模拟方程，确定所述储层的渗透率
Figure FDA0003864056580000041
其中，
Figure FDA0003864056580000042
为所述储层内的任意点的空间位置；以及t为时间。
7.根据权利要求5所述的确定储层损害程度的方法，其特征在于，所述特征参数为所述储层的表皮系数，
相应地，所述确定表征待诊断井的预设区域内的储层的损害程度的特征参数包括：
基于所述时空演化模拟方程，确定所述储层的渗透率
Figure FDA0003864056580000043
以及
基于所述储层的渗透率
Figure FDA0003864056580000044
及公式
Figure FDA0003864056580000045
确定所述储层的表皮系数
Figure FDA0003864056580000046
其中，
Figure FDA0003864056580000047
为所述储层的渗透率的初始值，rw为所述待诊断井的井筒半径，rsw为所述储层的损害半径，
Figure FDA0003864056580000048
为所述储层内的任意点的空间位置，以及t为时间。
8.一种贾敏效应损害储层的建模系统，其特征在于，所述建模系统包括：
速度确定装置，用于确定待诊断井的预设区域内的储层中的流体的达西表观速度，其中所述储层的渗透率低于预设渗透率；
第一建立装置，用于根据所述流体的达西表观速度与所述流体中的水分子的扩散系数，建立所述储层的水相运动方程；
第二建立装置，用于根据所述储层的孔隙的孔径分布特征及所述储层的预设渗透率模型，建立所述储层的渗透率的分布方程；以及
模拟方程确定装置，用于根据所述渗透率的分布方程及所述水相运动方程，确定贾敏效应损害储层的时空演化模拟方程，其中所述时空演化模拟方程用于模拟由贾敏效应引起的储层损害特征的四维时空演化过程，
其中，所述建立所述储层的渗透率的分布方程包括：
根据所述储层的孔隙的孔径分布特征，确定孔径为λ的孔隙的体积密度函数与所述储层的水相饱和度的孔径分布方程；以及
根据所述预设渗透率模型、所述孔径为λ的孔隙的体积密度函数及所述水相饱和度的孔径分布方程，建立所述储层的渗透率的分布方程，
其中，在所述储层的孔隙的孔径分布特征为所述储层中的孔径大于λ的孔隙的数量N(>λ)满足
Figure FDA0003864056580000051
的情况下，所述确定孔径为λ的孔隙的体积密度函数包括：
根据所述储层中的孔径大于λ的孔隙的数量N(>λ)，确定所述储层中的孔隙的总体积为
Figure FDA0003864056580000052
以及
根据所述储层中的孔隙的总体积Φmax及所述储层中的孔径小于λ的孔隙的体积
Figure FDA0003864056580000053
确定所述孔径为λ的孔隙的体积密度函数为
Figure FDA0003864056580000054
其中，D为所述孔隙的分形维数；以及λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)，
其中，确定所述储层的水相饱和度的孔径分布方程包括：
根据所述储层中的孔径大于λ的孔隙的数量N(>λ)，确定所述被非水相占据的孔隙的体积为
Figure FDA0003864056580000055
其中λpc为非水相受到所述最小毛细管阻力处的孔隙的特定孔径；以及
根据所述储层中的孔隙的总体积Φmax与被非水相占据的孔隙的体积Φnw(λ)，确定下式表示的所述水相饱和度的孔径分布方程，
Figure FDA0003864056580000061
其中，D为所述孔隙的分形维数；λ、λmin与λmax分别为所述孔隙的孔径、最小孔径与最大孔径；以及A＝αa'/(3-D)。
9.一种确定储层损害程度的系统，其特征在于，所述系统包括：
接收装置，用于接收基于根据权利要求8所述的贾敏效应损害储层的建模系统建立的时空演化模拟方程；以及
特征参数确定装置，用于基于所述时空演化模拟方程，确定表征待诊断井的预设区域内的储层的损害程度的特征参数。
10.一种机器可读存储介质，其特征在于，所述机器可读存储介质上存储有指令，该指令用于使得机器执行上述权利要求1-4中任一项所述的贾敏效应损害储层的建模方法和/或上述权利要求5-7中任一项所述的确定储层损害程度的方法。
CN202110991049.8A 2020-08-26 2021-08-26 贾敏损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统 Active CN113703069B (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2020108740813 2020-08-26
CN202010874081 2020-08-26

Publications (2)

Publication Number Publication Date
CN113703069A CN113703069A (zh) 2021-11-26
CN113703069B true CN113703069B (zh) 2022-12-02

Family

ID=78655493

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110991049.8A Active CN113703069B (zh) 2020-08-26 2021-08-26 贾敏损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统

Country Status (1)

Country Link
CN (1) CN113703069B (zh)

Family Cites Families (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5813798A (en) * 1997-03-28 1998-09-29 Whiffen; Greg Piecewise continuous control of groundwater remediation
US7657415B2 (en) * 2002-05-31 2010-02-02 Schlumberger Technology Corporation Subterranean formation treatment methods using a darcy scale and pore scale model
US10100638B2 (en) * 2014-11-20 2018-10-16 Schlumberger Technology Corporation Method for reservoir evaluation employing non-equilibrium asphaltene component
CN104504755B (zh) * 2014-12-30 2017-04-19 华中科技大学 一种山体中分布式地下设施温度场仿真方法
FR3041478A1 (fr) * 2015-09-21 2017-03-24 Commissariat Energie Atomique Determination d’une distribution spatiale de la permeabilite d’une electrode de cellule electrochimique
CA2932670A1 (en) * 2016-06-06 2017-12-06 Uti Limited Partnership Assessment of reservoir heterogeneity by using produced water chemistry
CN107506534B (zh) * 2017-08-04 2021-02-19 陕西延长石油（集团）有限责任公司 一种二氧化碳驱油封存中盖层封闭性评价方法
CN109441422B (zh) * 2018-12-03 2020-06-26 北京科技大学 一种页岩气井间距优化开采方法

Also Published As

Publication number Publication date
CN113703069A (zh) 2021-11-26

Similar Documents

Publication Publication Date Title
CN106979918B (zh) 一种获取致密油藏岩心的液体渗透率的方法及装置
Blunt Pore level modeling of the effects of wettability
JP5797262B2 (ja) 成熟炭化水素産地をシミュレートするための生産シミュレータ
CN111425191B (zh) 一种致密储层气-水相对渗透率计算方法及装置
US20100154514A1 (en) Method of determining the evolution of petrophysical properties of a rock during diagenesis
US11686190B2 (en) 4D quantitative and intelligent diagnosis method and system for spatio-temporal evolution of oil-gas reservoir damage types and extent
US7343275B2 (en) Method for modelling the production of hydrocarbons by a subsurface deposit which are subject to depletion
CN110598167A (zh) 低渗透油藏油水相对渗透率实验数据的处理方法
CN107525746A (zh) 一种表征致密油藏低速非达西渗流特征的方法及装置
CN114153007B (zh) 水锁效应损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN113703069B (zh) 贾敏损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
US8260594B2 (en) Method and system for predicting the apparent flow conductivity of a heterogeneous medium
CN113673111B (zh) 聚合物吸附损害油气层建模方法、损害程度时空演化4d定量与智能诊断技术
Hsu et al. Pore-scale analysis of the effects of contact angle hysteresis on blob mobilisation in a pore doublet
CN113705123B (zh) 外来颗粒损害油气层建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
Clerke et al. Spontaneous imbibition of water into oil saturated M_1 bimodal limestone
CN113657054B (zh) 乳化堵塞损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN113705122B (zh) 润湿反转损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN113655545B (zh) 应力敏感损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN113656982B (zh) 有机垢损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN112377179A (zh) 一种基于分形理论的致密砂岩储层束缚水饱和度预测方法
CA2821004A1 (en) Modeling immiscible two phase flow in a subterranean formation
CN113657053A (zh) 黏土膨胀损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统
CN112326512A (zh) 流体流动弥散的模拟方法
CN114139331B (zh) 无机沉淀损害油气层的建模方法、损害程度时空演化4d定量与智能诊断方法及其系统

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
CB03 Change of inventor or designer information

Inventor after: Jiang Guancheng

Inventor after: Ran Qihua

Inventor after: Liu Xiaobo

Inventor after: Cheng Rongchao

Inventor after: Wang Zenglin

Inventor after: Chen Gang

Inventor after: Dong Tengfei

Inventor after: Li Yizheng

Inventor after: He Kenbo

Inventor after: Yang Lili

Inventor after: Peng Chunyao

Inventor after: Luo Xiaohu

Inventor after: Luo Xuwu

Inventor after: Liang Xing

Inventor after: Tan Bin

Inventor before: Jiang Guancheng

Inventor before: Ran Qifa

Inventor before: Liu Xiaobo

Inventor before: Cheng Rongchao

Inventor before: Wang Zenglin

Inventor before: Chen Gang

Inventor before: Dong Tengfei

Inventor before: Li Yizheng

Inventor before: He Kenbo

Inventor before: Yang Lili

Inventor before: Peng Chunyao

Inventor before: Luo Xiaohu

Inventor before: Luo Xuwu

Inventor before: Liang Xing

Inventor before: Tan Bin

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant