[go: up one dir, main page]

CN106599442A - Recognition and evaluation method and device for physical properties of while-drilling reservoir based on comprehensive logging parameters - Google Patents

Recognition and evaluation method and device for physical properties of while-drilling reservoir based on comprehensive logging parameters Download PDF

Info

Publication number
CN106599442A
CN106599442A CN201611126856.9A CN201611126856A CN106599442A CN 106599442 A CN106599442 A CN 106599442A CN 201611126856 A CN201611126856 A CN 201611126856A CN 106599442 A CN106599442 A CN 106599442A
Authority
CN
China
Prior art keywords
drill bit
energy consumption
reservoir
work
judging
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
CN201611126856.9A
Other languages
Chinese (zh)
Other versions
CN106599442B (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.)
Daqing Drilling Engineering Co ltd
China National Petroleum Corp
Original Assignee
Daqing Petroleum Administration Bureau
China National Petroleum Corp
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 Daqing Petroleum Administration Bureau, China National Petroleum Corp filed Critical Daqing Petroleum Administration Bureau
Priority to CN201611126856.9A priority Critical patent/CN106599442B/en
Publication of CN106599442A publication Critical patent/CN106599442A/en
Application granted granted Critical
Publication of CN106599442B publication Critical patent/CN106599442B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Earth Drilling (AREA)

Abstract

本发明公开了一种基于综合录井参数的随钻储层物性识别评价方法和装置,涉及石油勘探开发领域,包括第一步,采集综合录井工程资料;第二步,计算钻头做功能耗E i ;第三步,计算钻头做功时垂向钻头做功能耗W H 和切向钻头做功能耗W L ;第四步,计算钻头做功能耗趋势值;第五步,计算钻头做功能耗比值W bi ;第六步,计算钻头做功能耗比值W bi 整体样本标准差σ;第七步,计算交汇面积S,判断储层孔隙度大小;第八步,储层进行物性分级;以及一种基于综合录井参数的随钻储层物性识别评价装置。解决钻后物性识别评价的时效性滞后、连续性差、准确性差等问题。

The invention discloses a method and device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive mud logging parameters, and relates to the field of petroleum exploration and development, including the first step of collecting comprehensive mud logging engineering data; the second step of calculating the power consumption of the drill E i ; the third step is to calculate the power consumption W H of the vertical drill bit and the power consumption W L of the tangential drill bit when the drill bit is working; the fourth step is to calculate the trend value of the power consumption of the drill bit ; The fifth step is to calculate the working power consumption ratio W bi of the drill bit; the sixth step is to calculate the overall sample standard deviation σ of the working power consumption ratio W bi of the drill bit; the seventh step is to calculate the intersection area S to judge the porosity of the reservoir; the eighth step The first step is to classify the physical properties of the reservoir; and a device for identifying and evaluating the physical properties of the reservoir while drilling based on comprehensive logging parameters. Solve the problems of time lag, poor continuity, and poor accuracy in the identification and evaluation of physical properties after drilling.

Description

基于综合录井参数的随钻储层物性识别评价方法和装置Reservoir physical property identification and evaluation method and device while drilling based on comprehensive logging parameters

技术领域technical field

本发明涉及石油勘探开发领域,具体说是一种基于综合录井参数的随钻储层物性识别评价方法和装置。The invention relates to the field of petroleum exploration and development, in particular to a method and device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive mud logging parameters.

背景技术Background technique

钻头破碎岩石所需能量与岩石强度有关,钻头做功实现破岩钻进;R.Teale提出了在钻进岩石过程中机械比能(钻头做功能耗)的概念,即钻头在钻压和扭矩作用下破碎单位体积岩石做的功(所需要的机械能量)。这一概念将破碎单位体积岩石所需能量它反映钻头的破岩效率,因此为钻井提供了一种评价钻井效率的方法,比能越大,说明钻头破岩效率越低,钻头与地层适应性越差,钻井参数有待优化。The energy required for the drill bit to break the rock is related to the strength of the rock. The drill bit works to achieve rock-breaking drilling; R.Teale proposed the concept of mechanical specific energy (drill bit work consumption) in the process of drilling rock, that is, the drill bit under the action of drilling pressure and torque The work done by crushing a unit volume of rock (the mechanical energy required). This concept reflects the energy required to break a unit volume of rock. It reflects the rock-breaking efficiency of the drill bit, so it provides a method to evaluate the drilling efficiency for drilling. The larger the specific energy, the lower the rock-breaking efficiency of the drill bit, and the adaptability of the drill bit to the formation. The worse it is, the drilling parameters need to be optimized.

钻头破碎岩石所需能量与岩石强度有关,因而钻头做功大小与地层物性好坏(孔隙度、渗透率)呈正相关,为此可为地层物性评价提供依据,钻头做功能耗越大,地层强度越大,物性条件越差,反之越好。现有的地层物性评价技术中由于都是基于钻后的物性评价,因此在物性评价的时效性和随钻方面处于空白,受钻井周期和其他工程因素的影响,尤其在碳酸盐岩地层、火山岩地层、以及泥岩裂缝型储层中,钻后评价存在着严重的滞后性。The energy required for the drill bit to break the rock is related to the strength of the rock, so the amount of work done by the drill bit is positively correlated with the physical properties of the formation (porosity, permeability), which can provide a basis for the evaluation of the physical properties of the formation. The greater the power consumption of the drill bit, the stronger the formation strength Larger, the worse the physical condition, and vice versa. Existing formation physical property evaluation technologies are all based on physical property evaluation after drilling, so there is a gap in the timeliness of physical property evaluation and while drilling, affected by drilling cycle and other engineering factors, especially in carbonate formations, In volcanic rock formations and fractured mudstone reservoirs, there is a serious lag in post-drilling evaluation.

发明内容Contents of the invention

有鉴于此,本发明提供一种基于综合录井参数的随钻储层物性识别评价方法和装置,可以在随钻时利用钻头做功能耗有效地评估出储层孔隙度大小以及对储层进行物性分级,以解决钻后物性识别评价的时效性滞后、连续性差、准确性差的问题。In view of this, the present invention provides a method and device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive logging parameters, which can effectively evaluate the porosity of reservoirs by using the power consumption of drill bits while drilling Classification of physical properties to solve the problems of time lag, poor continuity and poor accuracy of physical property identification and evaluation after drilling.

第一方面,本发明提供了一种基于综合录井参数的随钻储层物性识别评价方法,其特征在于,包括:In the first aspect, the present invention provides a method for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive logging parameters, which is characterized in that it includes:

第一步,采集综合录井工程资料;The first step is to collect comprehensive logging engineering data;

所述综合录井工程资料通过综合录井仪采集钻头每进尺0.1m的钻压Pi、钻头扭矩T、钻头转速n、机械钻速v、钻头直径dB和钻时Z数据;The comprehensive mud logging engineering data is collected through the comprehensive mud logging instrument, the drilling pressure Pi, the drill bit torque T, the drill bit speed n , the ROP v, the drill bit diameter d B and the drilling time Z data of the drill bit per footage of 0.1m;

第二步,计算钻头做功能耗EiThe second step is to calculate the power consumption E i of the drill bit;

所述钻头做功能耗Ei通过经典机械比能模型计算每0.1m的钻头做功能耗EiThe power consumption E i of the drill bit is calculated through the classical mechanical specific energy model to calculate the power consumption E i of the drill bit per 0.1m;

第三步,计算钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WLThe third step is to calculate the functional consumption W H of the vertical drill bit and the functional consumption W L of the tangential drill bit when the drill bit works;

第四步,计算钻头做功能耗趋势值 The fourth step is to calculate the power consumption trend value of the drill bit

