CN104929627A - 在rpm测井中计算井筒持水率和地层含水饱和度的方法 - Google Patents

在rpm测井中计算井筒持水率和地层含水饱和度的方法 Download PDF

Info

Publication number
CN104929627A
CN104929627A CN201510358153.8A CN201510358153A CN104929627A CN 104929627 A CN104929627 A CN 104929627A CN 201510358153 A CN201510358153 A CN 201510358153A CN 104929627 A CN104929627 A CN 104929627A
Authority
CN
China
Prior art keywords
oil
water
well
carbon
stratum
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.)
Granted
Application number
CN201510358153.8A
Other languages
English (en)
Other versions
CN104929627B (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.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC Research Institute 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 China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201510358153.8A priority Critical patent/CN104929627B/zh
Publication of CN104929627A publication Critical patent/CN104929627A/zh
Application granted granted Critical
Publication of CN104929627B publication Critical patent/CN104929627B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

本发明涉及一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，包括以下步骤：获得井眼和套管时RPM仪器测井中碳氧比随地层孔隙度变化规律图版；井眼和套管时RPM仪器测井中碳氧比计算公式；井眼充水地层充水和井眼充水地层充油时RPM仪器测井中碳氧比；井眼充油地层充油和井眼充油地层充水时RPM仪器测井中碳氧比；计算井眼里水的碳密度和氧密度、井眼里油的碳密度和氧密度；计算井眼、套管和砾石时RPM仪器测井中碳氧比；计算井眼、套管、砾石和泥质时RPM仪器测井中碳氧比；计算井眼、套管、砾石、泥质和钙质时RPM仪器测井中碳氧比；构建远近碳氧比交会图；根据远近碳氧比交会图获取井筒持水率和地层含水饱和度。

Description

在RPM测井中计算井筒持水率和地层含水饱和度的方法
技术领域
本发明涉及石油勘探开发领域，特别是关于一种在RPM测井中计算井筒持水率和地层含水饱和度的方法。
背景技术
RPM(储层性能监测仪)是Atlas公司生产的一种多功能脉冲中子测井仪，该仪器在2001年被引进于国内，已在渤海、上海、湛江等多个海上油气田进行115口井的测量，有效地解决了各油田所面临的诸多问题，如剩余油、油水界面、水淹层、低阻油气层流体性质和水淹方向确定等。目前中国石油大学(北京)的吴文圣利用Monte Carlo(蒙特卡罗)数值模拟的方法对小井径的双源距的RPM仪器测井中碳氧比测井的影响因素进行了校正，其中，影响因素包括井眼尺寸、井眼流体、套管尺寸、水泥环厚度、地层岩性、孔隙度、地层流体和油密度，并给出了能消除这些影响因素的数据处理方法。中国石油大学(华东)的张锋也通过MonteCarlo数值模拟的方法模拟了PNN(脉冲中子-中子)测井的影响因素，这些影响因素包括岩性、泥质含量、水泥环厚度、不同源距和仪器偏心，另外，张锋还模拟了RPM在测井过程中非弹和热中子俘获能谱情况下，影响因素孔隙度、井眼持水率、地层水矿化度、泥质含量和井眼套管对岩性判断的影响。二人都建立了模拟PNN测井在不同影响因素下的相关图版，但没有涉及海上稠油油藏RPM测井RPM仪器测井中碳氧比的解释方法。
目前，国内很直观的计算井筒持水率和地层含水饱和度的图版都不是在海上稠油油藏条件下所建立的理论图版，也没有井下砾石充填的情况，实际的井径也不是大井径，因此需要建立一套适合海上稠油油田RPM测井快速计算井筒持水率和地层含水饱和度的井下试验解释图版。
发明内容
针对上述问题，本发明的目的是提供一种能够准确、快速地在海上稠油油田RPM测井中计算井筒持水率和地层含水饱和度的方法。
为实现上述目的，本发明采取以下技术方案：一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，包括以下步骤：
1)利用蒙特卡罗模拟的方法，得到当仅考虑井眼尺寸和套管尺寸时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况下，纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版；
2)根据井眼尺寸和套管尺寸，并利用斯伦贝谢的郝佐格公式得到RPM仪器测井中碳氧比计算公式为：
式中，C/O为RPM仪器测井中碳氧比；A为碳和氧与中子反应平均截面的比值，是常数；a为每立方厘米油中碳原子的数目；b为每立方厘米砂岩中碳原子的数目；c为每立方厘米水中氧原子的数目；d为每立方厘米砂岩中氧原子的数目；为地层孔隙度；S0为地层含油饱和度，BC为井眼里碳密度的贡献；BO为井眼里氧密度的贡献；
3)当井眼充水地层充水时，地层中的含油饱和度S0＝0，RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充水地层充油时，地层中的含油饱和度S0＝1，RPM仪器测井中碳氧比计算公式(1)变形为：
式中，BC水为井眼里水的碳密度，BO水为井眼里水的氧密度；
4)当井眼充油地层充水时，地层中的含油饱和度S0＝0，RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充油地层充油时，地层中的含油饱和度S0＝1，RPM仪器测井中碳氧比计算公式(1)变形为：
式中，BC油为井眼里油的碳密度，BO油为井眼里油的氧密度；
5)在地层孔隙度其中为确定值时，分别计算近源距和远源距时所述步骤3)中的井眼里水的碳密度BC水、井眼里水的氧密度BO水和所述步骤4)中的井眼里油的碳密度BC油、井眼里油的氧密度BO油；
6)计算当影响因素包括井眼尺寸、套管尺寸和砾石，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
7)计算当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
8)计算当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
9)构建远近碳氧比交会图，在地层孔隙度时，以近源距RPM仪器测井中碳氧比为横坐标，远源距RPM仪器测井中碳氧比为纵坐标构建直角坐标系，将井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比标注在直角坐标系中形成四个基点，四个基点构成四边形，RPM仪器测井中碳氧比测量点都落在四边形区域内；
10)根据远近碳氧比交会图获取井筒持水率和地层含水饱和度。
所述步骤5)中，在地层孔隙度时，计算近源距和远源距时所述步骤3)中的井眼里水的碳密度BC水、井眼里水的氧密度BO水和所述步骤4)中的井眼里油的碳密度BC油、井眼里油的氧密度BO油的方法包括以下步骤：
(1)当井眼充水地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水水为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的近源距RPM仪器测井中碳氧比；bN水水为近源距时每立方厘米砂岩中碳原子的数目；cN水水为近源距时每立方厘米水中氧原子的数目；dN水水为近源距时每立方厘米砂岩中氧原子的数目；
(2)当井眼充水地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水油为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的近源距RPM仪器测井中碳氧比；aN水油为近源距时每立方厘米油中碳原子的数目；bN水油为近源距时每立方厘米砂岩中碳原子的数目；dN水油为近源距时每立方厘米砂岩中氧原子的数目；
(3)将所述步骤(1)和所述步骤(2)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值；
(4)当井眼充水地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水水为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的远源距RPM仪器测井中碳氧比；bF水水为远源距时每立方厘米砂岩中碳原子的数目；cF水水为远源距时每立方厘米水中氧原子的数目；dF水水为远源距时每立方厘米砂岩中氧原子的数目；
(5)当井眼充水地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水油为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的远源距RPM仪器测井中碳氧比；aF水油为远源距时每立方厘米油中碳原子的数目；bF水油为远源距时每立方厘米砂岩中碳原子的数目；dF水油为远源距时每立方厘米砂岩中氧原子的数目；
(6)将所述步骤(4)和所述步骤(5)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值；
(7)当井眼充油地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油水为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比；bN油水为近源距时每立方厘米砂岩中碳原子的数目；cN油水为近源距时每立方厘米水中氧原子的数目；dN油水为近源距时每立方厘米砂岩中氧原子的数目；
(8)当井眼充油地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油油为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比；aN油油为近源距时每立方厘米油中碳原子的数目；bN油油为近源距时每立方厘米砂岩中碳原子的数目；dN油油为近源距时每立方厘米砂岩中氧原子的数目；
(9)将所述步骤(7)和所述步骤(8)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油；
(10)当井眼充油地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油水为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的远源距RPM仪器测井中碳氧比；bF油水为远源距时每立方厘米砂岩中碳原子的数目；cF油水为远源距时每立方厘米水中氧原子的数目；dF油水为远源距时每立方厘米砂岩中氧原子的数目；
(11)当井眼充油地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油油为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的远源距RPM仪器测井中碳氧比；aF油油为远源距时每立方厘米油中碳原子的数目；bF油油为远源距时每立方厘米砂岩中碳原子的数目；dF油油为远源距时每立方厘米砂岩中氧原子的数目；
(12)将所述步骤(10)和所述步骤(11)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油。
所述步骤6)中计算当影响因素包括井眼尺寸、套管尺寸和砾石，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸和砾石时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，为砾石孔隙度；为总孔隙度，其中Vg为砾石含量；b1为每立方厘米砾石中碳原子的数目；d1为每立方厘米砾石中氧原子的数目；S01为砾石含油饱和度；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
所述步骤7)中计算当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸、砾石和泥质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vsh为泥质含量；b2为每立方厘米泥质中碳原子的数目；c2为每立方厘米泥质中氧原子的数目；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
所述步骤8)中，当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vca为钙质含量；b3为每立方厘米钙质中碳原子的数目；d3为每立方厘米钙质中氧原子的数目；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
所述步骤10)中，井筒持水率和地层含水饱和度与RPM仪器测井中碳氧比一一对应：井眼充水地层充水情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为1、地层含水饱和度为1；井眼充水地层充油情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为1、地层含水饱和度为0；井眼充油地层充水情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为0、地层含水饱和度为1；井眼充油地层充油情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为0、地层含水饱和度为0。
本发明由于采取以上技术方案，其具有以下优点：本发明由于采用根据纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，结合在井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况下RPM仪器测井中碳氧比计算公式，计算得到井眼里水的碳密度、井眼里水的氧密度、井眼里油的碳密度和井眼里油的氧密度，进而计算出在不同影响因素下的RPM仪器测井中碳氧比，由于井筒持水率和地层含水饱和度与RPM仪器测井中碳氧比一一对应，因此根据RPM仪器测井中碳氧比能够准确、快速地计算出海上稠油油田RPM测井中井筒持水率和地层含水饱和度。综上所述，本发明可以广泛应用于海上稠油油田RPM测井中井筒持水率和地层含水饱和度的计算中。
附图说明
图1是本发明的当影响因素包括井眼尺寸和套管尺寸时，纯砂岩地层近源距RPM仪器测井中碳氧比随孔隙度的变化规律图，其中，表示在井眼充水地层充水情况下纯砂岩地层近源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充水地层充油情况下纯砂岩地层近源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充油地层充水情况下纯砂岩地层近源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充油地层充油情况下纯砂岩地层近源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版；
图2是本发明的当影响因素包括井眼尺寸和套管尺寸时，纯砂岩地层远源距RPM仪器测井中碳氧比随孔隙度的变化规律图，其中，表示在井眼充水地层充水情况下纯砂岩地层远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充水地层充油情况下纯砂岩地层远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充油地层充水情况下纯砂岩地层远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版，表示在井眼充油地层充油情况下纯砂岩地层远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版；
图3是本发明的当影响因素包括井眼尺寸和套管尺寸，地层孔隙度时的远近碳氧比交会图；
图4是本发明的当影响因素包括井眼尺寸和套管尺寸，地层孔隙度时的远近碳氧比交会图；
图5是本发明的当影响因素包括井眼尺寸和套管尺寸，地层孔隙度时的远近碳氧比交会图；
图6是本发明的当影响因素包括井眼尺寸、套管尺寸和砾石，砾石含量Vg＝0.1，地层孔隙度时的远近碳氧比交会图；
图7是本发明的当影响因素包括井眼尺寸、套管尺寸和砾石，砾石含量Vg＝0.1，地层孔隙度时的远近碳氧比交会图；
图8是本发明的当影响因素包括井眼尺寸、套管尺寸和砾石，砾石含量Vg＝0.1，地层孔隙度时的远近碳氧比交会图；
图9是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，地层孔隙度时的远近碳氧比交会图；
图10是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，地层孔隙度时的远近碳氧比交会图；
图11是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，地层孔隙度时的远近碳氧比交会图；
图12是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，钙质含量Vca＝0.1，地层孔隙度时的远近碳氧比交会图；
图13是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，钙质含量Vca＝0.1，地层孔隙度时的远近碳氧比交会图；
图14是本发明的当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，钙质含量Vca＝0.1，地层孔隙度时的远近碳氧比交会图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示，本发明提供一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，包括以下步骤：
1)利用蒙特卡罗模拟的方法，得到当仅考虑井眼尺寸和套管尺寸时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况下，纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版。
2)根据井眼尺寸和套管尺寸，并利用斯伦贝谢的郝佐格(Hertzog)公式得到RPM仪器测井中碳氧比计算公式为：
式中，C/O为RPM仪器测井中碳氧比；A为碳和氧与中子反应平均截面的比值，是常数；a为每立方厘米油中碳原子的数目；b为每立方厘米砂岩中碳原子的数目；c为每立方厘米水中氧原子的数目；d为每立方厘米砂岩中氧原子的数目；为地层孔隙度。其中a,b,c,d为砂岩地层参数,为已知参数。
S0为地层含油饱和度，BC为井眼里碳密度的贡献；BO为井眼里氧密度的贡献。
3)当井眼充水地层充水时，地层中的含油饱和度S0＝0，井眼里碳密度的贡献BC为井眼里水的碳密度BC水，井眼里氧密度的贡献BO为井眼里水的氧密度BO水，此时RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充水地层充油时，地层中的含油饱和度S0＝1，井眼里碳密度的贡献BC为井眼里水的碳密度BC水，井眼里氧密度的贡献BO为井眼里水的氧密度BO水,此时RPM仪器测井中碳氧比计算公式(1)变形为：
4)当井眼充油地层充水时，地层中的含油饱和度S0＝0，井眼里碳密度的贡献BC为井眼里油的碳密度BC油，井眼里氧密度的贡献Bo为井眼里油的氧密度BO油，BO油，此时RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充油地层充油时，地层中的含油饱和度S0＝1，井眼里碳密度的贡献BC为井眼里油的碳密度BC油，井眼里氧密度的贡献Bo为井眼里油的氧密度BO油，BO油，此时RPM仪器测井中碳氧比计算公式(1)变形为：
5)在地层孔隙度其中为确定值时，分别计算近源距和远源距时步骤3)和步骤4)中的井眼里水的碳密度BC水、井眼里水的氧密度BO水、井眼里油的碳密度BC油和井眼里油的氧密度BO油，包括以下步骤：
(1)当井眼充水地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水水为从步骤1)中近源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的近源距RPM仪器测井中碳氧比；bN水水为近源距时每立方厘米砂岩中碳原子的数目；cN水水为近源距时每立方厘米水中氧原子的数目；dN水水为近源距时每立方厘米砂岩中氧原子的数目。
(2)当井眼充水地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水油为从步骤1)中近源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的近源距RPM仪器测井中碳氧比；aN水油为近源距时每立方厘米油中碳原子的数目；bN水油为近源距时每立方厘米砂岩中碳原子的数目；dN水油为近源距时每立方厘米砂岩中氧原子的数目。
(3)将步骤(1)和步骤(2)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值。
(4)当井眼充水地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水水为从步骤1)中远源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的远源距RPM仪器测井中碳氧比；bF水水为远源距时每立方厘米砂岩中碳原子的数目；cF水水为远源距时每立方厘米水中氧原子的数目；dF水水为远源距时每立方厘米砂岩中氧原子的数目。
(5)当井眼充水地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水油为从步骤1)中远源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的远源距RPM仪器测井中碳氧比；aF水油为远源距时每立方厘米油中碳原子的数目；bF水油为远源距时每立方厘米砂岩中碳原子的数目；dF水油为远源距时每立方厘米砂岩中氧原子的数目。
(6)将步骤(4)和步骤(5)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值。
(7)当井眼充油地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油水为从步骤1)中近源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比；bN油水为近源距时每立方厘米砂岩中碳原子的数目；cN油水为近源距时每立方厘米水中氧原子的数目；dN油水为近源距时每立方厘米砂岩中氧原子的数目。
(8)当井眼充油地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油油为从步骤1)中近源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比；aN油油为近源距时每立方厘米油中碳原子的数目；bN油油为近源距时每立方厘米砂岩中碳原子的数目；dN油油为近源距时每立方厘米砂岩中氧原子的数目。
(9)将步骤(7)和步骤(8)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油。
(10)当井眼充油地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油水为从步骤1)中远源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的远源距RPM仪器测井中碳氧比；bF油水为远源距时每立方厘米砂岩中碳原子的数目；cF油水为远源距时每立方厘米水中氧原子的数目；dF油水为远源距时每立方厘米砂岩中氧原子的数目。
(11)当井眼充油地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油油为从步骤1)中远源距RPM仪器测井中碳氧比随孔隙度的变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的远源距RPM仪器测井中碳氧比；aF油油为远源距时每立方厘米油中碳原子的数目；bF油油为远源距时每立方厘米砂岩中碳原子的数目；dF油油为远源距时每立方厘米砂岩中氧原子的数目。
(12)将步骤(10)和步骤(11)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油。
6)计算当影响因素包括井眼尺寸、套管尺寸和砾石，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸和砾石时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，为砾石孔隙度；为总孔隙度，其中Vg为砾石含量；b1为每立方厘米砾石中碳原子的数目；d1为每立方厘米砾石中氧原子的数目；S01为砾石含油饱和度。
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入公式(14)计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比C/ONCOR水水和远源距RPM仪器测井中碳氧比C/OFCOR水水。
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入公式(14)计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比C/ONCOR水油和远源距RPM仪器测井中碳氧比C/OFCOR水油。
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入公式(14)计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比C/ONCOR油水和远源距RPM仪器测井中碳氧比C/ONCOR油水。
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入公式(14)计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比C/ONCOR油油和远源距RPM仪器测井中碳氧比C/ONCOR油油。
7)计算当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸、砾石和泥质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vsh为泥质含量；b2为每立方厘米泥质中碳原子的数目；c2为每立方厘米泥质中氧原子的数目。
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入公式(15)计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比C/ONCOR水水和远源距RPM仪器测井中碳氧比C/OFCOR水水。
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入公式(15)计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比C/ONCOR水油和远源距RPM仪器测井中碳氧比C/OFCOR水油。
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入公式(15)计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比C/ONCOR油水和远源距RPM仪器测井中碳氧比C/ONCOR油水。
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入公式(15)计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比C/ONCOR油油和远源距RPM仪器测井中碳氧比C/ONCOR油油。
8)当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vca为钙质含量；b3为每立方厘米钙质中碳原子的数目；d3为每立方厘米钙质中氧原子的数目。
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入公式(16)计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比C/ONCOR水水和远源距RPM仪器测井中碳氧比C/OFCOR水水。
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入公式(16)计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比C/ONCOR水油和远源距RPM仪器测井中碳氧比C/OFCOR水油。
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入公式(16)计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比C/ONCOR油水和远源距RPM仪器测井中碳氧比C/ONCOR油水。
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入公式(16)计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比C/ONCOR油油和远源距RPM仪器测井中碳氧比C/ONCOR油油。
9)构建远近碳氧比交会图，在地层孔隙度时，以近源距RPM仪器测井中碳氧比为横坐标，远源距RPM仪器测井中碳氧比为纵坐标构建直角坐标系，将井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比标注在直角坐标系中形成四个基点，分别为WW(C/ONCOR水水，C/OFCOR水水)、WO(C/ONCOR水油，C/OFCOR水油)、OW(C/ONCOR油水，C/OFCOR油水)和OO(C/ONCOR油油，C/OFCOR油油)，四个基点构成四边形，所有RPM仪器测井中碳氧比测量点基本都落在四边形区域内。
10)根据远近碳氧比交会图获取井筒持水率YW和地层含水饱和度SW，井筒持水率YW和地层含水饱和度SW与RPM仪器测井中碳氧比一一对应，其中，点WW(C/ONCOR水水，C/OFCOR水水)对应井筒持水率YW＝1、地层含水饱和度SW＝1；点WO(C/ONCOR水油，C/OFCOR水油)对应井筒持水率YW＝1、地层含水饱和度SW＝0；点OW(C/ONCOR油水，C/OFCOR油水)对应井筒持水率YW＝0、地层含水饱和度SW＝1；点OO(C/ONCOR油油，C/OFCOR油油)对应井筒持水率YW＝0、地层含水饱和度SW＝0。
实施例
1)利用蒙特卡罗模拟的方法，当井眼尺寸为12.25in(英寸)和套管尺寸为9.625in时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况下，纯砂岩地层近源距RPM仪器测井中碳氧比随孔隙度的变化规律如图1所示，纯砂岩地层远源距RPM仪器测井中碳氧比随孔隙度的变化规律如图2所示。
2)影响因素包括井眼尺寸和套管尺寸，根据步骤1)中得到当地层孔隙度时，井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况的纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比(如表1所示)构建在时的远近碳氧比交会图。如图3、图4、图5所示，以近源距RPM仪器测井中碳氧比为横坐标，远源距RPM仪器测井中碳氧比为纵坐标构建直角坐标系，将井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比标注在直角坐标系中形成四个基点分别为WW(C/ONCOR水水，C/OFCOR水水)、WO(C/ONCOR水油，C/OFCOR水油)、OW(C/ONCOR油水，C/OFCOR油水)和OO(C/ONCOR油油，C/OFCOR油油)，四个基点构成四边形。
表1 影响因素包括井眼尺寸、套管尺寸时的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比
3)井筒持水率YW和地层含水饱和度SW与RPM仪器测井中碳氧比一一对应，根据远近碳氧比交会图获取井筒持水率YW和地层含水饱和度SW。其中点WW(C/ONCOR水水，C/OFCOR水水)对应井筒持水率YW＝1、地层含水饱和度SW＝1；点WO(C/ONCOR水油，C/OFCOR水油)对应井筒持水率YW＝1、地层含水饱和度SW＝0；点OW(C/ONCOR油水，C/OFCOR油水)对应井筒持水率YW＝0、地层含水饱和度SW＝1；点OO(C/ONCOR油油，C/OFCOR油油)对应井筒持水率YW＝0、地层含水饱和度SW＝0。
4)将表1所示的地层孔隙度时，井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况的纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，代入公式(2)、公式(3)、公式(4)和公式(5)计算得到在一定的地层孔隙度下近源距和远源距时的井眼里水的碳密度BC水、井眼里水的氧密度BO水、井眼里油的碳密度BC油和井眼里油的氧密度BO油，其中，A＝0.55。当地层孔隙度时，近源距和远源距时的井眼里水的碳密度BC水、井眼里水的氧密度BO水、井眼里油的碳密度BC油和井眼里油的氧密度BO油的计算结果如表2所示。
表2 当地层孔隙度时，近源距和远源距时的井眼里水的碳密度BC水、井眼里水的氧密度BO水、井眼里油的碳密度BC油和井眼里油的氧密度BO油
5)当影响因素包括井眼尺寸、套管尺寸和砾石，砾石含量Vg＝0.1时，令 S0＝0，根据公式(14)分别计算在井眼充水地层充水和井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比C/ONCOR水水和C/ONCOR油水，以及远源距RPM仪器测井中碳氧比C/OFCOR水水和C/OFCOR油水；令 S0＝1，根据公式(14)分别计算在井眼充水地层充油和井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比C/ONCOR水油和C/ONCOR油油，以及远源距RPM仪器测井中碳氧比C/OFCOR水油和C/OFCOR油油，计算结果如下(如表3所示)：
表3 影响因素包括井眼尺寸、套管尺寸和砾石时的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比
6)影响因素包括井眼尺寸、套管尺寸和砾石时，根据步骤5)中计算得到的当地层孔隙度时，井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比(如表3所示)构建在 时的远近碳氧比交会图(如图6、图7、图8所示)，根据远近碳氧比交会图获取井筒持水率YW和地层含水饱和度SW。
7)当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1时，令S0＝0，根据公式(15)分别计算在井眼充水地层充水和井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比C/ONCOR水水和C/ONCOR油水，以及远源距RPM仪器测井中碳氧比C/OFCOR水水和C/OFCOR油水；令S0＝1，根据公式(15)分别计算在井眼充水地层充油和井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比C/ONCOR水油和C/ONCOR油油，以及远源距RPM仪器测井中碳氧比C/OFCOR水油和C/OFCOR油油，计算结果如下(如表4所示)：
表4 影响因素包括井眼尺寸、套管尺寸、砾石和泥质时的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比
8)影响因素包括井眼尺寸、套管尺寸、砾石和泥质时，根据步骤7)中计算得到的当地层孔隙度时，井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比(如表4所示)构建在 时的远近碳氧比交会图(如图9、图10、图11所示)，根据远近碳氧比交会图获取井筒持水率YW和地层含水饱和度SW。
9)当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，砾石含量Vg＝0.1，泥质含量Vsh＝0.1，钙质含量Vca＝0.1时，令S0＝0，根据公式(16)分别计算在井眼充水地层充水和井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比C/ONCOR水水和C/ONCOR油水，以及远源距RPM仪器测井中碳氧比C/OFCOR水水和C/OFCOR油水；令S0＝1，根据公式(16)分别计算在井眼充水地层充油和井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比C/ONCOR水油和C/ONCOR油油，以及远源距RPM仪器测井中碳氧比C/OFCOR水油和C/OFCOR油油，计算结果如下(如表5所示)：
表5 影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质时的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比
10)影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质时，根据步骤9)中计算得到的当地层孔隙度时，井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比(如表5所示)构建在 时的远近碳氧比交会图(如图12、图13、图14所示)，根据远近碳氧比交会图获取井筒持水率YW和地层含水饱和度SW。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。

