CN114033362B - 一种基于录井和测井数据的页岩气声波曲线重建方法 - Google Patents

一种基于录井和测井数据的页岩气声波曲线重建方法 Download PDF

Info

Publication number
CN114033362B
CN114033362B CN202111260646.XA CN202111260646A CN114033362B CN 114033362 B CN114033362 B CN 114033362B CN 202111260646 A CN202111260646 A CN 202111260646A CN 114033362 B CN114033362 B CN 114033362B
Authority
CN
China
Prior art keywords
curve
coordinate
logging
coordinate value
logging data
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
CN202111260646.XA
Other languages
English (en)
Other versions
CN114033362A (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.)
Nie Lin
Original Assignee
Chengdu Yiyetulan Technology Co ltd
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 Chengdu Yiyetulan Technology Co ltd filed Critical Chengdu Yiyetulan Technology Co ltd
Priority to CN202111260646.XA priority Critical patent/CN114033362B/zh
Publication of CN114033362A publication Critical patent/CN114033362A/zh
Application granted granted Critical
Publication of CN114033362B publication Critical patent/CN114033362B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • E21B47/14Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • E21B43/25Methods for stimulating production
    • E21B43/26Methods for stimulating production by forming crevices or fractures
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Forestry; Mining
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Theoretical Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Geometry (AREA)
  • Acoustics & Sound (AREA)
  • Software Systems (AREA)
  • Computer Hardware Design (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • General Engineering & Computer Science (AREA)
  • Artificial Intelligence (AREA)
  • Remote Sensing (AREA)
  • Agronomy & Crop Science (AREA)
  • Animal Husbandry (AREA)
  • Marine Sciences & Fisheries (AREA)
  • Medical Informatics (AREA)
  • Health & Medical Sciences (AREA)
  • Economics (AREA)
  • General Health & Medical Sciences (AREA)
  • Human Resources & Organizations (AREA)
  • Marketing (AREA)
  • Primary Health Care (AREA)
  • Strategic Management (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于录井和测井数据的页岩气声波曲线重建方法,其包括:选择页岩气中的非烃段和含烃段;确定需要重建的构建曲线;根据构建曲线的右坐标初始值AIN建立构建曲线的最佳右坐标的选择范围;利用最佳右坐标的选择范围内的右坐标值计算对应的左坐标值INT(Ayi);建立最佳右坐标的选择范围内的右坐标值Ax、与其对应的最佳左坐标值和目标声波曲线的偏离度指标Di的关系模型;得到S个偏离度值;求得的构建曲线的最终左坐标Aright'和最终右坐标Aleft';建立转换方程,得到最终构建的新的声波曲线Bnew。本发明研发出了一种适用于页岩气水平井的测井声波重建技术,可以获得准确的岩石力学参数,为压裂方案设计、压后评估提供重要的数据支撑。

Description

一种基于录井和测井数据的页岩气声波曲线重建方法
技术领域
本发明涉及页岩气开发研究技术领域,具体涉及一种基于录井和测井数据的页岩气声波曲线重建方法。
背景技术
随着国内非常规油气资源勘探开发取得重要突破,储层特低孔渗的特征决定了体积压裂举足轻重的开发地位。压裂设计又是该技术体系中一个重要的技术环节,取得较准确的岩石力学参数是压裂设计的重要前提,同时也是完井后评估和产能评价的基础。
目前多采用完井测井的技术手段来取得岩石力学参数,但是由于测井仪器设计和解释的局限性,特别是测井采集的声波资料受到含气页岩和井筒含气的影响尤为突出,造成采集的测井信息失真,从而导致计算的岩石力学参数不准确,造成压裂设计分段选簇及施工参数等不合理,严重影响了压裂效果。
发明内容
针对现有技术中的上述不足,本发明提供了一种声波重建精度高的基于录井和测井数据的页岩气声波曲线重建方法。
为了达到上述发明目的,本发明所采用的技术方案为:
提供一种基于录井和测井数据的页岩气声波曲线重建方法,其包括以下步骤:
S1:根据页岩气钻井的录井和测井数据曲线,在数据曲线中选择页岩气中的非烃段和含烃段;
S2:选择与非烃段相关,且与含烃段不相关的曲线段作为需要重建的构建曲线;
S3:根据构建曲线的右坐标初始值AIN建立构建曲线的最佳右坐标的选择范围:(AIN-n,AIN-n+1,···,AIN,AIN+1,···,AIN+m);其中,n和m均为整数;
S4:利用最佳右坐标的选择范围内的右坐标值计算对应的左坐标值INT(Ayi);
S5:利用最佳左坐标值INT(Ayi)建立最佳左坐标的选择范围:{INT(Ayi)-n,INT(Ayi)-n+1,···,INT(Ayi),INT(Ayi)+1,···,INT(Ayi)+m};每个最佳左坐标的选择范围内的左坐标值INT(Ayi)对应最佳右坐标的选择范围的一个右坐标值Ax
S6:建立最佳右坐标的选择范围内的右坐标值Ax、与其对应的最佳左坐标值和目标声波曲线的偏离度指标Di的关系模型:
S7:利用关系模型计算构建曲线上每一点测井数据在左坐标值INT(Ayi)和右坐标值条件下Ax与目标声波曲线的偏离度D(A-B)i,得到S个偏离度值:
其中,s为构建曲线上测井数据的量;
S8:提取S个偏离度值中的最小值所对应的左坐标值和右坐标值,作为求得的构建曲线的最终左坐标Aright'和最终右坐标Aleft';
S9:根据目标声波曲线、最终左坐标Aright'和最终右坐标Aleft'建立转换方程,得到最终构建的新的声波曲线Bnew
其中,A为构建曲线。
进一步地,步骤S4包括:
S41:计算最佳右坐标的选择范围内每个右坐标值Ax对应的左坐标值Ay,得到m+n个左坐标值:
其中,RB为目标声波曲线的测井数据分别到右坐标和左坐标的距离之比,Ax为最佳右坐标的选择范围内的任意一个右坐标值,Ai为构建曲线中的第i个测井数据;
S42:根据均值定理计算m+n个左坐标值的理想值Ayi
S43:对理想值Ayi进行取整,得到最佳左坐标值INT(Ayi)。
进一步地,目标声波曲线的测井数据分别到右坐标和左坐标的距离之比RB的计算方法为:
其中,Lleft为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离,Lright为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离。
进一步地,目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离Lleft的计算方法为:
Lleft=Bi-Bleft
目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离Lright的计算方法为:
Lrightt=Bright-Bi
本发明的有益效果为:本发明研发出了一种适用于页岩气水平井的测井声波重建技术,可以获得准确的岩石力学参数,为压裂方案设计、压后评估提供重要的数据支撑。有效解决了由于测井仪器、气页岩和井筒含气的影响等方面造成采集的测井信息失真问题,实现精确的对出声波进行重建,消除误差。为后续各项岩石力学和物性参数的计算,提供了真实的声波数据,较现有技术而言,整个方法简单,效率高。
附图说明
图1为基于录井和测井数据的页岩气声波曲线重建方法的流程图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
本方案提供的基于录井和测井数据的页岩气声波曲线重建方法包括以下步骤:
S1:根据页岩气钻井的录井和测井数据曲线,在数据曲线中选择页岩气中的非烃段和含烃段;根据岩屑深颜、气测升高和钻时加快等录井特征,以及测井的低密度、高铀等测井表现,可以确定页岩非烃段和含烃段。
S2:选择与非烃段相关,且与含烃段不相关的曲线段作为需要重建的构建曲线;证明此构建曲线具有较大的误差,需要进行重建。
S3:根据构建曲线的右坐标初始值AIN建立构建曲线的最佳右坐标的选择范围:(AIN-n,AIN-n+1,···,AIN,AIN+1,···,AIN+m);其中,n和m均为整数;
S4:利用最佳右坐标的选择范围内的右坐标值计算对应的左坐标值INT(Ayi);具体包括如下步骤:
S41:计算最佳右坐标的选择范围内每个右坐标值Ax对应的左坐标值Ay,得到m+n个左坐标值:
其中,RB为目标声波曲线的测井数据分别到右坐标和左坐标的距离之比,Ax为最佳右坐标的选择范围内的任意一个右坐标值,Ai为构建曲线中的第i个测井数据;目标声波曲线的测井数据分别到右坐标和左坐标的距离之比RB的计算方法为:目标声波曲线是从实验室通过标准参数得到的标准曲线。
其中,Lleft为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离,Lright为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离。
目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离Lleft的计算方法为:
Lleft=Bi-Bleft
目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离Lright的计算方法为:
Lrightt=Bright-Bi
S42:根据均值定理计算m+n个左坐标值的理想值Ayi
S43:对理想值Ayi进行取整,得到最佳左坐标值INT(Ayi)。
S5:利用最佳左坐标值INT(Ayi)建立最佳左坐标的选择范围:{INT(Ayi)-n,INT(Ayi)-n+1,···,INT(Ayi),INT(Ayi)+1,···,INT(Ayi)+m};每个最佳左坐标的选择范围内的左坐标值INT(Ayi)对应最佳右坐标的选择范围的一个右坐标值Ax
S6:建立最佳右坐标的选择范围内的右坐标值Ax、与其对应的最佳左坐标值和目标声波曲线的偏离度指标Di的关系模型:
S7:利用关系模型计算构建曲线上每一点测井数据在左坐标值INT(Ayi)和右坐标值条件下Ax与目标声波曲线的偏离度D(A-B)i,得到S个偏离度值:
其中,s为构建曲线上测井数据的量;
S8:提取S个偏离度值中的最小值所对应的左坐标值和右坐标值,作为求得的构建曲线的最终左坐标Aright'和最终右坐标Aleft';
S9:根据目标声波曲线、最终左坐标Aright'和最终右坐标Aleft'建立转换方程,得到最终构建的新的声波曲线Bnew
其中,A为构建曲线。
本发明研发出了一种适用于页岩气水平井的测井声波重建技术,可以获得准确的岩石力学参数,为压裂方案设计、压后评估提供重要的数据支撑。有效解决了由于测井仪器、气页岩和井筒含气的影响等方面造成采集的测井信息失真问题,实现精确的对出声波进行重建,消除误差。为后续各项岩石力学和物性参数的计算,提供了真实的声波数据,较现有技术而言,整个方法简单,效率高。

Claims (4)

1.一种基于录井和测井数据的页岩气声波曲线重建方法,其特征在于,包括以下步骤:
S1:根据页岩气钻井的录井和测井数据曲线,在数据曲线中选择页岩气中的非烃段和含烃段;
S2:选择与非烃段相关,且与含烃段不相关的曲线段作为需要重建的构建曲线;
S3:根据构建曲线的右坐标初始值AIN建立构建曲线的最佳右坐标的选择范围:(AIN-n,AIN-n+1,···,AIN,AIN+1,···,AIN+m);其中,n和m均为整数;
S4:利用最佳右坐标的选择范围内的右坐标值计算对应的左坐标值INT(Ayi);
S5:利用最佳左坐标值INT(Ayi)建立最佳左坐标的选择范围:{INT(Ayi)-n,INT(Ayi)-n+1,···,INT(Ayi),INT(Ayi)+1,···,INT(Ayi)+m};每个最佳左坐标的选择范围内的左坐标值INT(Ayi)对应最佳右坐标的选择范围的一个右坐标值Ax
S6:建立最佳右坐标的选择范围内的右坐标值Ax、与其对应的最佳左坐标值和目标声波曲线的偏离度指标Di的关系模型:
S7:利用关系模型计算构建曲线上每一点测井数据在左坐标值INT(Ayi)和右坐标值条件下Ax与目标声波曲线的偏离度D(A-B)i,得到S个偏离度值:
其中,s为构建曲线上测井数据的量;
S8:提取S个偏离度值中的最小值所对应的左坐标值和右坐标值,作为求得的构建曲线的最终左坐标Aright'和最终右坐标Aleft';
S9:根据目标声波曲线、最终左坐标Aright'和最终右坐标Aleft'建立转换方程,得到最终构建的新的声波曲线Bnew
其中,A为构建曲线。
2.根据权利要求1所述的基于录井和测井数据的页岩气声波曲线重建方法,其特征在于,所述步骤S4包括:
S41:计算最佳右坐标的选择范围内每个右坐标值Ax对应的左坐标值Ay,得到m+n个左坐标值:
其中,RB为目标声波曲线的测井数据分别到右坐标和左坐标的距离之比,Ax为最佳右坐标的选择范围内的任意一个右坐标值,Ai为构建曲线中的第i个测井数据;
S42:根据均值定理计算m+n个左坐标值的理想值Ayi
S43:对理想值Ayi进行取整,得到最佳左坐标值INT(Ayi)。
3.根据权利要求2所述的基于录井和测井数据的页岩气声波曲线重建方法,其特征在于,所述目标声波曲线的测井数据分别到右坐标和左坐标的距离之比RB的计算方法为:
其中,Lleft为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离,Lright为目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离。
4.根据权利要求3所述的基于录井和测井数据的页岩气声波曲线重建方法,其特征在于,所述目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bleft的距离Lleft的计算方法为:
Lleft=Bi-Bleft
所述目标声波曲线中任意测井数据Bi到目标声波曲线的左坐标Bright的距离Lright的计算方法为:
Lrightt=Bright-Bi
CN202111260646.XA 2021-10-28 2021-10-28 一种基于录井和测井数据的页岩气声波曲线重建方法 Active CN114033362B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111260646.XA CN114033362B (zh) 2021-10-28 2021-10-28 一种基于录井和测井数据的页岩气声波曲线重建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111260646.XA CN114033362B (zh) 2021-10-28 2021-10-28 一种基于录井和测井数据的页岩气声波曲线重建方法

Publications (2)

Publication Number Publication Date
CN114033362A CN114033362A (zh) 2022-02-11
CN114033362B true CN114033362B (zh) 2024-04-26

Family

ID=80135650

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111260646.XA Active CN114033362B (zh) 2021-10-28 2021-10-28 一种基于录井和测井数据的页岩气声波曲线重建方法

Country Status (1)

Country Link
CN (1) CN114033362B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932836A (zh) * 2015-12-30 2017-07-07 中国石油化工股份有限公司 一种用于评价页岩气含气丰度的方法及系统
US10689954B1 (en) * 2019-04-24 2020-06-23 Dagang Oil Field Company Of Cnpc Research method of trajectory design and on-site tracking and adjustment of shale oil horizontal well
CN111650646A (zh) * 2020-05-26 2020-09-11 中国石油天然气集团有限公司 一种用于地震处理中的声波曲线校正方法
CN111927445A (zh) * 2020-09-06 2020-11-13 中国石油天然气集团有限公司 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106932836A (zh) * 2015-12-30 2017-07-07 中国石油化工股份有限公司 一种用于评价页岩气含气丰度的方法及系统
US10689954B1 (en) * 2019-04-24 2020-06-23 Dagang Oil Field Company Of Cnpc Research method of trajectory design and on-site tracking and adjustment of shale oil horizontal well
CN111650646A (zh) * 2020-05-26 2020-09-11 中国石油天然气集团有限公司 一种用于地震处理中的声波曲线校正方法
CN111927445A (zh) * 2020-09-06 2020-11-13 中国石油天然气集团有限公司 基于随钻岩矿数据拟合测井参数获取随钻地层压力的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
王瀚玮 ; 夏宏泉 ; 刘畅 ; 赵昊 ; .页岩储层脆性指数的随钻测井计算方法研究――以威远地区寒武系筇竹寺组为例.油气藏评价与开发.2018,(第03期),全文. *

Also Published As

Publication number Publication date
CN114033362A (zh) 2022-02-11

Similar Documents

Publication Publication Date Title
CN105781539B (zh) 一种致密油气储层饱和度测井计算方法
CN102353985B (zh) 基于非下采样Contourlet变换的拟声波曲线构建方法
CN108535773B (zh) 一种微地震事件检测方法及系统
CN102454399B (zh) 测井声波时差信号校正方法
CN104297785A (zh) 岩相约束储层物性参数反演方法及装置
CN107966732B (zh) 基于空间结构导向的地震属性变化率求取方法
CN106886046B (zh) 确定缝洞型气藏未投产区块可动用储量的方法
CN109900617B (zh) 一种基于声电成像测井图的裂缝性地层渗透率曲线计算方法
CN110955982B (zh) 一种变质岩储层渗透率计算方法、装置及计算机存储介质
CN108694264A (zh) 一种确定页岩气储层渗透率的方法
CN109829217A (zh) 压裂性裂缝油藏产能模拟方法及装置
CN109100793A (zh) 一种定量分析裂缝因素对储层影响的方法
CN107808068A (zh) 一种用于致密砂岩气藏高产富集评价的方法
CN111456709A (zh) 一种基于测井曲线的水平井多级压裂分段分簇方法
CN108871241A (zh) 一种定量评价粗糙裂缝酸刻蚀后粗糙程度变化的方法
CN108729912B (zh) 适用于油藏数值模拟的产量劈分方法
CN112796738A (zh) 一种阵列声波测井和常规测井结合的地层渗透率计算方法
CN114033362B (zh) 一种基于录井和测井数据的页岩气声波曲线重建方法
CN109582922B (zh) 一种水合物饱和度的现场快速判别方法及处理终端
CN110609328B (zh) 一种基于井震结合的含油饱和度预测方法
CN114861400A (zh) 一种绘制高精度水底地形与淤泥厚度双界面分布图的方法
CN106709992A (zh) 基于信息熵理论的海底地貌复杂度表示方法
CN112182992A (zh) 一种致密砂岩气藏水平井分段压裂裂缝位置优选方法
CN116341406B (zh) 一种确定倾斜油藏注采干扰时间的计算方法
CN106353837B (zh) 基于加权gardner公式的密度预测方法

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
TR01 Transfer of patent right

Effective date of registration: 20240520

Address after: No. 168, Section 1, Shiyou Road, Longchang City, Neijiang City, Sichuan Province, 642150, Building 8, Unit 3, 5th Floor, No.1

Patentee after: Nie Lin

Country or region after: China

Address before: No. 2420, floor 24, building 5, No. 46, erxianqiao East Road, Chenghua District, Chengdu, Sichuan 610000

Patentee before: Chengdu yiyetulan Technology Co.,Ltd.

Country or region before: China