所述计算钻头做功能耗趋势值是m个所述钻头做功能耗Ei平均值,m∈[3,10],j=1,2,3,4;The calculation of drill bit power consumption trend value is the average value of energy consumption E i of m drill bits, m∈[3,10], j=1, 2, 3, 4;

第五步,计算钻头做功能耗比值WbiThe fifth step is to calculate the power consumption ratio W bi of the drill bit;

所述钻头做功能耗比值Wbi为所述钻头做功能耗Ei和所述钻头做功能耗趋势值的比值,Wbi∈(0,1);The power consumption ratio W bi of the drill bit is the power consumption E i of the drill bit and the trend value of the power consumption of the drill bit The ratio of W bi ∈ (0,1);

第六步,计算钻头做功能耗比值Wbi整体样本标准差σ;The sixth step is to calculate the overall sample standard deviation σ of the power consumption ratio W bi of the drill bit;

第七步,计算交汇面积S,判断储层孔隙度大小;The seventh step is to calculate the intersection area S and determine the porosity of the reservoir;

所述交汇面积S通过利用所述垂向钻头做功能耗WH和所述切向钻头做功能耗WL对钻头深度h的积分得出:The intersection area S is obtained by integrating the work power consumption W H of the vertical drill bit and the work power consumption W L of the tangential drill bit to the depth h of the drill bit:

其中h为钻头深度,h2>h1Where h is the depth of the drill bit, h 2 >h 1 ;

第八步,储层进行物性分级;In the eighth step, the physical properties of the reservoir are classified;

所述储层进行物性分级的判断方式如下:The method for judging the physical property classification of the reservoir is as follows:

第一次判断:First judgment:

记所述钻头做功能耗比值Wbi为第一个样本Wbj,计算所述第一个样本Wbj的平均值,记为第一个平均值 Record the working power consumption ratio Wbi of the drill bit as the first sample Wbj , calculate the average value of the first sample Wbj , and record it as the first average value

所述且Wbj<1-σ,则判断为储层段,满足所述且Wbj<1-σ的Wbj记为第二个子样本Wcj,进入第二次判断,否则为非储层段,停止判断;said And W bj <1-σ, then it is judged as a reservoir section, which satisfies the And W bj with W bj <1-σ is recorded as the second sub-sample W cj , and enters the second judgment, otherwise it is a non-reservoir section, and the judgment is stopped;

第二次判断:Second judgment:

计算所述第二个子样本Wcj的平均值,记为第二个平均值 Calculate the average value of the second sub-sample W cj , denoted as the second average value

所述且1-2σ<Wcj<1-σ,满足且1-2σ<Wcj<1-σ的所有值Wcj记为第三个子样本Wdj,进入第三次判断,否则判断为次级段,停止判断;said And 1-2σ<W cj <1-σ, satisfy And all values W cj of 1-2σ<W cj <1-σ are recorded as the third sub-sample W dj , and enter the third judgment, otherwise it is judged as a secondary segment, and the judgment is stopped;

第三次判断:The third judgment:

计算第三个子样本Wdj的平均值,记为第三个平均值 Calculate the average value of the third sub-sample W dj , denoted as the third average

所述且1-3σ<Wdj<1-2σ,满足所述且1-3σ<Wdj<1-2σ的所有值Wdj记为第四个子样本Wej,进入第四次判断,否则判断为中级段,停止判断;said And 1-3σ<W dj <1-2σ, satisfying the And all values W dj of 1-3σ<W dj <1-2σ are recorded as the fourth sub-sample W ej , and enter the fourth judgment, otherwise it is judged as an intermediate stage, and the judgment is stopped;

第四次判断:Fourth judgment:

计算第四个子样本Wej的平均值,记为第四个平均值 Calculate the average value of the fourth sub-sample W ej , denoted as the fourth average

所述且Wej<1-3σ,则判断为优级段,停止判断,否则判断为良级段,停止判断。said And if W ej <1-3σ, it is judged as an excellent segment, and the judgment is stopped; otherwise, it is judged as a good segment, and the judgment is stopped.

优选地,所述的一种基于综合录井参数的随钻储层物性识别评价方法的经典机械比能模型为:Preferably, the classical mechanical specific energy model of the described method for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive logging parameters is:

优选地,所述的一种基于综合录井参数的随钻储层物性识别评价方法的钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算公式分别为;Preferably, the calculation formulas of the vertical drill bit power consumption W H and the tangential drill bit power consumption W L of the described method for identifying and evaluating the physical properties of the reservoir while drilling based on comprehensive mud logging parameters are as follows;

其中,a为同一钻头最大钻压的地层经验数据;b为同一钻头最大转速的地层经验数据。Among them, a is the formation experience data of the maximum WOB of the same drill bit; b is the formation experience data of the maximum speed of the same drill bit.

优选地,所述的一种基于综合录井参数的随钻储层物性识别评价方法的所述钻头做功能耗Ei采用五点钟形法滤波处理,所述经滤波处理公式为Ei *=β(Ei-2+Ei+2)+γ(Ei-1+Ei+1)+εEiPreferably, the power consumption E i of the drill bit in the method for identifying and evaluating the physical properties of reservoirs while drilling based on comprehensive mud logging parameters is filtered by the five-point bell method, and the filtered formula is E i * =β(E i-2 +E i+2 )+γ(E i-1 +E i+1 )+εE i ;

其中,所述β取值为0.11,γ取值为0.24,ε取值为0.3。Wherein, the value of β is 0.11, the value of γ is 0.24, and the value of ε is 0.3.

优选地,所述的一种基于综合录井参数的随钻储层物性识别评价方法的计算钻头做功能耗趋势值的方法采用滑动均值法;Preferably, the calculation of the drill bit power consumption trend value of the described method for identifying and evaluating the physical properties of the reservoir while drilling based on comprehensive logging parameters The method adopts the moving average method;

其中, in,

m∈[3,10],j=2,3,4。m∈[3,10],j=2,3,4.

优选地,所述的一种基于综合录井参数的随钻储层物性识别评价方法中m取值为5。Preferably, the value of m is 5 in the method for identifying and evaluating the physical properties of reservoirs while drilling based on comprehensive mud logging parameters.

第二方面,本发明提供了一种基于综合录井参数的随钻储层物性识别评价装置,包括:综合录井工程资料采集装置,所述综合录井工程资料采集装置分别与钻头做功能耗Ei计算装置和钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算装置连接;所述钻头做功能耗Ei计算装置与钻头做功能耗趋势值计算装置连接,所述钻头做功能耗趋势值计算装置与钻头做功能耗比值Wbi计算装置连接,所述钻头做功能耗比值Wbi计算装置与储层物性分级装置连接;所述钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算装置与交汇面积S计算与储层孔隙度大小判断装置连接;In a second aspect, the present invention provides a device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive mud logging parameters, including: a comprehensive mud logging engineering data acquisition device, and the comprehensive mud logging engineering data acquisition device is respectively connected to the drill bit for power consumption When the E i calculating device is connected with the drill bit working power consumption W H of the vertical drill bit and the calculating device of the tangential drilling bit working power consumption W L ; The computing device is connected to the trend value of the power consumption of the drill bit The calculation device is connected with the calculating device of the working power consumption ratio W bi of the drill bit, and the calculating device of the working power consumption ratio W bi of the drilling bit is connected with the reservoir physical property grading device; when the drilling bit works, the working power consumption W H of the vertical drill bit and the power consumption W H of the tangential drill bit The working power consumption W L calculation device is connected with the intersection area S calculation and the reservoir porosity judgment device;

所述综合录井工程资料采集装置采集钻头每进尺0.1m的综合录井工程资料,包括钻压Pi、钻头扭矩T、钻头转速n、机械钻速v、钻头直径dB和钻时Z;The integrated mud logging engineering data acquisition device collects the comprehensive mud logging engineering data of the drill bit per footage of 0.1m, including the pressure on bit Pi, the bit torque T, the bit speed n , the ROP v, the bit diameter d B and the drilling time Z;

