CN108918683A - 一种超临界二氧化碳压裂相变的声发射检测方法 - Google Patents

一种超临界二氧化碳压裂相变的声发射检测方法 Download PDF

Info

Publication number
CN108918683A
CN108918683A CN201810722666.6A CN201810722666A CN108918683A CN 108918683 A CN108918683 A CN 108918683A CN 201810722666 A CN201810722666 A CN 201810722666A CN 108918683 A CN108918683 A CN 108918683A
Authority
CN
China
Prior art keywords
carbon dioxide
supercritical carbon
test specimen
detection method
acoustic emission
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
CN201810722666.6A
Other languages
English (en)
Other versions
CN108918683B (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
Priority to CN201810722666.6A priority Critical patent/CN108918683B/zh
Publication of CN108918683A publication Critical patent/CN108918683A/zh
Application granted granted Critical
Publication of CN108918683B publication Critical patent/CN108918683B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
    • G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques

Landscapes

  • Physics & Mathematics (AREA)
  • Acoustics & Sound (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

本发明提供了一种超临界二氧化碳压裂相变的声发射检测方法。该检测方法包括：制备隔断试件；对隔断试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号的变化，得到超临界二氧化碳发生相变的声发射信号的标准曲线；制备样品试件；对样品试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳的声发射信号，根据标准曲线，获得超临界二氧化碳压裂过程的相变特征。本发明的超临界二氧化碳压裂相变的检测方法是一种超临界二氧化碳压裂相变发生情况的判定方法，可以获得超临界二氧化碳压裂过程中超临界二氧化碳的相变特征。

Description

一种超临界二氧化碳压裂相变的声发射检测方法
技术领域
本发明涉及一种超临界二氧化碳相变的检测方法，尤其涉及一种根据声发射信号判断超临界二氧化碳相变的检测方法，属于致密岩石油气资源开发技术领域。
背景技术
致密岩石油气藏在中国分布广泛，资源量巨大。超临界二氧化碳(SC-CO2)的密度接近液体，粘度接近气体，具有很强的溶剂化能力，扩散系数比液体大，具有良好的传质性，表面张力为零。正因为SC-CO2的这些优点，近些年，SC-CO2压裂受到广泛关注，在致密岩石油气藏开发领域有着广阔的应用前景。
然而，SC-CO2压裂过程复杂，是制约分析超临界二氧化碳压裂过程的关键难题。SC-CO2相态变化对于压裂效果具有积极地作用，有利于加快裂缝发育与提高压裂效果。
因此，探究SC-CO2压裂过程中SC-CO2的相态变化是本领域亟待解决的重要问题。
发明内容
为了解决上述技术问题，本发明的目的在于提供一种可以对超临界二氧化碳(SC-CO2)压裂过程进行声发射检测，并判断压裂过程中SC-CO2相变发生情况的方法。
为了实现上述目的，本发明提供了一种超临界二氧化碳压裂相变的声发射检测方法，该检测方法包括以下步骤：
制备隔断试件；
对隔断试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号的变化，得到超临界二氧化碳发生相变的声发射信号的标准曲线；
制备样品试件；
对样品试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳的声发射信号，根据标准曲线，获得超临界二氧化碳压裂过程的相变特征。
在本发明的一具体实施方式中，采用的隔断试件的内部具有隔断。
需要说明的是，在对隔断试件和试样试件进行超临界二氧化碳压裂实验试，压裂的实验参数和条件是相同的。
在本发明的一具体实施方式中，采用的隔断试件按照以下步骤制备：
将水泥与砂混合，制备试件，放置15天以上；
在试件上钻井筒孔以及探头孔；
在试件的表面切割声发射导线槽；
在井筒孔的底部预留长度10mm的空隙，空隙的上部通过水泥密封，水泥块的长度为30mm，粘接井筒，得到隔断试件。
在本发明的一具体实施方式中，将水泥与砂混合，制备得到具有一定形状的试件。对试件的形状没有特殊要求，一般采用的是正方体，比如试件的形状为长、宽、高均为300mm的正方体。
在本发明的一具体实施方式中，井筒孔的直径为20mm，长度为150mm。
在本发明的一具体实施方式中，探头孔的数量需要满足实验要求。具体地，可以钻12个探头孔；更进一步地，探头孔的直径为20mm，长度为25mm。
在本发明的一具体实施方式中，采用环氧树脂AB胶粘接井筒。
在本发明的一具体实施方式中，将水泥与砂按照模拟的样品的岩石比例进行混合。
在本发明的一具体实施方式中，采用的隔断试件中人为制造的空腔(隔断)，为超临界二氧化碳相变提供巨大的空间，保证超临界二氧化碳相变的发生。注入的超临界二氧化碳发生相变会对岩石内部造成动态载荷，动态载荷会在岩石表面产生冲击波。
在本发明的一具体实施方式中，采用的样品试件按照以下步骤制备：
将水泥与砂混合，制备试件，放置15天以上；
钻井筒孔以及探头孔；
在试件的表面切割声发射导线槽；
粘接井筒，得到样品试件。
在本发明的一具体实施方式中，将水泥与砂混合，制备得到具有一定形状的试件。对试件的形状没有特殊要求，一般采用的是正方体，比如试件的形状为长、宽、高均为300mm的正方体。
在本发明的一具体实施方式中，井筒孔的直径为20mm，长度为150mm。
在本发明的一具体实施方式中，探头孔的数量需要满足实验要求。具体地，可以钻12个探头孔；更进一步地，探头孔的直径为20mm，长度为25mm。
在本发明的一具体实施方式中，采用环氧树脂AB胶粘接井筒。
在本发明的一具体实施方式中，将水泥与砂按照模拟的样品的岩石比例进行混合。
在本发明的一具体实施方式中，进行超临界二氧化碳压裂实验时可以采用大型物理模拟压裂试验机。
在本发明的一具体实施方式中，根据地应力条件，确定三向地应力。
根据本发明的具体实施方式，采用本领域常规的大型物理模拟压裂试验机即可实现目的。在本发明中采用大型物理模拟压裂试验机进行真三轴试验，真三轴试验是使岩石试件处于三个主应力不相等(即σ1>σ2>σ3)的应力组合状态下的三轴压缩试验。
在本发明的一具体实施方式中，获得超临界二氧化碳发生相变的声发射信号的标准曲线时，是采用波形分析与三维损伤定位相结合，基于信号幅值、频率、能量与定位点位分布，提取相变信号特征。
在本发明的一具体实施方式中，声发射信号是通过声发射检测仪检测的。根据本发明的具体实施方式，采用本领域常规的声发射检测仪即可。
本发明的超临界二氧化碳压裂相变的声发射检测方法，具体包括以下步骤：
步骤一：制备隔断试件，包括以下步骤：
将水泥与砂按照模拟样品的岩石比例混合，制备规格为300mm×300mm×300mm的正方体试件，放置15天；
钻井筒孔，其中，井筒孔的直径为20mm，长度为150mm；
钻探头孔12个，其中，探头孔的直径为20mm，长度为25mm；
在试件的表面切割声发射导线槽；
在井筒孔底部留长度10mm的空隙，空隙的上部用水泥封住，水泥块的长度为30mm；
采用环氧树脂AB胶粘接井筒，得到隔断试件。
步骤二：对隔断试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号的变化，得到超临界二氧化碳发生相变的声发射信号的标准曲线，具体包括以下步骤：
将加工好的试件放入实验台上，连接声发射探头、声发射导路、注液管路；
启动超临界二氧化碳温度压力控制系统，保证注入的二氧化碳处于超临界状态；
开启声发射检测仪，设置声发射采样参数，波形图与参数表，测试调整探头检测精度；
对试件施加预选设定好的三向地应力，开启温度控制设备，保证试验机温度达到设定温度；
启动注入泵，向试件内部注入超临界二氧化碳，设定超临界二氧化碳的流量和流速，同时检测声发射信号，采集记录，确定相变过程中声发射信号特征；
压裂实验结束后，可以观察压裂裂缝破坏特征和方向，拍照记录。
步骤三：制备样品试件，包括以下步骤：
将水泥与砂按照模拟样品的岩石比例进行混合，制备规格为300mm×300mm×300mm的正方体试件，放置15天；
钻井筒孔，其中，井筒孔的直径为20mm，长度为150mm；
钻探头孔12个，其中，探头孔的直径为20mm，长度为25mm；
在试件的表面切割声发射导线槽；
采用环氧树脂AB胶粘接井筒，得到样品试件。
步骤四：对样品试件进行超临界二氧化碳压裂实验(实验步骤同步骤二)，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号，根据标准曲线，获得超临界二氧化碳压裂过程的相变特征。
本发明的超临界二氧化碳压裂相变的声发射检测方法考虑到裂缝发育过程中裂缝体积增加导致缝内压力降低，超临界二氧化碳向气态转变，此时压裂方式由静态向动态转变，最终导致岩石破坏。
本发明的超临界二氧化碳压裂相变的声发射检测方法，将无法观测的相变特征通过声发射检测手段转化为可以观测分析的声信号，进而获得超临界二氧化碳压裂过程中超临界二氧化碳的相变特征。
本发明的超临界二氧化碳压裂相变的声发射检测方法可以获得在真三轴条件下超临界二氧化碳压裂相态变化声发射特征与三维损伤定位特征。
本发明的超临界二氧化碳压裂相变的声发射检测方法可以模拟在真三轴条件下超临界二氧化碳的相变过程。
本发明的超临界二氧化碳压裂相变的声发射检测方法，能够通过声发射检测手段，获得表征相变的声发射信号特征，识别超临界二氧化碳压裂过程中的相变，为超临界二氧化碳压裂过程的相变分析提供参考，探究超临界二氧化碳的压裂致裂机理。
附图说明
图1为实施例1的超临界二氧化碳压裂相变的检测方法的工艺流程图。
图2为实施例1中的隔断试件的结构示意图。
图3为实施例1的压裂装置的结构示意图。
图4为实施例1中的注液压力和时间的曲线。
图5为实施例1中的注液压力和时间的曲线。
图6为实施例1中的注液压力和时间的曲线。
图7为实施例1中的注液压力和时间的曲线。
主要附图符号说明
1井筒 2环氧树脂AB胶 3水泥封块 4预制隔断 5水泥基质 6试件 7三向加压泵 8真空泵 9饱和液站 10千斤顶 11机架 12超临界二氧化碳温度压力控制设备 13二氧化碳罐 14电脑 15注入泵 16声发射检测仪
具体实施方式
为了对本发明的技术特征、目的和有益效果有更加清楚的理解，现对本发明的技术方案进行以下详细说明，但不能理解为对本发明的可实施范围的限定。
实施例1
本实施例提供了本发明提供了一种超临界二氧化碳压裂相变的声发射检测方法，其工艺如图1所示，该声发射检测方法具体包括以下步骤：
步骤一：制备隔断试件，其结构如图2所示，包括以下步骤：
将水泥与砂按1:1比例混合，制备规格为300mm×300mm×300mm正方体试件，放置15天，得到水泥基质5；
钻井筒孔，井筒孔用于使井筒1通过，其中，井筒孔的直径为20mm，长度为150mm；
钻探头孔12个，其中，碳井孔的直径为20mm，长为25mm；
在试件的表面切割声发射导线槽；
在井筒孔底部预留长度为10mm的空隙，空隙的上部用水泥封住，水泥封块3长度的30mm，得到预制隔断4；
采用环氧树脂AB胶2粘接井筒1，得到隔断试件，其结构如图2所示；
步骤二：采用如图3所示的装置(整个装置由机架11支撑)对隔断试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号的变化，得到隔断试件相变与声信号的标准曲线，具体包括以下步骤：
将加工好的试件6(隔断试件)放入实验台上，连接好声发射探头和声发射导路，连接好注液管路；
由二氧化碳罐13泵入二氧化碳，启动超临界二氧化碳温度压力控制设备12，保证注入的二氧化碳处于超临界状态；
开启声发射检测仪16，设置声发射采样参数，波形图与参数表(其中，采样频率为5000kHz，采样长度为3797，预采样为150微秒，采样间隔150微秒。波形图与参数表同步采集)，测试调整探头检测精度；
打开饱和液站9和真空泵8，为三向加压泵7提供动力；
通过三向加压泵7调节千斤顶10对试件6施加预选设定好的三向地应力(σ1＝30MPa,σ2＝20MPa,σ3＝15MPa)，开启温度控制系统，保证试验机温度达到设定温度(60℃)；
启动注入泵15，向试件6内部注入超临界二氧化碳，设定超临界二氧化碳的流量(18mL/min)，同时检测声发射信号，通过电脑14采集记录，确定相变过程中声发射信号特征；
压裂实验结束后，可以观察压裂裂缝破坏特征和方向，拍照记录。
步骤三：制备样品试件，包括以下步骤：
将水泥与砂按1:1比例混合，制备300mm×300mm×300mm的正方体试件，放置15天；
钻井筒孔，井筒孔的直径为20mm，长度为150mm；
钻探头孔12个，其中，探头孔的直径为20mm，长度为25mm；
在试件的表面切割声发射导线槽；
采用环氧树脂AB胶粘接井筒，得到样品试件。
步骤四：对样品试件进行超临界二氧化碳压裂实验(实验步骤同步骤二)，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号，根据标准曲线，获得超临界二氧化碳压裂过程的相变特征。得到的曲线如图4、图5、图6和图7所示。
由图5可以看出，注液压力过峰值点时，能量值呈现瞬间剧增。

Claims (10)

1.一种超临界二氧化碳压裂相变的声发射检测方法，其特征在于，该检测方法包括以下步骤：
制备隔断试件；
对所述隔断试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳发生相变时的声发射信号的变化，得到超临界二氧化碳发生相变的声发射信号的标准曲线；
制备样品试件；
对所述样品试件进行超临界二氧化碳压裂实验，并检测压裂过程中超临界二氧化碳的声发射信号，根据所述标准曲线，获得超临界二氧化碳压裂过程的相变特征。
2.根据权利要求1所述的检测方法，其特征在于，所述隔断试件的内部具有隔断。
3.根据权利要求1或2所述的检测方法，其特征在于，所述隔断试件按照以下步骤制备：
将水泥与砂混合，制备试件，放置15天以上；
在所述试件上钻井筒孔以及探头孔；
在所述试件的表面切割声发射导线槽；
在所述井筒孔的底部预留长度10mm的空隙，空隙的上部通过水泥密封，水泥块的长度为30mm，粘接井筒，得到所述隔断试件。
4.根据权利要求3所述的检测方法，其特征在于，所述试件的长、宽、高均为300mm；优选地，所述井筒孔的直径为20mm，长度为150mm。
5.根据权利要求3所述的检测方法，其特征在于，所述探头孔的数量为12个；优选地，所述探头孔的直径为20mm，长度为25mm。
6.根据权利要求1所述的检测方法，其特征在于，所述样品试件按照以下步骤制备：
将水泥与砂混合，制备试件，放置15天以上；
钻井筒孔以及探头孔；
在所述试件的表面割声发射导线槽；
粘接井筒，得到所述样品试件。
7.根据权利要求6所述的检测方法，其特征在于，所述试件的长、宽、高均为300mm；优选地，所述井筒孔的直径为20mm，长度为150mm。
8.根据权利要求6所述的检测方法，其特征在于，所述探头孔的数量为12个；优选地，所述探头孔的直径为20mm，长度为25mm。
9.根据权利要求1所述的检测方法，其特征在于，进行超临界二氧化碳压裂实验时根据地应力条件，确定三向地应力。
10.根据权利要求1所述的检测方法，其特征在于，获得所述超临界二氧化碳发生相变的声发射信号的标准曲线，是采用波形分析与三维损伤定位相结合，根据声信号的变化表征相变的特征。
CN201810722666.6A 2018-07-04 2018-07-04 一种超临界二氧化碳压裂相变的声发射检测方法 Active CN108918683B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810722666.6A CN108918683B (zh) 2018-07-04 2018-07-04 一种超临界二氧化碳压裂相变的声发射检测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810722666.6A CN108918683B (zh) 2018-07-04 2018-07-04 一种超临界二氧化碳压裂相变的声发射检测方法

Publications (2)

Publication Number Publication Date
CN108918683A true CN108918683A (zh) 2018-11-30
CN108918683B CN108918683B (zh) 2020-05-01

Family

ID=64424471

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810722666.6A Active CN108918683B (zh) 2018-07-04 2018-07-04 一种超临界二氧化碳压裂相变的声发射检测方法

Country Status (1)

Country Link
CN (1) CN108918683B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114018719A (zh) * 2021-11-04 2022-02-08 中国矿业大学 一种超临界二氧化碳压裂温压精准监测试验装置及方法
GB2613587A (en) * 2021-12-07 2023-06-14 Rolls Royce Plc Sensor assembly

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4095475A (en) * 1976-04-22 1978-06-20 Massachusetts Institute Of Technology Apparatus and method whereby wave energy is correlated with geometry of a manufactured part or the like or to positional relationships in a system
SU1363053A1 (ru) * 1986-07-01 1987-12-30 Предприятие П/Я В-2652 Устройство дл исследовани фазовых превращений в образцах
CN102606069A (zh) * 2012-04-01 2012-07-25 中国石油大学(华东) 超临界二氧化碳钻井井筒内相态的控制装置
CN104458918A (zh) * 2014-12-30 2015-03-25 重庆大学 超临界二氧化碳压裂页岩损伤定位监测装置及方法
CN105672974A (zh) * 2016-02-25 2016-06-15 重庆大学 三轴应力下超临界二氧化碳致裂页岩实验试件制作方法
CN106404549A (zh) * 2016-10-31 2017-02-15 中国石油大学(北京) 一种超临界二氧化碳压裂模拟实验装置

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4095475A (en) * 1976-04-22 1978-06-20 Massachusetts Institute Of Technology Apparatus and method whereby wave energy is correlated with geometry of a manufactured part or the like or to positional relationships in a system
SU1363053A1 (ru) * 1986-07-01 1987-12-30 Предприятие П/Я В-2652 Устройство дл исследовани фазовых превращений в образцах
CN102606069A (zh) * 2012-04-01 2012-07-25 中国石油大学(华东) 超临界二氧化碳钻井井筒内相态的控制装置
CN104458918A (zh) * 2014-12-30 2015-03-25 重庆大学 超临界二氧化碳压裂页岩损伤定位监测装置及方法
CN105672974A (zh) * 2016-02-25 2016-06-15 重庆大学 三轴应力下超临界二氧化碳致裂页岩实验试件制作方法
CN106404549A (zh) * 2016-10-31 2017-02-15 中国石油大学(北京) 一种超临界二氧化碳压裂模拟实验装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JIE KE ET AL.: "Detecting phase transitions in supercritical mixtures: an enabling tool for greener chemical reactions", 《PROC.R.SOC.A》 *

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114018719A (zh) * 2021-11-04 2022-02-08 中国矿业大学 一种超临界二氧化碳压裂温压精准监测试验装置及方法
CN114018719B (zh) * 2021-11-04 2024-01-26 中国矿业大学 一种超临界二氧化碳压裂温压精准监测试验装置及方法
GB2613587A (en) * 2021-12-07 2023-06-14 Rolls Royce Plc Sensor assembly
GB2613587B (en) * 2021-12-07 2023-12-06 Rolls Royce Plc Sensor assembly
US11927378B2 (en) 2021-12-07 2024-03-12 Rolls-Royce Plc Sensor assembly

Also Published As

Publication number Publication date
CN108918683B (zh) 2020-05-01

Similar Documents

Publication Publication Date Title
CN105510142B (zh) 一种煤岩多相不同流体三轴压裂试验装置与试验方法
CN107063963B (zh) 一种致密储层微裂缝扩展及渗流特征的测试装置和方法
CN104655495B (zh) 一种煤岩高温高压真三轴压裂渗流试验装置与试验方法
CN107515185A (zh) 一种煤体渗透率测试装置及其使用方法
CN111272576A (zh) 一种新的真三轴压裂渗流试验装置与方法
CN104535727B (zh) 一种水力加砂压裂系统
CN107905778A (zh) 超临界co2流体压裂增强型地热系统实验装置及方法
CN109298162A (zh) 不同相态二氧化碳致裂页岩装置及实验方法
CN107884288B (zh) 高温下含断续节理的岩石压缩剪切试样制作及试验方法
CN107727530B (zh) 基于温压回溯原理的一体化全程气密含气量测量仪的应用
CN104865124A (zh) 基于岩石应力-应变曲线和超声波纵波速度的页岩脆性指数测定方法
CN106092757B (zh) 一种静载与水压耦合作用破岩的试验装置及其使用方法
CN102735547A (zh) 真三轴状态下煤岩水压致裂试验方法
CN105352811A (zh) 一种小型水力压裂三轴评价测试加压装置及其方法
CN105403505A (zh) 固井界面胶结强度测试装置及方法
CN208239220U (zh) 原位孔内剪切测试仪
Zheng et al. Residual strength characteristics of CJPL marble under true triaxial compression
CN108918683A (zh) 一种超临界二氧化碳压裂相变的声发射检测方法
CN106499383A (zh) 一种评价固井屏障完整性的方法
Zhixi et al. Determination of rock fracture toughness and its relationship with acoustic velocity
CN107091623A (zh) 隧道围岩松动圈厚度计算方法
Ma et al. Comparisons of fracturing mechanism of tight sandstone using liquid CO2 and water
WO2021239149A1 (zh) 一种利用共振冲击波改造煤层气储层的系统及方法
CN105545290B (zh) 一种旋喷固化封孔模拟试验平台及试验方法
Yang et al. True triaxial hydraulic fracturing test and numerical simulation of limestone

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