考虑动态开裂的孔隙喉道网络模型的驱替模拟方法及装置制造方法及图纸

本申请实施例提供了一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法及装置，该方法包括：获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；驱替模拟所述二维孔隙喉道网络模型的渗流规律；在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率。本申请实施例可以降低非常规油气藏的驱替模拟过程与实际驱替过程的误差。

【技术实现步骤摘要】
考虑动态开裂的孔隙喉道网络模型的驱替模拟方法及装置
本申请涉及非常规油气藏开发研究
，尤其是涉及一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法及装置。
技术介绍
非常规油气藏的开发研究对于我国能源具有战略意义，它们的储层都具有双重孔隙喉道介质的特征，微孔隙分布与动态微孔隙开启、扩展对于渗流的影响研究是解决这些资源开发的关键。对于这些复杂的孔隙结构特征的非常规油气藏的开发问题，很难单纯从宏观的渗流理论去解决，因此，从微观结构特征以及微观渗流机理上开展研究，将是一个新的突破口。目前，在基于微观结构特征以及微观渗流机理的非常规油气藏开发模拟中，所用的孔隙喉道网络模型主要有等效连续体模型、离散网络模型、混合模型以及逾渗模型和孔隙喉道网络模型。其中，基于逾渗理论的孔隙喉道网络模型在常规油气藏开发机理的认识上已体现出一定的优势。然而，由于油藏开发是一个动态的过程，而这些孔隙喉道网络模型并没有考虑动态喉道开裂的过程，其驱替模拟过程与实际驱替过程有一定误差。
技术实现思路
本申请实施例的目的在于提供一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法及装置，以降低非常规油气藏的驱替模拟过程与实际驱替过程的误差。为达到上述目的，一方面，本申请实施例提供了一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，包括：获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；驱替模拟所述二维孔隙喉道网络模型的渗流规律；在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，所述根据所述特征参数构建二维孔隙喉道网络模型，包括：确定用于表征所述特征参数的分布规律的统计分布函数；按照所述统计分布函数生成孔隙喉道网络模型，并将所述特征参数赋值给所述孔隙喉道网络模型中的喉道单元和孔隙单元，从而形成二维孔隙喉道网络模型。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，所述确定用于表征所述特征参数的分布规律的统计分布函数，包括：基于公式Rmin＜R＜Rmax表征孔隙半径的分布规律；基于公式表征喉道半径的分布规律；其中，p(R)为孔隙半径的分布函数，f(x)为喉道半径的分布函数，R为孔隙半径，Rmax为孔隙半径的最大值，Rmin为孔隙半径的最大值，σ为分布函数的标准差，e为自然常数，μ为分布函数的期望值，x为喉道半径。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，所述驱替模拟所述二维孔隙喉道网络模型的渗流规律，包括：根据公式模拟所述二维孔隙喉道网络模型中的喉道单元的渗流规律；根据公式模拟所述二维孔隙喉道网络模型中的裂缝的渗流规律；其中，qij为二维孔隙喉道网络模型中第i行第j列位置处的流量，rij为二维孔隙喉道网络模型中第i行第j列位置处的喉道半径，为流体的动力粘滞系数，Lij为二维孔隙喉道网络模型中至第i行第j列位置处的裂缝长度，Pj为二维孔隙喉道网络模型中第j个节点处的压力，Pi为二维孔隙喉道网络模型中第i个节点处的压力，θ为接触角，σwn为润湿相和非润湿相间的界面张力，r为毛管半径，q为孔隙单元内的流量，e为自然常数，μ为孔隙内流体的动力粘滞系数，ΔP＝Pj-Pi-Pc。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，所述宏观参数包括：孔隙单元、喉道单元及裂缝的侵入状态、毛管压力曲线、相对渗透率曲线、流体饱和度和流体速度。另一方面，本申请实施例还提供了一种考虑动态开裂的孔隙喉道网络模型的驱替模拟装置，包括：模型构建模块，用于获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；渗流模拟模块，用于驱替模拟所述二维孔隙喉道网络模型的渗流规律；模型更新模块，用于在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；参数获取模块，用于在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟装置，所述根据所述特征参数构建二维孔隙喉道网络模型，包括：确定用于表征所述特征参数的分布规律的统计分布函数；按照所述统计分布函数生成孔隙喉道网络模型，并将所述特征参数赋值给所述孔隙喉道网络模型中的喉道单元和孔隙单元，从而形成二维孔隙喉道网络模型。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟装置，所述确定用于表征所述特征参数的分布规律的统计分布函数，包括：基于公式Rmin＜R＜Rmax表征孔隙半径的分布规律；基于公式表征喉道半径的分布规律；其中，p(R)为孔隙半径的分布函数，f(x)为喉道半径的分布函数，R为孔隙半径，Rmax为孔隙半径的最大值，Rmin为孔隙半径的最大值，σ为分布函数的标准差，e为自然常数，μ为分布函数的期望值，x为喉道半径。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟装置，所述驱替模拟所述二维孔隙喉道网络模型的渗流规律，包括：根据公式模拟所述二维孔隙喉道网络模型中的喉道单元的渗流规律；根据公式模拟所述二维孔隙喉道网络模型中的裂缝的渗流规律；其中，qij为二维孔隙喉道网络模型中第i行第j列位置处的流量，rij为二维孔隙喉道网络模型中第i行第j列位置处的喉道半径，为流体的动力粘滞系数，Lij为二维孔隙喉道网络模型中至第i行第j列位置处的裂缝长度，Pj为二维孔隙喉道网络模型中第j个节点处的压力，Pi为二维孔隙喉道网络模型中第i个节点处的压力，θ为接触角，σwn为润湿相和非润湿相间的界面张力，r为毛管半径，q为孔隙单元内的流量，e为自然常数，μ为孔隙内流体的动力粘滞系数，ΔP＝Pj-Pi-Pc。本申请实施例的考虑动态开裂的孔隙喉道网络模型的驱替模拟装置，所述宏观参数包括：孔隙单元、喉道单元及裂缝的侵入状态、毛管压力曲线、相对渗透率曲线、流体饱和度和流体速度。由以上本申请实施例提供的技术方案可见，本申请实施例首先获取孔隙喉道介质的特征参数，并根据特征参数构建二维孔隙喉道网络模型；其次驱替模拟二维孔隙喉道网络模型的渗流规律；然后在模拟渗流规律的过程中，当确认二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新二维孔隙喉道网络模型，然后在模拟渗流规律的过程中，可实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率，由于本申请实施例的孔隙喉道网络模型考虑了非常规油气藏的孔隙在驱替过程动态开裂，因此，本申请实施例的驱替过程模拟与非常规油气藏的实际驱替过程更为接近，从而降低了非常规油气藏的驱替模拟过程与实际驱替过程的误差，这种与非常规油气藏的实际驱替过程更为接近的驱替过程模拟，可以为后续制定非常规油气藏的驱替开发方案提供更为客观准确的参考依据。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面本文档来自技高网...