所述钻头做功能耗Ei计算装置,用于计算每0.1m的钻头做功能耗EiThe drill bit working power consumption E i calculation device is used to calculate the drill bit working power consumption E i per 0.1m;

所述钻头做功能耗趋势值计算装置用于计算m个所述钻头做功能耗Ei平均值得到钻头做功能耗趋势值m∈[3,10],j=1,2,3,4;The power consumption trend value of the drill bit The calculation device is used to calculate the average value of the power consumption E i of the m drill bits to obtain the trend value of the power consumption of the drill bit m∈[3,10],j=1,2,3,4;

所述钻头做功能耗比值Wbi计算装置用于计算所述钻头做功能耗Ei和所述钻头做功能耗趋势值的比值,得到钻头做功能耗比值WbThe working power consumption ratio W bi calculation device of the drill bit is used to calculate the working power consumption E i of the drill bit and the trend value of the working power consumption of the drill bit The ratio of the working power consumption of the drill bit is obtained W b ;

所述交汇面积S计算与储层孔隙度大小判断装置用于计算交汇面积S,根据所述交汇面积S判断储层孔隙度大小,其中所述交汇面积S通过利用所述垂向钻头做功能耗WH和所述切向钻头做功能耗WL对钻头深度h的积分得出;The intersection area S calculation and reservoir porosity size judging device is used to calculate the intersection area S, and judge the reservoir porosity size according to the intersection area S, wherein the intersection area S uses the vertical drill bit to perform power consumption The integration of W H and the power consumption W L of the tangential drill bit to the depth h of the drill bit is obtained;

所述储层物性分级装置对储层进行物性分级的判断方式如下:The method for judging the physical property classification of the reservoir by the reservoir physical property classification device is as follows:

第一次判断装置:The first judgment device:

记所述钻头做功能耗比值Wbi为第一个样本Wbj,计算所述第一个样本Wbj的平均值,记为第一个平均值 Record the working power consumption ratio Wbi of the drill bit as the first sample Wbj , calculate the average value of the first sample Wbj , and record it as the first average value

所述且Wbj<1-σ,则判断为储层段,满足所述且Wbj<1-σ的Wbj记为第二个子样本Wcj,进入第二次判断装置,否则为非储层段,停止判断;said And W bj <1-σ, then it is judged as a reservoir section, which satisfies the And W bj with W bj <1-σ is recorded as the second sub-sample W cj , and enters the second judgment device, otherwise it is a non-reservoir section, and the judgment is stopped;

第二次判断装置:The second judgment device:

计算所述第二个子样本Wcj的平均值,记为第二个平均值 Calculate the average value of the second sub-sample W cj , denoted as the second average value

所述且1-2σ<Wcj<1-σ,满足且1-2σ<Wcj<1-σ的所有值Wcj记为第三个子样本Wdj,进入第三次判断装置,否则判断为次级段,停止判断;said And 1-2σ<W cj <1-σ, satisfy And all values W cj of 1-2σ<W cj <1-σ are recorded as the third sub-sample W dj , and enter the third judgment device, otherwise it is judged as a secondary stage, and the judgment is stopped;

第三次判断装置:The third judgment device:

计算第三个子样本Wdj的平均值,记为第三个平均值 Calculate the average value of the third sub-sample W dj , denoted as the third average

所述且1-3σ<Wdj<1-2σ,满足所述且1-3σ<Wdj<1-2σ的所有值Wdj记为第四个子样本Wej,进入第四次判断装置,否则判断为中级段,停止判断;said And 1-3σ<W dj <1-2σ, satisfying the And all values W dj of 1-3σ<W dj <1-2σ are recorded as the fourth sub-sample W ej , and enter the fourth judgment device, otherwise it is judged as an intermediate stage, and the judgment is stopped;

第四次判断装置:The fourth judging device:

计算第四个子样本Wej的平均值,记为第四个平均值 Calculate the average value of the fourth sub-sample W ej , denoted as the fourth average

所述且Wej<1-3σ,则判断为优级段,停止判断,否则判断为良级段,停止判断。said And if W ej <1-3σ, it is judged as an excellent segment, and the judgment is stopped; otherwise, it is judged as a good segment, and the judgment is stopped.

优选地,所述所述钻头做功能耗Ei计算装置采用经典机械比能模型进行计算;Preferably, the calculation device for the energy consumption E i of the drill adopts a classical mechanical specific energy model for calculation;

其中,所述经典机械比能模型为:Wherein, the classical mechanical specific energy model is:

优选地,所述钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算装置的计算公式分别为:Preferably, when the drill bit is working, the calculation formulas of the vertical drill bit power consumption W H and the tangential drill bit power consumption W L calculation device are respectively:

其中,a为同一钻头最大钻压的地层经验数据;b为同一钻头最大转速的地层经验数据。Among them, a is the formation experience data of the maximum WOB of the same drill bit; b is the formation experience data of the maximum speed of the same drill bit.

优选地,所述钻头做功能耗Ei计算装置采用五点钟形法滤波处理,所述经滤波处理公式为Ei *=β(Ei-2+Ei+2)+γ(Ei-1+Ei+1)+εEiPreferably, the calculation device for the power consumption E i of the drill adopts a five-point bell method for filtering, and the filtered formula is E i * = β(E i-2 +E i+2 )+γ(E i -1 +E i+1 )+εE i ;

其中,所述β取值为0.11,γ取值为0.24,ε取值为0.3;Wherein, the value of β is 0.11, the value of γ is 0.24, and the value of ε is 0.3;

优选地,所述钻头做功能耗趋势点计算装置采用滑动均值法求取;Preferably, the drill bit acts as a power consumption trend point The calculation device adopts the sliding average method to calculate;

其中,m∈[3,10],j=1,2,3,4,…。in, m∈[3,10], j=1, 2, 3, 4, . . .

本发明具有如下有益效果:The present invention has following beneficial effects:

本发明提供一种基于综合录井参数的随钻储层物性识别评价方法和装置,可以在随钻时计算出垂向和切向钻头做功的交汇面积S,利用交汇面积S有效地评估出储层孔隙度大小,同时利用求取的钻头做功参数可以对储层进行物性分级,以解决测井物性评价时效性滞后、录井物性评价连续性差、准确性差等问题,为获取地层第一手物性资料提供可行性方案,弥补了测井在钻后评价中因井眼垮塌等因素影响而造成的物性评价失真的问题。The invention provides a method and device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive logging parameters, which can calculate the intersection area S of vertical and tangential drill bit work while drilling, and use the intersection area S to effectively evaluate the reservoir. At the same time, the physical properties of the reservoir can be classified by using the obtained drill bit work parameters, so as to solve the problems of time lag in physical property evaluation of well logging, poor continuity and poor accuracy of physical property evaluation of mud logging, etc. The data provide a feasible solution, which makes up for the distortion of physical property evaluation caused by factors such as borehole collapse in the post-drilling evaluation of well logging.

附图说明Description of drawings

通过以下参考附图对本发明实施例的描述,本发明的上述以及其它目的、特征和优点更为清楚,在图例中:Through the following description of the embodiments of the present invention with reference to the accompanying drawings, the above and other objects, features and advantages of the present invention are more clear, in the legend:

图1是本发明的储层物性识别评价方法的流程框图;Fig. 1 is the block flow diagram of the method for identifying and evaluating reservoir physical properties of the present invention;

图2是本发明判断储层孔隙度大小的流程示意图;Fig. 2 is the schematic flow chart of judging reservoir porosity size in the present invention;

图3是本发明判断储层进行物性分级的流程示意图;Fig. 3 is a schematic flow chart of judging the physical property classification of the reservoir in the present invention;

