CN109682853B - 一种基于fbg的冻土含冰量分布式原位测量方法及装置 - Google Patents
一种基于fbg的冻土含冰量分布式原位测量方法及装置 Download PDFInfo
- Publication number
- CN109682853B CN109682853B CN201910020807.4A CN201910020807A CN109682853B CN 109682853 B CN109682853 B CN 109682853B CN 201910020807 A CN201910020807 A CN 201910020807A CN 109682853 B CN109682853 B CN 109682853B
- Authority
- CN
- China
- Prior art keywords
- fbg
- frozen soil
- ice content
- sensor
- tubular sensor
- 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
Links
- 239000002689 soil Substances 0.000 title claims abstract description 170
- 238000000034 method Methods 0.000 title claims abstract description 50
- 238000012625 in-situ measurement Methods 0.000 title claims abstract description 34
- 238000010438 heat treatment Methods 0.000 claims abstract description 73
- 238000012545 processing Methods 0.000 claims abstract description 21
- 238000012544 monitoring process Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 13
- 238000005553 drilling Methods 0.000 claims abstract description 7
- 239000013307 optical fiber Substances 0.000 claims description 35
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 23
- 238000009792 diffusion process Methods 0.000 claims description 14
- 230000009471 action Effects 0.000 claims description 8
- 238000007405 data analysis Methods 0.000 claims description 8
- 238000005481 NMR spectroscopy Methods 0.000 claims description 7
- 206010070834 Sensitisation Diseases 0.000 claims description 6
- 230000004907 flux Effects 0.000 claims description 6
- 230000001681 protective effect Effects 0.000 claims description 6
- 230000008313 sensitization Effects 0.000 claims description 6
- 238000004458 analytical method Methods 0.000 claims description 5
- 230000036760 body temperature Effects 0.000 claims description 4
- 238000004806 packaging method and process Methods 0.000 claims description 4
- 229920000049 Carbon (fiber) Polymers 0.000 claims description 3
- 239000004917 carbon fiber Substances 0.000 claims description 3
- 238000001816 cooling Methods 0.000 claims description 3
- 239000004744 fabric Substances 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 2
- WABPQHHGFIMREM-UHFFFAOYSA-N lead(0) Chemical compound [Pb] WABPQHHGFIMREM-UHFFFAOYSA-N 0.000 claims 2
- 238000005259 measurement Methods 0.000 abstract description 12
- 230000008859 change Effects 0.000 abstract description 4
- 239000000835 fiber Substances 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 239000007789 gas Substances 0.000 description 7
- 230000008569 process Effects 0.000 description 6
- 239000002245 particle Substances 0.000 description 5
- 238000007710 freezing Methods 0.000 description 4
- 230000008014 freezing Effects 0.000 description 4
- 239000011148 porous material Substances 0.000 description 4
- 229920006335 epoxy glue Polymers 0.000 description 3
- 238000001566 impedance spectroscopy Methods 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 230000000087 stabilizing effect Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 230000007797 corrosion Effects 0.000 description 2
- 238000005260 corrosion Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005538 encapsulation Methods 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 238000012886 linear function Methods 0.000 description 2
- 239000007788 liquid Substances 0.000 description 2
- 238000001225 nuclear magnetic resonance method Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000004044 response Effects 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000012546 transfer Methods 0.000 description 2
- 241001391944 Commicarpus scandens Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000003139 buffering effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000004134 energy conservation Methods 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 230000017525 heat dissipation Effects 0.000 description 1
- 238000002513 implantation Methods 0.000 description 1
- 238000012613 in situ experiment Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 239000011810 insulating material Substances 0.000 description 1
- 238000009413 insulation Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052760 oxygen Inorganic materials 0.000 description 1
- 239000001301 oxygen Substances 0.000 description 1
- 230000000737 periodic effect Effects 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N25/00—Investigating or analyzing materials by the use of thermal means
- G01N25/20—Investigating or analyzing materials by the use of thermal means by investigating the development of heat, i.e. calorimetry, e.g. by measuring specific heat, by measuring thermal conductivity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
- G01K11/00—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00
- G01K11/32—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres
- G01K11/3206—Measuring temperature based upon physical or chemical changes not covered by groups G01K3/00, G01K5/00, G01K7/00 or G01K9/00 using changes in transmittance, scattering or luminescence in optical fibres at discrete locations in the fibre, e.g. using Bragg scattering
Landscapes
- Physics & Mathematics (AREA)
- General 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)
- Immunology (AREA)
- Pathology (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本发明公开了一种基于FBG的冻土含冰量分布式原位测量方法及装置,包括:加热电源、内置电阻丝和FBG的管状传感器、FBG解调仪和分析处理监测数据的计算机。将制作完成的管状传感器通过直埋或钻孔埋设入待测冻土中;连接加热电源对管状传感器内置的电阻丝进行短期加热,热量通过导热性能良好的管状传感器扩散到周围冻土中;管状传感器内置的FBG感测温度变化,并通过FBG解调仪采集并记录FBG的波长读数;将波长数据处理转化为温度变化,得到升温过程中的温度特征值;最后通过率定试验建立的冻土含冰量i与温度特征值ΔTt间的线性关系,得到冻土含冰量。本发明可以实现分布式、连续性测量冻土含冰量。
Description
技术领域
本发明公开了一种光纤布拉格光栅(FBG)的分布式原位测量冻土含冰量的方法及装置,涉及冻土含冰量测量技术领域。
背景技术
冻土是一种温度低于0℃且含有冰的土,是由土颗粒、冰、未冻水、气体组成的多相复杂体系。冻土对温度十分敏感且物理性质不稳定,其特性除了与土壤质地、容重和含水量有关外,还受含冰量的控制。不同于未冻土,冻土中的冰,使得冻土的性质既特殊又复杂。因此,测量冻土含冰量对于理论研究和工程实践都具有重要的意义。
目前,冻土含冰量测量方法有膨胀法、介电谱法、热脉冲法和核磁共振法(NMR)。
膨胀法基本原理是:将待测土壤样本放置在一个圆柱容器中润湿至饱和,在容器的一端插入一根管状的张力计并延至该容器的另一端。当土壤冻结时,由于土壤是饱和的,形成的冰的体积就会膨胀,为了保证土壤结构不被破坏,膨胀的体积会迫使部分液态水从土壤孔隙中流出并透过张力计进入预先标定好的毛细管中。通过测量毛细管中的液态水,再依据水冻结成冰的膨胀系数,冰的体积便可以计算出来。
介电谱法的基本原理是:通过测定冻土的介电常数来间接反映含冰量,通常通过确定两个频率下的介电谱并结合介电混合模型来推算饱和土壤的含冰量。
热脉冲法基本原理是:当冻土温度不变时,未冻水含量不再变化。借助热脉冲测定出冻土的热导率,便可计算出冻土含冰量。仅仅适用于温度较低的环境下,当温度接近冻结点时,测定结果不具可靠性。
核磁共振法(NMR)的测量原理是:一些原子的原子核就像一些小型的磁棒,在外加强磁场的作用下,这些原子核会定向排列。如果再加上无线电频波的作用,这些原子将吸收足够的能量从而在外加磁场内以另一个稳定的方向重新排列。将预先准备好的冻土土样放置于一个脉冲分析仪中,并对土样发射度的无线电脉冲。在无线电脉冲的作用下,围绕在土样周围的接收线圈将产生能够反映吸收能量原子数的电压,其大小与土样中固态冰的氧原子量成正比。分析仪通过检测上述电压便能测定出冻土样本中的含冰量。
上述四种方法均能够测量冻土含冰量,也都有各自的局限性:膨胀法仅仅局限在实验室环境下,容易破坏冻土原有结构且无法测定非饱和土壤的冻土含冰量;热脉冲法在温度较低的环境下可以测定出含冰量,但当温度接近冻结点时,测定的数据不够准确;介电谱法的使用受到土壤类型的制约;NMR法具有较高的测量精度,但分析仪体积庞大、价格昂贵,且安装与使用十分复杂,其应用范围仅仅局限于实验室环境,无法广泛使用。
发明内容
本发明的目的是提供一种基于FBG的分布式原位测量冻土含冰量的方法及装置,基于温度响应原理,运用FBG管状传感器,对内加热的管体进行温度感测,通过管体的升温曲线来确定温度特征值,从而测量冻土含冰量。解决了含冰量测量不准确、冻土原状结构易被破坏、难以对冻土含冰量进行原位测量等缺点。
为解决上述问题,本发明采用以下技术方案:一种基于FBG的冻土含冰量分布式原位测量方法,包括如下步骤:
步骤一、将制作、封装完备的管状传感器植入待测冻土的相应位置,所述管状传感器包括管体、光纤、加热电阻丝和FBG传感器;
步骤二、所述管状传感器管体在电流作用下开始升温;待扩散热流密度恒定后停止升温,管体开始降温;
步骤三、FBG解调仪采集并记录加热时间区间[t1,t2]内FBG的波长读数,所述的t1为扩散热流密度稳定后30s,t2为管体开始降温的时间;
步骤四、运用数据分析处理系统将波长数据转化为管体温度信息;计算管体的温度特征值,根据管体的温度特征值与冻土含冰量的线性关系:i=k1ΔTt+b1计算出冻土含冰量,其中i为冻土含冰量,ΔTt为管状传感器所测的温度特征值,k1、b1为常数,通过若干组冻土试样的率定试验确定;所述管状传感器温度特征值为通电加热后管状传感器周围扩散热流密度恒定后,选取一个特征时间区间[t1,t2]内的等时间间隔所测得的温度的算术平均值与初始温度的差值。
步骤一中是采用钻孔或直埋的方式将管状传感器竖向铺设于待测冻土土壤中。
步骤二所述管状传感器加热功率恒定,保证传感器稳态加热;步骤四中所述等时间间隔为10秒。
步骤三中所述时间区间[t1,t2]定义为特征时间区间,取值根据不同含水量、不同土壤条件而变化。
所述冻土含冰量i为冻土中冰的质量与全部水的质量的比值;冻土试样含冰量由核磁共振法确定。
一种用于所述的FBG的冻土含冰量分布式原位测量方法中的装置,包括加热电源、管状传感器、FBG解调仪和分析处理监测数据装置,所述管状传感器包括管体、光纤、加热电阻丝、FBG传感器,管体内贯穿有光纤和加热电阻丝,光纤上设置若干FBG传感器,加热电阻丝通过通电导线和加热电源连接,光纤通过光纤引线和FBG解调仪连接,用于采集并记录加热趋于稳定后的波长读数;所述分析处理监测数据装置连接FBG解调仪,并运用数据分析处理系统将波长数据转化为管体温度信息并计算管体的温度特征值。
所述管状传感器的管体材料采用高导热绝缘塑料,导热系数为1.13-1.20W/(m.k)。
所述管状传感器采用增敏性封装结构,由两个直径为R的半圆柱体组成,在其中一个半圆柱体截面中间位置开设直径为r的圆形小槽,光纤铺设在小槽中,两端呈自然松弛状态,不受拉力;往小槽内注入不固化导热膏;将两个半圆柱体用环氧胶固定封装,在管体外每隔间距D设置一个夹具固定管体。
所述管状传感器内相邻FBG传感器间距为d;管状传感器单独或串联使用;相邻管状传感器通过两端的螺口和螺母相连。
所述管状传感器外部光纤引线外设碳纤维布保护套,沿钻孔及管体传感器外壁连接FBG传感器。
光纤布拉格光栅(FBG)是在纤芯制造等周期光栅,当光纤沿轴向受拉或温度变化后,光纤沿轴向会发生形变,光纤折射率也随之变化导致输出信号光谱会发生变化,因此可以实现数值化测量。光纤光栅传感器是在光纤光栅的基础上,以布拉格条件为基础发展起来的一种本征波长调制型传感器,具有体积小(裸FBG传感器外径125μm)、抗电磁干扰能力强、性能不受湿度影响、稳定性好、耐腐蚀、灵敏度高的优点,近年来广泛应用于监测桥梁、大坝和岩土结构。
有益效果:
1.本发明基于冻土导热系数测量含冰量,在较低温度下测量可以降低水分迁移对测量结果的影响。
2.本发明可以直接测量原状冻土含冰量,对冻土扰动较小,避免了冻土在取样和运输及保存过程中的结构和成分变化。
3.本发明可以实现分布式、连续性测量冻土含冰量。
4.本发明可以实现对冻土含冰量变化的实时监测。
5.本发明经济安全,操作便捷,抗干扰能力强,精度可靠,稳定性高。
6.本法明可以应用于不同规模的实验研究,可以采用串联的方式提高空间分辨率。
附图说明
图1是本发明所述分布式原位测量传感器内部纵向结构示意图。
其中,1、螺口;2、不固化导热膏;3、FBG;4、光纤;5、光纤保护套;6、通电导线;7、电阻丝;8、高导热绝缘管体;
图2是本发明所述分布式原位测量传感器内部横截面示意图。
其中,2、不固化导热膏;4、光纤;7、电阻丝;8、高导热绝缘管体;9、环氧胶;r是增敏结构内部光纤穿出孔直径;R是增敏结构外部的圆形柱状管体的直径。
图3是本发明所述分布式原位测量传感器整体结构示意图。
其中,10、螺母;11、FBG管状传感器;12、夹具;13、光纤引出口;D是夹具间距。
图4是本发明所述冻土含冰量分布式原位测量系统示意图。
其中,14、加热电源;15、分析处理监测数据的计算机;16、FBG解调仪;17、光纤引线;18、通电导线;19、冻土;20、管状传感器;21、光纤光栅;22、电阻丝;23、钻孔。
图5是实施例1砂土中管状传感器温度特征值与冻土含冰量线性关系标定结果图。
图6是实施例1砂土含冰量测量结果图。
图7是实施例2中含冰量随冻土深度的变化曲线图。
具体实施方式
附图提供对本发明的进一步理解,并且构成说明书的一部分,与本发明的实施例一起用于解释本发明,并不构成对本发明的限制。下面结合附图和实施例对本发明做过更进一步的解释。
如图1及图2所示,一种基于FBG的冻土含冰量测量的内加热管状传感器,由内到外依次包括FBG传感器、光纤、加热电阻丝、不固化导热膏、带螺口的管体、光纤保护套、通电导线。
如图3所示,一种基于FBG的冻土含冰量测量的内加热管状传感器,传感器管体外部每隔间距D设置一个夹具固定管体,相邻管体之间通过螺母连接,螺母上预留光纤引线穿出孔口。
如图4所示,一种基于FBG的冻土含冰量分布式原位测量系统,包括加热电源1、分析处理监测数据的计算机2、FBG解调仪3、管状传感器7。所述加热电源1保持恒定功率,使电阻丝在稳定电流下工作;所述管状传感器7为具有内加热功能的FBG的管状传感器;所述FBG解调仪3连接带有FBG传感器的光纤,用于采集并记录加热趋于稳定后的波长读数;所述分析处理监测数据的计算机2连接FBG解调仪,并运用数据分析处理系统将波长数据转化为管体温度信息并计算管体的温度特征值。
FBG解调仪检测FBG反射的布拉格波长来反映温度信息。本试验中采用苏州南智传感科技有限公司生产的A-01FBG解调器对FBG进行解调,采集波长读数。解调仪采样记录间隔为10秒。
一种基于FBG的冻土含冰量分布式原位测量方法,包括如下步骤:
步骤一、将制作、封装完备的管状传感器植入待测冻土的相应位置,所述传感器为具有内加热功能的FBG管状传感器,由管体、光纤、内置电阻丝、FBG传感器组成;
步骤二、连接电源,对步骤一所述管状传感器通电加热,管体在电流作用下开始升温;待传感器周围扩散热流密度恒定后停止通电加热,管体开始降温;
步骤三、FBG解调仪采集并记录加热时间[t1,t2]区间内FBG的波长读数,t1为扩散热流密度稳定后30s,t2为管体开始降温的时间;
步骤四、运用数据分析处理系统将波长数据转化为管体温度信息;计算管体的温度特征值,根据管体的温度特征值与冻土含冰量的线性关系:i=k1ΔTt+b1计算出冻土含冰量,其中i为冻土含冰量,ΔTt为管状传感器所测的温度特征值,k1、b1为常数,通过若干组冻土试样的率定试验确定;所述管状传感器温度特征值为通电加热后管状传感器形成的温度场趋于稳定后,选取一个特征时间区间[t1,t2]内的等时间间隔所测得的温度的算术平均值与初始温度的差值。
进一步地,所述的基于FBG的冻土含冰量分布式原位测量方法,步骤一中将制作、封装完备的管状传感器植入待测冻土的相应位置,所述植入方法:采用钻孔或直埋的方式将传感器竖向铺设于待测冻土土壤中。
进一步地,所述的基于FBG的冻土含冰量分布式原位测量方法,步骤二所述管状传感器加热功率稳定,保证传感器稳态加热;步骤四中所述等时间间隔为10秒。
进一步地,所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,步骤三中所述时间区间[t1,t2]定义为特征时间区间,取值根据不同含水量、不同土壤条件而有所变化。
进一步地,所述冻土含冰量i为冻土中冰的质量与全部水的质量的比值。
进一步地,所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,所述率定试验包括如下步骤:
步骤一、若干组已知含冰量冻土土壤;
步骤二、将制作、封装完备的管状传感器植入待测冻土的相应位置,所述管状传感器为具有内加热功能的FBG管状传感器,由管体、光纤、加热电阻丝、FBG传感器组成;
步骤三、连接电源,对步骤二所述管状传感器通电加热,管体在电流作用下开始升温;待传感器周围扩散热流密度恒定后停止通电加热,管体开始降温;
步骤四、FBG解调仪采集并记录加热时间[t1,t2]区间内FBG的波长读数,t1为扩散热流密度稳定后30s,t2为管体开始降温的时间,将波长数据转化为管体温度信息;
步骤五、运用数据分析处理系统计算各管体的温度特征值,并拟合各管体的温度特征值与冻土含冰量线性关系:i=k1ΔTt+b1,其中i为冻土含冰量,ΔTt为管状传感器所测的温度特征值,k1、b1为常数;所述管状传感器温度特征值为通电加热后管状传感器周围扩散热流密度恒定后,选取一个特征时间区间[t1,t2]内的等时间间隔所测得的温度的算术平均值与初始温度的差值。
进一步地,所述的基于FBG的冻土含冰量分布式原位测量方法,所述若干组已知含冰量的冻土土壤为4~6组;冻土试样含冰量由核磁共振法确定。
一种基于FBG的冻土含冰量分布式原位测量装置,包括加热电源、管状传感器、FBG解调仪和分析处理监测数据的计算机。所述加热电源保持恒定功率,使电阻丝在稳定电流下工作;所述管状传感器为具有内加热功能的FBG管状传感器,由管体、光纤、加热电阻丝、FBG传感器组成;所述FBG解调仪连接带有FBG传感器的光纤,用于采集并记录加热趋于稳定后的波长读数;所述分析处理监测数据的计算机连接FBG解调仪,并运用数据分析处理系统将波长数据转化为管体温度信息并计算管体的温度特征值。
进一步地,所述冻土含冰量分布式原位测量装置,所述管状传感器的管体材料采购自昆山市奋发绝缘材料有限公司,导热系数为1.13-1.20W/(m.k)。所述高导热绝缘塑料具有导热性好、韧性好、比重轻、绝缘强、耐腐蚀、耐老化等显著特点。
进一步地,所述冻土含冰量分布式原位测量装置,所述FBG管状传感器管体内开设两条孔道,其中一条安装加热电阻丝,另一孔道安装带有FBG传感器的光纤。
进一步地,所述冻土含冰量分布式原位测量装置,基于FBG传感器对温度和应变同时响应的特性以及光纤易断的性质,所述FBG传感器通过增敏性封装提高温度灵敏系数:增敏结构由两个直径为R的半圆形柱体组成,在其中一个半圆形柱体截面中间位置开设直径为r的圆形小槽,光纤铺设在小槽中,两端呈自然松弛状态,不受拉力;往小槽内注入不固化导热膏以加快热传导的速度,同时对应变起到缓冲甚至消除作用;将两个半圆形柱体用环氧胶固定封装,考虑到分布式原位测量所需管体长度较大,在管体外每隔间距D设置一个夹具固定管体。
进一步地,所述冻土含冰量分布式原位测量装置,所述管状传感器内串联10个FBG传感器,相邻FBG传感器间距为d;多个管状传感器的串联实现冻土含冰量的准分布式原位测量;相邻管体通过两端的螺口和螺母相连。
进一步地,所述冻土含冰量分布式原位测量装置,所述管状传感器外部光纤引线通过碳纤维布保护套的保护性封装,沿钻孔及管体传感器外壁连接FBG传感器。
本发明原理:基于FBG的冻土含冰量分布式原位测量方法的基本原理是利用FBG传感器所测得的升温过程中温度特征值与冻土含冰量之间的线性关系,来测量冻土含冰量。进一步可以解释为:冻土的热传导性能由土壤颗粒、气体、未冻水、冰四者共同决定,其中,空气导热系数为0.024W/(m*K),水导热系数为0.60W/(m*K),冰导热系数为2.25W/(m*K)。土壤颗粒各项特征在测量过程中保持不变,气体导热系数远远小于未冻水和冰的导热系数,故可以忽略不计,土壤总含水率保持恒定,未冻水含量由含冰量决定,因此,冻土的热传导能力由含冰量决定。由于冰的导热系数是水的3~4倍,所以含冰量越高,冻土的热传导能力越强。将具有内加热功能的管体FBG传感器植入待测冻土中,通电后管体加热温度升高,与冻土之间形成温差,在含冰量越高的冻土中传热能力越强,而管体由于电源功率恒定产生的总能量固定,所以向冻土土体中扩散的能量越多,用于管体加热的能量就越少,管体的温度特征值就越低。因此,通过测量管体加热一定时间后得到的温度特征值,可以求得冻土含冰量。
实施例1
应用本发明所述方法和装置进行某冻土含冰量测量的室内实验。
管状传感器温度特征值与该冻土含冰量的线性关系的率定试验:
步骤一、通过钻孔采集某冻土地区不同深度冻土土壤试样5组,使用核磁共振法测量各深度冻土的含冰量,测量结果为:3%、4.6%、8.7%、9.3%、11.5%;
步骤二、用大体积环刀(环刀高度≥4cm,容积≥120cm3)取样,同时将制作、封装完备的管状传感器植入待测冻土的中心位置,所述管状传感器为具有内加热功能的FBG管状传感器,由管体、光纤、加热电阻丝、FBG传感器组成;
步骤三、连接电源,对步骤二所述管状传感器通电,管体在电流作用下开始升温;待管体周围扩散热流密度恒定后停止通电加热,管体开始降温;
步骤四、FBG解调仪每隔10秒采集并记录波长读数,将波长数据转化为管体温度信息;
步骤五、运用数据分析处理系统绘制管状传感器的升温曲线,选取的温度特征区间为加热后[14min,25min],计算该特征区间内各管体的温度特征值,并结合前面得到的冻土含冰量拟合得到各管体的温度特征值与冻土含冰量线性关系:i=-0.709ΔTt+0.1406,经线性拟合所得的标定系数R2=0.998,如图5所示。
按照上述率定试验相同的方法和步骤,通电加热、采集并解调信息、分析温度信息、计算温度特征值,将所求温度特征值代入率定试验所确定的一次函数i=-0.709ΔTt+0.1406中,计算冻土含冰量,如图6。从图6中可以看出来,使用本发明方法测得的冻土含冰量为8.2%。
本实施例1试验过程中所述冻土试样、FBG管状传感器均放置于与取样地同等温度的冷冻室中,加热电源、FBG解调仪、分析处理监测数据的计算机均放置于冷冻室外。
实施例2
应用本发明所述方法和装置进行分布式测量冻土含冰量的原位试验,对不同深度冻土土壤含冰量进行监测。试验地点选取青海某砂性土地区的冻土土壤。
步骤一、按照实施例1所述方法和步骤对该地区砂性土温度特征值与含冰量的线性关系进行率定试验,拟合得到温度特征值与含冰量的线性关系;
步骤二、通过钻孔将13个FBG管状传感器串联竖直埋设入待测冻土中,所述FBG管状传感器为10个FBG传感器串联组成,相邻传感器间距为20cm,可同时测量不同深度土壤的含冰量;
步骤三、将管状传感器连至稳定功率的电源和FBG传感器;连接电源,对步骤二所述FBG管状传感器通电加热,管体在电流作用下开始升温;待传感器周围扩散热流密度恒定后停止通电加热,管体开始降温;
步骤四、FBG解调仪每隔10秒自动采集并记录波长读数,将波长数据转化为管体温度信息;
步骤五、运用数据分析处理系统计算各FBG传感器感测的温度特征值,将所求温度特征值代入率定试验所确定的一次函数中,计算各深度的冻土含冰量,可以得到冻土含冰量随深度的变化曲线,如图7所示。
将本实施例中管状传感器温度特征值与冻土含冰量之间线性关系的理论推导过程阐述如下:
假设待测冻土具有均匀性,管体位于初始温度一致的无限土层中。将此类冻土中热量传递简化成一维问题。在管体表面取单位面积,根据欧姆定律,该单位面积单位时间内产生的能量为:
Q1=I2R (1)
式(1)中Q1为电阻丝单位面积上产生的能量,I为电流,R为电阻丝的电阻。I、R均为已知常数,Q1所以也为常数。
根据能量守恒,单位时间内用于管体加热的能量表示为:
Q2=Cm(T-T0)=CmΔTt (2)
式(2)中Q2为用于加热管体的能量;Cm为管体的比热容;T0为加热前管体的温度;T为加热后管体的实测温度,取值为管体周围扩散热流密度恒定后管体的平均温度;ΔTt定义为温度特征值。
管体单位面积在单位时间内散失的热量Q:
Q=Q1-Q2=I2R-CmΔTt (3)
将热源看做是无线长线热源,管体周边冻土的温度场满足:
式(5)中,r为以管体为圆心的半径;t为换热时间;T(r,t)为t时刻距管体中心半径为r处的岩土体温度;Q为管体换热量;L为管体长度;λ为岩土体导热系数;a为岩土体热扩散率;β为计算过程参数;T∞为距管体中心无穷远处岩土体温度(即地层初始温度);CS为岩土体单位体积热容,。
假设管体稳态导热,r取值管体半径rw,则加热后所测温度T:
基于无线长线热源模型,使用线性推导法,由(6)式可知T与时间t的对数呈线性关系。可以简化为:
T=klnt+m (7)
其中,T为管体内加热稳定后平均温度;P为管体热阻;γ为欧拉常数,取0.577216;k和m为FBG询问器测得的温度与时间的关系,基于最小二乘法拟合得到直线的斜率和截距。
结合温度响应的散热量Q可以得到导热率:
冻土由土骨架、气体、未冻水、冰组成,故冻土的导热率由四部分组成
λ=λss+λgg+λll+λii (11)
其中,s、g为冻土中土壤颗粒、气体质量与冻土总质量的比值,l、i为冻土中未冻水、冰的质量与全部水的质量的比值。气体含量g极小可以忽略,土骨架、未冻水、冰的导热系数λs、λl、λs以及土壤颗粒的含量s可以通过相关实验和资料得到,未冻水含量l=1-i,故冻土导热率是关于含冰量的线性函数,可以表示为:
λ=ai+b (12)
式(12)中,a=λi,b=λss+λgg+λll,a、b都是常数。
联立式(3)(10)、(12)可以得到含冰量和温度特征值ΔTt的线性关系:
可以进一步简化得:
i=k1ΔTt+b1 (14)
其中k1、b1都是常数。
由式(14)可以看出管状传感器的温度特征值与冻土含冰量之间成线性关系。因此可以通过监测管状传感器的温度来计算冻土含冰量。
以上所述仅是本发明的优选实施方式,应当指出,对于本技术领域的普通技术人员,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应包含在本发明的保护范围之内。
Claims (8)
1.一种基于FBG的冻土含冰量分布式原位测量方法,其特征在于,包括如下步骤
步骤一、将制作、封装完备的管状传感器植入待测冻土的相应位置,所述管状传感器包括管体、光纤、加热电阻丝和FBG传感器;
步骤二、所述管状传感器管体在电流作用下开始升温;待扩散热流密度恒定后停止升温,管体开始降温;
步骤三、FBG解调仪采集并记录加热特征时间区间[t 1,t 2]内FBG的波长读数,所述的t 1为扩散热流密度稳定后30s,t 2为管体开始降温的时间;
步骤四、运用数据分析处理系统将波长数据转化为管体温度信息;计算管体的温度特征值,根据管体的温度特征值与冻土含冰量的线性关系:计算出冻土含冰量,其中/>为冻土含冰量,/>为管状传感器所测的温度特征值,/>、/>为常数,通过若干组冻土试样的率定试验确定;所述管状传感器温度特征值为通电加热后管状传感器周围扩散热流密度恒定后,选取一个特征时间区间[t1,t2]内的等时间间隔所测得的温度的算术平均值与初始温度的差值。
2.根据权利要求1所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,步骤一中是采用钻孔或直埋的方式将管状传感器竖向铺设于待测冻土土壤中。
3.根据权利要求1所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,步骤二所述管状传感器加热功率恒定,保证传感器稳态加热;步骤四中所述等时间间隔为10秒。
4.根据权利要求1所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,步骤三中所述时间区间[t1,t2]定义为特征时间区间,取值根据不同含水量、不同土壤条件而变化。
5.根据权利要求1所述的基于FBG的冻土含冰量分布式原位测量方法,其特征在于,所述冻土含冰量为冻土中冰的质量与全部水的质量的比值;冻土试样含冰量由核磁共振法确定。
6.一种用于权利要求1~5任一所述的FBG的冻土含冰量分布式原位测量方法中的装置,其特征在于,包括加热电源、管状传感器、FBG解调仪和分析处理监测数据装置,所述管状传感器包括管体、光纤、加热电阻丝、FBG传感器,管体结构采用增敏性封装,由两个直径为R的半圆柱体组成,在其中一个半圆柱体截面中间位置开设直径为r的圆形小槽,光纤铺设在小槽中,两端呈自然松弛状态,小槽内注入不固化导热膏;管体材料采用高导热绝缘塑料,管体内贯穿光纤和加热电阻丝,光纤上设置若干FBG传感器,加热电阻丝通过通电导线和加热电源连接,光纤通过光纤引线和FBG解调仪连接,用于采集并记录加热趋于稳定后的波长读数;所述分析处理监测数据装置连接FBG解调仪,并运用数据分析处理系统将波长数据转化为管体温度信息并计算管体的温度特征值。
7.根据权利要求6所述用于FBG的冻土含冰量分布式原位测量方法中的装置,其特征在于,所述管状传感器内相邻FBG传感器间距为d;管状传感器单独或串联使用;相邻管状传感器通过两端的螺口和螺母相连。
8.根据权利要求6所述用于FBG的冻土含冰量分布式原位测量方法中的装置,其特征在于,所述管状传感器外部光纤引线外设碳纤维布保护套,沿钻孔及管体传感器外壁连接FBG传感器。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910020807.4A CN109682853B (zh) | 2019-01-09 | 2019-01-09 | 一种基于fbg的冻土含冰量分布式原位测量方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910020807.4A CN109682853B (zh) | 2019-01-09 | 2019-01-09 | 一种基于fbg的冻土含冰量分布式原位测量方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109682853A CN109682853A (zh) | 2019-04-26 |
CN109682853B true CN109682853B (zh) | 2024-02-13 |
Family
ID=66192852
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910020807.4A Active CN109682853B (zh) | 2019-01-09 | 2019-01-09 | 一种基于fbg的冻土含冰量分布式原位测量方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109682853B (zh) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110793940B (zh) * | 2019-11-13 | 2022-04-05 | 南京大学 | 一种基于光纤光栅的土体基质吸力准分布式原位测量方法及装置 |
CN111879794A (zh) * | 2020-07-28 | 2020-11-03 | 中国科学院西北生态环境资源研究院 | 一种区域表层土壤含冰量的监测系统 |
CN112268926B (zh) * | 2020-11-03 | 2021-08-24 | 北京交通大学 | 一种冻土含冰量测量方法及装置 |
CN118348062A (zh) * | 2024-04-17 | 2024-07-16 | 河海大学 | 基于fbg热脉冲双探针的土体含水率测量方法与装置 |
CN118392795A (zh) * | 2024-04-23 | 2024-07-26 | 河海大学 | 一种土体-传感器接触热阻测定与土体含水率校正方法 |
Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU767275A1 (ru) * | 1978-04-11 | 1980-09-30 | За витель | Способ измерени льдистости мерзлых грунтов |
EP0405229A2 (de) * | 1989-06-24 | 1991-01-02 | Forschungszentrum Karlsruhe GmbH | Verfahren und eine Vorrichtung zur Messung des volumetrischen Wasseranteils mineralischer und/oder organischer Gemische |
US6003620A (en) * | 1996-07-26 | 1999-12-21 | Advanced Coring Technology, Inc. | Downhole in-situ measurement of physical and or chemical properties including fluid saturations of cores while coring |
CN102135513A (zh) * | 2011-01-04 | 2011-07-27 | 刘波 | 冻土未冻水含量的测试方法 |
WO2012002016A1 (ja) * | 2010-06-28 | 2012-01-05 | 株式会社フジクラ | 超電導線材の常電導転移の検出方法 |
CN202706066U (zh) * | 2012-04-25 | 2013-01-30 | 青海省电力公司检修公司 | 一种高海拔冻土层输电线路铁塔基础在线监测系统 |
CN103175742A (zh) * | 2013-03-06 | 2013-06-26 | 北京交通大学 | 一种冻土动载荷直剪仪 |
CN103364320A (zh) * | 2013-07-19 | 2013-10-23 | 河海大学 | 多孔介质结构体渗流的分布式光纤测试方法 |
CN103383345A (zh) * | 2012-05-04 | 2013-11-06 | 中国石油天然气股份有限公司 | 一种冻土区油气管道周围含水量监测系统 |
CN103454309A (zh) * | 2013-09-04 | 2013-12-18 | 南京大学 | 一种土壤含水率分布式测量方法及系统 |
WO2015165138A1 (zh) * | 2014-04-30 | 2015-11-05 | 河海大学 | 一种透明冻土及其制备方法和应用 |
CN105092824A (zh) * | 2015-09-09 | 2015-11-25 | 中国科学院寒区旱区环境与工程研究所 | 一种自记式冻土含冰量测量装置 |
CN105241920A (zh) * | 2015-11-06 | 2016-01-13 | 天津城建大学 | 利用比热计算确定土冻结过程中未冻水含量的方法 |
CN105651812A (zh) * | 2015-12-30 | 2016-06-08 | 湖北工业大学 | 一种基于dts检测灌注桩完整性的检测系统设计方法 |
RU2592915C1 (ru) * | 2015-06-20 | 2016-07-27 | Федеральное государственное бюджетное Учреждение науки - Институт мерзлотоведения им.П.И.Мельникова Сибирского отделения Российской академии наук | Способ определения количества незамерзшей воды в мерзлых грунтах |
CN106525280A (zh) * | 2016-11-25 | 2017-03-22 | 中交第公路勘察设计研究院有限公司 | 一种高海拔多年冻土区分布式高精度温度监测系统及方法 |
CN106546557A (zh) * | 2016-10-28 | 2017-03-29 | 中南大学 | 一种土体未冻水含量测量装置及测量方法 |
CN106979791A (zh) * | 2017-04-20 | 2017-07-25 | 苏州南智传感科技有限公司 | 内加热fbg传感器及其封装方法 |
JP2017128893A (ja) * | 2016-01-19 | 2017-07-27 | 東京電力ホールディングス株式会社 | 凍土方式遮水壁の造成方法 |
CN107132172A (zh) * | 2017-04-20 | 2017-09-05 | 苏州南智传感科技有限公司 | 基于ihat‑fbg的岩土体渗流速率与含水率监测系统及方法 |
CN107272073A (zh) * | 2017-04-26 | 2017-10-20 | 北京交通大学 | 一种应用探地雷达计算冻土相对含水量的方法 |
RU2654832C1 (ru) * | 2017-06-21 | 2018-05-22 | Автономная некоммерческая образовательная организация высшего образования "Сколковский институт науки и технологий" | Способ определения содержания незамерзшей воды в мерзлых грунтах |
CN108279211A (zh) * | 2018-03-25 | 2018-07-13 | 南京大学 | 一种测量土体含水率的准分布式系统、传感装置和方法 |
WO2018227281A1 (en) * | 2017-06-12 | 2018-12-20 | Advanced Opto-Mechanical Systems And Technologies Inc. | Multi-parameter distributed fiber optic sensor system and methods of sensor manufacturing |
CN208297333U (zh) * | 2018-03-25 | 2018-12-28 | 南京大学 | 一种测量土体含水率的准分布式系统和传感装置 |
CN109141269A (zh) * | 2018-10-09 | 2019-01-04 | 中国地震局地壳应力研究所 | 分布式光纤光栅孔壁应变计 |
-
2019
- 2019-01-09 CN CN201910020807.4A patent/CN109682853B/zh active Active
Patent Citations (26)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU767275A1 (ru) * | 1978-04-11 | 1980-09-30 | За витель | Способ измерени льдистости мерзлых грунтов |
EP0405229A2 (de) * | 1989-06-24 | 1991-01-02 | Forschungszentrum Karlsruhe GmbH | Verfahren und eine Vorrichtung zur Messung des volumetrischen Wasseranteils mineralischer und/oder organischer Gemische |
US6003620A (en) * | 1996-07-26 | 1999-12-21 | Advanced Coring Technology, Inc. | Downhole in-situ measurement of physical and or chemical properties including fluid saturations of cores while coring |
WO2012002016A1 (ja) * | 2010-06-28 | 2012-01-05 | 株式会社フジクラ | 超電導線材の常電導転移の検出方法 |
CN102135513A (zh) * | 2011-01-04 | 2011-07-27 | 刘波 | 冻土未冻水含量的测试方法 |
CN202706066U (zh) * | 2012-04-25 | 2013-01-30 | 青海省电力公司检修公司 | 一种高海拔冻土层输电线路铁塔基础在线监测系统 |
CN103383345A (zh) * | 2012-05-04 | 2013-11-06 | 中国石油天然气股份有限公司 | 一种冻土区油气管道周围含水量监测系统 |
CN103175742A (zh) * | 2013-03-06 | 2013-06-26 | 北京交通大学 | 一种冻土动载荷直剪仪 |
CN103364320A (zh) * | 2013-07-19 | 2013-10-23 | 河海大学 | 多孔介质结构体渗流的分布式光纤测试方法 |
CN103454309A (zh) * | 2013-09-04 | 2013-12-18 | 南京大学 | 一种土壤含水率分布式测量方法及系统 |
WO2015165138A1 (zh) * | 2014-04-30 | 2015-11-05 | 河海大学 | 一种透明冻土及其制备方法和应用 |
RU2592915C1 (ru) * | 2015-06-20 | 2016-07-27 | Федеральное государственное бюджетное Учреждение науки - Институт мерзлотоведения им.П.И.Мельникова Сибирского отделения Российской академии наук | Способ определения количества незамерзшей воды в мерзлых грунтах |
CN105092824A (zh) * | 2015-09-09 | 2015-11-25 | 中国科学院寒区旱区环境与工程研究所 | 一种自记式冻土含冰量测量装置 |
CN105241920A (zh) * | 2015-11-06 | 2016-01-13 | 天津城建大学 | 利用比热计算确定土冻结过程中未冻水含量的方法 |
CN105651812A (zh) * | 2015-12-30 | 2016-06-08 | 湖北工业大学 | 一种基于dts检测灌注桩完整性的检测系统设计方法 |
JP2017128893A (ja) * | 2016-01-19 | 2017-07-27 | 東京電力ホールディングス株式会社 | 凍土方式遮水壁の造成方法 |
CN106546557A (zh) * | 2016-10-28 | 2017-03-29 | 中南大学 | 一种土体未冻水含量测量装置及测量方法 |
CN106525280A (zh) * | 2016-11-25 | 2017-03-22 | 中交第公路勘察设计研究院有限公司 | 一种高海拔多年冻土区分布式高精度温度监测系统及方法 |
CN106979791A (zh) * | 2017-04-20 | 2017-07-25 | 苏州南智传感科技有限公司 | 内加热fbg传感器及其封装方法 |
CN107132172A (zh) * | 2017-04-20 | 2017-09-05 | 苏州南智传感科技有限公司 | 基于ihat‑fbg的岩土体渗流速率与含水率监测系统及方法 |
CN107272073A (zh) * | 2017-04-26 | 2017-10-20 | 北京交通大学 | 一种应用探地雷达计算冻土相对含水量的方法 |
WO2018227281A1 (en) * | 2017-06-12 | 2018-12-20 | Advanced Opto-Mechanical Systems And Technologies Inc. | Multi-parameter distributed fiber optic sensor system and methods of sensor manufacturing |
RU2654832C1 (ru) * | 2017-06-21 | 2018-05-22 | Автономная некоммерческая образовательная организация высшего образования "Сколковский институт науки и технологий" | Способ определения содержания незамерзшей воды в мерзлых грунтах |
CN108279211A (zh) * | 2018-03-25 | 2018-07-13 | 南京大学 | 一种测量土体含水率的准分布式系统、传感装置和方法 |
CN208297333U (zh) * | 2018-03-25 | 2018-12-28 | 南京大学 | 一种测量土体含水率的准分布式系统和传感装置 |
CN109141269A (zh) * | 2018-10-09 | 2019-01-04 | 中国地震局地壳应力研究所 | 分布式光纤光栅孔壁应变计 |
Non-Patent Citations (2)
Title |
---|
缩短核磁共振测定冻土未冻水含量实验时间的方法;李东阳;刘波;刘念;马永君;王莉;;冰川冻土(第06期);第1502-1507页 * |
青藏高原北麓河地区原状多年冻土导热系数的试验研究;温智, 盛煜, 马巍, 邓友生, 吴基春;冰川冻土(第02期);第182-187页 * |
Also Published As
Publication number | Publication date |
---|---|
CN109682853A (zh) | 2019-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109682853B (zh) | 一种基于fbg的冻土含冰量分布式原位测量方法及装置 | |
He et al. | Development and application of the heat pulse method for soil physical measurements | |
Sayde et al. | Feasibility of soil moisture monitoring with heated fiber optics | |
Ciocca et al. | Heated optical fiber for distributed soil‐moisture measurements: A lysimeter experiment | |
Tarara et al. | Measuring soil water content in the laboratory and field with dual‐probe heat‐capacity sensors | |
Cao et al. | A soil moisture estimation method using actively heated fiber Bragg grating sensors | |
Cao et al. | An improved distributed sensing method for monitoring soil moisture profile using heated carbon fibers | |
Jarrell et al. | Soil water and temperature status | |
CN110793940B (zh) | 一种基于光纤光栅的土体基质吸力准分布式原位测量方法及装置 | |
Manohar et al. | Measurement of apparent thermal conductivity by the thermal probe method | |
Cheng et al. | In situ measuring soil ice content with a combined use of dielectric tube sensor and neutron moisture meter in a common access tube | |
CN110658123A (zh) | 一种基于光纤主动变温的非饱和土渗透系数的原位测试方法 | |
Wang et al. | A modified soil water content measurement technique using actively heated fiber optic sensor | |
AU2020102910A4 (en) | A measuring instrument for soil thermal conductivity profile characteristics | |
Zhu et al. | Monitoring infiltration of capillary barrier with actively heated fibre Bragg gratings | |
US5121993A (en) | Triaxial thermopile array geo-heat-flow sensor | |
CN209542497U (zh) | 一种基于fbg的冻土含冰量分布式原位测量装置 | |
Zhao et al. | Sandy soil moisture content measurement method based on heated fiber Bragg grating | |
CN114858720A (zh) | 原位冻土冰水组分辨识与动态示踪的光电测试系统和方法 | |
JP2001281071A (ja) | 長距離地盤熱物性量計測装置 | |
CN112268926B (zh) | 一种冻土含冰量测量方法及装置 | |
CN214310255U (zh) | 一种原位监测土壤非饱和导水率的装置 | |
Liu et al. | In situ soil moisture and thermal properties estimated using a dual-probe heat-pulse | |
Bassi et al. | New sensors for moisture monitoring in historic walls: preliminary results | |
Valente et al. | On-chip integrated silicon bulk-micromachined soil moisture sensor based on the DPHP method |
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 |