【技术保护点】
1.一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，包括：获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；驱替模拟所述二维孔隙喉道网络模型的渗流规律；在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率。

【技术特征摘要】
1.一种考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，包括：获取孔隙喉道介质的特征参数，并根据所述特征参数构建二维孔隙喉道网络模型；驱替模拟所述二维孔隙喉道网络模型的渗流规律；在模拟所述渗流规律的过程中，当确认所述二维孔隙喉道网络模型内喉道单元的顶点压力达到临界压力值时，使对应的喉道单元张开而形成裂缝，以更新所述二维孔隙喉道网络模型；在模拟所述渗流规律的过程中，实时获取不同阶段的二维孔隙喉道网络模型的宏观参数及驱油效率。2.如权利要求1所述的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，所述根据所述特征参数构建二维孔隙喉道网络模型，包括：确定用于表征所述特征参数的分布规律的统计分布函数；按照所述统计分布函数生成孔隙喉道网络模型，并将所述特征参数赋值给所述孔隙喉道网络模型中的喉道单元和孔隙单元，从而形成二维孔隙喉道网络模型。3.如权利要求2所述的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，所述确定用于表征所述特征参数的分布规律的统计分布函数，包括：基于公式表征孔隙半径的分布规律；基于公式表征喉道半径的分布规律；其中，p(R)为孔隙半径的分布函数，f(x)为喉道半径的分布函数，R为孔隙半径，Rmax为孔隙半径的最大值，Rmin为孔隙半径的最大值，σ为分布函数的标准差，e为自然常数，μ为分布函数的期望值，x为喉道半径。4.如权利要求1所述的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，所述驱替模拟所述二维孔隙喉道网络模型的渗流规律，包括：根据公式模拟所述二维孔隙喉道网络模型中的喉道单元的渗流规律；根据公式模拟所述二维孔隙喉道网络模型中的裂缝的渗流规律；其中，qij为二维孔隙喉道网络模型中第i行第j列位置处的流量，rij为二维孔隙喉道网络模型中第i行第j列位置处的喉道半径，为流体的动力粘滞系数，Lij为二维孔隙喉道网络模型中至第i行第j列位置处的裂缝长度，Pj为二维孔隙喉道网络模型中第j个节点处的压力，Pi为二维孔隙喉道网络模型中第i个节点处的压力，θ为接触角，σwn为润湿相和非润湿相间的界面张力，r为毛管半径，q为孔隙单元内的流量，e为自然常数，μ为孔隙内流体的动力粘滞系数，ΔP＝Pj-Pi-Pc。5.如权利要求1所述的考虑动态开裂的孔隙喉道网络模型的驱替模拟方法，其特征在于，所述宏观参数包括：孔隙单元、喉道单元及裂缝的侵入状态、毛管压力曲线、相对渗透率曲线、流体饱和度和流体速度。6.一种考虑动态开裂的孔隙喉道网络模...

【专利技术属性】
技术研发人员：吕伟峰，贾宁洪，王智刚，李彤，刘庆杰，杨济如，杨胜建，
申请(专利权)人：中国石油天然气股份有限公司，
类型：发明
国别省市：北京,11