图4是本发明的储层物性识别评价装置的示意框图;Fig. 4 is a schematic block diagram of the reservoir physical property identification and evaluation device of the present invention;

图5是本发明判断储层进行物性分级装置的流程示意图。Fig. 5 is a schematic flowchart of the device for judging reservoirs and classifying physical properties of the present invention.

具体实施方式detailed description

以下基于实例对本发明进行描述,但是值得说明的是,本发明并不限于这些实施例。在下文对本发明的细节描述中,详尽描述了一些特定的细节部分。然而,对于没有详尽描述的部分,本领域技术人员也完全可以理解本发明。The present invention is described below based on examples, but it should be noted that the present invention is not limited to these examples. In the following detailed description of the invention, some specific details are set forth in detail. However, the present invention can be fully understood by those skilled in the art even for parts not described in detail.

此外,本领域普通技术人员应当理解,所提供的附图只是为说明本发明的目的、特征和优点,附图并不是实际按照比例绘制的。In addition, those of ordinary skill in the art should understand that the provided drawings are only for illustrating the objects, features and advantages of the present invention, and the drawings are not actually drawn to scale.

同时,除非上下文明确要求,否则整个说明书和权利要求书中的“包括”、“包含”等类似词语应当解释为包含的含义而不是排他或穷举的含义;也就是说,是“包含但不限于”的含义。At the same time, unless the context clearly requires, the words "include", "include" and other similar words in the entire specification and claims should be interpreted as an inclusive meaning rather than an exclusive or exhaustive meaning; that is, "include but not limited to the meaning of ".

图1是本发明的储层物性识别评价方法的流程框图。如图1所示,第一步101,采集综合录井工程资料,通过综合录井仪采集钻头每进尺0.1m的钻压Pi、钻头扭矩T、钻头转速n、机械钻速v、钻头直径dB和钻时Z参数。Fig. 1 is a flow chart of the reservoir physical property identification and evaluation method of the present invention. As shown in Fig. 1, the first step 101 is to collect comprehensive mud logging engineering data, and collect the drilling pressure P i , drill bit torque T, drill bit speed n, ROP v, and drill bit diameter per 0.1m of the drill bit through the comprehensive mud logging instrument d B and Z parameters while drilling.

第二步102,计算钻头做功能耗Ei,钻头做功能耗Ei通过经典机械比能模型计算每0.1m的钻头做功能耗EiThe second step 102 is to calculate the power consumption E i of the drill bit , and calculate the power consumption E i of the drill bit per 0.1m through the classical mechanical specific energy model;

所述经典机械比能模型为:The classical mechanical specific energy model is:

对所述钻头做功能耗Ei采用五点钟形法滤波处理,经滤波处理得到Ei=0.11(Ei-2+Ei+2)+0.24(Ei-1+Ei+1)+0.3EiThe power consumption E i of the drill bit is filtered by the five-point bell method, and E i =0.11(E i-2 +E i+2 )+0.24(E i-1 +E i+1 ) is obtained after filtering +0.3E i ;

第三步103,计算钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL,钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL通过第一步101所述综合录井工程资料、垂向钻头做功能耗公式和切向钻头做功能耗公式计算得出;The third step 103 is to calculate the functional consumption W H of the vertical drill and the functional consumption W L of the tangential drill when the drill is working. The functional consumption W H of the vertical drill and the functional consumption W L of the tangential drill are passed through the first step The comprehensive mud logging engineering data described in step 101, the vertical drill bit energy consumption formula and the tangential drill bit energy consumption formula are calculated;

其中,垂向钻头做功能耗公式 Among them, the energy consumption formula of the vertical drill

其中,切向钻头做功能耗公式 Among them, the energy consumption formula of tangential drill

公式中,a为同一钻头最大钻压的地层经验数据;b为同一钻头最大转速的地层经验数据。In the formula, a is the formation experience data of the maximum drilling pressure of the same drill bit; b is the formation experience data of the maximum speed of the same drill bit.

第四步104,计算钻头做功能耗趋势值所述计算钻头做功能耗趋势值采用滑动均值法计算,选取5个所述钻头做功能耗Ei作为一组,依次计算钻头做功能耗趋势值 The fourth step 104 is to calculate the power consumption trend value of the drill bit The calculation of drill bit power consumption trend value Using the sliding average method to calculate, select 5 said drill bits as a group, and calculate the trend value of the drill bit’s energy consumption in turn

其中,j=2,3,4。Among them, j=2,3,4.

第五步105,计算钻头做功能耗比值Wbi,钻头做功能耗比值Wbi为所述钻头做功能耗Ei和所述钻头做功能耗趋势值的比值, The fifth step 105 is to calculate the power consumption ratio of the drill bit W bi , the power consumption ratio W bi of the drill bit is the power consumption E i of the drill bit and the trend value of the power consumption of the drill bit ratio of

第六步106,计算钻头做功能耗比值Wbi整体样本标准差σ;整体样本标准差为:The sixth step 106 is to calculate the overall sample standard deviation σ of the working power consumption ratio W bi of the drill bit; the overall sample standard deviation is:

其中,μ为Wbi的数学期望值,样本总量为N;Among them, μ is the mathematical expectation value of W bi , and the total sample size is N;

第七步107,计算交汇面积S,所述交汇面积S通过利用所述垂向钻头做功能耗WH和所述切向钻头做功能耗WL对钻头深度h的积分得出:In the seventh step 107, the intersection area S is calculated, and the intersection area S is obtained by integrating the energy consumption W H of the vertical drill bit and the energy consumption W L of the tangential drill bit on the depth h of the drill bit:

其中h为钻头深度,h2>h1;交汇面积S具体判断储层孔隙度大小的方法在图2中进行详细说明。Where h is the depth of the drill bit, h 2 >h 1 ; the method of judging the porosity of the reservoir by the intersection area S is described in detail in Fig. 2 .

第八步108,对储层进行物性分级,根据样本标准差σ、钻头做功能耗比值Wbi均值及Wbi子样本的均值对储层进行物性分级,对储层进行物性分级的判断方式在图3中详细描述。The eighth step 108 is to classify the physical properties of the reservoir, and classify the physical properties of the reservoir according to the sample standard deviation σ, the average value of the drill bit power consumption ratio W bi and the average value of the W bi sub-samples, and the judgment method for the physical property classification of the reservoir is Described in detail in Figure 3.

图2是本发明判断储层间空隙大小的流程示意图。如图2所示,在图1中第七步107计算交汇面积S,通过步骤201对交汇面积S是否大于0进行判断,交汇面积S>0时,为正向交汇202,说明储层孔隙度小,交汇面积S≤0时,为负向交汇204,储层孔隙度大;正向交汇202时,交汇面积S越大表示储层孔隙越不发育,物性越差;负向交汇204时,交汇面积S越大表示储层孔隙越发育,物性越好。Fig. 2 is a schematic flow chart of judging the size of inter-reservoir gaps in the present invention. As shown in Figure 2, the seventh step 107 in Figure 1 calculates the intersection area S, and judges whether the intersection area S is greater than 0 through step 201. When the intersection area S>0, it is a positive intersection 202, indicating the porosity of the reservoir Small, when the intersection area S≤0, it is a negative intersection 204, and the porosity of the reservoir is large; when the positive intersection 202, the larger the intersection area S indicates that the reservoir pores are less developed and the physical properties are worse; when the negative intersection 204, The larger the intersection area S, the more developed the reservoir pores and the better the physical properties.

图3是本发明判断储层进行物性分级的流程示意图。如图3所示:Fig. 3 is a schematic flow chart of judging the physical property classification of the reservoir in the present invention. As shown in Figure 3:

步骤301,记所述钻头做功能耗比值Wbi为第一个样本Wbj,,进入步骤302;Step 301, record the working power consumption ratio W bi of the drill bit as the first sample W bj, , and enter step 302;

