CN112649256A - 一种基于人工合成双矿物岩石求取矿物弹性模量的方法 - Google Patents

一种基于人工合成双矿物岩石求取矿物弹性模量的方法 Download PDF

Info

Publication number
CN112649256A
CN112649256A CN201910960276.7A CN201910960276A CN112649256A CN 112649256 A CN112649256 A CN 112649256A CN 201910960276 A CN201910960276 A CN 201910960276A CN 112649256 A CN112649256 A CN 112649256A
Authority
CN
China
Prior art keywords
mineral
elastic modulus
rock
modulus
artificially synthesized
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.)
Pending
Application number
CN201910960276.7A
Other languages
English (en)
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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201910960276.7A priority Critical patent/CN112649256A/zh
Publication of CN112649256A publication Critical patent/CN112649256A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • G01N1/286Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q involving mechanical work, e.g. chopping, disintegrating, compacting, homogenising
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N29/00Investigating 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/04Analysing solids
    • G01N29/07Analysing solids by measuring propagation velocity or propagation time of acoustic waves
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N9/00Investigating density or specific gravity of materials; Analysing materials by determining density or specific gravity
    • G01N9/36Analysing materials by measuring the density or specific gravity, e.g. determining quantity of moisture
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/01Indexing codes associated with the measuring variable
    • G01N2291/011Velocity or travel time
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N2291/00Indexing codes associated with group G01N29/00
    • G01N2291/02Indexing codes associated with the analysed material
    • G01N2291/023Solids
    • G01N2291/0232Glass, ceramics, concrete or stone

