CN103645514B - 多同性源电极阵列电阻率的地下工程超前探测方法及系统 - Google Patents

多同性源电极阵列电阻率的地下工程超前探测方法及系统 Download PDF

Info

Publication number
CN103645514B
CN103645514B CN201310729095.6A CN201310729095A CN103645514B CN 103645514 B CN103645514 B CN 103645514B CN 201310729095 A CN201310729095 A CN 201310729095A CN 103645514 B CN103645514 B CN 103645514B
Authority
CN
China
Prior art keywords
electrode
workplace
array
potential electrode
electrodes
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
CN201310729095.6A
Other languages
English (en)
Other versions
CN103645514A (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.)
Shandong Bai20 Huitong Engineering Technology Co ltd
Original Assignee
Shandong University
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 Shandong University filed Critical Shandong University
Priority to CN201310729095.6A priority Critical patent/CN103645514B/zh
Publication of CN103645514A publication Critical patent/CN103645514A/zh
Application granted granted Critical
Publication of CN103645514B publication Critical patent/CN103645514B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种基于多同性源电极阵列电阻率的地下工程超前探测方法，适用于隧道或坑道等地下工程钻爆法或TBM施工过程中的不良地质探测：（1）在工作面上构建阵列式测量电极系；（2）在边墙上组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方；（3）利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取工作面及各断面中各测量电极M与测量电极N间的电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；（4）采用基于光滑约束的最小二乘线性反演迭代处理，将测得的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。

Description