步骤302,计算第一个样本Wbj的平均值,记为第一个平均值 进入步骤303;Step 302, calculate the average value of the first sample W bj , recorded as the first average value Go to step 303;

步骤303,第一个样本Wbj<第一个平均值且第一个样本Wbj<1-σ,则判断为储层段B,进入步骤304,否则为非储层段A,停止判断;Step 303, the first sample W bj < the first average value And if the first sample W bj <1-σ, then it is judged as reservoir section B, and enter step 304, otherwise, it is non-reservoir section A, and stop judging;

步骤304,记录满足且Wbj<1-σ的Wbj,记为第二个子样本Wcj,进入步骤305;Step 304, the record satisfies W bj with W bj <1-σ is recorded as the second sub-sample W cj , and enters step 305;

步骤305,计算所述第二个子样本Wcj的平均值,记为第二个平均值 Step 305, calculating the average value of the second sub-sample W cj , which is recorded as the second average value

其中,m为Wcj的个数,进入步骤306;which is Wherein, m is the number of Wcj, enter step 306;

步骤306,第二个子样本Wcj<第二个平均值且1-2σ<第二个子样本Wcj<1-σ,进入步骤307,否则判断为次级段C,停止判断;Step 306, the second sub-sample W cj < the second mean value And 1-2σ<the second sub-sample W cj <1-σ, go to step 307, otherwise it is judged as sub-section C, stop judging;

步骤307,记录满足且1-2σ<Wcj<1-σ的所有值Wcj,记为第三个子样本Wdj,进入步骤308;Step 307, the record satisfies And all values W cj of 1-2σ<W cj <1-σ are recorded as the third sub-sample W dj , and enter step 308;

步骤308,计算第三个子样本Wdj的平均值,记为第三个平均值其中,n为Wdj的个数,进入步骤309;Step 308, calculate the average value of the third sub-sample W dj , recorded as the third average value which is Wherein, n is the number of Wdj, enter step 309;

步骤309,第三个子样本Wdj<第三个平均值且1-3σ<第三个子样本Wdj<1-2σ,则进入步骤310,否则判断为中级段D,停止判断;Step 309, the third sub-sample W dj < the third average value And 1-3σ<the third sub-sample W dj <1-2σ, then go to step 310, otherwise it is judged as intermediate segment D, stop judging;

步骤310,记录满足的所有值Wdj,记为第四个子样本Wej,进入步骤311;Step 310, the record satisfies All values W dj of are recorded as the fourth sub-sample W ej , and enter step 311;

步骤311,计算第四个子样本的平均值,记为第四个平均值其中,p为Wej的个数,进入步骤312;Step 311, calculate the average value of the fourth sub-sample, which is recorded as the fourth average value which is Wherein, p is the number of Wej, enter step 312;

步骤312,第四个子样本Wej<四个平均值且第四个子样本Wej<1-3σ,则判断为优级段F,停止判断,否则判断为良级段E,停止判断;Step 312, the fourth sub-sample W ej < four mean values And if the fourth sub-sample W ej <1-3σ, then it is judged to be an excellent grade segment F, and the judgment is stopped; otherwise, it is judged to be a good grade segment E, and the judgment is stopped;

其中,步骤301、步骤302和步骤303为第一次判断;步骤304、步骤305和步骤306为第二次判断;步骤307、步骤308和步骤309为第三次判断;步骤310、步骤311和步骤312为第四次判断。Wherein, step 301, step 302 and step 303 are judged for the first time; Step 304, step 305 and step 306 are judged for the second time; Step 307, step 308 and step 309 are judged for the third time; step 310, step 311 and Step 312 is the fourth judgment.

图4是本发明的储层物性识别评价装置的示意框图。如图4所示,一种基于综合录井参数的随钻储层物性识别评价装置包括:综合录井工程资料采集装置401,综合录井工程资料采集装置分别与钻头做功能耗Ei计算装置402和钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算装置403连接;钻头做功能耗Ei计算装置402与钻头做功能耗趋势值计算装置404连接,钻头做功能耗趋势值计算装置404与钻头做功能耗比值Wbi计算装置405连接,钻头做功能耗比值Wbi计算装置405与钻头做功能耗比值整体样本标准差σ计算装置406连接,钻头做功能耗比值整体样本标准差σ计算装置406与储层物性分级装置408连接;Fig. 4 is a schematic block diagram of the reservoir physical property identification and evaluation device of the present invention. As shown in Figure 4, a device for identifying and evaluating physical properties of reservoirs while drilling based on comprehensive mud logging parameters includes: a comprehensive mud logging engineering data acquisition device 401, and the comprehensive mud logging engineering data acquisition device is respectively used as a calculation device for the energy consumption E i of the drill bit 402 is connected with the calculation device 403 of the vertical drill power consumption W H and the tangential drill power consumption W L when the drill bit is working; the drill power consumption E i calculation device 402 is connected with the drill power consumption trend value The calculation device 404 is connected, and the trend value of the power consumption of the drill bit is The calculation device 404 is connected with the calculation device 405 for the power consumption ratio Wbi of the drill bit, and the calculation device 405 for the power consumption ratio W bi of the drill bit is connected with the standard deviation σ calculation device 406 for the overall sample standard deviation of the power consumption ratio of the drill bit, and the overall sample standard for the power consumption ratio of the drill bit is The difference σ calculation device 406 is connected with the reservoir physical property classification device 408;

钻头做功能耗比值Wbi计算装置405与储层物性分级装置408连接;钻头做功时垂向钻头做功能耗WH和切向钻头做功能耗WL计算装置403与交汇面积S计算与储层孔隙度大小判断装置407连接。The calculation device 405 for the power consumption ratio W bi of the drill bit is connected to the reservoir physical property classification device 408; when the drill bit works, the power consumption W H of the vertical drill bit and the power consumption W L of the tangential drill bit are calculated with the calculation device 403 and the intersection area S calculation and reservoir The porosity judgment device 407 is connected.

综合录井工程资料采集装置401采集钻头每进尺0.1m的综合录井工程资料,包括钻压Pi、钻头扭矩T、钻头转速n、机械钻速v、钻头直径dB和钻时Z。The comprehensive mud logging engineering data acquisition device 401 collects comprehensive mud logging engineering data per 0.1m of drill bit, including bit pressure P i , bit torque T, bit rotational speed n, ROP v, bit diameter d B and drilling time Z.

钻头做功能耗Ei计算装置402,用于计算每0.1m的钻头做功能耗EiThe drilling power consumption E i calculation device 402 is used to calculate the drilling power consumption E i per 0.1 m.

钻头做功能耗趋势值计算装置404用于计算m个所述钻头做功能耗Ei平均值得到钻头做功能耗趋势值m∈[3,10],j=1,2,3,4。Trend value of power consumption of drill bits The calculation device 404 is used to calculate the average value of the power consumption E i of the m drill bits to obtain the trend value of the power consumption of the drill bits m∈[3,10],j=1,2,3,4.

所述钻头做功能耗比值Wbi计算装置405用于计算所述钻头做功能耗Ei和所述钻头做功能耗趋势值的比值,得到钻头做功能耗比值WbThe working power consumption ratio W bi calculation device 405 of the drill bit is used to calculate the working power consumption E i of the drill bit and the trend value of the working power consumption of the drill bit The ratio of the bit to get the power consumption ratio of the drill bit W b .

交汇面积S计算与储层孔隙度大小判断装置407用于计算交汇面积S,根据所述交汇面积S判断储层孔隙度大小,其中所述交汇面积S通过利用所述垂向钻头做功能耗WH和所述切向钻头做功能耗WL对钻头深度h的积分得出。The intersecting area S calculation and reservoir porosity judgment device 407 is used to calculate the intersecting area S, and judge the porosity of the reservoir according to the intersecting area S, wherein the intersecting area S uses the vertical drill bit to perform power consumption W H and the work power consumption W L of the tangential drill bit are integral to the drill bit depth h.

