CN109557280A - 一种储层应力敏感性的确定方法 - Google Patents

一种储层应力敏感性的确定方法 Download PDF

Info

Publication number
CN109557280A
CN109557280A CN201710891799.1A CN201710891799A CN109557280A CN 109557280 A CN109557280 A CN 109557280A CN 201710891799 A CN201710891799 A CN 201710891799A CN 109557280 A CN109557280 A CN 109557280A
Authority
CN
China
Prior art keywords
reservoir
pressure
net impact
permeability
stress
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.)
Pending
Application number
CN201710891799.1A
Other languages
English (en)
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 Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
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 Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201710891799.1A priority Critical patent/CN109557280A/zh
Publication of CN109557280A publication Critical patent/CN109557280A/zh
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials

Abstract

本发明涉及一种储层应力敏感性的确定方法,属于油气田开发技术领域。本发明根据不同油气储层应力变化方式采用对应的实验方法来确定岩石样品渗透率与净应力关系,并将渗透率与净应力关系曲线沿净应力增加的方向进行外延至当前净应力点,并以实际储层原始净应力为起点来进行确定实际储层的应力敏感性,使得到的应力敏感性更加符合储层的实际情况,得到的储层应力敏感性更加精准,为储层保护,油田开发、开发规划等提供更有价值的数值。

Description