多同性源电极阵列电阻率的地下工程超前探测方法及系统
技术领域
本发明涉及隧道或者坑道等地下工程钻爆法或TBM施工过程中不良地质的探测方法，尤其涉及一种基于多同性源电极阵列电阻率的地下工程超前探测方法及系统。
背景技术
随着水利水电工程和交通工程的建设重点逐步向地形地质条件极端复杂的西部山区和岩溶地区转移，所遇到的工程地质条件异常复杂，施工中面临着突水、突泥等诸多高风险地质灾害。突涌水灾害给施工安全带来了重大灾难和生命财产损失，而预先探明隧道掌子面前方的地质灾害赋存源是规避和防治重大灾害事故的有效途径，因此开展隧道超前预报工作研究显得尤为迫切和重要。
作为一种重要的工程地球物理探测技术，直流电阻率法以其适应性强、方便快捷、成本低等诸多优点被广泛应用于岩土工程勘察、滑坡地质灾害监测、环境岩土工程评价等领域中。与地面直流电阻率法相比，由于隧道掌子面空间狭小，导致地面的很多探测方法无法应用于隧道突涌水灾害源超前预报。目前直流电阻率法在隧道（巷道）超前预报中大都采用三极法的探测模式，即将供电点A置于掌子面，供电点B置于隧道后方“无穷远”处，测量电极M、N沿着隧道轴向移动测量。
这一方法主要面临以下问题：①由于测量电极MN在掌子面后方移动测量，易受到掌子面后方异常体的干扰；②由于隧道的特殊环境，隧道三极法的探测模式测线垂直于异常体，导致对掌子面近处的异常体响应敏感，对远处的异常体响应相对不敏感，且供电与接收属于一维的探测模式，分辨率不高。
发明内容
本发明的目的是为了克服上述现有技术的不足，提出了一种动电源点供电、阵列电极接收的基于多同性源电极阵列电阻率的地下工程超前探测方法，适用于隧道或坑道等地下工程钻爆法或TBM施工过程中的不良地质探测。利用多同性源供电的方式来有效压制地下工程工作面后方异常体的干扰，并抑制工作面近处异常体的敏感性，从而相对提高了工作面远处异常体的敏感性，阵列式测量电极的接收模式提高了对异常体的空间分辨率。
本发明采用的技术方案为：
一种多同性源电极阵列电阻率的地下工程超前探测方法，它包括以下步骤：
（1）在地下工程工作面上布置多条平行的测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
（2）沿地下工程工作面轮廓线、在阵列式测量电极系的外围均匀布多个供电电极A，然后以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后远方；
（3）利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取各测量电极M与测量电极N间的电势差；然后依次选择工作面后方下一工作断面上的供电电极A提供相同大小的正电流，并继续测量阵列式测量电极系的电势，上位机计算相应电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；
（4）上位机采用基于光滑约束的最小二乘线性反演迭代处理，将步骤（3）中测得的阵列式测量电极系的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。
所述步骤（4）中，反演目标函数如下：
Φ＝(Δd-AΔm)T(Δd-AΔm)+λ(CΔm)T(CΔm)    （1）
式中，Φ为反演目标函数，Δd为观测数据的差值，A为偏导数矩阵，Δm为模型参数修正量，C为光滑度矩阵，λ为拉格朗日常数，决定了光滑约束的权重；
对公式（1）进行变分之后得到具有光滑约束的三维电阻率反演方程：
(ATA+λCTC)Δm＝ATΔd    （2）
m(k+1)＝m(k)+Δm    （3）
式中，m(k)为第k次迭代中的模型参数向量，m(k+1)为下一次反演迭代中的模型参数向量。
进行反演成像时，首先给模型参数赋初值并进行正演计算，将正演得到的理论观测数据与实际测量数据进行比较，如果方差小于用户的设定值，则输出模型参数；
否则，计算偏导数矩阵A、光滑度矩阵C、以及理论观测数据与实际测量数据的差值Δd，代入方程（2）求得每次反演迭代中的模型参数修正量Δm，再代入方程（3）得到下一次迭代中的模型参数m(k+1)，直至正演的理论观测数据与测量数据的方差小于设定值，则输出模型参数，从而通过反演得到工作面前方地质体的三维电阻率图像。
一种多同性源电极阵列电阻率的地下工程超前探测系统，它包括布置在地下工程工作面上的多条平行测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
在地下工程工作面轮廓线以及阵列式测量电极系的外围均匀布多个供电电极A，以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；
同时，将另一个供电电极B与测量电极N都置于工作面后远方；
所述各供电电极均与电源连接，所述各测量电极则与上位机连接。
所述平行测线有3条，分别为L1、L2和L3，它们布置在距离底部和的高度上，其中H为拱顶距离底部的高度。
所述供电电极B与测量电极N置于工作面后方距离工作面无穷远处，例如100-150m处，此处的100-150m，是无穷远位置，大于5倍的AM即可，100m和150m都是大于的，一般认为结果是没有差异的。
本发明具有以下有益效果：
（1）利用多同性源供电可以有效地压制地下工程工作面后方异常体的干扰，并能够降低掌子面近处异常体的敏感性，从而使得工作面远处异常体的敏感性相对提高；
（2）采用阵列式测量电极的接收模式，有效地提高了对异常体的空间分辨率，增大了探测距离。
（3）采用多个同性源供电电极，与传统三极法相比，降低了掌子面近处异常比，提高了掌子面原处的异常，如表1所示。对于掌子面近处的异常体（4m），同性源异常比13.1%小于传统定点源异常比16.1%，降低了掌子面近处异常体的敏感性；对于掌子面远处的异常体同性源异常比均大于传统定点源异常，从而提高了对远处异常体的探测能力。
表1不同探测距离异常对比
附图说明
图1是本发明的探测方法流程图。
图2是本发明的探测示意图。
图3是本发明的三维电阻率反演流程。
图4是本发明的三维电阻率反演成像图。
图5是本发明实施例3的三维成像结果图。
图中，1是地下工程工作面，2是阵列式测量电极系，3是供电电极系。
具体实施方式
下面结合附图与实施例对本发明做进一步说明。
本发明探测方法如图1中，流程包括阵列式测量电极系和供电电极系的布设、数据采集、三维电阻率反演等流程，可以通过反演得到工作面前方地质体的三维电阻率图像。
实施例1：
以某隧道为例，隧道断面为拱形，跨度12m高8m，里程ZK97+853附近出现涌水，前方可能存在含水构造，为探明前方不良地质情况，为施工的顺利进行提供保障，采用基于多同性源电极阵列的电阻率方法进行超前探测。具体工作方式如图2所示：
（1）在隧道掌子面上布置3条平行测线L1、L2和L3，分别布置在距离底板3m、6m和9m的高度上，每条测线上的测量电极M的数目为13个，组成阵列式测量电极系；
（2）沿掌子面轮廓线、在阵列式测量电极系的外围均匀布置4个供电电极A，然后从这4个供电电极A出发在边墙上分别布置4条测线，每条测线上的供电电极A的数目为5个，距离工作面的距离依次为0m、2m、6m、14m和30m，，组成供电电极系；
（3）另一个供电电极B置于掌子面后方150m的底板上，另一个测量电极N置于掌子面后方120m的底板上；
（4）探测时，首先掌子面上的4个供电电极A供入相同大小的正电流0.1A，供电电极B为负极，阵列式测量电极系中的每个测量电极M分别测量电势，并记录下来每个测量电极M与测量电极N间的电势差UMN，然后依次选择工作面后方下一断面上的4个供电电极A供入相同大小的正电流0.1A，并继续测量阵列式测量电极系的电势，不断重复该步骤，直至供电电极系移动到预定位置，探测结束；
（5）采用常规的基于光滑约束的最小二乘线性方法对采集到的电势差数据UMN进行反演，基本流程如图3所示，首先设定模型参数初值，本次初始电阻率值设为1000Ω·m，进行正演计算，将正演得到的理论观测数据与实际测量数据进行比较，如果方差小于用户的设定值（本实例相对值为3%），则输出模型参数。
否则，计算偏导数矩阵A、光滑度矩阵C、以及理论观测数据与实际测量数据的差值Δd，代入方程（2）求得每次反演迭代中的模型参数修正量Δm，再代入方程（3）得到下一次迭代中的模型参数m(k+1)，直至正演的理论观测数据与测量数据的方差小于设定值3%，则输出模型参数，从而通过反演得到掌子面前方地质体的三维电阻率图像。
本次探测的三维成像结果如图3所示，在ZK97+879～ZK97+883段落处围岩富水性较强。
实施例2：
以某隧道为例，隧道断面为拱形，跨度12m高8m，里程ZK97+853附近出现涌水，重复试例步骤（1）～（5），本次探测三维成像结果如图4所示，在ZK97+871～ZK97+875段落处围岩富水性较强。
实施例3：
以某隧道为例，隧道断面为拱形，跨度12m高8m，里程ZK97+853附近出现涌水，重复试例步骤（1）～（5），本次探测三维成像结果如图5所示，在ZK97+874～ZK97+875段落处围岩富水性较强。