进一步地,交汇面积S>0时,为正向交汇,说明储层孔隙度小,交汇面积S≤0时,为负向交汇,储层孔隙度大;正向交汇时,交汇面积S越大表示储层孔隙越不发育,物性越差;负向交汇时,交汇面积S越大表示储层孔隙越发育,物性越好。Furthermore, when the intersection area S>0, it is positive intersection, indicating that the porosity of the reservoir is small; when the intersection area S≤0, it is negative intersection, and the porosity of the reservoir is large; when the intersection area is positive, the larger the intersection area S is It means that the less developed the reservoir pores, the worse the physical properties; when the negative intersection, the larger the intersection area S means the more developed the reservoir pores, the better the physical properties.

所述储层物性分级装置408是利用所述钻头做功能耗比值Wbi计算装置405的钻头做功能耗比值Wbi均值及Wbi中子样本的均值和钻头做功能耗比值整体样本标准差σ计算装置406计算的样本标准差σ对储层进行物性分级,对储层进行物性分级的判断方式如图5所示。The reservoir physical property grading device 408 utilizes the drill bit power consumption ratio W bi average value of the drill bit power consumption ratio W bi calculation device 405, the mean value of W bi neutron samples, and the drill bit power consumption ratio overall sample standard deviation σ The sample standard deviation σ calculated by the calculation device 406 classifies the physical properties of the reservoir, and the judgment method for classifying the physical properties of the reservoir is shown in FIG. 5 .

图5是本发明判断储层进行物性分级装置的流程示意图,如图5所示:Fig. 5 is a flow diagram of the present invention judging the physical property grading device of the reservoir, as shown in Fig. 5:

第一次判断装置:The first judgment device:

记所述钻头做功能耗比值Wbi为第一个样本Wbj,计算所述第一个样本Wbj的平均值,记为第一个平均值 Record the working power consumption ratio Wbi of the drill bit as the first sample Wbj , calculate the average value of the first sample Wbj , and record it as the first average value

所述且Wbj<1-σ,则判断为储层段,满足所述且Wbj<1-σ的Wbj记为第二个子样本Wcj,进入第二次判断装置,否则为非储层段,停止判断;said And W bj <1-σ, then it is judged as a reservoir section, which satisfies the And W bj with W bj <1-σ is recorded as the second sub-sample W cj , and enters the second judgment device, otherwise it is a non-reservoir section, and the judgment is stopped;

第二次判断装置:The second judgment device:

计算所述第二个子样本Wcj的平均值,记为第二个平均值 Calculate the average value of the second sub-sample W cj , denoted as the second average value

所述且1-2σ<Wcj<1-σ,满足且1-2σ<Wcj<1-σ的所有值Wcj记为第三个子样本Wdj,进入第三次判断装置,否则判断为次级段,停止判断;said And 1-2σ<W cj <1-σ, satisfy And all values W cj of 1-2σ<W cj <1-σ are recorded as the third sub-sample W dj , and enter the third judgment device, otherwise it is judged as a secondary stage, and the judgment is stopped;

第三次判断装置:The third judgment device:

计算第三个子样本Wdj的平均值,记为第三个平均值 Calculate the average value of the third sub-sample W dj , denoted as the third average

所述且1-3σ<Wdj<1-2σ,满足所述且1-3σ<Wdj<1-2σ的所有值Wdj记为第四个子样本Wej,进入第四次判断装置,否则判断为中级段,停止判断;said And 1-3σ<W dj <1-2σ, satisfying the And all values W dj of 1-3σ<W dj <1-2σ are recorded as the fourth sub-sample W ej , and enter the fourth judgment device, otherwise it is judged as an intermediate stage, and the judgment is stopped;

第四次判断装置:The fourth judging device:

计算第四个子样本Wej的平均值,记为第四个平均值 Calculate the average value of the fourth sub-sample W ej , denoted as the fourth average

所述且Wej<1-3σ,则判断为优级段,停止判断,否则判断为良级段,停止判断。said And if W ej <1-3σ, it is judged as an excellent segment, and the judgment is stopped; otherwise, it is judged as a good segment, and the judgment is stopped.

显然,本领域的技术人员应该明白,上述的本发明的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本发明不限制于任何特定的硬件和软件结合。Obviously, those skilled in the art should understand that each module or each step of the above-mentioned present invention can be realized by a general-purpose computing device, and they can be concentrated on a single computing device, or distributed in a network formed by multiple computing devices Optionally, they can be implemented with program codes executable by a computing device, so that they can be stored in a storage device and executed by a computing device, or they can be made into individual integrated circuit modules, or they can be integrated into Multiple modules or steps are fabricated into a single integrated circuit module to realize. As such, the present invention is not limited to any specific combination of hardware and software.

以上所述实施例仅为表达本发明的实施方式,其描述较为具体和详细,但并不能因此而理解为对发明专利范围的限制。应当指出的是,对本领域的普通技术人员来说,在不脱离本发明的前提下,还可以做出若干变形、同等替换、改进等,这些都属于本发明的保护范围。因此,本发明专利的保护范围应以所附权利要求为准。The above-mentioned examples are only to express the implementation of the present invention, and the description thereof is relatively specific and detailed, but should not be construed as limiting the patent scope of the invention. It should be noted that those skilled in the art can make several modifications, equivalent replacements, improvements, etc. without departing from the present invention, and these all belong to the protection scope of the present invention. Therefore, the protection scope of the patent for the present invention should be based on the appended claims.

Claims (10)