一种储层应力敏感性的确定方法
技术领域
本发明涉及一种储层应力敏感性的确定方法,属于油气田开发技术领域。
背景技术
世界上所有的低渗透储层几乎都存在应力敏感性。应力敏感性是油气藏开采过程中,由于储集层的孔隙压力降低而导致多孔介质固体颗粒原有的受力平衡状态被打破后建立新的压力平衡状态的过程,也是岩石变形和流体渗流耦合作用的过程。发生应力敏感后的储集层孔隙空间被压缩变形,渗流效果发生改变,宏观表现为渗透率的降低,进而影响油气的产出,导致产量降低。因此,搞清储层应力敏感性对于储层保护、开发规划、产能预测等具有重要意义。
目前对于储层的应力敏感性评价方法没有统一的标准,问题主要集中在两个方面:第一,测定应力敏感性曲线都一刀切的采用围压应力敏感测定方法。目前室内测定应力敏感性的方法主要包括两类:围压应力敏感和回压应力敏感,但前者操作较为方便,多被试验时采用,后者实验操作较为复杂,因为复杂室内测定都直接采用围压应力敏感曲线测定方法。实际的油气藏储层应力变化分为两类:一类是适宜采用注水开发的储层,能够及时注水补充地层能量;另一类是注水困难、注水压力高的储层,由于地层损失能量很难得到及时补充,而不得不采用衰竭式开采,如深层高压低渗储层等,气藏一般采用的也是衰竭式开采,因此也归为该类;两类储层应力变化机理是不同的。而目前不管是哪类油气藏储层应力变化均采用采用围压应力敏感测定方法,实际上,这种方法并不适于所有的油气藏储层应力敏感性评价。
第二,应力敏感性评价方法没有统一的标准,在使用室内测得应力敏感性曲线进行储层敏感性评价时存在较多误区。很多学者直接应用室内的应力敏感曲线进行储层敏感性评价,也有学者提出应该将应力敏感曲线按照曲线趋势进行分段评价,但两者都不完善,应力敏感性一定要与和地层应力条件相结合,否则就没了意义,那么如何应用室内测得的应力敏感曲线进行储层敏感性评价呢?为此,提出岩石应力敏感性和储层应力敏感性的概念,岩石应力敏感性就是脱离了实际储层条件的单纯的岩心的应力敏感性,即采用室内测得的完整的敏感性曲线进行评价的敏感性,属于纯科研范围,而不是储层的实际敏感性;储层应力敏感性,就要还原到实际储层的应力条件,但是由于实际储层的应力条件一般都较为苛刻,对于室内试验设备的压力要求都较高,很多室内压力设备很难满足要求,那么如何应用室内已测得的岩石应力敏感性曲线进行储层应力敏感性评价呢?
综上可知,使用目前的室内测定方法所测得应力敏感曲线并不能直接用来进行实际储层的应力敏感性评价,这也是目前室内试验设备的一大局限,如何科学高效地利用目前室内测定的结果解决实际问题呢,即在保证室内试验设备满足实验条件下的情况下,要使研究方法更科学,评价的结果更符合实际,将具有十分重要的意义。
发明内容
本发明的目的是提供一种储层应力敏感性的确定方法,以解决目前确定的储层应力敏感性的不实性和片面型。
本发明为解决上述技术问题而提供一种储层应力敏感性的确定方法,该应力敏感性确定方法的方案一:该方法包括以下步骤:
1)根据目标储层的实际储层应力变化状况选择对应的实验方法进行实验,得到各压力点对应的净应力和渗透率值,并将得到的净应力和渗透率值进行回归拟合,确定净应力与渗透率的关系;
2)根据目标储层的深度、岩石密度、原始地层压力和当前地层压力确定目标储层的原始净应力和当前净应力;
3)根据步骤1)确定的净应力和与渗透率的关系计算目标储层的原始渗透率和当前渗透率,将净应力与渗透率的关系曲线沿净应力增加的方向进行外延至当前净应力点,并以原始净应力为起点,获取一组净应力点,分别求取各净应力点对应的渗透率值,以得到各净应力下的渗透率数值;
4)以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率与起点的渗透率的比值,该比值为相应的净应力点对应的渗透率保持率,绘制渗透率保持率与净应力关系曲线,该曲线即为储层应力敏感曲线。
应力敏感性确定方法的方案二:在应力敏感性确定方法的方案一的基础上,所述步骤1)中实际的油气藏储层应力变化分为两类:一类是采用注水开发的储层,该类储层采用围压应力敏感评价实验;一类是采用衰竭式开采的储层,该类储层采用回压应力敏感性评价实验。
应力敏感性确定方法的方案三:在应力敏感性确定方法的方案二的基础上,所述回压应力敏感性评价实验是按照恒定围压、恒定注入压力、降低回压的方法进行,升压过程中,保持围压、注入压力和回压一同分级上升,且保持围压大于注入压力,注入压力等于回压,并在围压升至设定值时,保持围压和注入压力恒定,逐级降低回压,分别测定不同回压下流体渗流稳定的流量,直至测定完成所有设定的回压点。
应力敏感性确定方法的方案四:在应力敏感性确定方法的方案二的基础上,所述围压应力敏感评价实验是按照恒定注入压力和回压值的方法进行,并始终保持注入压力比回压高出恒定值,逐级升高围压开展实验,分别测定各个围压下的稳定流量,直至测定完成所有设定的围压点。
应力敏感性确定方法的方案五:在应力敏感性确定方法的方案一的基础上,所述步骤1)中在将净应力和渗透率值进行回归拟合时可分段进行。
应力敏感性确定方法的方案六:在应力敏感性确定方法的方案一的基础上,所述步骤2)中目标储层的原始净应力为目标储层上覆压力与目标储层的原始地层压力的差值,所述目标储层的当前净应力为目标储层上覆压力与目标储层的当前地层压力的差值,所述目标储层上覆压力是根据目标储层的深度和岩石密度计算得到。
应力敏感性确定方法的方案七:在应力敏感性确定方法的方案一的基础上,所述步骤3)在进行外延时,采用设定间隔的方式确定外延的净应力点。
应力敏感性确定方法的方案八:在应力敏感性确定方法的方案一的基础上,该方法还包括应力敏感对目标储层的损害程度的确定:以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率变化率,选取渗透率变化率的最大值作为目标储层的储层应力敏感损害率,所述的储层应力敏感损害率即为应力敏感对目标储层的损害程度。
有益效果:本发明根据不同油气储层应力变化方式采用对应的实验方法来确定岩石样品渗透率与净应力关系,并将渗透率与净应力关系曲线沿净应力增加的方向进行外延至当前净应力点,并以实际储层原始净应力为起点来进行确定实际储层的应力敏感性,使得到的储层应力敏感性更加符合储层的实际情况,得到的储层应力敏感性更加精准,为储层保护,油田开发、开发规划等提供更有价值的数值。
附图说明
图1是本发明实施例中测定的岩石渗透率与净应力关系曲线图;
图2是本发明实施例中测定的岩石应力敏感曲线图;
图3是本发明实施例中确定的岩石渗透率与净应力拟合关系曲线图;
图4是本发明实施例中净应力外延后的渗透率与净应力关系曲线图;
图5是本发明实施例中确定的储层应力敏感曲线图;
图6是本发明实施例中确定的渗透变化率与净应力关系曲线图。
具体实施方式
下面结合附图对本发明的具体实施方式做进一步的说明。
本发明储层应力敏感性确定方法的实施例
本实施例的储层应力敏感性确定方法根据不同油气储层应力变化方式采用对应的实验方法,将实验测定的应力敏感曲线进行内插和外延,以目标储层的原始净应力为起点来计算目标储层的应力敏感性,使确定的目标储层应力敏感性更加明确,符合储层实际情况。下面结合具体储层数据为例进行详细说明。
1.根据目标储层的实际储层应力变化状况选择对应的实验方法进行实验,得到各压力点对应的净应力和渗透率值。
实际的油气藏储层应力变化分为两类:一类是适宜采用注水开发的储层,能够及时注水补充地层能量;一类是注水困难、注水压力高的储层,由于地层损失能量很难得到及时补充,而不得不采用衰竭式开采,如深层高压低渗储层等,气藏一般采用的也是衰竭式开采,此也归为该类。这两类储层应力变化机理是不同的,因此本发明针对这两类储层应力变化分别选取围压应力敏感评价实验和回压应力敏感性评价实验。
回压应力敏感性评价实验是按照恒定围压、恒定注入压力、降低回压的方法进行,实验升压过程中,保持围压、注入压力、回压一同分级缓慢上升,保持围压比注入压力高3MPa,注入压力和回压相等,并考虑到实验设备压力要求,将围压升至设定值(本实施例为40MPa)为止,保持围压和注入压力恒定,逐级降低回压,分别测定不同回压下流体渗流稳定的流量,直至测定完成所有设定的回压点,结束实验。
围压应力敏感评价实验是按照恒定注入压力和回压值,并始终保持注入压力比回压高出恒定值,逐级升高围压开展实验,分别测定个各围压下的稳定的流量,直至测定完成所有设定的围压点,结束实验。
本实施例中选取的目标储层是衰竭式开采油藏,因此选取回压应力敏感性实验方法来进行实验。
根据得到的实验数据计算净应力和渗透率值,净应力(或称净围压)是围压与孔隙内压的差值,但是考虑岩石前后端受力不均,孔隙内压为注入压力与回压的平均值;按照达西定律计算上述所有压力点对应的渗透率值,定义渗透率保持率某回压下的渗透率与初始回压下渗透率的比值的百分数。本实施例通过所有实验点的计算,得到的该块岩石净应力与渗透率的关系曲线如图1所示,净应力与渗透率保持率关系曲线如图2所示。
2.将得到的净应力和渗透率值进行回归拟合,确定净应力与渗透率的关系。
将实验测定得到的净应力和渗透率数值进行回归拟合,当实验数据拟合比较困难时,可将净应力与渗透率的关系曲线进行分段拟合,例如,当得到的岩石净应力与渗透率的关系曲线明显无法用一种关系式来精确表示,根据岩石净应力与渗透率的关系曲线将其进行分段处理,以得到精确的拟合关系式。本实施例通过拟合方式得到净应力和渗透率关系如图3所示,其中R是实验数据点拟合时的相关系数,是表征变量之间线性相关程度的量,采用EXCEL趋势线拟合时,R是自动计算出的,R值越大,说明所拟合的趋势线与试验点符合程度越好。
3.根据目标储层的深度、岩石密度、原始地层压力和当前地层压力确定目标储层的原始净应力和当前净应力。
首先采集目标储层的深度和岩石密度参数,按照下式计算目标储层实际的上覆压力数值:
其中PC为上覆层压力,单位为MPa;ρr为上覆层平均密度,一般为2.3g/cm3;H为储层深度,单位为m;
然后结合原始地层压力和当前地层压力的数值计算原始地层和当前地层条件下的净应力的数值,储层的净应力为储层上覆压力与地层压力的差值。
本实施例选取的目标储层深度为3750m,目标储层的原始地层压力为40MPa,当前地层压力为29MPa,计算的上覆压力为74.35MPa。则对应的原始净应力为34.35MPa,当前储层的净应力为45.35MPa。
4.根据步骤2确定的净应力与渗透率的关系计算目标储层的原始渗透率和当前渗透率,通过采用插值法在原始净应力和当前净应力之间以2MPa为间隔进行内插,或以当前净应力的基础上按照上述插值的方法进行外延至未来规划的净应力点,获取一系列净应力点,分别求取各点的渗透率值,以得到各净应力下的渗透率数值。
本实施例在34.35~45.35MPa之间以2MPa为间隔分别插入一些净应力点,并同样以2MPa为间隔向后延长净应力点至60MPa,按照图3拟合的公式分别求取各组岩石净应力对应的渗透率值,并绘制成渗透率与净应力关系曲线,如图4所示。
5.以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率与起点的渗透率的比值,该比值为相应的净应力点对应的渗透率保持率。
以原始净应力对应的渗透率为起点,第n个净应力点对应的渗透率与起点渗透率的比值的百分数即为该点的渗透率保持率,依次计算所有净应力点对应的渗透率保持率,数据处理结束。
本实施例将图4中的渗透率数值以34.35MPa对应的渗透率数值为起点,后面净应力对应的渗透率与34.35MPa对应渗透率的比值的百分数,即为该净应力下渗透率保持率,依次计算后续各净应力点对应的渗透率保持率。将上述计算出的净应力与渗透率保持率数值绘制成以净应力为横坐标、以渗透率保持率为纵坐标的关系曲线即为储层的应力敏感曲线,如图5所示。
此外,本发明还能够根据得到的储层的应力敏感曲线确定应力敏感损害程度,确定过程为:以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率变化率,渗透率变化率为渗透率保持率减去1,选取渗透率变化率的最大值作为目标储层的应力敏感损害率,应力敏感损害率即为应力敏感对目标储层的损害程度。
本实施例将图4中的渗透率数值以34.35MPa对应的渗透率数值为起点,计算三块样品的渗透率变化率,计算结果如图6所示,在当前地层压力下对应三块样品的渗透率变化率分别为90.78%、88.74%和87.20%,从图6数值变化趋势来看,当前压力下对应的渗透率变化率数值即为从原始净应力条件降至目前净应力条件下的渗透率变化率最大值,分别为9.22%、11.26%和12.80%,该数值依照目前在用的储层敏感性流动实验评价方法,敏感损害程度为“弱”。