Claims (3)

1.一种多同性源电极阵列电阻率的地下工程超前探测方法，其特征是，它包括以下步骤：
(1)在地下工程工作面上布置多条平行的测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
(2)沿地下工程工作面轮廓线、在阵列式测量电极系的外围均匀布多个供电电极A，然后以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后方，距离工作面100-150米处；
(3)利用位于工作面的供电电极A以及远方的供电电极B提供正负电流，利用阵列式测量电极系中的每个测量电极M分别测量电势，并获取各测量电极M与测量电极N间的电势差；然后依次选择工作面后方下一工作断面上的供电电极A提供相同大小的正电流，并继续测量阵列式测量电极系的电势，上位机计算相应电势差；不断重复本步骤，直至供电电极系移动到预定位置，探测结束；
(4)上位机采用基于光滑约束的最小二乘线性反演迭代处理，将步骤(3)中测得的阵列式测量电极系的电势差数据进行反演，得到工作面前方地质体的三维电阻率图像。
2.一种多同性源电极阵列电阻率的地下工程超前探测系统，其特征是，它包括布置在地下工程工作面上的多条平行测线，在各测线上设置数量相同、间距为d的测量电极M，构成阵列式测量电极系；
在地下工程工作面轮廓线以及阵列式测量电极系的外围均匀布多个供电电极A，以这些供电电极A为出发点在边墙上布置供电测线，并在每条供电测线上再布置若干间隔相同的位于同一工作断面上的供电电极A，从而组成供电电极系；同时，将另一个供电电极B与测量电极N都置于工作面后方，距离工作面100-150米处；
所述供电电极A和供电电极B均与电源连接，所述测量电极M和测量电极N则与上位机连接。
3.如权利要求2所述的基于多同性源电极阵列电阻率的地下工程超前探测系统，其特征是，所述平行测线有3条，分别为L1、L2和L3，它们布置在距离底部和的高度上，其中H为拱顶距离底部的高度。
CN201310729095.6A 2013-12-25 2013-12-25 多同性源电极阵列电阻率的地下工程超前探测方法及系统 Active CN103645514B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310729095.6A CN103645514B (zh) 2013-12-25 2013-12-25 多同性源电极阵列电阻率的地下工程超前探测方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310729095.6A CN103645514B (zh) 2013-12-25 2013-12-25 多同性源电极阵列电阻率的地下工程超前探测方法及系统

Publications (2)

Publication Number Publication Date
CN103645514A CN103645514A (zh) 2014-03-19
CN103645514B true CN103645514B (zh) 2015-01-07

Family