Claims (6)

1.一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，包括以下步骤：
1)利用蒙特卡罗模拟的方法，得到当仅考虑井眼尺寸和套管尺寸时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充油和井眼充油地层充水四种情况下，纯砂岩地层近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比随地层孔隙度的变化规律图版；
2)根据井眼尺寸和套管尺寸，并利用斯伦贝谢的郝佐格公式得到RPM仪器测井中碳氧比计算公式为：
式中，C/O为RPM仪器测井中碳氧比；A为碳和氧与中子反应平均截面的比值，是常数；a为每立方厘米油中碳原子的数目；b为每立方厘米砂岩中碳原子的数目；c为每立方厘米水中氧原子的数目；d为每立方厘米砂岩中氧原子的数目；为地层孔隙度；S0为地层含油饱和度，BC为井眼里碳密度的贡献；BO为井眼里氧密度的贡献；
3)当井眼充水地层充水时，地层中的含油饱和度S0＝0，RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充水地层充油时，地层中的含油饱和度S0＝1，RPM仪器测井中碳氧比计算公式(1)变形为：
式中，BC水为井眼里水的碳密度，BO水为井眼里水的氧密度；
4)当井眼充油地层充水时，地层中的含油饱和度S0＝0，RPM仪器测井中碳氧比计算公式(1)变形为：
当井眼充油地层充油时，地层中的含油饱和度S0＝1，RPM仪器测井中碳氧比计算公式(1)变形为：
式中，BC油为井眼里油的碳密度，BO油为井眼里油的氧密度；
5)在地层孔隙度其中为确定值时，分别计算近源距和远源距时所述步骤3)中的井眼里水的碳密度BC水、井眼里水的氧密度BO水和所述步骤4)中的井眼里油的碳密度BC油、井眼里油的氧密度BO油；
6)计算当影响因素包括井眼尺寸、套管尺寸和砾石，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
7)计算当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
8)计算当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
9)构建远近碳氧比交会图，在地层孔隙度时，以近源距RPM仪器测井中碳氧比为横坐标，远源距RPM仪器测井中碳氧比为纵坐标构建直角坐标系，将井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比标注在直角坐标系中形成四个基点，四个基点构成四边形，RPM仪器测井中碳氧比测量点都落在四边形区域内；
10)根据远近碳氧比交会图获取井筒持水率和地层含水饱和度。
2.如权利要求1所述的一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，其特征在于：所述步骤5)中，在地层孔隙度时，计算近源距和远源距时所述步骤3)中的井眼里水的碳密度BC水、井眼里水的氧密度BO水和所述步骤4)中的井眼里油的碳密度BC油、井眼里油的氧密度BO油的方法包括以下步骤：
(1)当井眼充水地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水水为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的近源距RPM仪器测井中碳氧比；bN水水为近源距时每立方厘米砂岩中碳原子的数目；cN水水为近源距时每立方厘米水中氧原子的数目；dN水水为近源距时每立方厘米砂岩中氧原子的数目；
(2)当井眼充水地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR水油为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的近源距RPM仪器测井中碳氧比；aN水油为近源距时每立方厘米油中碳原子的数目；bN水油为近源距时每立方厘米砂岩中碳原子的数目；dN水油为近源距时每立方厘米砂岩中氧原子的数目；
(3)将所述步骤(1)和所述步骤(2)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值；
(4)当井眼充水地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水水为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充水情况下的远源距RPM仪器测井中碳氧比；bF水水为远源距时每立方厘米砂岩中碳原子的数目；cF水水为远源距时每立方厘米水中氧原子的数目；dF水水为远源距时每立方厘米砂岩中氧原子的数目；
(5)当井眼充水地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR水油为变化规律图版中得到的当地层的孔隙度时，在井眼充水地层充油情况下的远源距RPM仪器测井中碳氧比；aF水油为远源距时每立方厘米油中碳原子的数目；bF水油为远源距时每立方厘米砂岩中碳原子的数目；dF水油为远源距时每立方厘米砂岩中氧原子的数目；
(6)将所述步骤(4)和所述步骤(5)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里水的碳密度BC水和井眼里水的氧密度BO水的值；
(7)当井眼充油地层充水时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油水为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的近源距RPM仪器测井中碳氧比；bN油水为近源距时每立方厘米砂岩中碳原子的数目；cN油水为近源距时每立方厘米水中氧原子的数目；dN油水为近源距时每立方厘米砂岩中氧原子的数目；
(8)当井眼充油地层充油时，近源距RPM仪器测井中碳氧比的计算式如下：
式中，C/ONCOR油油为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的近源距RPM仪器测井中碳氧比；aN油油为近源距时每立方厘米油中碳原子的数目；bN油油为近源距时每立方厘米砂岩中碳原子的数目；dN油油为近源距时每立方厘米砂岩中氧原子的数目；
(9)将所述步骤(7)和所述步骤(8)中的计算式联立方程组求解得到在近源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油；
(10)当井眼充油地层充水时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油水为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充水情况下的远源距RPM仪器测井中碳氧比；bF油水为远源距时每立方厘米砂岩中碳原子的数目；cF油水为远源距时每立方厘米水中氧原子的数目；dF油水为远源距时每立方厘米砂岩中氧原子的数目；
(11)当井眼充油地层充油时，远源距RPM仪器测井中碳氧比的计算式如下：
式中，C/OFCOR油油为从变化规律图版中得到的当地层的孔隙度时，在井眼充油地层充油情况下的远源距RPM仪器测井中碳氧比；aF油油为远源距时每立方厘米油中碳原子的数目；bF油油为远源距时每立方厘米砂岩中碳原子的数目；dF油油为远源距时每立方厘米砂岩中氧原子的数目；
(12)将所述步骤(10)和所述步骤(11)中的计算式联立方程组求解得到在远源距时，当地层的孔隙度时井眼里油的碳密度BC油和井眼里油的氧密度BO油。
3.如权利要求1所述的一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，其特征在于：所述步骤6)中计算当影响因素包括井眼尺寸、套管尺寸和砾石，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸和砾石时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，为砾石孔隙度；为总孔隙度，其中Vg为砾石含量；b1为每立方厘米砾石中碳原子的数目；d1为每立方厘米砾石中氧原子的数目；S01为砾石含油饱和度；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸和砾石的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
4.如权利要求1所述的一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，其特征在于：所述步骤7)中计算当影响因素包括井眼尺寸、套管尺寸、砾石和泥质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)影响因素包括井眼尺寸、套管尺寸、砾石和泥质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vsh为泥质含量；b2为每立方厘米泥质中碳原子的数目；c2为每立方厘米泥质中氧原子的数目；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石和泥质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
5.如权利要求1所述的一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，其特征在于：所述步骤8)中，当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质，时，在井眼充水地层充水、井眼充水地层充油、井眼充油地层充水和井眼充油地层充油四种情况下的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比，包括以下步骤：
(1)当影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质时，利用斯伦贝谢的郝佐格公式计算RPM仪器测井中碳氧比的计算公式为：
式中，Vca为钙质含量；b3为每立方厘米钙质中碳原子的数目；d3为每立方厘米钙质中氧原子的数目；
(2)在井眼充水地层充水时，S0＝0，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(3)在井眼充水地层充油时，S0＝1，BC＝BC水，BO＝BO水，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(4)在井眼充油地层充水时，S0＝0，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比；
(5)在井眼充油地层充油时，S0＝1，BC＝BC油，BO＝BO油，代入所述步骤(1)中的计算公式计算影响因素包括井眼尺寸、套管尺寸、砾石、泥质和钙质的近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比。
6.如权利要求1所述的一种在RPM测井中计算井筒持水率和地层含水饱和度的方法，其特征在于：所述步骤10)中，井筒持水率和地层含水饱和度与RPM仪器测井中碳氧比一一对应：井眼充水地层充水情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为1、地层含水饱和度为1；井眼充水地层充油情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为1、地层含水饱和度为0；井眼充油地层充水情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为0、地层含水饱和度为1；井眼充油地层充油情况下近源距RPM仪器测井中碳氧比和远源距RPM仪器测井中碳氧比形成的基点对应井筒持水率为0、地层含水饱和度为0。
CN201510358153.8A 2015-06-25 2015-06-25 在rpm测井中计算井筒持水率和地层含水饱和度的方法 Active CN104929627B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510358153.8A CN104929627B (zh) 2015-06-25 2015-06-25 在rpm测井中计算井筒持水率和地层含水饱和度的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201510358153.8A CN104929627B (zh) 2015-06-25 2015-06-25 在rpm测井中计算井筒持水率和地层含水饱和度的方法