Claims (8)

1.一种储层应力敏感性的确定方法,其特征在于,该方法包括以下步骤:
1)根据目标储层的实际储层应力变化状况选择对应的实验方法进行实验,得到各压力点对应的净应力和渗透率值,并将得到的净应力和渗透率值进行回归拟合,确定净应力与渗透率的关系;
2)根据目标储层的深度、岩石密度、原始地层压力和当前地层压力确定目标储层的原始净应力和当前净应力;
3)根据步骤1)确定的净应力和与渗透率的关系计算目标储层的原始渗透率和当前渗透率,将净应力与渗透率的关系曲线沿净应力增加的方向进行外延至当前净应力点,并以原始净应力为起点,获取一组净应力点,分别求取各净应力点对应的渗透率值,以得到各净应力下的渗透率数值;
4)以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率与起点的渗透率的比值,该比值为相应的净应力点对应的渗透率保持率,绘制渗透率保持率与净应力关系曲线,该曲线即为储层应力敏感曲线。
2.根据权利要求1所述的储层应力敏感性的确定方法,其特征在于,所述步骤1)中实际的油气藏储层应力变化分为两类:一类是采用注水开发的储层,该类储层采用围压应力敏感评价实验;一类是采用衰竭式开采的储层,该类储层采用回压应力敏感性评价实验。
3.根据权利要求2所述的储层应力敏感性的确定方法,其特征在于,所述回压应力敏感性评价实验是按照恒定围压、恒定注入压力、降低回压的方法进行,升压过程中,保持围压、注入压力和回压一同分级上升,且保持围压大于注入压力,注入压力等于回压,并在围压升至设定值时,保持围压和注入压力恒定,逐级降低回压,分别测定不同回压下流体渗流稳定的流量,直至测定完成所有设定的回压点。
4.根据权利要求2所述的储层应力敏感性的确定方法,其特征在于,所述围压应力敏感评价实验是按照恒定注入压力和回压值的方法进行,并始终保持注入压力比回压高出恒定值,逐级升高围压开展实验,分别测定各个围压下的稳定流量,直至测定完成所有设定的围压点。
5.根据权利要求1所述的储层应力敏感性的确定方法,其特征在于,所述步骤1)中在将净应力和渗透率值进行回归拟合时可分段进行。
6.根据权利要求1所述的储层应力敏感性的确定方法,其特征在于,所述步骤2)中目标储层的原始净应力为目标储层上覆压力与目标储层的原始地层压力的差值,所述目标储层的当前净应力为目标储层上覆压力与目标储层的当前地层压力的差值,所述目标储层上覆压力是根据目标储层的深度和岩石密度计算得到。
7.根据权利要求1所述的储层应力敏感性的确定方法,其特征在于,所述步骤3)在进行外延时,采用设定间隔的方式确定外延的净应力点。
8.根据权利要求1所述的储层应力敏感性的确定方法,其特征在于,该方法还包括应力敏感对目标储层的损害程度的确定:以目标储层原始净应力对应的渗透率为起点,计算各净应力点对应的渗透率变化率,选取渗透率变化率的最大值作为目标储层的储层应力敏感损害率,所述的储层应力敏感损害率即为应力敏感对目标储层的损害程度。
CN201710891799.1A 2017-09-27 2017-09-27 一种储层应力敏感性的确定方法 Pending CN109557280A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710891799.1A CN109557280A (zh) 2017-09-27 2017-09-27 一种储层应力敏感性的确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710891799.1A CN109557280A (zh) 2017-09-27 2017-09-27 一种储层应力敏感性的确定方法