ID=50250762

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310729095.6A Active CN103645514B (zh) 2013-12-25 2013-12-25 多同性源电极阵列电阻率的地下工程超前探测方法及系统

Country Status (1)

Country Link
CN (1) CN103645514B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114088782A (zh) * 2021-10-18 2022-02-25 中国矿业大学 一种应力与渗流作用下煤岩体突水危险区域电位判识方法

Families Citing this family (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104678443A (zh) * 2015-01-22 2015-06-03 深圳市市政设计研究院有限公司 一种孤石探测仪及其控制方法
CN108845361B (zh) * 2018-06-06 2019-07-23 中南大学 一种多个供电点同时供电的激电扫面方法
CN108845362A (zh) * 2018-06-06 2018-11-20 中南大学 一种多个供电点同时供电的激电扫面方法
CN109001816A (zh) * 2018-06-06 2018-12-14 中南大学 一种多对供电点同时供电的激电扫面方法
CN108828678B (zh) * 2018-08-25 2020-05-29 安徽省公路工程检测中心 一种隧道施工超前地质探测系统
CN110824568B (zh) * 2019-11-14 2021-10-08 山东大学 一种搭载于盾构机刀盘-内置聚焦式电法探水系统及方法
CN111880233B (zh) * 2020-07-17 2023-06-13 中煤科工集团西安研究院有限公司 一种煤矿井下水平钻孔中直流电超前透视探测装置及方法
CN111650652B (zh) * 2020-07-17 2022-09-30 安徽理工大学 坑道三方向视电阻率超前探测方法
CN112505790B (zh) * 2020-09-27 2021-12-24 山东大学 高空间分辨率复合式电阻率探测计算方法
CN112415602B (zh) * 2020-10-15 2021-11-19 山东大学 基于深度分辨率的隧道电阻率超前探测优化方法及系统
CN112415605B (zh) * 2020-10-16 2022-05-13 山东大学 多路屏蔽型电缆、激发极化探测装置及方法
CN112505786B (zh) * 2020-11-20 2023-04-14 安徽理工大学 坑道三方向视极化率超前探测方法
CN112666614B (zh) * 2021-01-15 2022-09-06 西南科技大学 基于电法勘探和数字高程模型泥石流物源静储量计算方法
CN113311487B (zh) * 2021-05-18 2022-08-09 山东大学 隧道掌子面发射的频率域激发极化超前探水方法与装置
CN114778948A (zh) * 2022-06-17 2022-07-22 中铁大桥科学研究院有限公司 动水隧道岩体电阻率监测方法及相关设备

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4321540A (en) * 1977-09-01 1982-03-23 Compagnie Generale De Geophysique Electric prospecting of the subsoil with lineal electrodes
CN102508303A (zh) * 2011-11-23 2012-06-20 山东大学 地下工程聚焦层析激发极化超前探测方法
CN102749653A (zh) * 2012-07-30 2012-10-24 山东大学 用于地下工程高分辨率三维电阻率ct成像的电极递送装置
CN103278857A (zh) * 2013-05-13 2013-09-04 江苏大学 一种井下直流超前探测装置的设计方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2001166061A (ja) * 1999-12-10 2001-06-22 Ohbayashi Corp トンネル切羽前方の探査方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4321540A (en) * 1977-09-01 1982-03-23 Compagnie Generale De Geophysique Electric prospecting of the subsoil with lineal electrodes
CN102508303A (zh) * 2011-11-23 2012-06-20 山东大学 地下工程聚焦层析激发极化超前探测方法
CN102749653A (zh) * 2012-07-30 2012-10-24 山东大学 用于地下工程高分辨率三维电阻率ct成像的电极递送装置
CN103278857A (zh) * 2013-05-13 2013-09-04 江苏大学 一种井下直流超前探测装置的设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
基于激发极化法的隧道含水地质构造超前探测研究;李术才 等;《岩石力学与工程学报》;20110731;第30卷(第7期);第1297-1309页 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114088782A (zh) * 2021-10-18 2022-02-25 中国矿业大学 一种应力与渗流作用下煤岩体突水危险区域电位判识方法
CN114088782B (zh) * 2021-10-18 2022-05-17 中国矿业大学 一种应力与渗流作用下煤岩体突水危险区域电位判识方法

Also Published As

Publication number Publication date
CN103645514A (zh) 2014-03-19

Similar Documents

Publication Publication Date Title
CN103645514B (zh) 多同性源电极阵列电阻率的地下工程超前探测方法及系统
CN102508303B (zh) 地下工程聚焦层析激发极化超前探测方法
US10613245B2 (en) Mine TEM three-component detection method
CN101706585B (zh) 一种用于地下掘进工程中的隐患电法超前预报方法
CN103076635B (zh) Tbm施工隧道前向三维激发极化法超前探测装置系统及方法
CN104330838B (zh) 地下工程突涌水通道精细探查及其注浆封堵效果评价方法
CN106772621B (zh) 一种近全方位电阻率隧道超前地质预报方法
CN108957563B (zh) 一种隧道施工超前地质探测系统以及探测方法
CN109143378B (zh) 一种在煤矿巷道内顺层超前探测含水构造的二次时差方法
US20230184983A1 (en) Vector-resistivity-based real-time advanced detection method for water-bearing hazard body
CN108873073B (zh) 一种基于网络并行电法的三维跨孔电阻率层析成像方法
CN106869913B (zh) 一种利用井地电位技术探测油田注水井水驱前缘的方法
CN109343131B (zh) 一种侧向供电的地井直流电法探测方法与装置
CN111650652B (zh) 坑道三方向视电阻率超前探测方法
Kirsch et al. Geoelectrical methods
CN104820247B (zh) 矿井分布式三维电法仪及其观测方法
CN105044789B (zh) 一种矿井瞬变电磁超前探测的巷道效应消除方法和装置
CN104502986B (zh) 物探激电测深数据层析法处理方法
CN110080740B (zh) 一种检测井下煤层气水力压裂效果的方法及系统
CN106405672A (zh) 一种tbm施工隧道三维激发极化超前预报干扰去除方法
RU2340918C2 (ru) Способ геоэлектроразведки
CN205958772U (zh) 基于网络并行电法技术的方位电测井装置
CN104749645A (zh) 一种煤层基底奥灰顶界面起伏形态探测的方法和装置
CN111856594B (zh) 坑道三方向视频散率超前探测方法
Yang et al. Dynamic monitoring of mining destruction on coal seam floor with constrained time-lapse resistivity imaging inversion

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20140319

Assignee: Shandong Zhiyan Exploration Technology Co.,Ltd.

Assignor: Shandong University

Contract record no.: 2019370000006

Denomination of invention: Underground construction advanced detection method and system based on resistivity of multiple same source electrode arrays

Granted publication date: 20150107

License type: Exclusive License

Record date: 20190124

EE01 Entry into force of recordation of patent licensing contract
TR01 Transfer of patent right

Effective date of registration: 20220425

Address after: 250100 room 1002, floor 10, building a3-1, Hanyu Golden Valley, No. 7000, jingshidong Road, high tech Zone, Jinan, Shandong Province

Patentee after: Shandong Zhiyan Exploration Technology Co.,Ltd.

Address before: 250061, No. ten, No. 17923, Lixia District, Ji'nan City, Shandong Province

Patentee before: SHANDONG University

TR01 Transfer of patent right
EC01 Cancellation of recordation of patent licensing contract

Assignee: Shandong Zhiyan Exploration Technology Co.,Ltd.

Assignor: SHANDONG University

Contract record no.: 2019370000006

Date of cancellation: 20220507

EC01 Cancellation of recordation of patent licensing contract
TR01 Transfer of patent right

Effective date of registration: 20220513

Address after: 250100 Room 601, 6 / F, building a7-5, Hanyu Golden Valley, No. 7000, Jingshi Road, Jinan area, China (Shandong) pilot Free Trade Zone, Jinan City, Shandong Province

Patentee after: Shandong bai20 Huitong Engineering Technology Co.,Ltd.

Address before: 250100 room 1002, floor 10, building a3-1, Hanyu Golden Valley, No. 7000, jingshidong Road, high tech Zone, Jinan, Shandong Province

Patentee before: Shandong Zhiyan Exploration Technology Co.,Ltd.

TR01 Transfer of patent right
PE01 Entry into force of the registration of the contract for pledge of patent right

Denomination of invention: Underground Engineering Advance Detection Method and System for Resistivity of Multiple Source Electrode Arrays

Effective date of registration: 20230830

Granted publication date: 20150107

Pledgee: Qilu Bank Co.,Ltd. Jinan Kaiyuan Branch

Pledgor: Shandong bai20 Huitong Engineering Technology Co.,Ltd.

Registration number: Y2023980054639

PE01 Entry into force of the registration of the contract for pledge of patent right