CN106950347B - 一种评价泥页岩各组分体积的方法 - Google Patents
一种评价泥页岩各组分体积的方法 Download PDFInfo
- Publication number
- CN106950347B CN106950347B CN201710124289.1A CN201710124289A CN106950347B CN 106950347 B CN106950347 B CN 106950347B CN 201710124289 A CN201710124289 A CN 201710124289A CN 106950347 B CN106950347 B CN 106950347B
- Authority
- CN
- China
- Prior art keywords
- volume
- mud shale
- mineral
- kerogen
- model
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 239000011707 mineral Substances 0.000 claims abstract description 92
- 229910052500 inorganic mineral Inorganic materials 0.000 claims abstract description 87
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims abstract description 50
- 229910052799 carbon Inorganic materials 0.000 claims abstract description 50
- 239000011435 rock Substances 0.000 claims abstract description 43
- 239000011148 porous material Substances 0.000 claims abstract description 29
- 239000000470 constituent Substances 0.000 claims abstract description 26
- 238000005457 optimization Methods 0.000 claims abstract description 22
- 238000004458 analytical method Methods 0.000 claims abstract description 18
- 238000003062 neural network model Methods 0.000 claims abstract description 16
- 238000013528 artificial neural network Methods 0.000 claims abstract description 12
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 238000002474 experimental method Methods 0.000 claims abstract description 6
- 238000002790 cross-validation Methods 0.000 claims abstract description 5
- 238000012549 training Methods 0.000 claims description 19
- 238000002360 preparation method Methods 0.000 claims description 14
- 229910052683 pyrite Inorganic materials 0.000 claims description 11
- 239000011028 pyrite Substances 0.000 claims description 11
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 11
- 150000004760 silicates Chemical class 0.000 claims description 8
- 238000002441 X-ray diffraction Methods 0.000 claims description 7
- 239000002734 clay mineral Substances 0.000 claims description 7
- 238000001514 detection method Methods 0.000 claims description 7
- 238000012546 transfer Methods 0.000 claims description 7
- HEDRZPFGACZZDS-UHFFFAOYSA-N Chloroform Chemical compound ClC(Cl)Cl HEDRZPFGACZZDS-UHFFFAOYSA-N 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 6
- 238000005259 measurement Methods 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 6
- BVKZGUZCCUSVTD-UHFFFAOYSA-L Carbonate Chemical compound [O-]C([O-])=O BVKZGUZCCUSVTD-UHFFFAOYSA-L 0.000 claims description 5
- 229910001748 carbonate mineral Inorganic materials 0.000 claims description 5
- 239000004927 clay Substances 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 238000013461 design Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims description 2
- 239000004575 stone Substances 0.000 claims description 2
- XUIMIQQOPSSXEZ-UHFFFAOYSA-N Silicon Chemical compound [Si] XUIMIQQOPSSXEZ-UHFFFAOYSA-N 0.000 claims 1
- 229940125717 barbiturate Drugs 0.000 claims 1
- 239000003921 oil Substances 0.000 claims 1
- 229910052710 silicon Inorganic materials 0.000 claims 1
- 239000010703 silicon Substances 0.000 claims 1
- 230000004044 response Effects 0.000 abstract description 12
- 230000000694 effects Effects 0.000 description 12
- 239000003079 shale oil Substances 0.000 description 10
- 229930195733 hydrocarbon Natural products 0.000 description 7
- 150000002430 hydrocarbons Chemical class 0.000 description 7
- 238000005481 NMR spectroscopy Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 5
- 239000000203 mixture Substances 0.000 description 5
- 239000004215 Carbon black (E152) Substances 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000011161 development Methods 0.000 description 4
- 239000005416 organic matter Substances 0.000 description 4
- 235000013351 cheese Nutrition 0.000 description 3
- 238000009826 distribution Methods 0.000 description 3
- 235000021185 dessert Nutrition 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 229910052900 illite Inorganic materials 0.000 description 2
- VGIBGUSAECPPNB-UHFFFAOYSA-L nonaaluminum;magnesium;tripotassium;1,3-dioxido-2,4,5-trioxa-1,3-disilabicyclo[1.1.1]pentane;iron(2+);oxygen(2-);fluoride;hydroxide Chemical compound [OH-].[O-2].[O-2].[O-2].[O-2].[O-2].[F-].[Mg+2].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[Al+3].[K+].[K+].[K+].[Fe+2].O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2.O1[Si]2([O-])O[Si]1([O-])O2 VGIBGUSAECPPNB-UHFFFAOYSA-L 0.000 description 2
- 238000005325 percolation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- GUJOJGAPFQRJSV-UHFFFAOYSA-N dialuminum;dioxosilane;oxygen(2-);hydrate Chemical compound O.[O-2].[O-2].[O-2].[Al+3].[Al+3].O=[Si]=O.O=[Si]=O.O=[Si]=O.O=[Si]=O GUJOJGAPFQRJSV-UHFFFAOYSA-N 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 210000004884 grey matter Anatomy 0.000 description 1
- 239000011799 hole material Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- NLYAJNPCOHFWQQ-UHFFFAOYSA-N kaolin Chemical compound O.O.O=[Al]O[Si](=O)O[Si](=O)O[Al]=O NLYAJNPCOHFWQQ-UHFFFAOYSA-N 0.000 description 1
- 229910052622 kaolinite Inorganic materials 0.000 description 1
- 210000000867 larynx Anatomy 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 229910052901 montmorillonite Inorganic materials 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000011158 quantitative evaluation Methods 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 229910021646 siderite Inorganic materials 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 229910021647 smectite Inorganic materials 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
- G01N33/24—Earth materials
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/08—Investigating permeability, pore-volume, or surface area of porous materials
- G01N15/088—Investigating volume, surface area, size or distribution of pores; Porosimetry
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
- G06N3/084—Backpropagation, e.g. using gradient descent
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Physics (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Artificial Intelligence (AREA)
- Evolutionary Computation (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Biomedical Technology (AREA)
- Software Systems (AREA)
- Dispersion Chemistry (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Food Science & Technology (AREA)
- Medicinal Chemistry (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本发明属于泥页岩组分评价技术领域，公开了一种利用常规测井曲线评价泥页岩各组分体积的方法，包括：基于抽提后泥页岩的有机碳分析、孔隙度测试及全岩鉴定的实验，结合泥页岩各组分密度，对泥页岩各组分体积进行标定，建立泥页岩组分体积模型；在ΔlogR法评价总有机碳含量的基础上，结合抽提前后有机碳的关系，计算干酪根体积，采用交叉验证的方法优化计算各矿物组分和孔隙体积的BP神经网络模型。本发明在保证泥页岩各组分体积之和为1的前提下，不仅发挥了BP神经网络多输入、多输出的优势，而且解决了泥页岩各组分与测井响应之间复杂的非线性问题。
Description
技术领域
本发明属于泥页岩组分评价技术领域，尤其涉及一种评价泥页岩各组分体积的方法。
背景技术
近年来，页岩油气作为非常规油气领域的重要发展方向之一，因其资源量巨大而受到广泛关注。页岩油气是指已经生成但滞留于富有机质泥页岩的微-纳米级孔隙中的烃类，富有机质泥页岩既是烃源岩，又是储集岩，其源储一体的性质决定了页岩油气能否被有效开采主要取决于烃类的富集程度和渗流能力。泥页岩中含烃量的大小主要与总有机碳含量(TOC)和储层孔隙度有关，而烃类的渗流能力主要受储集空间(孔、喉分布及连通性)的影响。泥页岩总有机碳是评价岩石生烃潜力的重要参数，在岩石中所占的体积很小，以干酪根和残留烃的形式存在。此外，相对于常规油气藏，页岩油气储集空间较为致密，且渗透率极低，一般无自然流动能力，需要大规模水力压裂才能形成工业产能，其可压裂性受矿物组成的影响。因此，在页岩油气勘探开发前期，对泥页岩各组分(干酪根、孔隙、矿物)体积的评价显得格外重要。
泥页岩干酪根、孔隙以及各矿物体积分别可以通过实验室内Leco/Rock-Eval分析、覆压孔渗以及全岩XRD分析等实验手段直接或间接得到，且精度较高，但是受取样成本和实验经费限制，难以连续定量评价岩石各组分的含量。国外多数石油公司逐渐利用伽马能谱测井(NGS)、核磁共振测井(NMR)、元素俘获测井(ECS)等特殊的测井技术对泥页岩地层干酪根体积、各矿物含量以及孔隙度等进行解释，并取得了较大的成功。但是因价格昂贵等问题，上述特殊测井技术在国内并未得到广泛应用，有这些特殊测井资料的井位相对较少，因此，亟需一种利用常规测井曲线对泥页岩干酪根、各矿物组分及孔隙体积预测的方法。
现有技术一：廖东良(2014)利用ECS测井解释的干酪根和各矿物含量作为已知条件，基于线性全体积模型，建立页岩各组分(干酪根、骨架矿物、油气)的多个测井响应方程并对其进行求解，以此计算页岩地层各矿物和干酪根含量(专利申请号：201410318700.5和201410319217.9))。
现有技术一的缺点：
受氧化物闭合模型的精度及元素质量转化为矿物质量的系数的影响，ECS测井获得的地层矿物含量和岩心全岩分析(XRD)实测值之间存在着一定的差异，而其建立的测井响应方程求取的为ECS测井解释的结果，这无疑存在着误差；此外，该体积模型采用的是线性全体积模型，对于非均质性较强的泥页岩地层，各组分的分布形式存在着较大差异使得对测井的响应不是简单的线性叠加。
现有技术二的技术方案：
刘欢(2016)在获得待测样品的标准俘获伽马能谱的基础上，计算目的页岩地层的矿物质量含量，并利用线性体积模型构建多个测井响应方程，以此确定页岩的矿物体积含量。
现有技术二的缺点：
该方法是建立在准确的获取俘获伽马能谱，而这一技术采用的是非常规测井范畴，对于没有非常规测井资料的工区难以使用；此外，构建的测井响应方程是线性模型。
现有技术三的技术方案：
钟广法等(测井资料反演南海北部陆坡渐新统的矿物组分，2006)根据实测的岩芯分析资料，将地层简化为陆源碎屑、碳酸盐岩、粘土矿物和孔隙四个组分，并选用与孔隙度关系密切的测井曲线，建立测井响应方程组，根据实测分析值与测井之间的关系，反算地层各组分的测井响应参数，以此计算地层各组分含量。
现有技术三的缺点：
该方法针对的对象是常规砂岩储层，而对于泥页岩储层来说，除了含有上述四个组分外，其有机质和黄铁矿含量较为发育，且非均质性较强，各组分分布形式存在较大差异，反算的地层各组分的测井响应参数出现负值等不合理现象，因此，该方法已不能适用于泥页岩储层各组分预测。
现有技术四的技术方案：
张晋言等(利用测井资料评价泥页岩油气“五性”指标，2012)采用单矿物组分和测井曲线响应之间的关系，分别建立了泥页岩中的泥质含量、砂质含量、灰质含量与各测井曲线的关系模型。
现有技术四的缺点：
该方法采用测井曲线逐一拟合泥页岩单矿物组分，但最终求得的各组分体积之和不等于1；此外，该方法所采用的测井曲线计算泥页岩单矿物组分多数为经验模型且地域性较强，不宜推广。
现有技术五的技术方案：
Jacobi等(2008)利用密度测井和核磁共振测井确定的孔隙度的差值计算干酪根体积。
现有技术五的缺点
该方法利用到核磁共振测井技术，属于非常规测井范畴，对于没有核磁共振测井资料的井位难以推广应用。
现有技术六的技术方案：
Lewis等(2004)根据总有机碳含量、岩石密度、干酪根密度与干酪根体积之间的关系，实现了在评价总有机碳含量基础上，结合密度测井曲线计算干酪根体积。
现有技术六的缺点：
该方案在计算干酪根体积时采用的是总有机碳含量，而总有机碳并非全部来自干酪根，还有部分来自岩石中残留的油气的贡献，因此，其评价结果偏高。
综上所述，现有技术存在的问题是：第一，泥页岩非均质性较强，且各组分分布形式复杂，其测井响应不是简单的线性叠加，全体积线性模型已不再适用；第二，根据总有机碳含量评价干酪根体积时，未考虑到残留油气中有机碳对总有机碳的贡献，导致评价的干酪根体积偏高。
发明内容
针对现有技术存在的问题，本发明提供了一种评价泥页岩各组分体积的方法。
本发明是这样实现的，一种评价泥页岩各组分体积的方法，所述评价泥页岩各组分体积的方法包括：
基于抽提后泥页岩的有机碳分析、孔隙度测试及全岩鉴定的实验，结合泥页岩各组分密度，对泥页岩各组分体积进行标定，建立泥页岩组分体积模型；
在ΔlogR法评价总有机碳含量的基础上，结合抽提前后有机碳的关系，计算干酪根体积，并结合测井曲线一并作为BP神经网络模型的输入数据，各矿物组分和孔隙体积作为期望输出数据；采用交叉验证的方法优化各矿物组分和孔隙体积的BP神经网络预测模型。进一步，所述计算干酪根体积的方法包括：
采用改进的ΔlogR法，基于计算TOC和实测TOC之间的误差最小化，自动选取基线，优化叠合系数，把TOC含量背景值作为待定系数，其TOC计算模型为：
TOC＝A×ΔlogR+B (1)
TOC为泥页岩总有机碳含量；ΔlogR为在算术坐标下，电阻率曲线和声波时差测井曲线在细粒非生油岩处叠合后，两条测井曲线在对数电阻率坐标上的间距；A和B为模型计算系数；
岩石内干酪根中有机碳含量TOCk，通过氯仿抽提后岩样的有机碳分析得到，且该值与岩石总有机碳含量TOC一般呈现线性关系，由TOC计算得到，即：
TOCK＝C×TOC+D (2)
C和D为计算系数，通过泥页岩抽提前后有机碳分析实验结果拟合得到；
根据干酪根中有机碳含量TOCk，岩石密度ρb和干酪根密度ρK计算得到干酪根体积Vk：
式中，Kvr为干酪根与有机碳之间的转化系数，一般取值为1.2；
因此，联立公式(1)(2)(3)得干酪根体积的测井评价模型为：
进一步，各矿物组分和孔隙体积的BP神经网络模型的建立方法，包括：数据的准备和网络模型参数的优化；
所述数据的准备包括期望输出数据的准备、输入数据的优选及数据预处理；
所述网络模型参数优化的对象为隐含层节点数S，节点传递函数。
进一步，所述期望输出数据的准备包括：
建立泥页岩组分模型：根据各矿物的化学成分、密度属性的差异，把泥页岩的矿物类型划分为4类：粘土类、硅酸盐类、碳酸盐类和黄铁矿；基于泥页岩矿物类型划分，并结合干酪根和孔隙，把泥页岩组成划分为6个组分，即粘土类矿物、硅酸盐类矿物、碳酸盐类矿物、黄铁矿、干酪根和孔隙。
进一步，所述期望输出数据的准备还包括进行各组分体积的标定：根据干酪根及各矿物组分的密度，结合泥页岩组分模型对泥页岩各组分的体积进行标定；其中，干酪根体积VK根据公式(3)计算得到；
孔隙体积VP为岩心分析测试的总孔隙度φ，即：
实验室内根据全岩分析XRD得到的矿物含量Mi(XRD)均为质量百分数，因此，结合各矿物的密度ρi即求出各矿物的体积Vi，其计算公式为：
基于泥页岩组分全体积模型，所有矿物体积、干酪根体积与孔隙体积之和为1；但值得注意的是，实验室内XRD分析检测不到干酪根的含量，即XRD分析得到的各矿物质量比例为各矿物质量与总矿物质量之比，而非各矿物质量与岩石质量之比，对公式(6)计算的各矿物的体积Vi进行校正，其校正公式为：
式中，Vmi为校正后的各矿物体积；
根据公式(5)～(7)分别求得实测的泥页岩各矿物组分和孔隙体积，并把该部分作为模型的期望输出数据。
进一步，输入数据的优选方法包括：
对于建立的泥页岩组分模型，选用与各矿物和孔隙体积相关性较高的测井曲线作为模型输入变量；各组分体积与测井曲线的相关性的判别方法通过公式(8)实现，优选皮尔森相关系数在0.01水平上显著相关的测井曲线；此外，泥页岩矿物组分及孔隙体积受到干酪根体积的约束，把优选的测井曲线和干酪根体积一并作为模型的输入数据；
式中，r为皮尔森相关系数；xi和yi均为变量；N为个数。
进一步，数据的预处理包括：
根据输入的测井曲线的量纲不同及网络收敛速度，把数据归一化到-1和1之间，归一化计算公式为：
式中，x为输入变量；z为x经过归一化后的变量；xmax和xmin分别为输入变量的最大值和最小值。
进一步，网络模型参数的优化包括：
BP神经网络模型采用单隐含层神经网络；参与BP神经网络模型优化的数据包括期望输出数据和输入数据，把参与BP神经网络模型优化的数据随机分成训练样本、验证样本、检测样本三个部分，训练样本和验证样本参与网络训练，检测样本不参与网络训练，仅用来检测网络模型的预测性能；
采用训练样本和验证样本交叉验证的方法，对BP网络模型参数进行优化，并根据检测样本对优化的网络进行检测，基于训练样本、验证样本、检测样本的输出值与期望值之间误差之和的最小化，自动调整隐含层节点数S、节点传递函数TF，直到模型精度满足要求为止。
本发明的优点及积极效果为：
本发明在泥页岩样品的有机碳分析、孔隙度测试以及矿物含量检测结果，建立泥页岩组分模型，并提出一种利用常规测井曲线预测泥页岩组分(矿物、干酪根和孔隙)的评价方法，该方法结合BP神经网络和ΔlogR技术，BP神经网络输入参数包括干酪根体积和测井曲线，矿物组分(粘土类矿物、硅酸盐类矿物、碳酸盐类矿物、黄铁矿)和孔隙体积为输出结果。干酪根体积是在改进的ΔlogR模型评价TOC基础上，结合泥页岩抽提前后有机碳实验结果，通过干酪根体积和有机碳转换公式获得；各矿物组分和孔隙体积是根据优化后的BP神经网络求得。与前人采用的线性全体积模型和单组分多测井拟合法相比，该方法不仅保证了估算的各组分体积之和为1，同时解决了泥页岩各组分体积与测井响应之间复杂的非线性问题，此外，该方法在计算干酪根体积时同时去除了可溶有机质中残留碳的影响。
本发明以大民屯凹陷古近系沙河街组泥页岩为例，根据提出的方法分别对泥页岩中干酪根、各矿物组分及孔隙的体积进行了应用，并分别与实测值进行了对比。
如图3所示，本发明计算的泥页岩各组分体积与实测值均分布在y＝x附近，其中，干酪根、粘土类矿物、硅酸盐类矿物、碳酸盐类矿物及孔隙度的计算值和实测值相关系数(R2)均在75％以上，且各组分体积之和为100％，效果较好。此外，检测数据点的在y＝x两侧的均匀分布保证了该模型的预测性能。但是对黄铁矿的预测效果不是很好，可能与其含量较低有关。与实测值相比，本发明的预测结果表现出较好的匹配效果，精度较高，能够适用于泥页岩各组分体积的预测。
泥页岩组成(矿物、干酪根和孔隙)评价对于页岩油气的富集和压裂研究有重要意义，针对国内元素俘获等特殊测井资料有限的情况，能否利用常规测井资料预测泥页岩各组分体积关系到下一步对页岩油甜点的预测，因此本发明对页岩油勘探开发具有重要意义。
附图说明
图1是本发明实施例提供的评价泥页岩各组分体积的方法流程图。
图2是本发明实施例提供的泥页岩组分模型示意图。
图3是本发明实施例提供的评价泥页岩各组分体积的效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
本发明中定义为：TOC：总有机碳含量；TOCk：岩石内干酪根中有机碳含量；ECS：元素俘获测井；XRD：全岩分析。
下面结合附图对本发明应用原理作详细描述。
如图1所示，本发明实施例提供的评价泥页岩各组分体积的方法，包括：
基于抽提后泥页岩的有机碳分析、孔隙度测试及全岩鉴定的实验，结合泥页岩各组分密度，对泥页岩各组分体积进行标定，建立泥页岩组分体积模型；
在ΔlogR法评价总有机碳含量的基础上，结合抽提前后有机碳的关系，计算干酪根体积，并结合测井曲线作为BP神经网络模型的输入数据，各矿物组分和孔隙体积作为期望输出数据；采用交叉验证的方法优化各矿物组分和孔隙体积的BP神经网络预测模型。
计算干酪根体积的方法包括：
ΔlogR法计算泥页岩总有机碳含量TOC，主要是将算术坐标下的声波时差(AC)曲线和对数坐标下电阻率(RT)曲线在细粒非生油岩处重合，并确定为基线位置，两条曲线在对数电阻率坐标上的间距即为ΔlogR，即：
式中，R和Δt分别是电阻率和声波时差值；Rbseline和Δtbseline分别为细粒非生油岩段的电阻率和声波时差基线值；K为叠合系数。
ΔlogR一般与岩石总有机碳含量呈正相关性，但是受人为选取基线、叠合系数采用定值0.02、工区的TOC含量背景值不确定等影响，通常该方法测井计算的TOC与实测TOC之间的相关性很难达到预期目标。为此，采用改进的ΔlogR法，基于计算TOC和实测TOC之间的误差最小化，自动选取基线，优化叠合系数，把TOC含量背景值作为待定系数，其TOC计算模型为：
TOC＝A×ΔlogR+B (2)
TOC为泥页岩总有机碳含量；ΔlogR为在算术坐标下，电阻率曲线和声波时差测井曲线在细粒非生油岩处叠合后，两条测井曲线在对数电阻率坐标上的间距；A和B为模型计算系数；
岩石内干酪根中有机碳含量(TOCk)，可通过氯仿抽提后岩样的有机碳分析得到，且该值与岩石总有机碳含量TOC一般呈现较好的线性关系，可由TOC计算得到，即：
TOCK＝C×TOC+D (3)
根据干酪根中有机碳含量(TOCk)，岩石密度ρb和干酪根密度ρK可计算得到干酪根体积(Vk)：
式中，Kvr为干酪根与有机碳之间的转化系数，一般取值为1.2。
因此，联立公式(2)(3)(4)可得干酪根体积的测井评价模型为：
各矿物组分及孔隙体积评价包括：
关于各矿物组分和孔隙体积的BP神经网络模型的建立，其步骤总体包括两个部分：数据的准备和网络模型参数的优化。其中，数据的准备包括期望输出数据的准备(各矿物组分和孔隙体积)、输入数据(测井曲线)的优选及数据预处理等，网络模型参数优化的对象主要为隐含层节点数S，节点传递函数等。
期望输出数据的准备：
(1)泥页岩组分模型
对于泥页岩来说，其一般发育伊利石、绿泥石、高岭石、蒙脱石、伊蒙混层、石英、长石、方解石、白云石、菱铁矿以及黄铁矿等无机矿物，根据有限的测井资料难以评价出上述所有矿物，因此，需要简化矿物类型。根据各矿物的化学成分、密度属性的差异，把泥页岩的矿物类型划分为4类：粘土类、硅酸盐类、碳酸盐类和黄铁矿(表1)。
如图2所示，基于上述泥页岩矿物类型划分方案，并结合到干酪根和孔隙，把泥页岩组成细为6个组分模型，即粘土类矿物、硅酸盐类矿物、碳酸盐类矿物、黄铁矿、干酪根和孔隙。
表1 泥页岩组成细分方案
(2)各组分体积的标定
考虑到干酪根、各矿物组分等的密度存在显著的差异，结合组分模型对泥页岩各组分的体积进行标定。其中，干酪根体积VK可根据公式(4)计算得到。
孔隙体积VP为岩心分析测试的总孔隙度φ，即：
实验室内根据全岩分析(XRD)得到的矿物含量Mi(XRD)均为质量百分数，因此，结合各矿物的密度ρi(各矿物的密度详见表1)即可求出各矿物的体积Vi，其计算公式为：
基于泥页岩组分全体积模型，所有矿物体积、干酪根体积与孔隙体积之和为1。但值得注意的是，实验室内XRD分析检测不到干酪根的含量，即XRD分析得到的各矿物质量比例为各矿物质量与总矿物质量之比，而非各矿物质量与岩石质量之比，因此需要对公式(7)计算的各矿物的体积Vi进行校正，其校正公式为：
式中，Vmi为校正后的各矿物体积。
因此，根据泥页岩样品公式(6)～(8)可分别求得实测的泥页岩各矿物组分和孔隙体积，并把该部分作为模型的期望输出数据。
输入数据的优选：
对于上述建立的泥页岩组分模型，选用与各矿物和孔隙体积相关性较高的测井曲线作为模型输入变量，其预测效果愈好。各组分体积与测井曲线的相关性的判别方法见公式(9)，优选皮尔森相关系数在0.01水平上显著相关的测井曲线。此外，泥页岩矿物组分及孔隙体积受到干酪根体积的约束，因此，把优选的测井曲线和干酪根体积一并作为模型的输入数据。
式中，r为皮尔森相关系数；xi和yi均为变量；N为个数。
数据的预处理：
考虑到输入的测井曲线的量纲不同及网络收敛速度，把数据归一化到-1和1之间，归一化计算公式为：
式中，x为输入变量；z为x经过归一化后的变量；xmax和xmin分别为输入变量的最大值和最小值。
模型参数的优化：
单隐含层神经网络可以有效的逼近任意连续函数，为了加快程序运行速度，BP神经网络模型采用单隐含层神经网络。
参与BP神经网络模型优化的数据包括期望输出数据和输入数据，考虑到训练样本在神经网络建立过程中占有重要的地位，样本是否具有代表性，直接影响着网络模型的效果以及预测性能。因此，把参与BP神经网络模型优化的数据随机分成训练样本、验证样本、检测样本三个部分，训练样本和验证样本参与网络训练，检测样本不参与网络训练，仅用来检测网络模型的预测性能。
BP神经网络模型参数优化的主要对象是隐含层节点数S及节点传递函数TF。采用训练样本和验证样本交叉验证(cross-validation)的方法，对BP网络模型参数进行优化，并根据检测样本对优化的网络进行检测，基于训练样本、验证样本、检测样本的输出值与期望值之间误差之和的最小化，自动调整隐含层节点数S、节点传递函数TF等，直到模型精度满足要求为止。
下面结合积极效果对本发明应用原理作进一步描述。
本发明以大民屯凹陷古近系沙河街组泥页岩为例，根据提出的方法分别对泥页岩中干酪根、各矿物组分及孔隙的体积进行了应用，并分别与实测值进行了对比。
如图3所示，本发明计算的泥页岩各组分体积与实测值均分布在y＝x附近，其中，干酪根、粘土类矿物、硅酸盐类矿物、碳酸盐类矿物及孔隙度的计算值和实测值相关系数(R2)均在75％以上，且各组分体积之和为100％，效果较好。此外，检测数据点的在y＝x两侧的均匀分布保证了该模型的预测性能。但是对黄铁矿的预测效果不是很好，可能与其含量较低有关。与实测值相比，本发明的预测结果表现出较好的匹配效果，精度较高，能够适用于泥页岩各组分体积的预测。
泥页岩组成(矿物、干酪根和孔隙)评价对于页岩油气的富集和压裂研究有重要意义，针对国内元素俘获等特殊测井资料有限的情况，能否利用常规测井资料预测泥页岩各组分体积关系到下一步对页岩油甜点的预测，因此本发明对页岩油勘探开发具有重要意义。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。
Claims (4)
1.一种评价泥页岩各组分体积的方法，其特征在于，所述评价泥页岩各组分体积的方法包括：
基于抽提后泥页岩的有机碳分析、孔隙度测试及全岩鉴定的实验，结合泥页岩各组分密度，对泥页岩各组分体积进行标定，建立泥页岩组分体积模型；
在ΔlogR法评价总有机碳含量的基础上，结合抽提前后有机碳的关系，计算干酪根体积，并把干酪根体积和测井曲线一并作为BP神经网络模型的输入数据，各矿物组分和孔隙体积作为期望输出数据；采用交叉验证的方法优化各矿物组分和孔隙体积的BP神经网络预测模型；
所述计算干酪根体积的方法包括：
采用改进的ΔlogR法，基于计算TOC和实测TOC之间的误差最小化，自动选取基线，优化叠合系数，把TOC含量背景值作为待定系数，其TOC计算模型为：
TOC＝A×ΔlogR+B；
TOC为泥页岩总有机碳含量；ΔlogR为在算术坐标下，电阻率曲线和声波时差测井曲线在细粒非生油岩处叠合后，两条测井曲线在对数电阻率坐标上的间距；A和B为模型计算系数；
岩石内干酪根中有机碳含量TOCk，通过氯仿抽提后岩样的有机碳分析得到，且该值与岩石总有机碳含量TOC一般呈现线性关系，由TOC计算得到，即：
TOCK＝C×TOC+D；
C和D为计算系数，通过泥页岩抽提前后有机碳分析实验结果拟合得到；
根据干酪根中有机碳含量TOCk，岩石密度ρb和干酪根密度ρK计算得到干酪根体积Vk：
式中，Kvr为干酪根与有机碳之间的转化系数，取值为1.2；
因此，联立TOC计算模型公式、岩石内干酪根中有机碳含量TOCk公式、干酪根体积Vk公式得干酪根体积的测井评价模型为：
各矿物组分和孔隙体积的BP神经网络模型的建立方法，包括：数据的准备和网络模型参数的优化；
所述数据的准备包括期望输出数据的准备、输入数据的优选及数据预处理；
所述网络模型参数优化的对象为隐含层节点数S，节点传递函数；
所述期望输出数据的准备包括：
建立泥页岩组分模型：
根据各矿物的化学成分、密度属性的差异，把泥页岩的矿物类型划分为4类：粘土类、硅酸盐类、碳酸盐类和黄铁矿；基于泥页岩矿物类型划分，并结合干酪根和孔隙，把泥页岩组成划分为6个组分，即粘土类矿物、硅酸盐类矿物、碳酸盐类矿物、黄铁矿、干酪根和孔隙；
所述期望输出数据的准备还包括进行各组分体积的标定：根据干酪根及各矿物组分的密度，结合泥页岩组分模型对泥页岩各组分的体积进行标定；其中，干酪根体积VK根据公式计算得到；
孔隙体积VP为岩心分析测试的总孔隙度φ，即：
实验室内根据全岩分析XRD得到的矿物含量Mi(XRD)均为质量百分数，因此，结合各矿物的密度ρi即求出各矿物的体积Vi，其计算公式为：
基于泥页岩组分全体积模型，所有矿物体积、干酪根体积与孔隙体积之和为1；但值得注意的是，实验室内XRD分析检测不到干酪根的含量，即XRD分析得到的各矿物质量比例为各矿物质量与总矿物质量之比，而非各矿物质量与岩石质量之比，对公式计算的各矿物的体积Vi进行校正，其校正公式为：
式中，Vmi为校正后的各矿物体积；
根据公式分别求得实测的泥页岩各矿物组分和孔隙体积，并把该部分作为模型的期望输出数据。
2.如权利要求1所述的评价泥页岩各组分体积的方法，其特征在于，
输入数据的优选方法包括：
对于建立的泥页岩组分模型，选用与各矿物和孔隙体积相关性较高的测井曲线作为模型输入变量；各组分体积与测井曲线的相关性的判别方法通过公式实现，式中，r为皮尔森相关系数；xi和yi均为变量；N为个数；
优选皮尔森相关系数在0.01水平上显著相关的测井曲线；此外，由于泥页岩矿物组分及孔隙体积受到干酪根体积的约束，因此，优选的测井曲线和干酪根体积一并作为模型的输入数据。
3.如权利要求1所述的评价泥页岩各组分体积的方法，其特征在于，数据的预处理包括：
根据输入的测井曲线的量纲不同及网络收敛速度，把数据归一化到-1和1之间，归一化计算公式为：
式中，x为输入变量；z为x经过归一化后的变量；xmax和xmin分别为输入变量的最大值和最小值。
4.如权利要求1所述的评价泥页岩各组分体积的方法，其特征在于，网络模型参数的优化包括：
BP神经网络模型采用单隐含层神经网络；参与BP神经网络模型优化的数据包括期望输出数据和输入数据，把参与BP神经网络模型优化的数据随机分成训练样本、验证样本、检测样本三个部分，训练样本和验证样本参与网络训练，检测样本不参与网络训练，仅用来检测网络模型的预测性能；
采用训练样本和验证样本交叉验证的方法，对BP网络模型参数进行优化，并根据检测样本对优化的网络进行检测，基于训练样本、验证样本、检测样本的输出值与期望值之间误差之和的最小化，自动调整隐含层节点数S、节点传递函数TF，直到模型精度满足要求为止。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710124289.1A CN106950347B (zh) | 2017-03-03 | 2017-03-03 | 一种评价泥页岩各组分体积的方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710124289.1A CN106950347B (zh) | 2017-03-03 | 2017-03-03 | 一种评价泥页岩各组分体积的方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106950347A CN106950347A (zh) | 2017-07-14 |
| CN106950347B true CN106950347B (zh) | 2019-03-05 |
Family
ID=59466612
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710124289.1A Expired - Fee Related CN106950347B (zh) | 2017-03-03 | 2017-03-03 | 一种评价泥页岩各组分体积的方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106950347B (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110211144A (zh) * | 2019-06-10 | 2019-09-06 | 中国石油大学(华东) | 一种泥页岩扫描电镜图片孔隙提取方法及系统 |
Families Citing this family (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107560994B (zh) * | 2017-09-05 | 2019-08-02 | 中国石油大学(华东) | 一种评价泥页岩有机质、粘土和其它矿物中孔径分布方法 |
| CN109989746A (zh) * | 2017-12-29 | 2019-07-09 | 中国石油天然气股份有限公司 | 碳酸盐岩储层评价的方法和装置 |
| CN108458960B (zh) * | 2018-03-27 | 2019-10-29 | 中国石油大学(华东) | 富有机质泥页岩含氢组分、孔隙度及孔径的评价方法 |
| CN109188557B (zh) * | 2018-09-29 | 2020-07-10 | 中国石油天然气股份有限公司 | 一种应用生烃增压模型计算排烃效率的方法 |
| CN109612794B (zh) * | 2018-12-06 | 2021-05-11 | 同济大学 | 一种烃源岩中不同赋存态有机质的分离与定量方法 |
| CN112746835B (zh) * | 2019-10-31 | 2024-10-01 | 中国石油化工集团有限公司 | 一种优化的深层页岩气地质甜点测井综合评价方法 |
| CN111022040B (zh) * | 2019-11-29 | 2022-07-12 | 中国石油集团川庆钻探工程有限公司 | 一种页岩气储层岩心有机化合物体积的计算方法 |
| CN111058840A (zh) * | 2019-12-18 | 2020-04-24 | 延安大学 | 一种基于高阶神经网络的有机碳含量(toc)评价方法 |
| CN111028095A (zh) * | 2019-12-19 | 2020-04-17 | 中国地质大学（武汉） | 一种基于测井曲线的定量识别页岩岩相的方法 |
| CN111749677A (zh) * | 2020-06-30 | 2020-10-09 | 宋立才 | 测井信息处理方法和装置、电子设备及存储介质 |
| CN111912961B (zh) * | 2020-08-24 | 2022-12-09 | 东北石油大学 | 一种泥页岩中干酪根吸附油定量评价方法及系统 |
| CN111946336A (zh) * | 2020-09-06 | 2020-11-17 | 中国石油天然气集团有限公司 | 基于矿物分析数据获取页岩油储层核磁孔隙度的方法 |
| CN112489736A (zh) * | 2020-12-09 | 2021-03-12 | 中国石油大学(北京) | 一种矿物含量分析方法、装置、设备及存储介质 |
| CN112987125B (zh) * | 2021-02-22 | 2021-12-17 | 中国地质大学（北京） | 一种基于测井数据的页岩脆性指数预测方法 |
| CN114035227B (zh) * | 2021-11-11 | 2023-07-07 | 中国海洋石油集团有限公司 | 基于随钻xrd全岩录井的变质岩潜山储层孔隙度预测方法 |
| CN115684554B (zh) * | 2022-12-06 | 2023-08-11 | 东北石油大学三亚海洋油气研究院 | 富有机质泥页岩的四端元岩相划分方法、系统及介质 |
| CN118518538A (zh) * | 2024-05-27 | 2024-08-20 | 招商局检测认证(重庆)有限公司 | 一种通过密度体积对abc磷酸二氢铵干粉灭火剂进行定性的判定方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
| CN103822866A (zh) * | 2014-03-21 | 2014-05-28 | 中国石油大学(华东) | 一种评价目的层段泥页岩孔隙度方法 |
| CN105317429A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种确定页岩地层干酪根含量的方法 |
| CN105927218A (zh) * | 2016-05-20 | 2016-09-07 | 中国石油大学(北京) | 一种陆相泥页岩储层含气量预测方法及装置 |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
-
2017
- 2017-03-03 CN CN201710124289.1A patent/CN106950347B/zh not_active Expired - Fee Related
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103792581A (zh) * | 2013-12-31 | 2014-05-14 | 刘致水 | 基于岩石物理模型的页岩气储层地应力测井预测方法 |
| CN103822866A (zh) * | 2014-03-21 | 2014-05-28 | 中国石油大学(华东) | 一种评价目的层段泥页岩孔隙度方法 |
| CN105317429A (zh) * | 2014-07-04 | 2016-02-10 | 中国石油化工股份有限公司 | 一种确定页岩地层干酪根含量的方法 |
| CN105986815A (zh) * | 2015-02-27 | 2016-10-05 | 中国石油化工股份有限公司 | 一种用于识别页岩地层地质甜点的方法 |
| CN105927218A (zh) * | 2016-05-20 | 2016-09-07 | 中国石油大学(北京) | 一种陆相泥页岩储层含气量预测方法及装置 |
Non-Patent Citations (2)
| Title |
|---|
| 页岩气储层Ramirez体积物理模型的改进;赵培强 等;《石油学报》;20140531;第35卷(第3期);第480-485页 |
| 页岩气储层有机碳含量与测井参数的关系及预测模型;孟召平 等;《煤炭学报》;20150228;第40卷(第2期);第247-253页 |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110211144A (zh) * | 2019-06-10 | 2019-09-06 | 中国石油大学(华东) | 一种泥页岩扫描电镜图片孔隙提取方法及系统 |
| CN110211144B (zh) * | 2019-06-10 | 2021-06-22 | 中国石油大学(华东) | 一种泥页岩扫描电镜图片孔隙提取方法及系统 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106950347A (zh) | 2017-07-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN106950347B (zh) | 一种评价泥页岩各组分体积的方法 | |
| Wood | Water distribution in the Montney tight gas play of the western Canadian sedimentary basin: significance for resource evaluation | |
| CN102175832B (zh) | 一种确定典型储层最佳饱和度计算模型的方法 | |
| CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
| CN104806232B (zh) | 一种确定碳酸盐岩储层孔隙度下限的方法 | |
| CN105842416B (zh) | 一种基于兰格缪尔模型的页岩吸附气含量动态计算方法 | |
| CN107795320B (zh) | 一种水平井碳酸盐岩储层参数的计算方法 | |
| CN110554441A (zh) | 一种基于钙质含量校正的致密砂岩含气饱和度计算方法 | |
| CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
| CN107092036A (zh) | 一种基于储层真电阻率反演的火山岩地层流体识别方法及系统 | |
| Feng et al. | Accurate determination of water saturation in tight sandstone gas reservoirs based on optimized Gaussian process regression | |
| Coats et al. | Modeling conformance as dispersion | |
| CN110688781B (zh) | 一种低渗非均质气藏储层测井解释方法 | |
| CN105350959A (zh) | 由测井岩性密度确定页岩气储层含气饱和度的方法 | |
| CN111625750A (zh) | 一种针对致密砂岩储层气层与气水同层的识别方法 | |
| CN109667576A (zh) | 一种高矿化度成因低阻油层测井识别方法 | |
| CN115793094B (zh) | 一种曲线叠合重构识别复杂页岩层系岩性的方法及应用 | |
| CN108805158A (zh) | 一种致密油储层成岩相划分方法 | |
| CN112784404A (zh) | 一种基于常规测井资料的砂砾岩束缚水饱和度计算方法 | |
| Li et al. | Evaluation of irreducible water saturation by electrical imaging logging based on capillary pressure approximation theory | |
| CN109655394B (zh) | 一种孔喉结构参数约束下的核磁共振t2谱渗透率计算方法 | |
| CN109343121A (zh) | 一种确定热接触变质岩储层中岩浆热液规模的方法 | |
| Coats et al. | Modeling conformance as dispersion | |
| Tillero | Stepping forward: An automated rock type index and a new predictive capillary pressure function for better estimation of permeability and water saturation. Case study, Urdaneta-01 heavy oil reservoir | |
| Olson et al. | Micropaleontologic proxies for sea-level change and stratigraphic discontinuities |
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 | ||
| CP02 | Change in the address of a patent holder | ||
| CP02 | Change in the address of a patent holder |
Address after: 266580 Qingdao Changjiang Road, Huangdao District, Shandong, No. 66 Patentee after: China University of Petroleum (East China) Address before: 266580 No. two, 271 North Road, Shandong, Dongying Patentee before: China University of Petroleum (East China) |
|
| CF01 | Termination of patent right due to non-payment of annual fee | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190305 |