Publications (1)

Publication Number Publication Date
CN109557280A true CN109557280A (zh) 2019-04-02

Family

ID=65863990

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710891799.1A Pending CN109557280A (zh) 2017-09-27 2017-09-27 一种储层应力敏感性的确定方法

Country Status (1)

Country Link
CN (1) CN109557280A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110470581A (zh) * 2019-08-02 2019-11-19 中国石油大学(北京) 确定储层应力敏感程度的方法、装置及存储介质
CN111999227A (zh) * 2020-08-19 2020-11-27 中国石油大学(北京) 一种应力敏感性测试方法和装置
CN114166726A (zh) * 2022-02-14 2022-03-11 西南石油大学 一种岩心渗透率张量敏感性损伤测试设备与评价方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106153518A (zh) * 2016-06-22 2016-11-23 西南石油大学 致密砂岩气储层压裂液损害实验评价方法
CN107132170A (zh) * 2017-04-10 2017-09-05 中国石油天然气股份有限公司 一种储层应力敏感性的确定方法和装置

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106153518A (zh) * 2016-06-22 2016-11-23 西南石油大学 致密砂岩气储层压裂液损害实验评价方法
CN107132170A (zh) * 2017-04-10 2017-09-05 中国石油天然气股份有限公司 一种储层应力敏感性的确定方法和装置

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
朱华银 等: "大庆深层火山岩储层应力敏感性研究", 《天然气地质学》 *
杨满平 等: "考虑储层初始有效应力的岩石应力敏感性分析", 《天然气地球科学》 *
林琳: "室内试验研究储层应力敏感性", 《石油化工应用》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110470581A (zh) * 2019-08-02 2019-11-19 中国石油大学(北京) 确定储层应力敏感程度的方法、装置及存储介质
CN110470581B (zh) * 2019-08-02 2021-05-11 中国石油大学(北京) 确定储层应力敏感程度的方法、装置及存储介质
CN111999227A (zh) * 2020-08-19 2020-11-27 中国石油大学(北京) 一种应力敏感性测试方法和装置
CN114166726A (zh) * 2022-02-14 2022-03-11 西南石油大学 一种岩心渗透率张量敏感性损伤测试设备与评价方法
US11604132B1 (en) 2022-02-14 2023-03-14 Southwest Petroleum University Testing device and evaluation method for sensitivity damage of core permeability tensor

