CN114927176A - 考虑时间效应的含水合物沉积物蠕变本构模型构建方法 - Google Patents

考虑时间效应的含水合物沉积物蠕变本构模型构建方法 Download PDF

Info

Publication number
CN114927176A
CN114927176A CN202210596664.3A CN202210596664A CN114927176A CN 114927176 A CN114927176 A CN 114927176A CN 202210596664 A CN202210596664 A CN 202210596664A CN 114927176 A CN114927176 A CN 114927176A
Authority
CN
China
Prior art keywords
creep
hydrate
stage
model
loading
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.)
Granted
Application number
CN202210596664.3A
Other languages
English (en)
Other versions
CN114927176B (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.)
Qingdao Institute of Marine Geology
Original Assignee
Qingdao Institute of Marine Geology
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 Qingdao Institute of Marine Geology filed Critical Qingdao Institute of Marine Geology
Priority to CN202210596664.3A priority Critical patent/CN114927176B/zh
Publication of CN114927176A publication Critical patent/CN114927176A/zh
Application granted granted Critical
Publication of CN114927176B publication Critical patent/CN114927176B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
    • 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/02—Details
    • G01N3/06—Special adaptations of indicating or recording means
    • 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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
    • 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/0058—Kind of property studied
    • G01N2203/0069—Fatigue, creep, strain-stress relations or elastic constants
    • G01N2203/0071—Creep
    • 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/0658—Indicating or recording means; Sensing means using acoustic or ultrasonic detectors
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computing Systems (AREA)
  • Pure & Applied Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Software Systems (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • General Engineering & Computer Science (AREA)
  • Operations Research (AREA)
  • Algebra (AREA)
  • Chemical & Material Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Pathology (AREA)
  • Immunology (AREA)
  • Biochemistry (AREA)
  • Analytical Chemistry (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Computational Linguistics (AREA)
  • Evolutionary Computation (AREA)
  • Molecular Biology (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Artificial Intelligence (AREA)
  • Probability & Statistics with Applications (AREA)
  • Evolutionary Biology (AREA)
  • Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)

Abstract

本发明公开一种考虑时间效应的含水合物沉积物蠕变本构模型构建方法，通过含水合物沉积物分级加载蠕变实验，获取试样蠕变应变时程曲线和超声波速变化数据，以超声波速变化表征加载过程的蠕变损伤，并考虑加载过程中有效轴向载荷变化。形成一种含水合物沉积物蠕变本构模型构建方法。本发明方案考虑时间效应对含水合物沉积物储层变形的影响，真实反映含水合物沉积物储层实际失稳特征，有效提高水合物试采井壁稳定分析计算精度，为现场水合物试采作业提供理论支撑。

Description

考虑时间效应的含水合物沉积物蠕变本构模型构建方法
技术领域
本发明属于海洋天然气水合物沉积物本构模型技术领域，具体涉及一种考虑时间效应的含水合物沉积物蠕变本构模型构建方法。
背景技术
天然气属于低碳清洁化石能源，大力发展天然气产业，提高天然气在一次能源消费结构中的占比是实现“碳达峰、碳中和”的重要举措。随着天然气消费量的日益增加，面对巨大的天然气能源的增长需求，常规天然气资源开发已不能满足，大力发展非常规天然气资源具有重大意义。非常规天然气资源包括致密气、页岩气、煤层气、天然气水合物等。
海洋水合物资源开发过程中，需要预先在水合物储层建立井眼，再采用降压法、热激发、CO2置换和固态流化等方法进行开采。钻井与试采过程中，井壁附近水合物沉积物储层受到外力作用产生变形，随外力的持续作用岩石的变形不断增加，最终将导致井壁失稳，引起地层垮塌，海底滑坡等地质灾害，保证井壁稳定与地层稳定对水合物安全开发和环境保护具有重要意义。
对于水合物储层井壁稳定分析，通常采用理论推导或数值模拟等方法进行，这些方法以水合物沉积物本构模型为基础，耦合井壁围岩应力，进行储层变形计算。目前研究中水合物沉积物本构模型均采用弹性或弹塑性模型，忽略长期载荷作用下蠕变效应的影响。实际工程中水合物储层在恒定应力或应变的长期作用下，会产生蠕变破坏。建立水合物沉积物蠕变本构模型，进行井壁稳定计算，能够更真实准确反映岩石破坏的实际特征。因此，亟需一种考虑时间效应，能反应储层蠕变影响的含水合物沉积物本构模型，进行井壁稳定，地层稳定分析计算，保障水合物试采安全。
发明内容
本发明为解决现有技术中存在的缺陷，提出一种考虑时间效应的含水合物沉积物蠕变本构模型构建方法，考虑时间效应对含水合物沉积物储层变形的影响，真实反映含水合物沉积物储层实际失稳特征，以提高水合物试采井壁稳定分析计算精度，为现场水合物试采作业提供理论支撑。
本发明是采用以下的技术方案实现的：一种考虑时间效应的含水合物沉积物蠕变本构模型构建方法，包括以下步骤：
步骤A、对含水合物沉积物开展分级加载蠕变实验，获取分级加载应变时程曲线和超声波速变化曲线，，并将蠕变应变时程曲线分为稳态蠕变阶段、长期强度阶段和非稳态蠕变阶段；
步骤B、基于声波速度变化和加载过程中轴向有效应变力变化，建立含水合物沉积物分级加载蠕变本构模型；
步骤B1、构建稳态蠕变阶段的蠕变方程：
基于广义开尔文体模型，构建稳态蠕变阶段的本构方程，广义开尔文体模型由胡克体和开尔文体串联构成，对稳态蠕变阶段的本构方程进行求解，得到稳态蠕变阶段的蠕变方程；
步骤B2、稳态长期强度阶段的蠕变方程
(1)对长期强度阶段应变时程曲线描述，在广义开尔文体模型的基础上串联一个理想粘塑体模型，变形为西原模型，构建长期强度阶段的本构方程，并对其进行求解得到初步的长期强度阶段的蠕变方程；
(2)蠕变实验稳定加载过程中，轴向有效载荷随轴向位移的增加逐渐减小，考虑该因素对长期强度阶段西原模型进行改进，将其中理想粘塑体模型中摩擦片改为圣维南体，构建考虑轴向有效载荷变化的长期强度阶段本构方程，对初步的长期强度阶段的蠕变方程进行修正；
步骤B3、对于非稳态蠕变阶段，在考虑轴向有效载荷变化影响的西原模型的基础上，再串联一个非线性粘塑体模型，进而得到非稳态蠕变阶段的蠕变方程；
步骤B4、基于超声速度变化，获得声波修正系数对步骤B1、B2和B3三个阶段的蠕变方程修正，得到考虑声速变化和实验中有效载荷变化的含天然气水合物沉积物全阶段蠕变本构模型。
进一步的，所述步骤A中，含水合物沉积物分级加载蠕变实验具体包括以下步骤：
步骤A1、水合物生成：
(1)含水合物沉积物试样制备：含水合物沉积物基质为砂或泥质粉砂，与水均匀混合装填进反应釜腔室；
(2)含水合物沉积物蠕变设备组装：反应釜腔室两端安装超声传感器，封闭反应釜，连接管线并测压检漏，进行数据采集连接；
(3)含水合物沉积物生成：加载预设围压和孔压，调节至设定温度，水合物随即开始生成；
步骤A2、含水合物沉积物分级加载蠕变实验：
(1)预加载：施加轴向初始载荷，使轴向载荷与有效围压达到平衡，持续一定时间，待压力/位移数据稳定后再进行下一步加载；
(2)加载：加载至设置载荷，保持围压孔压稳定，开始进行分级加载蠕变实验，记录轴向位移数据并采集超声速度数据；当前一级蠕变变形趋于稳定，进行下一级加载，直至试样破坏，实验结束；
(3)实验过程中，获取分级加载蠕变应变时程曲线和超声速度随加载进行变化曲线。
进一步的，所述步骤A1中，实验管路中的循环介质采用去离子水或纯水，并加入0.03％的十二烷基硫酸钠促进水合物生成，水合物生成时间大于100h。
进一步的，所述步骤A2中，每级加载时间大于24h，当蠕变速度小于0.0015mm/h时，进行下一级加载。
进一步的，所述步骤B1中，稳态阶段的本构方程为：
Figure BDA0003668212300000031
对上述本构方程进行求解，得到
Figure BDA0003668212300000032
其中，E0、E1为胡克体弹性模量，σ为轴向应力，σy为长期强度，ε为轴向应变，η1为牛顿粘性系数。
进一步的，所述步骤B2中，长期强度阶段的本构方程为：
Figure BDA0003668212300000033
对上式进行求解，进而得到初步的长期强度阶段的蠕变方程：
Figure BDA0003668212300000034
其中，σy为长期强度，η2为牛顿粘性系数；t为加载时间。
进一步的，所述步骤B2中，考虑轴向有效载荷变化影响，修正后的长期强度阶段蠕变方程为：
Figure BDA0003668212300000035
进一步的，所述步骤B3中，考虑实验中轴向有效载荷变化，得到非稳态蠕变阶段的蠕变方程如下：
Figure BDA0003668212300000036
其中，σf为破坏强度，n为非线性参数。
进一步的，所述步骤B4中，得到的全阶段蠕变本构模型为：
Figure BDA0003668212300000041
声波修正系数D表达式为：
Figure BDA0003668212300000042
其中，V初为加载初始声速；V末为加载结束声速。
进一步的，对于所得到的全阶段蠕变本构模型，采用最小二乘法、回归分析法或粒子群优化算法对含水合物沉积物分级加载实验蠕变应变时程曲线进行拟合，通过不断调整模型参数，使模型拟合曲线与试验曲线尽可能重合，将拟合得到的最优参数作为含水合物沉积物稳态蠕变阶段基础参数。
与现有技术相比，本发明的优点和积极效果在于：
(1)本方案以西原模型为基础，将其中理想粘塑体中胡克体改为圣维南体，并增加了非线性粘塑体模型，建立改进的西原模型，构建考虑时间效应影响的含水合物沉积物全阶段蠕变本构模型，更为真实的反映含水合物沉积物的实际失稳特征；
(2)根据声速变化，引入声速修正系数，对蠕变本构模型中胡克元件弹性模量进行修正，使模型可以同时兼顾考虑材料弹性模量增强与材料损伤的影响，反映加载过程中蠕变弹性参数变化。声速变化与弹性模量变化直接相关，故采用声速对弾性模量进行修正。
(3)本发明中考虑加载过程中应变引起的有效轴向载荷变化的影响，对蠕变模型参数进行修正，更为准确的表征全阶段变形特征，提高蠕变模型准确性。
附图说明
图1是本发明实施例的考虑时间效应的含水合物沉积物蠕变本构模型构建方法流程示意图；
图2是本发明实施例的含水合物沉积物蠕变曲线各变形阶段示意图；
图3是本发明实施例的常规西原模型示意图；
图4是本发明实施例的改进的西原模型示意图；
图5是本发明实施例的含水合物沉积物分级加载实验蠕变应变时程曲线；
图6是本发明实施例的含水合物沉积物蠕变本构模型计算值与实验值对比示意图；
其中：6、胡克体；7、开尔文体；8、理想粘塑体；9、圣维南体；10、非线性粘塑体。
具体实施方式
为了能够更加清楚地理解本发明的上述目的、特征和优点，下面结合附图及实施例对本发明做进一步说明。在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用不同于在此描述的其他方式来实施，因此，本发明并不限于下面公开的具体实施例。
本发明实施例提出一种考虑时间效应的含水合物沉积物蠕变本构模型构建方法，流程示意图如图1所示，包括以下步骤：
步骤A、对含水合物沉积物开展分级加载蠕变实验，获取分级加载应变时程曲线和超声波速变化曲线；
步骤B、基于声波速度变化和加载过程中轴向有效应变力变化，建立含水合物沉积物分级加载蠕变本构模型。
本实施例中，在步骤A中，含水合物沉积物分级加载蠕变实验可以采用三轴应力条件下含水合物沉积物渗透率测量装置及其方法(ZL2017110703083)进行，并在该装置反应釜腔室两端增设超声传感器，测量蠕变过程中声速变化。含水合物沉积物分级加载蠕变实验包括两部分，首先进行水合物生成，再进行分级加载蠕变实验，具体包括以下步骤：
步骤A1、水合物生成：
(1)含水合物沉积物试样制备：含水合物沉积物基质为砂或泥质粉砂，去离子水中加入0.03％的十二烷基硫酸钠(SDS)促进水合物生成，二者均匀混合，装填进反应釜腔室；
(2)含水合物沉积物蠕变设备组装：反应釜腔室两端安装超声传感器，封闭反应釜，反应釜放置于固定支架上，连接管线并测压检漏，进行数据采集连接；
(3)含水合物沉积物生成：加载预设围压和孔压，调节至设定温度，水合物随即开始生成。
其中，需要注意的是：水合物生成时间大于100h，保证试样中水完全生成水合物，蠕变实验过程中水合物饱和度保持稳定，孔压保持稳定；实验管路中需使用去离子水或纯水，避免管路堵塞。
步骤A2、含水合物沉积物分级加载蠕变实验：
(1)预加载：施加轴向初始载荷，使轴向载荷与有效围压达到平衡，加载持续24h以上，待压力/位移数据稳定后再进行下一步加载；
(2)加载：加载至设置载荷，保持围压孔压稳定，开始进行分级加载蠕变实验，记录轴向位移数据并采集超声速度数据。当前一级蠕变变形趋于稳定，进行下一级加载，直至试样破坏，实验结束；
其中，需要注意的是，为保证充分考虑时间效应的影响，每级加载需大于24h，当蠕变速度小于0.0015mm/h，进行下一级加载。
(3)实验过程中，获取分级加载蠕变应变时程曲线和超声速度随加载进行变化曲线。
本实施例采用稳态蠕变应变速率法进行含水合物沉积物长期强度判断，以长期强度为分界，将蠕变应变时程曲线分为三部分，轴向载荷小于长期强度为稳态蠕变阶段，轴向载荷大于长期强度小于破坏强度为长期强度阶段，轴向载荷大于破坏强度为非稳态蠕变阶段。如图2所示，为含水合物沉积物蠕变曲线各变形阶段示意图，曲线ABC表示稳态蠕变阶段，曲线DEFG为长期强度阶段，曲线HIJK为非稳态蠕变阶段。稳态蠕变阶段随时间增加曲线增长逐渐减小，直至达到稳定应变；长期强度阶段曲线斜率先较小，后保持一段直线变化，最终，由于轴向有效载荷小于长期强度曲线趋于稳定，回归稳态蠕变阶段；非稳态蠕变阶段，加载一段时间后，曲线斜率逐渐增大，试样宏观破坏。
步骤B、基于步骤A所获得的分级加载应变时程曲线和超声波速变化曲线，建立含水合物沉积物分级加载蠕变本构模型；
步骤B1、对稳态蠕变阶段应变时程曲线进行拟合：
基于广义开尔文体模型，构建本构方程，广义开尔文体模型由胡克体和开尔文体串联构成，本构方程为：
Figure BDA0003668212300000061
式中：E0、E1为胡克体弹性模量，单位MPa；η1为牛顿粘性系数，单位MPa；σ为轴向应力，单位MPa；ε为轴向应变，量纲为1，σy为长期强度；单位MPa。
对上式进行求解，得到稳态蠕变阶段蠕变方程表达式为：
Figure BDA0003668212300000062
步骤B2、对长期强度阶段应变时程曲线进行拟合：
长期强度阶段应变时程曲线描述，在广义开尔文体模型的基础上串联一个理想粘塑体模型，变形为西原模型，本构方程为：
Figure BDA0003668212300000063
对上式进行求解，得到初步的长期强度阶段蠕变方程表达式为：
Figure BDA0003668212300000064
其中，σy为长期强度，η2为牛顿粘性系数，单位MPa；t为加载时间，单位h；
如图3所示，为西原模型示意图，西原模型由广义开尔文体模型串联一个理想粘塑体模型组成，可以描述试样稳态蠕变和长期强度阶段中的蠕变应变特征。
考虑实验中轴向有效载荷变化，加载过程中应变引起的有效轴向载荷变化对水合物沉积物蠕变本构模型表征过程中会带来很大的误差。蠕变实验稳定加载过程中，轴向有效载荷随轴向位移的增加逐渐减小，考虑这种变化的影响，对长期强度阶段西原模型进行改进，将其中理想粘塑体模型中摩擦片改为圣维南体，构建考虑轴向有效载荷变化的长期强度阶段本构方程，对初步的长期强度阶段的蠕变方程进行修正，修正后的长期强度阶段的蠕变方程表达式变为：
Figure BDA0003668212300000071
步骤B3、对非稳态蠕变阶段应变时程曲线进行拟合
对于非稳态蠕变阶段，在考虑轴向有效载荷变化影响的西原模型的基础上，再串联一个非线性粘塑体模型，得到非稳态蠕变阶段蠕变方程表达式为：
Figure BDA0003668212300000072
式中：σf为破坏强度；单位MPa；n为非线性参数。
图4为改进西原模型示意图，模型在西原模型的基础上，将理想粘塑体中摩擦片改变为圣维南体，用于表征轴向有效载荷变化产生的影响，增加一个非线性粘塑体，用于表征非稳态蠕变阶段蠕变应变特征。
步骤B4、得到含天然气水合物沉积物全阶段蠕变本构模型
蠕变加载过程中轴向载荷对试样内部弹性参数变化会产生影响，常规岩土材料本构模型建立时采用声发射数进行表征，但是声发射数仅能表征材料损伤，不能表征材料属性增强。本发明基于超声速度变化对上述三个阶段的胡克体元件弹性模量进行修正，修正后的考虑声速变化和实验中有效载荷变化的含天然气水合物沉积物全阶段蠕变本构模型表达式为：
Figure BDA0003668212300000073
声波修正系数D计算表达式为：
Figure BDA0003668212300000074
式中：V初为加载初始声速；单位为m/s；V末为加载结束声速；单位为m/s。
基于上述模型(式7)，采用最小二乘法、回归分析法或粒子群优化算法对含水合物沉积物分级加载实验蠕变应变时程曲线进行拟合；通过不断调整模型参数，使模型拟合曲线与试验曲线尽可能重合，将拟合得到的最优参数作为含水合物沉积物稳态蠕变阶段基础参数。具体为：
本实施例中，弹性模量E0，弹性模量E1，粘性系数η1，通过稳态蠕变阶段曲线拟合获得；σy为考虑轴向有效载荷变化的长期强度，即长期强度阶段中近直线段终点位置对应的轴向有效载荷。弹性模量E2，粘性系数η2，通过长期强度阶段曲线拟合获得；粘性系数η3和系数n，通过非稳态蠕变阶段曲线拟合获得。
本实施例中，采用0.15～0.3mm粒径的砂(孔隙度为38.33％)作为松散沉积物，水合物生成所需气体是纯度为99.9％的甲烷气体。本实验中含水合物沉积物饱和度为30％，围压为7.1MPa，孔压6.1MPa。分级加载路径为2.4，2.6，2.8，3.0。长期强度阶段直线上升拐点位置，有效轴向载荷为2.732MPa，蠕变模型中长期强度σy为2.732MPa，非稳态蠕变曲线拐点作为破坏强度为2.81MPa，各阶段超声波波速与对应修正系数如表1所示。
表1各阶段超声波波速与对应声波修正系数表
轴向载荷 V<sub>初</sub>/(m/s) V<sub>末</sub>/(m/s) D
2.4 2776 2841 0.0234
2.6 2841 2847.4 0.00225
2.8 2847.4 2832 0.00544
3.0 2832 2785.4 0.0167
结合图5和图6所示，本发明实施例中对含水合物沉积物分级加载实验蠕变应变时程曲线进行拟合，得到弹性模量，粘性系数分别如表2所示。
表2弹性模量和粘性系数表
轴向载荷 E<sub>0</sub>/MPa E<sub>1</sub>/MPa E<sub>2</sub>/MPa η<sub>1</sub>/(MPa.h) η<sub>2</sub>/(MPa.h) η<sub>3</sub>/(MPa.h) n R<sup>2</sup>/％
2.4 514 8.7 / 209.8 / / / 96.02
2.6 514 31.8 / 321.8 / / / 98.06
2.8 514 56.7 4.5 432.8 4.7 / / 98.48
3.0 514 81.6 4.5 545.9 9.9 9.5e5 4.2 95.69
如图6所示，本发明实施例中含水合物沉积物蠕变本构模型计算值与实验值对比示意图，图中散点为实验值，实线为模型计算值，各阶段曲线拟合误差系数均大于95％，所提出的含水合物沉积物分级加载蠕变本构模型构建方法合理准确。
以上所述，仅是本发明的较佳实施例而已，并非是对本发明作其它形式的限制，任何熟悉本专业的技术人员可能利用上述揭示的技术内容加以变更或改型为等同变化的等效实施例应用于其它领域，但是凡是未脱离本发明技术方案内容，依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化与改型，仍属于本发明技术方案的保护范围。

Claims (10)

1.考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于，包括以下步骤：
步骤A、对含水合物沉积物开展分级加载蠕变实验，获取分级加载应变时程曲线和超声波速变化曲线，并将蠕变应变时程曲线分为稳态蠕变阶段、长期强度阶段和非稳态蠕变阶段；
步骤B、基于声波速度变化和加载过程中轴向有效应变力变化，建立含水合物沉积物分级加载蠕变本构模型；
步骤B1、构建稳态蠕变阶段的蠕变方程：
基于广义开尔文体模型，构建稳态蠕变阶段的本构方程，广义开尔文体模型由胡克体和开尔文体串联构成，对稳态蠕变阶段的本构方程进行求解，得到稳态蠕变阶段的蠕变方程；
步骤B2、构建长期强度阶段的蠕变方程：
(1)对长期强度阶段应变时程曲线描述，在广义开尔文体模型的基础上串联一个理想粘塑体模型，变形为西原模型，构建长期强度阶段的本构方程，并对其进行求解得到初步的长期强度阶段的蠕变方程；
(2)蠕变实验稳定加载过程中，轴向有效载荷随轴向位移的增加逐渐减小，考虑该因素对长期强度阶段的西原模型进行改进，将其中理想粘塑体模型中摩擦片改为圣维南体，构建考虑轴向有效载荷变化的长期强度阶段本构方程，对初步的长期强度阶段的蠕变方程进行修正；
步骤B3、构建非稳态蠕变阶段的蠕变方程：
对于非稳态蠕变阶段，在考虑轴向有效载荷变化影响的西原模型的基础上，再串联一个非线性粘塑体模型，进而得到非稳态蠕变阶段的蠕变方程；
步骤B4、基于超声速度变化，获得声波修正系数对步骤B1、B2和B3三个阶段的蠕变方程修正，得到考虑声速变化和实验中有效载荷变化的含天然气水合物沉积物全阶段蠕变本构模型。
2.根据权利要求1所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤A中，含水合物沉积物分级加载蠕变实验具体包括以下步骤：
步骤A1、水合物生成：
(1)含水合物沉积物试样制备：含水合物沉积物基质为砂或泥质粉砂，与循环介质均匀混合装填进反应釜腔室；
(2)含水合物沉积物蠕变设备组装：反应釜腔室两端安装超声传感器，封闭反应釜，连接管线并测压检漏，进行数据采集连接；
(3)含水合物沉积物生成：加载预设围压和孔压，调节至设定温度，水合物随即开始生成；
步骤A2、含水合物沉积物分级加载蠕变实验：
(1)预加载：施加轴向初始载荷，使轴向载荷与有效围压达到平衡，持续一定时间，待压力/位移数据稳定后再进行下一步加载；
(2)加载：加载至设置载荷，保持围压孔压稳定，开始进行分级加载蠕变实验，记录轴向位移数据并采集超声速度数据；当前一级蠕变变形趋于稳定，进行下一级加载，直至试样破坏，实验结束；
(3)实验过程中，获取分级加载蠕变应变时程曲线和超声速度随加载进行变化曲线。
3.根据权利要求2所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤A1中，实验管路中的循环介质采用去离子水或纯水，并加入0.03％的十二烷基硫酸钠促进水合物生成，水合物生成时间大于100h。
4.根据权利要求2所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤A2中，每级加载时间大于24h，当蠕变速度小于0.0015mm/h时，进行下一级加载。
5.根据权利要求1所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤B1中，稳态阶段的本构方程为：
Figure FDA0003668212290000021
对上述本构方程进行求解，得到
Figure FDA0003668212290000022
其中，E0、E1为胡克体弹性模量，σ为轴向应力，σy为长期强度，ε为轴向应变，η1为牛顿粘性系数。
6.根据权利要求5所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤B2中，长期强度阶段的本构方程为：
Figure FDA0003668212290000023
对上式进行求解，进而得到初步的长期强度阶段的蠕变方程：
Figure FDA0003668212290000024
其中，σy为长期强度，η2为牛顿粘性系数；t为加载时间。
7.根据权利要求6所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤B2中，考虑轴向有效载荷变化影响的，对初步的长期强度阶段的蠕变方程进行修正，修正后的长期强度阶段蠕变方程为：
Figure FDA0003668212290000031
8.根据权利要求7所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤B3中，考虑实验中轴向有效载荷变化，得到非稳态蠕变阶段的蠕变方程如下：
Figure FDA0003668212290000032
其中，σf为破坏强度，n为非线性参数。
9.根据权利要求8所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：所述步骤B4中，得到的全阶段蠕变本构模型为：
Figure FDA0003668212290000033
声波修正系数D表达式为：
Figure FDA0003668212290000034
其中，V初为加载初始声速；V末为加载结束声速。
10.根据权利要求9所述的考虑时间效应的含水合物沉积物蠕变本构模型构建方法，其特征在于：对于所得到的全阶段蠕变本构模型，采用最小二乘法、回归分析法或粒子群优化算法对含水合物沉积物分级加载实验蠕变应变时程曲线进行拟合，通过不断调整模型参数，使模型拟合曲线与试验曲线尽可能重合，将拟合得到的最优参数作为含水合物沉积物稳态蠕变阶段基础参数。
CN202210596664.3A 2022-05-30 2022-05-30 考虑时间效应的含水合物沉积物蠕变本构模型构建方法 Active CN114927176B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210596664.3A CN114927176B (zh) 2022-05-30 2022-05-30 考虑时间效应的含水合物沉积物蠕变本构模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210596664.3A CN114927176B (zh) 2022-05-30 2022-05-30 考虑时间效应的含水合物沉积物蠕变本构模型构建方法

Publications (2)

Publication Number Publication Date
CN114927176A true CN114927176A (zh) 2022-08-19
CN114927176B CN114927176B (zh) 2024-04-16

Family

ID=82812059

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210596664.3A Active CN114927176B (zh) 2022-05-30 2022-05-30 考虑时间效应的含水合物沉积物蠕变本构模型构建方法

Country Status (1)

Country Link
CN (1) CN114927176B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106370526A (zh) * 2016-10-26 2017-02-01 华北科技学院 岩石流变扰动效应三轴实验设备及其实验方法
CN112834357A (zh) * 2021-01-08 2021-05-25 青岛海洋地质研究所 海底天然气水合物沉积物储层旁压蠕变试验系统及方法
CN113919148A (zh) * 2021-09-29 2022-01-11 中南大学 一种岩石非线性蠕变模型建立方法
CN114547908A (zh) * 2022-03-02 2022-05-27 大连海事大学 一种热力耦合石灰岩分数阶蠕变损伤本构模型构建方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106370526A (zh) * 2016-10-26 2017-02-01 华北科技学院 岩石流变扰动效应三轴实验设备及其实验方法
CN112834357A (zh) * 2021-01-08 2021-05-25 青岛海洋地质研究所 海底天然气水合物沉积物储层旁压蠕变试验系统及方法
CN113919148A (zh) * 2021-09-29 2022-01-11 中南大学 一种岩石非线性蠕变模型建立方法
CN114547908A (zh) * 2022-03-02 2022-05-27 大连海事大学 一种热力耦合石灰岩分数阶蠕变损伤本构模型构建方法

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
MINGTAO CHEN ET AL.: "Review on the Test Methods and Devices for Mechanical Properties of Hydrate-Bearing Sediments", 《RECENT ADVANCES IN MARINE ENGINEERING GEOLOGY》, vol. 14, no. 10, 20 May 2022 (2022-05-20) *
吴能友等: "记笔记 海洋天然气水合物储层蠕变行为的主控因素与研究展望", 《海洋地质与第四纪地质》, vol. 41, no. 05, 28 October 2021 (2021-10-28), pages 3 - 11 *
杨秀荣;姜谙男;江宗斌;: "含水状态下软岩蠕变试验及损伤模型研究", 岩土力学, no. 1, 10 July 2018 (2018-07-10), pages 167 - 174 *
武东生;孟陆波;李天斌;赖林;: "灰岩三轴高温后效流变特性及长期强度研究", 岩土力学, no. 1, 10 June 2016 (2016-06-10), pages 183 - 191 *
牛双建;党元恒;冯文林;杨大方;: "损伤破裂砂岩单轴蠕变特性试验研究", 岩土力学, no. 05, 10 May 2016 (2016-05-10), pages 1249 - 1258 *

Also Published As

Publication number Publication date
CN114927176B (zh) 2024-04-16

Similar Documents

Publication Publication Date Title
CN111502652B (zh) 一种三孔介质气藏水平井产量递减及生产动态预测方法
CN112036098A (zh) 一种深层油气藏水力裂缝扩展数值模拟的方法
CN109580388B (zh) 一种岩土材料剪切屈服面与体积屈服面的测定方法
CN102509179A (zh) 气井动态产能预测方法
CN108843303B (zh) 一种基于泥岩蠕变模型的油水井套管损坏预测方法
CN110929462A (zh) 一种使用生产数据计算低渗气藏真实压力及气藏储量的方法
CN114152554B (zh) 水力剪切刺激干热岩储层增透模拟实验系统与实验方法
CN116401966A (zh) 一种多孔裂隙岩体时变浆液扩散过程的模拟方法
CN104374684A (zh) 用于测试采动过程卸荷煤岩体渗透率的系统及其应用
CN115587674A (zh) 油藏改建储气库扩容达产过程气井动态产能预测方法
CN114927176A (zh) 考虑时间效应的含水合物沉积物蠕变本构模型构建方法
CN113919201A (zh) 一种水力压裂裂缝多尺度扩展的网格自适应方法
Gordeliy et al. Characteristic features of salt-cavern behavior
CN113123783B (zh) 一种精细评价地层体积含量的方法
CN113294147B (zh) 一种考虑重力因素影响的单洞型断溶体储层试井解释方法
CN111666665A (zh) 一种提高碳酸岩动静态杨氏模量转换模型精度的方法
CN116205163A (zh) 一种天然气水合物藏数值模拟中的地层稳定性分析方法
CN114611306A (zh) 一种基于压力差转换系数的单井控制储量计算方法
CN105389439A (zh) 天然气关井井筒建模模拟实现方法
CN115452681A (zh) 一种储层油气相对渗透率确定方法
CN108074282B (zh) 基于三维地质建模软件的物性参数的厚度校正方法
CN110671102B (zh) 一种气井临界出砂压差的确定方法及系统
Ng et al. A state-dependent constitutive model for gas hydrate-bearing sediment incorporating phase change
CN115822533B (zh) 一种底水砂岩油藏驱油方法
CN117784239B (zh) 过去甲烷水合物稳定带底界面位置的判断方法

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