CN108732021B - 一种co2压裂过程的应力场变化测量方法 - Google Patents

一种co2压裂过程的应力场变化测量方法 Download PDF

Info

Publication number
CN108732021B
CN108732021B CN201710252945.6A CN201710252945A CN108732021B CN 108732021 B CN108732021 B CN 108732021B CN 201710252945 A CN201710252945 A CN 201710252945A CN 108732021 B CN108732021 B CN 108732021B
Authority
CN
China
Prior art keywords
dimensional
dimensional physical
fracturing
physical model
states
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
CN201710252945.6A
Other languages
English (en)
Other versions
CN108732021A (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 Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
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 Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN201710252945.6A priority Critical patent/CN108732021B/zh
Priority to US16/317,815 priority patent/US10564080B2/en
Priority to PCT/CN2018/083339 priority patent/WO2018192481A1/zh
Publication of CN108732021A publication Critical patent/CN108732021A/zh
Application granted granted Critical
Publication of CN108732021B publication Critical patent/CN108732021B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N3/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N3/08—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces
    • G01N3/10—Investigating strength properties of solid materials by application of mechanical stress by applying steady tensile or compressive forces generated by pneumatic or hydraulic pressure
    • G01N3/12—Pressure testing
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0001—Type of application of the stress
    • G01N2203/0003—Steady
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0014—Type of force applied
    • G01N2203/0016—Tensile or compressive
    • G01N2203/0019—Compressive
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/003—Generation of the force
    • G01N2203/0042—Pneumatic or hydraulic means
    • G01N2203/0044—Pneumatic means
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/0058—Kind of property studied
    • G01N2203/006—Crack, flaws, fracture or rupture
    • G01N2203/0062—Crack or flaws
    • G01N2203/0066—Propagation of crack
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0676—Force, weight, load, energy, speed or acceleration
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2203/00—Investigating strength properties of solid materials by application of mechanical stress
    • G01N2203/02—Details not specific for a particular testing method
    • G01N2203/06—Indicating or recording means; Sensing means
    • G01N2203/067—Parameter measured for estimating the property
    • G01N2203/0682—Spatial dimension, e.g. length, area, angle

Landscapes

  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提供一种CO2压裂过程的应力场变化测量方法，区别于现有利用数值模拟显示压裂应力场变化的方法，它综合采用CT扫描、数字重构、3D打印、CO2压裂实验、应力冻结及光弹性量测技术，不仅可以透明显示三维物理模型压裂时内部裂缝的空间分布与扩展形态，而且可以获得压裂裂缝扩展过程中内部三维应力相位图，实现固体CO2压裂过程中三维应力场及其演化规律的透明显示与定量表征；同时利用多个相同透明的三维物理模型实现上述显示与表征，确保了测量结果的准确性和可靠性，解决了现有数值模拟方法中存在的计算精度低及准确性和可靠性难以保证的问题。

Description

一种CO2压裂过程的应力场变化测量方法
技术领域
本发明涉及内部应力场测量技术领域，特别涉及一种CO2压裂过程的应力场变化测量方法。
背景技术
由于非常规油气水力压裂过程中不仅对水资源的消耗量巨大，而且压裂液中存在的大量有害添加剂容易造成生态环境的污染，因此，采用超临界CO2压裂便是当前重点发展的一种无水压裂技术。
由于CO2压裂过程的复杂性，水力压裂中的一些基础理论和方法并不完全适用于超临界CO2的压裂研究，而对于超临界CO2的现有研究大都处于对于CO2压裂后形成缝网结构的定性分析上，无法实现其压裂过程的观测和表征。
现有技术中也存在采用数值模拟的方法，来定量分析裂缝扩展过程中应力场的分布演化规律的方案，但是数值模拟的方法计算精度容易受外界条件影响，使得数值模拟结果的准确性和可靠性无法保证。
发明内容
本发明提供一种CO2压裂过程的应力场变化测量方法，以解决现有技术中由于计算精度易受外界影响而导致的数值模拟结果的准确性和可靠性无法保证的问题。
为实现上述目的，本申请提供的技术方案如下：
一种CO2压裂过程的应力场变化测量方法，包括：
CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验；
计算机根据CO2压裂实验后的三维物理模型的CT扫描结果，进行数字重构，获得完整缝网的三维数字模型，并在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型；
3D打印机根据所述不同裂缝状态下的三维数字模型进行打印，得到多个不同裂缝状态下的三维物理模型；
CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验；
温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理；
切片机对进行降温处理之后的各个三维物理模型进行切片，得到各个三维物理模型的三个正交平面二维切片；
计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力；所述应力条纹分布为三个满足实验要求的正交平面二维切片进行二维光弹实验后得到的。
优选的，所述不同的裂缝状态包括：裂缝起裂前、裂缝扩展到预设长度及裂缝停止扩展。
优选的，所述CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验，包括：
CO2泵送***向裂缝起裂前的三维物理模型中注入超临界CO2，直至注入压力较破裂压力低预设压差时，停止注入超临界CO2并维持当前注入压力；
CO2泵送***向裂缝扩展到预设长度的三维物理模型中注入超临界CO2，直至注入压力较形成当前裂缝的压力低预设压差时，停止注入超临界CO2并维持当前注入压力。
优选的，在所述CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验之前，还包括：
计算机对所述真实储层岩心的扫描结果进行数字重构，获得所述真实储层岩心的三维数字模型；
3D打印机根据所述真实储层岩心的三维数字模型进行打印，得到多个相同且透明的所述真实储层岩心的三维物理模型；
温箱对多个所述真实储层岩心的三维物理模型进行处理，使多个所述真实储层岩心的三维物理模型的性质稳定；
三轴加载装置对多个所述真实储层岩心的三维物理模型进行三向伺服加载，直至达到目标值。
优选的，在所述CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验之前，还包括：
温箱对所述不同裂缝状态下的三维物理模型进行处理，使所述不同裂缝状态下的三维物理模型的性质稳定；
三轴加载装置对所述不同裂缝状态下的三维物理模型进行三向伺服加载，直至达到目标值。
优选的，所述温箱对多个所述真实储层岩心的三维物理模型进行处理，使多个所述真实储层岩心的三维物理模型的性质稳定，包括：
温箱对多个所述真实储层岩心的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使多个所述真实储层岩心的三维物理模型的性质稳定。
优选的，所述温箱对所述不同裂缝状态下的三维物理模型进行处理，使所述不同裂缝状态下的三维物理模型的性质稳定，包括：
温箱对所述不同裂缝状态下的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使所述不同裂缝状态下的三维物理模型的性质稳定。
优选的，所述CO2泵送***对真实储层岩心的三维物理模型进行CO2压裂实验之前，还包括：
CO2泵送***的温度控制***调节CO2的注入压力和温度。
优选的，所述温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理，包括：
温箱以2℃/h的速度降温至室温。
优选的，所述计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力，包括：
计算机根据所述应力条纹分布进行相移法提取，得到相应的应力场；
计算机根据公式
Figure BDA0001272595000000031
计算得到载荷作用下三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx；
计算机根据三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx，计算得到三维物理模型内各个点的三维最大剪应力。
本发明提供的所述CO2压裂过程的应力场变化测量方法，通过CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验；然后由计算机根据CO2压裂实验后的三维物理模型的CT扫描结果，进行数字重构，获得完整缝网的三维数字模型，并在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型；再由3D打印机根据所述不同裂缝状态下的三维数字模型进行打印，得到多个不同裂缝状态下的三维物理模型；然后CO2泵送***再对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验；通过温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理；由切片机对进行降温处理之后的各个三维物理模型进行切片，得到各个三维物理模型的三个正交平面二维切片；最后通过计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力。也即，区别于现有技术中的数值模拟的方法，本发明所述CO2压裂过程的应力场变化测量方法，结合数字重构、3D打印技术、CO2压裂实验、温度控制和二维光弹性实验，不仅可以直观显示三维物理模型的不同裂缝状态，而且定量化的记录三维物理模型内部应力相位图，从而实现三维应力场分布的可视化与定量表征；并且，真实反映了非连续结构体的复杂结构的同时，通过多个相同且透明的三维物理模型同时进行上述实验，保证了测量结果的准确性和可靠性，避免了现有技术中由于计算精度易受外界影响而导致的数值模拟结果的准确性和可靠性无法保证的问题。
附图说明
为了更清楚地说明本发明实施例或现有技术内的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述内的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1是本发明实施例提供的CO2压裂过程的应力场变化测量方法的流程图；
图2是本发明另一实施例提供的CO2压裂过程的应力场变化测量方法的另一流程图；
图3是本发明另一实施例提供的CO2压裂过程的应力场变化测量方法的另一流程图；
图4是本发明另一实施例提供的热光曲线图；
图5是本发明另一实施例提供的CO2压裂过程的应力场变化测量方法的另一流程图。
具体实施方式
下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本申请一部分实施例，而不是全部的实施例。基于本申请中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本申请保护的范围。
本发明提供一种CO2压裂过程的应力场变化测量方法，以解决现有技术中由于计算精度易受外界影响而导致的数值模拟结果的准确性和可靠性无法保证的问题。
具体的，该CO2压裂过程的应力场变化测量方法，参见图1，包括：
S101、CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验；
具体的，可以在CO2注入前通过调节气体的压力和温度改变CO2气体的相态，保证压裂流体处于特定的气态或者超临界态。
也即，优选的，CO2泵送***包括温度控制***，用于调节CO2的注入压力和温度。
在CO2压裂实验之前，需要三轴加载装置分别对真实储层岩心的各个三维物理模型进行固定，在真实储层岩心的各个三维物理模型被压裂后，停止超临界CO2的注入，卸载三轴加载装置的外部三向应力，将真实储层岩心的各个三维物理模型从温箱中取出。
S102、计算机根据CO2压裂实验后的三维物理模型的CT扫描结果，进行数字重构，获得完整缝网的三维数字模型，并在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型；
在具体的实际应用中，可以将完整的压裂过程分成三个不同的裂缝扩展阶段，即该不同的裂缝状态包括：裂缝起裂前、裂缝扩展到预设长度及裂缝停止扩展；然后在重构得到的完整缝网模型中截取不同扩展阶段的裂缝形态。
优选的，所述在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型，包括：
基于图像处理技术，对所述CT扫描结果的图像进行区域阈值分割、截取及重构操作，得到裂缝扩展到预设长度的三维数字模型。
该预设长度可以根据具体的应用环境而定，此处不做具体限定，均在本申请的保护范围内。
S103、3D打印机根据所述不同裂缝状态下的三维数字模型进行打印，得到多个不同裂缝状态下的三维物理模型；
根据上述不同裂缝状态的设置，利用3D技术重新制作物理模型，这样就得到裂缝扩展过程中三个不同阶段的物理模型，即没有初始裂缝的完整模型、带有一定长度初始裂缝的模型、裂缝停止扩展的物理模型。
S104、CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验；
值得说明的是，步骤S104中的第二次CO2压裂试验，并不是重复步骤S101中的第一次CO2压裂实验。步骤S101中的第一次CO2压裂实验，其目的是得到真实储层岩心的多个相同且透明的三维物理模型在特定的温度和应力条件下最终的压裂缝网形态。
根据上述内容可知，CO2压裂实际上经历了不同的起裂和扩展过程，才最终形成完整的压裂缝网。而裂缝的起裂和扩展过程是一个高速的动态过程，无法在第一次CO2压裂实验过程中开展冻结实验。因此，为了研究裂缝的不同扩展阶段、压裂裂缝尖端应力场的演化，就需要通过步骤S102，人为的将第一次CO2压裂实验后得到的完整裂缝网络划分成不同的扩展阶段，即将动态的裂缝扩展过程分解为不同的准静态裂缝扩展阶段；并基于图像处理技术，截取重构得到带有一定裂缝长度的不同阶段扩展阶段的三维数字模型。再通过步骤S103利用3D打印技术，制作得到相应的三维物理模型，进而可以通过步骤S104对各个阶段的三维物理模型开展第二次CO2压裂实验，使其各自维持相应的注入压力后进行应力冻结实验，即可得到CO2压裂过程中三维应力场及其演化规律的透明显示与定量表征。
具体的，对于不同裂缝状态下的三维物理模型，可以设置不同的目标压力，优选的，选择如下压力设置：
CO2泵送***向裂缝起裂前的三维物理模型中注入超临界CO2，直至注入压力较破裂压力低预设压差时，停止注入超临界CO2并维持当前注入压力；
CO2泵送***向裂缝扩展到预设长度的三维物理模型中注入超临界CO2，直至注入压力较形成当前裂缝的压力低预设压差时，停止注入超临界CO2并维持当前注入压力。
该预设压差可以根据具体的应用环境进行设定，此处不做具体限定，均在本申请的保护范围内。
压裂过程中的破裂压力，是指裂缝开始起裂时的压力。也就是一次完整的压裂过程中，得到的注入压力-时间曲线的峰值点。一旦注入压力超过破裂压力，裂缝就开始起裂扩展，试件(即进行实验的三维物理模型)发生破坏。本步骤中对裂缝起裂前的三维物理模型的注入压力要低于第一次压裂得到的破裂压力，是因为要保证其处于即将起裂的状态，以便于进行应力冻结实验。
除第一次CO2压裂实验外，后续开展的CO2压裂实验是为了实现应力冻结，因此要保证注入压力和外部围压维持在一定定值，且不会使裂缝继续扩展。因此，后续CO2压裂实验的注入压力以第一次CO2压裂实验得到压力-时间曲线中的应力值作为参考，保证施加的注入压力小于能够使裂缝继续扩展所需的驱动力。
S105、温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理；
优选的，可以设置温箱以2℃/h的速度降温至室温。
然后，即可卸载三轴加载装置的外部三向应力，将进行相应目标压力下的CO2压裂实验后的各个三维物理模型从温箱中取出。
S106、切片机对进行降温处理之后的各个三维物理模型进行切片，得到各个三维物理模型的三个正交平面二维切片；
具体的，可以选取三个进行降温处理之后的三维物理模型，分别沿三个正交平面x-y、y-z、z-x进行切片，得到三维物理模型的三个正交平面二维切片，切片厚度为1mm。
S107、计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力；
所述应力条纹分布为三个满足实验要求的正交平面二维切片进行二维光弹实验后得到的。
在具体的实际应用中，可以先通过对三个正交平面二维切片进行磨平抛光，使其成为满足实验要求的正交平面二维切片，再对其进行二维光弹实验，得到应力条纹分布，然后通过步骤S107得到三维物理模型内各个点的三维最大剪应力。
本实施例提供的所述CO2压裂过程的应力场变化测量方法，通过上述过程，区别于现有的利用数值模拟显示压裂应力场变化的方法，它综合采用CT扫描、数字重构、3D打印、CO2压裂实验、应力冻结及光弹性量测技术，不仅可以透明显示三维物理模型压裂时内部裂缝的空间分布与扩展形态，而且可以获得压裂裂缝扩展过程中内部三维应力相位图，实现固体CO2压裂过程中三维应力场及其演化规律的透明显示与定量表征。并且，本实施例不仅可以真实地反映固体内部复杂的非连续结构，同时，利用多个相同透明的三维物理光弹模型实现了对内部应力场及其演化的透明分析与表征，确保了测量结果的准确性和可靠性，解决了现有数值模拟方法及传统应力实验技术中普遍存在的裂纹扩展计算精度和准确性难以保证、复杂结构模型难以制备和压裂扩展应力场无法连续及透明表征的难题。
本发明另一实施例提供了另外一个具体的CO2压裂过程的应力场变化测量方法，在上述实施例的基础之上，在图1的步骤S101之前，参见图2，还包括：
S201、计算机对所述真实储层岩心的扫描结果进行数字重构，获得所述真实储层岩心的三维数字模型；
S202、3D打印机根据所述真实储层岩心的三维数字模型进行打印，得到多个相同且透明的所述真实储层岩心的三维物理模型；
具体可以采用具有光弹特性的Vero Clear材料打印基质部分，采用Fullcure705材料打印孔隙或裂隙部分。
S203、温箱对多个所述真实储层岩心的三维物理模型进行处理，使多个所述真实储层岩心的三维物理模型的性质稳定；
优选的，温箱对多个所述真实储层岩心的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使多个所述真实储层岩心的三维物理模型的性质稳定。
S204、三轴加载装置对多个所述真实储层岩心的三维物理模型进行三向伺服加载，直至达到目标值。
达到预设的冻结温度后，保持恒温1小时，保证三维物理模型内外均达到统一的冻结温度，然后即可通过步骤S204对三维物理模型进行三向伺服加载达到目标值。
该目标值可以视其具体应用环境而定，此处不做具体限定，均在本申请的保护范围内。
另外，在对所述不同裂缝状态下的三维物理模型进行应力冻结实验时，需要保证在试件达到冻结温度时，所受到的外部应力作用维持恒定，且试件处于线弹性阶段，不能发生破坏，因此，优选的，在图1的步骤S104之前，参见图3，还包括：
S301、温箱对所述不同裂缝状态下的三维物理模型进行处理，使所述不同裂缝状态下的三维物理模型的性质稳定；
优选的，温箱对所述不同裂缝状态下的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使所述不同裂缝状态下的三维物理模型的性质稳定。
S302、三轴加载装置对所述不同裂缝状态下的三维物理模型进行三向伺服加载，直至达到目标值。
达到预设的冻结温度后，保持恒温1小时，保证三维物理模型内外均达到统一的冻结温度，然后即可通过步骤S204对三维物理模型进行三向伺服加载达到目标值。
该目标值可以视其具体应用环境而定，此处不做具体限定，均在本申请的保护范围内。
为使三维物理模型的性质稳定，温箱对上述三维物理模型进行处理，采用60℃作为各个三维物理模型的冻结温度，当然也可以采用其他温度，此处不做具体限定，具体取值可以参见图4所示的三维物理模型的热光曲线，该热光曲线一般通过记录对径受压圆盘在不同温度下的条纹级数变化来得到。
具体可以针对实验中使用的光弹性材料Vero Clear，通过设计不同温度点的反射式光弹实验进行详细测试。利用3D打印技术打印得到
Figure BDA0001272595000000101
的圆盘试件，为消除实验中的反光，可以先对圆盘试件的一侧进行喷涂。然后在高温箱中利用自行设计的连杆装置对圆盘进行对径加载，高温箱可以设置不同的控温曲线程序并留有拍摄窗。实验中按照10℃/h的速度升到目标温度，恒温15分钟左右，利用摄像机拍摄记录下圆盘中心处的条纹级数，根据测得不同温度点的条纹数据绘制Vero Clear的热光曲线，如图4所示。
然后选取热光曲线过渡态阶段的各个温度作为目标温度，测试打印模型处于目标温度时的物理力学性质，本实施例综合热光曲线中不同温度点的变化趋势以及材料物理力学性质，选定过渡态温度60℃作为冻结温度，当然也可以选择其他温度，此处不做具体限定，均在本申请的保护范围内。
本发明另一实施例提供了另外一个具体的CO2压裂过程的应力场变化测量方法，在上述实施例及图1至图3的基础之上，步骤S107如图5所示，包括：
S401、计算机根据所述应力条纹分布进行相移法提取，得到相应的应力场；
S402、计算机根据公式
Figure BDA0001272595000000111
计算得到载荷作用下三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx；
S403、计算机根据三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx，计算得到三维物理模型内各个点的三维最大剪应力。
具体的工作原理与上述实施例相同，此处不再一一赘述。
本发明中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
以上所述，仅是本发明的较佳实施例而已，并非对本发明作任何形式上的限制。虽然本发明已以较佳实施例揭露如上，然而并非用以限定本发明。任何熟悉本领域的技术人员，在不脱离本发明技术方案范围情况下，都可利用上述揭示的方法和技术内容对本发明技术方案做出许多可能的变动和修饰，或修改为等同变化的等效实施例。因此，凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所做的任何简单修改、等同变化及修饰，均仍属于本发明技术方案保护的范围内。

Claims (8)

1.一种CO2压裂过程的应力场变化测量方法，其特征在于，包括：
CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验；
计算机根据CO2压裂实验后的三维物理模型的CT扫描结果，进行数字重构，获得完整缝网的三维数字模型，并在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型；
3D打印机根据所述不同裂缝状态下的三维数字模型进行打印，得到多个不同裂缝状态下的三维物理模型；
CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验；
温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理；
切片机对进行降温处理之后的各个三维物理模型进行切片，得到各个三维物理模型的三个正交平面二维切片；
计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力；所述应力条纹分布为三个满足实验要求的正交平面二维切片进行二维光弹实验后得到的；
所述不同的裂缝状态包括：裂缝起裂前、裂缝扩展到预设长度及裂缝停止扩展；
所述在所述完整缝网的三维数字模型中截取不同的裂缝状态，生成不同裂缝状态下的三维数字模型，包括：
基于图像处理技术，对所述CT扫描结果的图像进行区域阈值分割、截取及重构操作，得到裂缝扩展到预设长度的三维数字模型；
所述CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验，包括：
CO2泵送***向裂缝起裂前的三维物理模型中注入超临界CO2，在注入压力与破裂压力之差达到预设压差时，停止注入超临界CO2并维持当前注入压力；
CO2泵送***向裂缝扩展到预设长度的三维物理模型中注入超临界CO2，在注入压力与形成当前裂缝的压力之差达到预设压差时，停止注入超临界CO2并维持当前注入压力。
2.根据权利要求1所述的CO2压裂过程的应力场变化测量方法，其特征在于，在所述CO2泵送***对真实储层岩心的多个相同且透明的三维物理模型进行CO2压裂实验之前，还包括：
计算机对所述真实储层岩心的扫描结果进行数字重构，获得所述真实储层岩心的三维数字模型；
3D打印机根据所述真实储层岩心的三维数字模型进行打印，得到多个相同且透明的所述真实储层岩心的三维物理模型；
温箱对多个所述真实储层岩心的三维物理模型进行处理，使多个所述真实储层岩心的三维物理模型的性质稳定；
三轴加载装置对多个所述真实储层岩心的三维物理模型进行三向伺服加载，直至达到目标值。
3.根据权利要求1所述的CO2压裂过程的应力场变化测量方法，其特征在于，在所述CO2泵送***对所述不同裂缝状态下的三维物理模型进行相应目标压力下的CO2压裂实验之前，还包括：
温箱对所述不同裂缝状态下的三维物理模型进行处理，使所述不同裂缝状态下的三维物理模型的性质稳定；
三轴加载装置对所述不同裂缝状态下的三维物理模型进行三向伺服加载，直至达到目标值。
4.根据权利要求2所述的CO2压裂过程的应力场变化测量方法，其特征在于，所述温箱对多个所述真实储层岩心的三维物理模型进行处理，使多个所述真实储层岩心的三维物理模型的性质稳定，包括：
温箱对多个所述真实储层岩心的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使多个所述真实储层岩心的三维物理模型的性质稳定。
5.根据权利要求3所述的CO2压裂过程的应力场变化测量方法，其特征在于，所述温箱对所述不同裂缝状态下的三维物理模型进行处理，使所述不同裂缝状态下的三维物理模型的性质稳定，包括：
温箱对所述不同裂缝状态下的三维物理模型进行缓慢升温，直至升温至60℃时，保持恒温1小时，使所述不同裂缝状态下的三维物理模型的性质稳定。
6.根据权利要求1所述的CO2压裂过程的应力场变化测量方法，其特征在于，所述CO2泵送***对真实储层岩心的三维物理模型进行CO2压裂实验之前，还包括：
CO2泵送***的温度控制***调节CO2的注入压力和温度。
7.根据权利要求1所述的CO2压裂过程的应力场变化测量方法，其特征在于，所述温箱对进行相应目标压力下的CO2压裂实验后的各个三维物理模型进行降温处理，包括：
温箱以2℃/h的速度降温至室温。
8.根据权利要求1所述的CO2压裂过程的应力场变化测量方法，其特征在于，所述计算机对应力条纹分布进行处理，得到各个三维物理模型内各个点的三维最大剪应力，包括：
计算机根据所述应力条纹分布进行相移法提取，得到相应的应力场；
计算机根据公式
Figure FDA0004250296550000031
计算得到载荷作用下三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx；
计算机根据三个正交平面内的最大剪应力分布(τmax)xy，(τmax)yz，(τmax)zx，计算得到三维物理模型内各个点的三维最大剪应力。
CN201710252945.6A 2017-04-18 2017-04-18 一种co2压裂过程的应力场变化测量方法 Active CN108732021B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN201710252945.6A CN108732021B (zh) 2017-04-18 2017-04-18 一种co2压裂过程的应力场变化测量方法
US16/317,815 US10564080B2 (en) 2017-04-18 2018-04-17 Method for measuring stress field variations during CO2 fracturing process
PCT/CN2018/083339 WO2018192481A1 (zh) 2017-04-18 2018-04-17 一种co2压裂过程的应力场变化测量方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710252945.6A CN108732021B (zh) 2017-04-18 2017-04-18 一种co2压裂过程的应力场变化测量方法

Publications (2)

Publication Number Publication Date
CN108732021A CN108732021A (zh) 2018-11-02
CN108732021B true CN108732021B (zh) 2023-06-30

Family

ID=63924103

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710252945.6A Active CN108732021B (zh) 2017-04-18 2017-04-18 一种co2压裂过程的应力场变化测量方法

Country Status (1)

Country Link
CN (1) CN108732021B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110243689B (zh) * 2019-06-11 2022-04-08 河海大学 一种硬岩真三轴co2爆燃压裂实验装置及其方法
CN110426405A (zh) * 2019-07-30 2019-11-08 北京建筑大学 一种页岩超临界co2压裂试验分析方法
CN110500090A (zh) * 2019-08-21 2019-11-26 中国石油大学(北京) 立体式体积压裂模拟方法及装置
CN110924933A (zh) * 2019-11-18 2020-03-27 中国石油集团川庆钻探工程有限公司 一种动态模拟页岩压裂缝网可视化实验方法
WO2021134260A1 (zh) * 2019-12-30 2021-07-08 中国矿业大学(北京) 用于识别平面复杂结构模型塑性区的实验***及实验方法
CN111502627A (zh) * 2020-04-20 2020-08-07 华美孚泰油气增产技术服务有限责任公司 一种水力压裂t型裂缝实验装置及其实验方法
GB2595062B (en) * 2020-04-24 2022-05-04 Univ China Mining Experimental method and system for simulating evolution of reservoir fracture stress field

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0209476A2 (fr) * 1985-07-19 1987-01-21 ETAT FRANCAIS représenté par Le Ministère de l'Urbanisme du Logement et des Transports C.E.T.E. Procédé et dispositif de détection d'apparition et de quantification d'évolution de fissures à la surface d'un matériau
US5277062A (en) * 1992-06-11 1994-01-11 Halliburton Company Measuring in situ stress, induced fracture orientation, fracture distribution and spacial orientation of planar rock fabric features using computer tomography imagery of oriented core
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103592186A (zh) * 2013-11-29 2014-02-19 重庆大学 真三轴受力条件下页岩水力压裂损伤演化装置与实验方法
CN104949890A (zh) * 2015-05-26 2015-09-30 中国矿业大学 一种模拟断层滑移的实验方法
CN106289988A (zh) * 2015-05-29 2017-01-04 中国科学院地质与地球物理研究所 超临界二氧化碳岩石压裂试验***
CN106501090A (zh) * 2016-09-26 2017-03-15 中国石油天然气股份有限公司 用于水力压裂模拟实验的裂缝表征方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8412500B2 (en) * 2007-01-29 2013-04-02 Schlumberger Technology Corporation Simulations for hydraulic fracturing treatments and methods of fracturing naturally fractured formation
US20140067269A1 (en) * 2012-08-30 2014-03-06 Schlumberger Technology Corporation Undisturbed fracture orientation
US10788604B2 (en) * 2014-06-25 2020-09-29 Schlumberger Technology Corporation Fracturing and reactivated fracture volumes
US10197489B2 (en) * 2015-07-28 2019-02-05 Chevron U.S.A. Inc. Processes and systems for characterizing and optimizing fracturing fluids

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0209476A2 (fr) * 1985-07-19 1987-01-21 ETAT FRANCAIS représenté par Le Ministère de l'Urbanisme du Logement et des Transports C.E.T.E. Procédé et dispositif de détection d'apparition et de quantification d'évolution de fissures à la surface d'un matériau
US5277062A (en) * 1992-06-11 1994-01-11 Halliburton Company Measuring in situ stress, induced fracture orientation, fracture distribution and spacial orientation of planar rock fabric features using computer tomography imagery of oriented core
CN102852516A (zh) * 2012-04-19 2013-01-02 北京大学 用于油气藏开发的全缝长三维压裂数值模拟的方法和装置
CN103592186A (zh) * 2013-11-29 2014-02-19 重庆大学 真三轴受力条件下页岩水力压裂损伤演化装置与实验方法
CN104949890A (zh) * 2015-05-26 2015-09-30 中国矿业大学 一种模拟断层滑移的实验方法
CN106289988A (zh) * 2015-05-29 2017-01-04 中国科学院地质与地球物理研究所 超临界二氧化碳岩石压裂试验***
CN106501090A (zh) * 2016-09-26 2017-03-15 中国石油天然气股份有限公司 用于水力压裂模拟实验的裂缝表征方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Visual representation and characterization of three-dimensional hydrofracturing cracks within heterogeneous rock through 3D printing and transparent models;Peng Liu;Int J Coal Sci Technol;284-294 *
储层射孔压裂裂缝起裂与扩展的数值分析;赵熙;采矿与安全工程学报;544-550 *
基于三维重构技术的裂缝扩展规律研究;肖雯;;复杂油气藏(04);76-79 *

Also Published As

Publication number Publication date
CN108732021A (zh) 2018-11-02

Similar Documents

Publication Publication Date Title
CN108732021B (zh) 一种co2压裂过程的应力场变化测量方法
CN108072467B (zh) 一种非连续结构体内部应力场的测量方法
Liu et al. Experimental research on the compressive fracture toughness of wing fracture of frozen soil
Wang et al. Co-effects of bedding planes and parallel flaws on fracture evolution in anisotropic rocks
CN104406867B (zh) 基于复型与小时间尺度寿命预测的疲劳裂纹扩展测试方法
Trimonova et al. The Study of the Unstable Fracure Propagation in the Injection Well: Numerical and Laboratory Modeling
CN104297030A (zh) 一种模拟冻土自然形成过程的制样装置及制样方法
CN108872394B (zh) 一种确定室内岩石加卸载过程中内部视应力场的方法及***
Cao et al. Experimental investigation of plane shear fracture characteristics of sandstone after cyclic freeze–thaw treatments
CN107885930A (zh) 继电器弹性金属材料加速贮存退化失效机理变化判别方法
CN106368688A (zh) 研究岩石受热破裂条件和裂缝分布的实验方法
CN105259036A (zh) 地层岩石力学参数的测量方法
McGuire et al. Experimental measurements of stress and chemical controls on the evolution of fracture permeability
CN105181529A (zh) 一种测定凝胶颗粒粘弹性的方法
CN108716395A (zh) 基于冰块压裂的可视化试验方法
CN112649282A (zh) 一种模拟应力扰动条件下诱发断层/裂缝激活的方法
Kim et al. Laboratory measurement of microproppant placement quality using split core plug permeability under stress
Chen et al. 3D morphology and formation mechanism of fractures developed by true triaxial stress
Isaka et al. Investigation of the aperture-dependent flow characteristics of a supercritical carbon dioxide-induced fracture under high-temperature and high-pressure conditions: a numerical study
Ju et al. Characterization of stress field evolution during 3D internal fracture propagation using additively printed models and frozen stress techniques
Fan et al. Spatial failure mode analysis of frozen sandstone under uniaxial compression based on CT technology
CN104076125A (zh) 可重复使用模拟裂缝性漏失地层的钢铁岩心及应用方法
Zhang Fractal dimension of fracture surface in rock material after high temperature
CN108414324A (zh) 一种用于纯ⅱ型动态断裂的试验件及实验方法
Baptista-Pereira et al. Effects of fluid diffusivity on hydraulic fracturing processes using visual analysis

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