Publications (2)

Publication Number Publication Date
CN104929627A true CN104929627A (zh) 2015-09-23
CN104929627B CN104929627B (zh) 2017-07-28

Family

ID=54117020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510358153.8A Active CN104929627B (zh) 2015-06-25 2015-06-25 在rpm测井中计算井筒持水率和地层含水饱和度的方法

Country Status (1)

Country Link
CN (1) CN104929627B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243220A (zh) * 2015-10-26 2016-01-13 中国海洋石油总公司 Rpm-c/o测井消除环空内流体影响的方法
CN105370273A (zh) * 2015-12-10 2016-03-02 中国海洋石油总公司 一种碳氧比测井环境影响因素校正的方法
CN106958443A (zh) * 2017-05-02 2017-07-18 长江大学 一种碳氧比测井钻井液侵入影响因素的校正方法及系统
US11460602B2 (en) 2019-09-05 2022-10-04 Saudi Arabian Oil Company Systems and methods for saturation logging of hydrocarbon wells

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4492864A (en) * 1981-07-31 1985-01-08 Texaco Inc. Neutron well logging
CN101487390A (zh) * 2009-02-23 2009-07-22 大庆油田有限责任公司 一种确定油层原始含油饱和度的阿尔奇模式方法
CN101906963A (zh) * 2010-07-23 2010-12-08 中国石油化工集团公司 利用c/o和地层宏观俘获截面交会技术确定饱和度方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4492864A (en) * 1981-07-31 1985-01-08 Texaco Inc. Neutron well logging
CN101487390A (zh) * 2009-02-23 2009-07-22 大庆油田有限责任公司 一种确定油层原始含油饱和度的阿尔奇模式方法
CN101906963A (zh) * 2010-07-23 2010-12-08 中国石油化工集团公司 利用c/o和地层宏观俘获截面交会技术确定饱和度方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
吴进波: "海上稠油油藏RPM测井解释方法研究", 《中国优秀硕士学位论文全文数据库 基础科学辑》 *
吴进波等: "海上稠油油藏RPM校正方法研究", 《石油工业技术监督》 *
汪忠浩等: "RPM 测井在海上M 油田动态监测中的应用研究", 《2014油气藏监测与管理国际会议（2014 ICRSM)论文集》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105243220A (zh) * 2015-10-26 2016-01-13 中国海洋石油总公司 Rpm-c/o测井消除环空内流体影响的方法
CN105243220B (zh) * 2015-10-26 2018-06-19 中国海洋石油总公司 Rpm-c/o测井消除环空内流体影响的方法
CN105370273A (zh) * 2015-12-10 2016-03-02 中国海洋石油总公司 一种碳氧比测井环境影响因素校正的方法
CN105370273B (zh) * 2015-12-10 2018-01-26 中国海洋石油集团有限公司 一种碳氧比测井环境影响因素校正的方法
CN106958443A (zh) * 2017-05-02 2017-07-18 长江大学 一种碳氧比测井钻井液侵入影响因素的校正方法及系统
CN106958443B (zh) * 2017-05-02 2020-10-09 长江大学 一种碳氧比测井钻井液侵入影响因素的校正方法及系统
US11460602B2 (en) 2019-09-05 2022-10-04 Saudi Arabian Oil Company Systems and methods for saturation logging of hydrocarbon wells