Similar Documents

Publication Publication Date Title
CN102778421B (zh) 地质历史时期砂岩储层渗透率演化恢复方法
CN109490507B (zh) 页岩气井体积压裂后压裂液返排率的实验测试装置及方法
CN108376188B (zh) 黄土湿陷性的评价计算方法
CN107462936B (zh) 利用压力监测资料反演低渗透储层非达西渗流规律的方法
CN110147561A (zh) 一种含天然裂缝致密油气储层体积压裂缝网预测方法
CN106932324B (zh) 一种确定高含水砂岩油藏储层渗透率变化规律的方法
CN106295095B (zh) 基于常规测井资料预测低渗透砂岩储层产能的方法
CN106599449A (zh) 一种溶洞体积计算的试井解释方法
CN204125898U (zh) 一种模拟地震作用下坝基渗流的试验装置
CN105447762B (zh) 一种流体替换的低渗透油藏水淹信息的计算方法
CN109557280A (zh) 一种储层应力敏感性的确定方法
CN104634804A (zh) 一种利用核磁共振t2谱确定储层相对渗透率的方法
CN110259444A (zh) 基于流场诊断的水驱油藏渗流场可视化表征与评价方法
CN109577929A (zh) 一种超低渗致密油藏水平井建立有效驱替的定量评价方法
CN110208487A (zh) 一种基于ct扫描的页岩水化损伤测试方法
CN104747149B (zh) 二元复合驱剂窜通道的检测方法
CN109031436A (zh) 一种基于储层分类下的低阻低渗疏松砂岩储层含水饱和度评价方法
Wei et al. A quick analysis approach for early-time well test data of vertically fractured wells: Applications in Changqing oilfield, China
CN110374562A (zh) 一种油藏堵调方法
CN106320255A (zh) 一种漫滩水流滩槽过流量的计算方法
CN106545335B (zh) 致密油藏衰竭式开发合理生产压差的确定方法
CN107992690A (zh) 一种诱导应力干扰下多裂缝扩展形态均衡程度的评价方法
CN104131813B (zh) 一种非常规气藏地层参数的获取方法及系统
CN108343431B (zh) 一种基岩裂隙地下水系统流网探究方法
CN105606490A (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20190402

RJ01 Rejection of invention patent application after publication