1. A method for identifying and evaluating physical properties of a reservoir while drilling based on comprehensive logging parameters is characterized by comprising the following steps:
firstly, collecting comprehensive logging engineering data;
the comprehensive logging engineering data is used for acquiring the bit pressure P of 0.1m per advancing rule of the drill bit through a comprehensive logging instrumentiBit torque T, bit speed n, rate of penetration v, bit diameter dBAnd time-of-drilling Z data;
secondly, calculating the working energy consumption E of the drill biti
The drill bitWork application energy consumption EiCalculating the working energy consumption E of the drill bit of every 0.1m through a classical mechanical specific energy modeli
Thirdly, calculating the working energy consumption W of the vertical drill when the drill worksHAnd the working energy consumption W of the tangential drill bitL
Fourthly, calculating the trend value of the work-doing energy consumption of the drill bit
Calculating the trend value of the work-doing energy consumption of the drill bitM said drill bits do work and consume energy EiAverage value, m ∈ [3,10 ]],j=1,2,3,4,…;
Fifthly, calculating the working energy consumption ratio W of the drill bitbi
The working energy consumption ratio W of the drill bitbiEnergy consumption E for applying work to the drill bitiAnd the trend value of the work-doing energy consumption of the drill bitRatio of (A) to (B), Wbi∈(0,1);
Sixthly, calculating the working energy consumption ratio W of the drill bitbiThe whole sample standard deviation σ;
step seven, calculating the intersection area S and judging the porosity of the reservoir;
the intersection area S is obtained by utilizing the vertical drill bit to do work and consuming energy WHAnd the working energy consumption W of the tangential bitLIntegration of the bit depth h yields:
S = &Integral; h 1 h 2 ( W H - W L ) d h ,
wherein h is the depth of the drill bit, h2>h1
Eighthly, grading the physical properties of the reservoir;
the physical property grading judgment mode of the reservoir is as follows:
judging for the first time:
recording the working energy consumption ratio W of the drill bitbiIs the first sample WbjCalculating said first sample WbjIs taken as the first average value
The above-mentionedAnd WbjIf the sum is less than 1-sigma, the reservoir section is judged to be satisfiedAnd WbjW < 1-sigmabjIs recorded as a second subsample WcjEntering a second judgment, otherwise, judging as a non-reservoir section, and stopping the judgment;
and (4) judging for the second time:
calculating said second subsample WcjIs taken as the second average value
The above-mentionedAnd 1-2 sigma<Wcj<1-sigma, satisfyAnd 1-2 sigma<Wcj<All values W of 1-sigmacjIs recorded as the thirdSub-sample WdjEntering a third judgment, otherwise, judging as a secondary section, and stopping the judgment;
and (3) judging for the third time:
calculate the third subsample WdjIs taken as the third average value
The above-mentionedAnd 1-3 sigma<Wdj<1-2 σ, satisfy theAnd 1-3 sigma<Wdj<All values W of 1-2 sigmadjIs recorded as the fourth subsample WejEntering the fourth judgment, otherwise, judging the section as a middle stage, and stopping the judgment;
fourth judgment:
calculate the fourth subsample WejIs taken as the fourth average value
The above-mentionedAnd WejIf the value is less than 1-3 sigma, judging as a superior segment and stopping judging, otherwise, judging as a superior segment and stopping judging.
2. The method for identifying and evaluating while-drilling reservoir physical properties based on comprehensive logging parameters as recited in claim 1, wherein the classical mechanical specific energy model is:
E i = 4 P i &pi;d B 2 + 480 n T d B 2 &upsi; , i = 1 , 2 , 3 ... .
3. the method for identifying and evaluating physical properties of reservoir while drilling based on comprehensive logging parameters as recited in claim 1 or 2, wherein the energy consumption W for vertical drill bit to do work when the drill bit does work is WHAnd the working energy consumption W of the tangential drill bitLThe calculation formulas are respectively;
W H = ( W i + W i W i a &times; n b 3 ) &times; n &times; Z ;
W L = &pi; &times; d B 4 &times; T &times; n &times; Z ;
wherein a is the stratum empirical data of the maximum bit pressure of the same drill bit; b is the formation experience data of the maximum rotation speed of the same drill bit.
4. The method for identifying and evaluating physical properties of reservoir while drilling based on comprehensive logging parameters as claimed in claim 3, wherein the energy consumption E of work done by the outgoing drill bitiFiltering by adopting a five-point bell-shaped method, wherein the filtered processing formula is Ei=β(Ei-2+Ei+2)+γ(Ei-1+Ei+1)+Ei
Wherein the beta value is 0.11, the gamma value is 0.24, and the gamma value is 0.3.
5. The method for identifying and evaluating physical properties of reservoir while drilling based on comprehensive logging parameters as claimed in any one of claims 1 to 4, wherein the trend value of energy consumption of bit doing work is calculatedThe method (1) adopts a sliding mean method;
wherein,
6. the method for identifying and evaluating while-drilling reservoir physical properties based on comprehensive logging parameters as recited in claim 5, wherein the value of m is 5.
7. A device for identifying and evaluating physical properties of a reservoir while drilling based on comprehensive logging parameters is characterized by comprising:
a comprehensive logging engineering data acquisition device (401) which is respectively connected with the working energy consumption E of the drill bitiCalculating device (402) and energy consumption W for doing work of vertical drill bit when drill bit does workHAnd the working energy consumption W of the tangential drill bitLA computing device (403) connected; the working energy consumption E of the drill bitiCalculating a trend value of energy consumption of the drill bit and the computing device (402)The computing device (404) is connected, and the trend value of the work energy consumption of the drill bitCalculating a ratio W of work done by the device (404) to energy consumed by the drill bitbiA calculating device (405) is connected, and the ratio W of the work energy consumption of the drill bitbiThe calculating device (405) is connected with the drill bit work energy consumption ratio integral sample standard deviation sigma calculating device (406), and the drill bit work energy consumption ratio integral sample standard deviation sigma calculating device (406) is connected with the reservoir physical property grading device (408); energy consumption W for vertical drill bit doing work when drill bit does workHAnd the working energy consumption W of the tangential drill bitLThe calculation device (403) is connected with the intersection area S calculation and reservoir porosity judgment device (407);
the comprehensive logging engineering data acquisition device (401) acquires comprehensive logging engineering data of 0.1m per advancing rule of a drill bit, including bit pressure PiBit torque T, bit speed n, rate of penetration v, bit diameter dBAnd time of drilling Z;
the working energy consumption E of the drill bitiCalculating means (402) for calculating the energy consumption E per 0.1m of bit worki
The trend value of the energy consumption of the drill bit doing workThe computing device (404) is used for computing the work energy consumption E of the m drill bitsiThe average value obtains the trend value of the energy consumption of the drill bit doing workm∈[3,10],j=1,2,3,4;
The working energy consumption ratio W of the drill bitbiA calculation device (405) for calculating the energy consumption E of the bit workiAnd the trend value of the work-doing energy consumption of the drill bitTo obtain the ratio W of the work-doing energy consumption of the drill bitbi
The drill bit work energy consumption ratio integral sample standard deviation sigma calculating device (406) is used for calculating the drill bit work energy consumption ratio WbiThe standard deviation σ of (a);
the intersection area S calculation and reservoir porosity judgment device (407) is used for calculating an intersection area S and judging the reservoir porosity according to the intersection area S, wherein the intersection area S is obtained by utilizing the vertical drill bit to do work and consuming energy WHAnd the working energy consumption W of the tangential bitLObtaining the integral of the depth h of the drill bit;
the determination method of the reservoir physical property grading device (408) for grading the physical property of the reservoir is as follows:
first judging means:
recording the working energy consumption ratio W of the drill bitbiIs the first sample WbjCalculating said first sample WbjIs taken as the first average value
The above-mentionedAnd WbjIf the sum is less than 1-sigma, the reservoir section is judged to be satisfiedAnd WbjW < 1-sigmabjIs recorded as a second subsample WcjEntering a second judging device, otherwise, stopping judging for a non-reservoir section;
the second judging device:
calculating said second subsample WcjIs taken as the second average value
The above-mentionedAnd 1-2 sigma<Wcj<1-sigma, satisfyAnd 1-2 sigma<Wcj<All values W of 1-sigmacjIs recorded as a third subsample WdjEntering a third judging device, otherwise judging as a secondary section, and stopping judging;
the third judging device:
calculate the third subsample WdjIs taken as the third average value
The above-mentionedAnd 1-3 sigma<Wdj<1-2 σ, satisfy theAnd 1-3 sigma<Wdj<All values W of 1-2 sigmadjIs recorded as the fourth subsample WejEntering a fourth judging device, otherwise, judging the section as a middle stage, and stopping judging;
fourth judging means:
calculate the fourth subsample WejIs taken as the fourth average value
The above-mentionedAnd WejIf the value is less than 1-3 sigma, judging as a superior segment and stopping judging, otherwise, judging as a superior segment and stopping judging.
8. The device for identifying and evaluating physical properties of reservoir while drilling based on comprehensive logging parameters as claimed in claim 7, wherein the energy consumption E of the drill bit for doing work isiThe calculating device (402) adopts a classical mechanical specific energy model for calculation;
wherein the classical mechanical specific energy model is:
E i = 4 P i &pi;d B 2 + 480 n T d B 2 &upsi; , i = 1 , 2 , 3 ... .
9. the device for identifying and evaluating physical properties of reservoir while drilling based on comprehensive logging parameters as recited in claim 7 or 8, wherein the energy consumption W for vertical drill bit to do work when the drill bit does work is WHAnd the working energy consumption W of the tangential drill bitLA calculation formula score of the calculation means (403)Respectively, the following steps:
W H = ( W i + W i W i a &times; n b 3 ) &times; n &times; Z ; W L = &pi; &times; d B 4 &times; T &times; n &times; Z ;
wherein a is the stratum empirical data of the maximum bit pressure of the same drill bit; b is the formation experience data of the maximum rotation speed of the same drill bit.
10. The device for identifying and evaluating reservoir physical properties while drilling based on comprehensive logging parameters as claimed in claim 9, wherein:
the working energy consumption E of the drill bitiThe computing device (402) adopts a five-point bell-shaped filtering process, and the filtered processing formula is Ei=β(Ei-2+Ei+2)+γ(Ei-1+Ei+1)+Ei
Wherein the beta value is 0.11, the gamma value is 0.24, and the gamma value is 0.3;
the working energy consumption trend point of the drill bitThe calculating device (404) adopts a sliding mean method to obtain;
wherein,
CN201611126856.9A 2016-12-09 2016-12-09 Based on comprehensive logging parameters with bore reservoir properties identification and evaluation method and apparatus Active CN106599442B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201611126856.9A CN106599442B (en) 2016-12-09 2016-12-09 Based on comprehensive logging parameters with bore reservoir properties identification and evaluation method and apparatus

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201611126856.9A CN106599442B (en) 2016-12-09 2016-12-09 Based on comprehensive logging parameters with bore reservoir properties identification and evaluation method and apparatus