Landscapes

  • 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)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Acoustics & Sound (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

本发明提出了一种基于人工合成双矿物岩石求取矿物弹性模量的方法,该方法包括:制备第一矿物和第二矿物,其中第一矿物的弹性模量是已知的,第二矿物的弹性模量是待测的;将第一矿物和第二矿物混合在一起压制成型,制成人工合成双矿物岩石;通过测量人工合成双矿物岩石的纵横波速度和密度,求出第二矿物的弹性模量。采用本发明方法求取干酪根的弹性模量,经过测量计算得到干酪根的弹性模量,将结果与其它两种方法(冷压法和热压法)求得的干酪根模量进行对比,发现本方法求取干酪根的弹性模量与岩石物理手册中的数据最为接近,说明本发明在求取矿物弹性模量是具有较高的准确性。

Description

一种基于人工合成双矿物岩石求取矿物弹性模量的方法
技术领域
本发明涉及地震勘探岩石物理技术,是一种基于实验求取矿物弹性模量的方法,更具体地,涉及一种基于人工合成双矿物岩石求取矿物弹性模量的方法。
背景技术
随着岩石物理技术的发展和应用,其在油气勘探和开发中也起到极大的作用。在之前的应用中,岩石物理主要集中在纵横波速度,随着油气勘探的复杂化,岩石物理建模也不只需要速度,还需要弹性模量的数据,而野外采集的岩石绝大部分都是多矿物组成的岩石,为了精确的建模,必须知道每一种矿物的弹性模量,故而求取每一种矿物的弹性模量显得很重要。
一种较常用的方法是将单矿物加压(有时还需要加温)成型,测试其密度和速度就可以换算出该种矿物的弹性模量。但是,这种方法存在缺陷,其一是许多矿物质能难找到单晶固体;其二是该方法的理论基础是假设矿物本身的弹性模型与单晶体聚合之后的弹性模量一致,但这对于氧化物和硅酸盐等矿物是难以达到要求的。
Brace在1969年提出用间接的方法求取矿物质的压缩系数,他将两种压缩系数相近的矿物混合到一起压制成型,通过测量混合物和压缩系数和已知的另一种矿物的压缩系数便可计算出所求矿物的压缩系数,本发明主要是以该方法为基础,将其改进并引入地震勘探领域用以求取未知矿物的弹性模量。
发明内容
本发明是一种求取矿物弹性模量的方法,应用一种已知其弹性模量的矿物和所求矿物混合到一起并压制成型,通过测量人工岩石的纵横波速度和密度,结合已知矿物的弹性模量,便可以求出未知矿物的弹性模量。
根据本发明的一个发明,提供一种基于人工合成双矿物岩石求取矿物弹性模量的方法,该方法包括:
制备第一矿物和第二矿物,其中第一矿物的弹性模量是已知的,第二矿物的弹性模量是待测的;
将第一矿物和第二矿物混合在一起压制成型,制成人工合成双矿物岩石;
通过测量人工合成双矿物岩石的纵横波速度和密度,求出第二矿物的弹性模量。
进一步地,利用Voigt和Reuss公式,根据人工合成双矿物岩石中各矿物的体积含量和弹性模量来预测人工合成双矿物岩石的等效弹性模量的上下限。
具体地,人工合成双矿物岩石的等效弹性模量的Voigt上限MV是:
MV=f1M1+f2M2
等效弹性模量的Reuss下限MR是:
Figure BDA0002228685820000021
其中f1和f2分别表示第一矿物和第二矿物的体积含量,且f1+f2=1,M1和M2分别表示第一矿物和第二矿物的弹性模量。
进一步地,通过测量人工合成双矿物岩石的纵横波速度和密度,计算获得平均弹性模量MVRH的值。
具体地,平均弹性模量MVRH表示为:
Figure BDA0002228685820000022
可选地,所述第一矿物和第二矿物为粉末状,且所述第一矿物和第二矿物的粒径相同。
优选地,所述第一矿物的弹性模量接近第二矿物的弹性模量。
进一步地,将第一矿物和第二矿物混合并加入胶结剂充分混合在一起压制成型,其中胶结剂含量低于矿物粉末总重量的5%,胶结剂可以是工业所用的液态混合胶。
进一步地,先将两种矿物的混合物在5-10MPa的压力下保压30-60分钟,再将压力提高至15-20MPa保压24-48小时。
本发明针对求取矿物弹性模量提出了一种新的方法,该方法是将两种矿物粉末进行混合,其中一种矿物为弹性模量已知的矿物,另一种矿物是待测矿物,将二者充分混合加入交联剂,并在高压下压制成型,通过测量样品的速度和密度从而得到样品的弹性模量,将该弹性模量认为是用V-R-H平均法得到的模量,将其结果代入Voigt-Reuss公式便可求出未知矿物的弹性模量,通过与其它两种方法的对比,结果显示该方法求得的矿物弹性模量更为准确。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述,本公开的上述以及其它目的、特征和优势将变得更加明显,其中,在本公开示例性实施方式中,相同的参考标号通常代表相同部件。
图1示出了本发明的基于人工合成双矿物岩石求取矿物弹性模量的方法流程图。
图2示出了根据本发明实施例的压制成型的双矿物人工岩石图片。
图3示出了采用冷压法和热压法制作的人工岩石样品图片。
图4示出了根据本发明实施例的双矿物法求取矿物弹性模量的误差图。
图5示出了根据本发明实施例的制作双矿物人工岩石的流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式,然而应该理解,可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反,提供这些实施方式是为了使本公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明的理论基础是Brace与1969年提出的理论:如果聚合物(人工岩石)的密度和矿物的密度接近,那么聚合物的弹性模量与矿物的弹性模量一定具有密切的相关关系。据此本发明设计将两种矿物混合在一起压制成型,其中一种矿物的弹性模量是已知的(以下称为矿物1),另一种矿物是待测的(以下称为矿物2),通过测量成型的人工岩石的纵横波速度和密度便可利用Voigt-Reuss-Hill平均公式求出矿物2的弹性模量。
如图1所示,提供一种基于人工合成双矿物岩石求取矿物弹性模量的方法,该方法包括:
制备第一矿物和第二矿物,其中第一矿物的弹性模量是已知的,第二矿物的弹性模量是待测的;
将第一矿物和第二矿物混合在一起压制成型,制成人工合成双矿物岩石;
通过测量人工合成双矿物岩石的纵横波速度和密度,求出第二矿物的弹性模量。
进一步地,Voigt和Reuss公式是根据岩石中各矿物的体积含量和弹性模量来预测岩石的等效弹性模量的上下限。两种矿物成分的等效弹性模量的Voigt上限MV是:
MV=f1M1+f2M2
等效弹性模量的Reuss下限MR是:
Figure BDA0002228685820000041
其中f1和f2分别表示矿物的体积含量,有f1+f2=1,M1和M2分别表示两种矿物的弹性模量。
Voigt-Reuss-Hill平均是Voigt上限和Reuss下限的算术平均,表示为:
Figure BDA0002228685820000042
经过测量人工合成双矿物岩石的纵横波速度和密度,可以得到MVRH的值,而f1、f2和M1均为已知,联立上述公式便可求得第二矿物的弹性模量。
本发明针对求取矿物弹性模量提出了一种新的方法,该方法是将两种矿物粉末进行混合,其中一种矿物为弹性模量已知的矿物,另一种矿物是待测矿物,将二者充分混合加入交联剂,并在高压下压制成型,通过测量样品的速度和密度从而得到样品的弹性模量,将该弹性模量认为是用V-R-H平均法得到的模量,将其结果代入Voigt-Reuss公式便可求出未知矿物的弹性模量,通过与其它两种方法的对比,结果显示该方法求得的矿物弹性模量更为准确。
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
接下来,参照附图2-5描述本发明的实施例。图5示出了根据本发明实施例的制作双矿物人工岩石的流程图。
首先,如图5所示,制作双矿物人工岩石。在本实施例中,选取碳粉作为第一矿物,碳粉的弹性模量已知,体积模量为2.9GPa,剪切模量为2.6GPa。干酪根是页岩的重要组成部分,干酪根的物性参数对研究页岩的性质具有重要的意义,而对于干酪根的弹性模量却始终无法准确求得。选取干酪根作为第二矿物,利用本发明的方法求取干酪根的弹性模量。
首先,将干酪根粉末和碳粉按照1:1进行充分混合,搅拌均匀获得混合粉末。在混合粉末中加入少量胶结剂搅拌均匀。
清洁模具,在模具内表面涂抹凡士林。
将混合粉末逐层放入模具。
先将混合物在5MPa的压力下保压30分钟,再将压力提高至15MPa保压24小时,制作得到的双矿物人工岩石见附图2。
应用超声测试系统测量双矿物人工岩石样品的纵横波速度,结合其密度便可以用如下公式计算双矿物人工岩石样品的体积模量K和剪切模量模量μ:
Figure BDA0002228685820000051
μ=ρVS 2
其中,Vp为纵波速度,Vs为横波速度,ρ为样品密度。双矿物人工岩石样品的MVRH值是样品的弹性模量(体积模量K或剪切模量模量μ)。
Voigt-Reuss-Hill平均是Voigt上限和Reuss下限的算术平均,表示为:
Figure BDA0002228685820000061
两种矿物成分的等效弹性模量的Voigt上限MV是:
MV=f1M1+f2M2
等效弹性模量的Reuss下限MR是:
Figure BDA0002228685820000062
其中f1和f2分别表示矿物的体积含量,有f1+f2=1,M1和M2分别表示两种矿物的弹性模量,f1、f2和M1均为已知,联立上述公式便可求得干酪根的弹性模量。
对比例
为了对比本发明的优劣,同时还设计了两组对照试验。一种是将干酪根在常温下以15MPa的压力下压制成型并测量其速度和密度,进而换算出其弹性模量;另一种是将将干酪根在120度的高温下以15MPa的压力下压制成型并测量其速度和密度,进而换算出其弹性模量,制作的样品见附图3。图3示出了采用冷压法和热压法制作的样品图片。
最终的计算结果见下表,由表可以看出,利用本发明求得的干酪根模量更为准确。下表给出了三种实验方法求得干酪根弹性模量的对比。
Figure BDA0002228685820000063
Figure BDA0002228685820000071
用本发明用于求取干酪根的弹性模量。首先选取已知弹性模量的碳粉作为第一矿物,将未知弹性模量的干酪根作为第二矿物,并将两者混合压制成型,经过测量计算得到干酪根的弹性模量,将结果与其它两种方法(冷压法和热压法)求得的干酪根模量进行对比,发现本方法求取干酪根的弹性模量与岩石物理手册中的数据最为接近,说明本发明在求取矿物弹性模量是具有较高的准确性。
此外,本发明还分析了该方法的误差,设相对误差为C,则有:
Figure BDA0002228685820000072
相应的误差图见附图4,图4示出了根据本发明实施例的双矿物法求取矿物弹性模量的误差图。
由图4可以看出当第二矿物的模量在2.9GPa附近时,求得结果误差较小,而且越接近误差越小。这就要求本发明在使用时要使得第一矿物的弹性模量越接近第二矿物的弹性模量,求得结果越准确。
以上已经描述了本公开的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择,旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进,或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (10)

1.一种基于人工合成双矿物岩石求取矿物弹性模量的方法,其特征在于,该方法包括:
制备第一矿物和第二矿物,其中第一矿物的弹性模量是已知的,第二矿物的弹性模量是待测的;
将第一矿物和第二矿物混合在一起压制成型,制成人工合成双矿物岩石;
通过测量人工合成双矿物岩石的纵横波速度和密度,求出第二矿物的弹性模量。
2.根据权利要求1所述的方法,其特征在于,利用Voigt和Reuss公式,根据人工合成双矿物岩石中各矿物的体积含量和弹性模量来预测人工合成双矿物岩石的等效弹性模量的上下限。
3.根据权利要求2所述的方法,其特征在于,人工合成双矿物岩石的等效弹性模量的Voigt上限MV是:
MV=f1M1+f2M2
等效弹性模量的Reuss下限MR是:
Figure FDA0002228685810000011
其中f1和f2分别表示第一矿物和第二矿物的体积含量,且f1+f2=1,M1和M2分别表示第一矿物和第二矿物的弹性模量。
4.根据权利要求3所述的方法,其特征在于,通过测量人工合成双矿物岩石的纵横波速度和密度,计算获得平均弹性模量MVRH的值。
5.根据权利要求4所述的方法,其特征在于,平均弹性模量MVRH表示为:
Figure FDA0002228685810000021
6.根据权利要求1所述的方法,其特征在于,所述第一矿物和第二矿物为粉末状,且所述第一矿物和第二矿物的粒径相同。
7.根据权利要求1所述的方法,其特征在于,所述第一矿物的弹性模量接近第二矿物的弹性模量。
8.根据权利要求1所述的方法,其特征在于,将第一矿物和第二矿物混合并加入胶结剂充分混合在一起压制成型。
9.根据权利要求8所述的方法,其特征在于,胶结剂含量低于矿物粉末总量的5%。
10.根据权利要求1所述的方法,其特征在于,先将两种矿物的混合物在5-10MPa的压力下保压30-60分钟,再将压力提高至15-20MPa保压24-48小时。
CN201910960276.7A 2019-10-10 2019-10-10 一种基于人工合成双矿物岩石求取矿物弹性模量的方法 Pending CN112649256A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910960276.7A CN112649256A (zh) 2019-10-10 2019-10-10 一种基于人工合成双矿物岩石求取矿物弹性模量的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910960276.7A CN112649256A (zh) 2019-10-10 2019-10-10 一种基于人工合成双矿物岩石求取矿物弹性模量的方法

Publications (1)

Publication Number Publication Date
CN112649256A true CN112649256A (zh) 2021-04-13

Family

ID=75343492

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910960276.7A Pending CN112649256A (zh) 2019-10-10 2019-10-10 一种基于人工合成双矿物岩石求取矿物弹性模量的方法

Country Status (1)

Country Link
CN (1) CN112649256A (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102757212A (zh) * 2012-07-07 2012-10-31 中国石油天然气股份有限公司 一种砂岩储层低渗透人造岩心的制作方法
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN105988136A (zh) * 2015-01-30 2016-10-05 中石化石油工程技术服务有限公司 一种利用纵横波速度信息进行含气性分析的方法
CN106842364A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 不整合生物灰岩储层地震岩石物理模型建立方法
CN109490953A (zh) * 2017-09-13 2019-03-19 中国石油化工股份有限公司 致密砂岩储层岩石物理建模方法及系统
CN110275202A (zh) * 2019-05-30 2019-09-24 河海大学 一种用于致密油储层脆性预测的方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102757212A (zh) * 2012-07-07 2012-10-31 中国石油天然气股份有限公司 一种砂岩储层低渗透人造岩心的制作方法
WO2016041189A1 (zh) * 2014-09-19 2016-03-24 杨顺伟 一种评价页岩气储层及寻找甜点区的方法
CN105988136A (zh) * 2015-01-30 2016-10-05 中石化石油工程技术服务有限公司 一种利用纵横波速度信息进行含气性分析的方法
CN106842364A (zh) * 2015-12-04 2017-06-13 中国石油化工股份有限公司 不整合生物灰岩储层地震岩石物理模型建立方法
CN109490953A (zh) * 2017-09-13 2019-03-19 中国石油化工股份有限公司 致密砂岩储层岩石物理建模方法及系统
CN110275202A (zh) * 2019-05-30 2019-09-24 河海大学 一种用于致密油储层脆性预测的方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
W.F BRACE,ET AL.: "Isothermal Compressibility of Kyanite, Andalustie, and Sillimanite form Synthetic Aggregates", JOURNAL OF GEOPHYSICAL RESEARCH, vol. 74, no. 8, 15 April 1969 (1969-04-15), pages 2089 - 2098 *
刘建;张林海;陶谦;刘伟;周仕明;王立双;: "多组分油井水泥石弹性模量预测模型", 钻井液与完井液, no. 02 *
石玉梅;谢韬;: "地下固体介质弹性模量反演", 地球物理学报, no. 12 *

Similar Documents

Publication Publication Date Title
Chawre Correlations between ultrasonic pulse wave velocities and rock properties of quartz-mica schist
Luan et al. Laboratory measurements of brittleness anisotropy in synthetic shale with different cementation
CN107203005B (zh) 一种定量化计算裂缝描述参数的方法
CN108399270B (zh) 一种确定页岩地层中各向异性泥质比例的方法
Nguyen et al. A composite sphere assemblage model for porous oolitic rocks
Wollny et al. Coupling of microstructural and macrostructural computational approaches for asphalt pavements under rolling tire load
Tavassoti-Kheiry et al. Application of ultrasonic pulse velocity testing of asphalt concrete mixtures to improve the prediction accuracy of dynamic modulus master curve
Shilko et al. The determining influence of the competition between pore volume change and fluid filtration on the strength of permeable brittle solids
Gong et al. Ultrasonic velocity and mechanical anisotropy of synthetic shale with different types of clay minerals
Robeyst et al. Relating ultrasonic measurements on fresh concrete with mineral additions to the microstructure development simulated by Cemhyd3D
Romanowski et al. Investigations on oilwell cement strength response to ultrasonic measurements in the presence of additives
CN112649256A (zh) 一种基于人工合成双矿物岩石求取矿物弹性模量的方法
Trusler et al. Second and third interaction virial coefficients of the (methane+ propane) system determined from the speed of sound
Zhou et al. A two-dimensional miniature cone penetration test system for centrifuge modelling
Jutten et al. Relationship between cement slurry composition, mechanical properties, and cement-bond-log output
Hernández et al. Application of micromechanics to the characterization of mortar by ultrasound
CN111896202A (zh) 一种用于地震动力边坡模型试验相似材料配比的确定方法
Zhang et al. Effects of Asphalt Source, Asphalt Grade, and Inclusion of Additives on Asphalt Foaming Characteristics
Harari et al. Pore-compressibility study of Arabian carbonate reservoir rocks
Zeng et al. Measurement of G max and Estimation of K of Saturated Clay Using Bender Elements in an Oedometer
CN110954949A (zh) 一种致密砂岩软孔隙度分布反演方法
Moon et al. Acoustic method for determining the static gel strength of slurries
Tu et al. Influence of infiltration velocity on the measurement of the surface energy components of asphalt binders using the Wilhelmy plate method
Kelessidis et al. Comprehensive assessment of additive and class g cement properties affecting rheology, fluid loss, setting time and long term characteristics of elastic cements
Moriwaki et al. Method for determining the coefficient of permeability of clays

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