Also Published As

Publication number Publication date
CN104929627B (zh) 2017-07-28

Similar Documents

Publication Publication Date Title
Senger et al. Effects of igneous intrusions on the petroleum system: a review
Feng et al. Production performance of gas hydrate accumulation at the GMGS2-Site 16 of the Pearl River Mouth Basin in the South China Sea
CN104453834A (zh) 一种井组注采关系优化调整方法
CN104929627A (zh) 在rpm测井中计算井筒持水率和地层含水饱和度的方法
Zhou et al. Experimental and numerical study on spontaneous imbibition of fracturing fluids in the horn river shale gas formation
Li et al. Research on sequence stratigraphy, hydrogeological units and commingled drainage associated with coalbed methane production: a case study in Zhuzang syncline of Guizhou province, China
Haizlip et al. The impact of high non condensable gas concentrations on well performance—Kizildere Geothermal Reservoir, Turkey
Wei et al. A state-of-the-art review and prospect of gas hydrate reservoir drilling techniques
Chao et al. Sedimentary characteristics, evolution and controlling factors of the alluvial fans, the north of Qinghai Lake
Sun et al. Numerical simulation of produced water reinjection technology for water-soluble gas recovery
Temizel et al. An analysis of scale buildup in seawater injection of waterflooding operations
Long et al. Numerical investigation on gas hydrate production by depressurization in hydrate-bearing reservoir
Moridis et al. Numerical simulation studies of gas production scenarios from hydrate accumulations at the Mallik Site, McKenzie Delta, Canada
Haeckel et al. 3-D numerical modelling of gas hydrate accumulations at the Alaska North Slope
Bacon et al. Validation of CO2 injection simulations with monitoring well data
Wach et al. Well placement, cost reduction, and increased production using reservoir models based on outcrop, core, well-log, seismic data, and modern analogs: Onshore and offshore Western Trinidad
CN110469321B (zh) 一种确定地层破裂压力梯度的录井方法
Garven Fault-Related CO2 Degassing, Geothermics, and Fluid Flow in Southern California Basins--Physiochemical Evidence and Modeling
Tian et al. Potential evaluation of CO^ sub 2^ storage and enhanced oil recovery of tight oil reservoir in the Ordos Basin, China
Haklıdır et al. Geochemistry of fluids in new wells of Kizildere geothermal field in Turkey
Boles et al. Fault-related CO2 degassing, geothermics, and fluid flow in southern California basins---Physiochemical evidence and modeling
Akhmetzhanova et al. Water exploration works for reservoir pressure maintenance at Brown Field X
Xu et al. Analysis of residual oil recovery based on fluvial reservoir models with architecture interface
Hai-min et al. Development technique and practice of horizontal wells for complex fault-block reservoirs in Jidong oil field
Adu‐Gyamfi et al. Simulation study of chemo‐mechanical impacts of CO2 injection in morrow b sandstone reservoir

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CP01 Change in the name or title of a patent holder

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC Research Institute

Patentee before: China National Offshore Oil Corporation

CP01 Change in the name or title of a patent holder