Publications (2)

Publication Number Publication Date
CN106599442A true CN106599442A (en) 2017-04-26
CN106599442B CN106599442B (en) 2019-10-29

Family

ID=58598757

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201611126856.9A Active CN106599442B (en) 2016-12-09 2016-12-09 Based on comprehensive logging parameters with bore reservoir properties identification and evaluation method and apparatus

Country Status (1)

Country Link
CN (1) CN106599442B (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108533254A (en) * 2018-03-20 2018-09-14 盘锦中录油气技术服务有限公司 Reservoir hydrocarbons water layer logging parameters Dominated Factors weight analysis method
CN112031749A (en) * 2019-05-16 2020-12-04 中国石油集团工程技术研究院有限公司 Comprehensive performance evaluation method of drill bit for oil and gas drilling
CN113638689A (en) * 2021-07-26 2021-11-12 山东大学 A quantitative drilling device and method
CN113863847A (en) * 2021-09-17 2021-12-31 山东大学 Constant-speed drilling device and stratum evaluation method
CN116861610A (en) * 2022-03-24 2023-10-10 大庆石油管理局有限公司 Method for determining porosity while drilling based on comprehensive logging parameters
CN117386359A (en) * 2023-12-13 2024-01-12 中国石油集团川庆钻探工程有限公司 Fracture-cavity type carbonate reservoir porosity while-drilling quantitative evaluation method

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1341803A (en) * 2000-08-28 2002-03-27 霍利贝顿能源服务公司 Method for predicting drilling system performance for given formation and its system
US20040238219A1 (en) * 2003-05-30 2004-12-02 Nichols Richard A. Drilling string torsional energy control assembly and method
CN104153768A (en) * 2014-07-04 2014-11-19 中国海洋石油总公司 Granite reservoir stratum reservoir performance evaluation method
CN104504182A (en) * 2014-12-09 2015-04-08 中国石油集团川庆钻探工程有限公司 Method for rapidly identifying unconventional gas reservoir
WO2016066177A2 (en) * 2014-10-30 2016-05-06 El-Neiri Mohamed Hussien Mohamed Dynamic underbalanced drilling technique

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1341803A (en) * 2000-08-28 2002-03-27 霍利贝顿能源服务公司 Method for predicting drilling system performance for given formation and its system
US20040238219A1 (en) * 2003-05-30 2004-12-02 Nichols Richard A. Drilling string torsional energy control assembly and method
CN104153768A (en) * 2014-07-04 2014-11-19 中国海洋石油总公司 Granite reservoir stratum reservoir performance evaluation method
WO2016066177A2 (en) * 2014-10-30 2016-05-06 El-Neiri Mohamed Hussien Mohamed Dynamic underbalanced drilling technique
CN104504182A (en) * 2014-12-09 2015-04-08 中国石油集团川庆钻探工程有限公司 Method for rapidly identifying unconventional gas reservoir

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108533254A (en) * 2018-03-20 2018-09-14 盘锦中录油气技术服务有限公司 Reservoir hydrocarbons water layer logging parameters Dominated Factors weight analysis method
CN112031749A (en) * 2019-05-16 2020-12-04 中国石油集团工程技术研究院有限公司 Comprehensive performance evaluation method of drill bit for oil and gas drilling
CN113638689A (en) * 2021-07-26 2021-11-12 山东大学 A quantitative drilling device and method
CN113863847A (en) * 2021-09-17 2021-12-31 山东大学 Constant-speed drilling device and stratum evaluation method
CN116861610A (en) * 2022-03-24 2023-10-10 大庆石油管理局有限公司 Method for determining porosity while drilling based on comprehensive logging parameters
CN117386359A (en) * 2023-12-13 2024-01-12 中国石油集团川庆钻探工程有限公司 Fracture-cavity type carbonate reservoir porosity while-drilling quantitative evaluation method
CN117386359B (en) * 2023-12-13 2024-03-08 中国石油集团川庆钻探工程有限公司 Fracture-cavity type carbonate reservoir porosity while-drilling quantitative evaluation method

Also Published As

Publication number Publication date
CN106599442B (en) 2019-10-29

Similar Documents

Publication Publication Date Title
CN106599442B (en) Based on comprehensive logging parameters with bore reservoir properties identification and evaluation method and apparatus
Pittman Relationship of porosity and permeability to various parameters derived from mercury injection-capillary pressure curves for sandstone
CN104695950B (en) Volcanic Reservoir PRODUCTION FORECASTING METHODS
CN110114552A (en) For determining the method and system of well drilling detritus depth
CN104806232B (en) A kind of method for determining porosity lower limit of fracture
CN104350233B (en) It is evaluated using the formation volume of the differential data of standardization
CN110130883A (en) Method and device for determining rock formation parameters
CN110397402B (en) Drilling method and device
CN112664174B (en) Tunnel surrounding rock grade determination method and system based on multiple drill holes
CN105913332B (en) Method and device for determining dense gas reservoir development index
CN104500054B (en) The determination method and device of formation pore pressure
CN103180548A (en) System and method for hydrocarbon gas pay zone characterization in a subterranean reservoir
CN106223941B (en) Method is determined based on the shale gas reservoir organic porosity of well-log information
CN104412130A (en) Formation compositional evaluation using normalized differential data
CN110596757A (en) A Correction Method of P-wave and S-wave Velocity in Shale Formation
CN106157181A (en) A kind of based on boring the parameter method to soil property country rock quick real-time grading
CN107143330A (en) Shale gas reservoir quality surveys mud logging evaluation method
CN104948176A (en) Method for identifying carbonate reservoir fractures based on permeability increasing rate
CN105257252A (en) Method for optimizing shale gas horizontal well clustering perforation well section by using logging information
CN105204081B (en) A kind of method predicting shale gas Clay Mineral and constituent content thereof
CN109555515A (en) Formation collapsed pressure determines method and apparatus
CN104153768A (en) Granite reservoir stratum reservoir performance evaluation method
CN114818451A (en) Mechanical drilling rate prediction method, device, storage medium and equipment
CN112381938B (en) Stratum identification method based on trenchless parameter while drilling machine learning
CN113073959B (en) Drilling well loss mechanism prediction and auxiliary decision method based on fuzzy decision tree

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: 20250125

Address after: No. 12 Aiguo Road, Ranghu Road District, Daqing City, Heilongjiang Province, China 163453

Patentee after: Daqing Drilling Engineering Co.,Ltd.

Country or region after: China

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Address before: 100007 Oil Mansion, 9 North Avenue, Dongcheng District, Beijing, Dongzhimen

Patentee before: China National Petroleum Corp.

Country or region before: China

Patentee before: DA QING PETROLEUM ADMINISTRATION BUREAU

TR01 Transfer of patent right