CN105931125A - Tight oil segmented multi-cluster volume fractured horizontal well production prediction method - Google Patents

Abstract

The invention provides a tight oil segmented multi-cluster volume fractured horizontal well production prediction method. According to the method, the change of horizontal segment transverse oiliness is reflected by the oil reservoir transverse heterogeneous coefficient based on the analysis of tight oil field development test production affecting factors, the size of oiliness is reflected by horizontal segment along oil reservoir standard gas hydrocarbon value envelope surface area acquired through logging, the change of reservoir permeability is reflected by logging interpretation permeability, and the artificial fracture development degree is comprehensively reflected by the logging brittleness index and the artificial fracture liquid quantity into the ground so that the objectives that great prediction of tight oil segmented multi-cluster volume fractured horizontal well production can be realized, the investment risk of the tight oil scale development can be reduced and the production building benefit can be enhanced can be realized.

Description

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation

Technical field

The invention belongs to oil field development technical field, be specifically related to a kind of fine and close oil many bunches of volume fracturing horizontal wells of segmentation and produce Amount Forecasting Methodology.

Background technology

Fine and close oil refers to be clipped in or be close in the compact reservoir of high quality source rock series of strata, without extensive long-distance migration And the oil-gas accumulation formed, the most pressure break is not without natural production capacity, and the many bunches of volume fracturing exploitations of current horizontal well in segments have become cause The Main Means of close oil exploitation；Fine and close oil level well development is single with the most significant difference of low-permeability oil deposit horizontal well development in the past It is longer that the product cycle is invested greatly and built to well, for reducing investment risk in fine and close oily scale development, reaches raising and produces the target building benefit, Many bunches of volume fracturing horizontal well production predictions of segmentation have become as the hot issue of fine and close oil exploitation, are also that a difficult point is asked simultaneously Topic, is mainly reflected in following 2 aspects.

1, the fine and close oil horizontal oiliness of reservoir changes greatly, and the change of horizontal well production and horizontal oiliness has pass closely System, existing Horizontal Well Log Interpretation model and method of testing are difficult to reflect the change of horizontal oiliness anisotropism.This is Because horizontal well is different from straight well, the nearly horizontal penetrating ground of horizontal well, and the medium of wellbore is the most symmetrical, is subject to simultaneously The impact of gravity factor, instrument well logging state is typically eccentric, and the eccentric measurement to various loggers all has in various degree Impact, most notable of which be transformation scale quite in the case of, the dependency of well log interpretation oil saturation and yield Poor, need the influence factor that screening is sensitive to yield.

2, compact reservoir rock brittleness index is relatively big, and reservoir reconstruction is larger, and the fracture network of formation is more complicated, retouches The technology of pressure-break length, width and the height of stating artificial seam net is the most immature, current conventional man-made fracture detection technique-- The result of down-hole micro-seismic monitoring method test is only microseism signal parameter, it is impossible to represent real fracture pattern.

From the point of view of literature survey frac water horizontal well production prediction method situation, existing Forecasting Methodology can be summarized as two classes: one Class is method for numerical simulation, and another kind of is economics analysis model prediction method.But whether method for numerical simulation, or theoretical Analytical model Forecasting Methodology, when carrying out production forecast, due to many bunches of volume fracturing horizontal well man-made fracture ginsengs of fine and close oil segmentation Number description technique is the most immature, and the well log interpretation oil saturation of existence routine can not reflect the horizontal oiliness of horizontal well The problem of change, causes existing method precision when filed application poor.Need from the fine and close many bunches of bodies of profit horizontal well segmentation of impact Overstocking and split the oiliness of yield, permeability and three aspects of development degree of micro cracks in oil, screening becomes with fine and close profit horizontal well actual production Change sensitive parameter.

Summary of the invention

It is an object of the invention to overcome method for numerical simulation and economics analysis method to be difficult to preferably reflect that fine and close oil reservoir is horizontal Change greatly to oiliness and describe difficult problem with man-made fracture development degree.

To this end, the invention provides a kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation, including with Lower step:

Step 1) after horizontal well drilling terminates, obtain oil reservoir horizontal coefficient of heterogeneity HE_f；

Step 2) according to horizontal segment in well logging result along the permeability of journey oil reservoir, obtain horizontal segment along the well logging infiltration of journey oil reservoir Rate enveloping surface area K'；

Step 3) according to well logging result real-time during horizontal well drilling, it is eliminated when boring and entrance mud flow rate Standard total hydrocarbon in gas logging value Q_T, application definite integral obtains net horizontal section along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T'；

Step 4) with Barnett shale brittleness index formula, by compressional wave time difference, shear wave slowness and rock volume density value Calculate rock brittleness index；

Step 5) according to calculated rock brittleness index, apply definite integral principle, obtain horizontal segment crisp along journey oil reservoir Sex index enveloping surface area B I'；

Step 6) according to during many bunches of volume fracturings of segmentation every section enter ground liquid measure add up read group total individual well always enter ground liquid Amount L_v；

Step 7) obtain many bunches of volume fracturing horizontal well production Q predictor formulas of segmentation:

Q=1.628*Ln (HE_f×K'×Q_T'×BI'×L_V)-50.717。

Step 1) described oil reservoir horizontal coefficient of heterogeneity HE_fHorizontal segment oil reservoir total length L is met for boring_oWith discontinuous oil layer section The ratio of number N:

${\mathrm{HE}}_f=\frac{L_o}{N}=\frac{\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{L}_i}{N}$

In formula, HE_fFor the horizontal coefficient of heterogeneity of oil reservoir, m/ section；L_oHorizontal segment oil reservoir total length, m is met for boring；N meets water for boring In flat section, oil layer section number, individual；L_iI-th oil layer section length in horizontal segment, m is met for boring.

Step 2) horizontal segment along journey oil reservoir well logging permeability enveloping surface area K' be calculated by following formula:

$K^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{k}_{ij}\×\mathrm{\Δ}l)$

Wherein, K' is horizontal segment to log well permeability enveloping surface area along journey oil reservoir, mD m；

k_ijFor the permeability of jth sample point, mD in i-th oil layer section in brill chance horizontal segment；

Δ l is well logging sampling interval, m.

Step 3) net horizontal section is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T' it is calculated by following formula:

${Q_T}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{Q}_{Tij}\×\mathrm{\Δ}l)$

$Q_{Tij}=\frac{2}{\mathrm{\π}}\×100\×arctg\[\frac{{\mathrm{\ν}}_{ij}{t}_{ij}}{{\mathrm{\ν}}_o{t}_o}tg(\frac{\mathrm{\π}}{2}{q}_{tij}\×\frac{1}{100})\]$

Wherein, Q_T' for horizontal segment along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area, % m；Q_TijFor boring chance level The standard total hydrocarbon in gas logging value that in Duan, i-th oil layer section jth sample point is corresponding, %；

Step 4) rock brittleness index is calculated by following formula:

${\mathrm{BI}}_{ij}=\frac{{\mathrm{\ΔE}}_{ij}+{\mathrm{\Δ\μ}}_{ij}}{2}\×100$

Wherein,

$E_{ij}=\frac{{\mathrm{\ρ}}_{bij}}{{{\mathrm{\Δt}}_{sij}}^2}\×\left(\frac{3{{\mathrm{\Δt}}_{sij}}^2-4{{\mathrm{\Δt}}_{p_{ij}}}^2}{{{\mathrm{\Δt}}_{sij}}^2-{{\mathrm{\Δt}}_{p_{ij}}}^2}\right)\×{10}^5,{\mathrm{\μ}}_{ij}=\frac{{{\mathrm{\Δt}}_{sij}}^2-2{{\mathrm{\Δt}}_{p_{ij}}}^2}{2({{\mathrm{\Δt}}_{sij}}^2-{{\mathrm{\Δt}}_{p_{ij}}}^2)};$

${\mathrm{\Δt}}_{sij}=\frac{{\mathrm{\Δt}}_{p_{ij}}}{{\[1-1.15\frac{1/{\mathrm{\ρ}}_{bij}+{(1/{\mathrm{\ρ}}_{bij})}^3}{e^{1/{\mathrm{\ρ}}_{b_{ij}}}}\]}^{1.5}}$

In formula, BI_ijThe well logging brittleness index that in horizontal segment, i-th oil layer section jth sample point is corresponding, % is met for boring； ΔE_ijFor Young's modulus after the normalization that i-th oil layer section jth sample point in brill chance horizontal segment is corresponding, dimensionless；Δμ_ij For boring the Poisson's ratio after meeting the normalization that in horizontal segment, i-th oil layer section jth sample point is corresponding, dimensionless；E_ijMeet for boring The Young's modulus that in horizontal segment, i-th oil layer section jth sample point is corresponding, 10⁴MPa；μ_ijI-th oil in horizontal segment is met for boring The Poisson's ratio that interval jth sample point is corresponding, dimensionless；ρ_bijI-th oil layer section jth sampled point in horizontal segment is met for boring The rock volume density that place is corresponding, g/cm³；Δt_sijCorresponding for boring i-th oil layer section jth sample point in chance horizontal segment The shear wave slowness of rock, μ s/m；Δt_pijThe vertical of i-th oil layer section jth sample point is corresponding in horizontal segment rock is met for boring The ripple time difference, μ s/m.

Step 5) horizontal segment obtained along journey oil reservoir brittleness index enveloping surface area B I' by following formula:

${\mathrm{BI}}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{BI}}_{ij}\×\mathrm{\Δ}l).$

Step 6) described in individual well always enter ground liquid measureWherein, L_vGround liquid measure, m is always entered for individual well³；L_vi Ground liquid measure, m is entered for many bunches of volume fracturings of segmentation i-th section³。

Described standard total hydrocarbon in gas logging value Q_TBy antitrigonometric function method, well logging real-time during horizontal well drilling is obtained Total hydrocarbon in gas logging value q_tIt is corrected obtaining:

$Q_{Tij}=\frac{2}{\mathrm{\π}}\×100\×arctg\[\frac{{\mathrm{\ν}}_{ij}{t}_{ij}}{{\mathrm{\ν}}_o{t}_o}tg(\frac{\mathrm{\π}}{2}{q}_{tij}\×\frac{1}{100})\]$

In formula, ν_ijI-th oil layer section jth sample point actual well drilled flow quantity in horizontal segment, m is met for boring³/min；t_ij When meeting the actual brill of i-th oil layer section jth sample point in horizontal segment for boring, min/m；ν_oFor block standard well liquid entrance, m³/min；t_oWhen boring for block standard, min/m；q_tijCorresponding for boring i-th oil layer section jth sample point in chance horizontal segment With boring total hydrocarbon in gas logging value, %.

The shear wave slowness Δ t of described rock_sijCalculating formula when being used for logging well only to compressional wave time difference Δ t_pijTest, do not have There is the situation of test shear wave slowness.

The invention has the beneficial effects as follows: this method that the present invention provides fine and close oil field field development experiments yield effect because of On the basis of element is analyzed, propose, with the change of the oil reservoir horizontal coefficient of heterogeneity reflection horizontal oiliness of horizontal segment, to obtain with well logging Horizontal segment along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area reflection oiliness size, with well log interpretation permeability reflect The change of reservoir permeability, to log well brittleness index and man-made fracture enters ground liquid measure concentrated expression man-made fracture development degree, from And realize the preferable prediction of fine and close oil many bunches of volume fracturing horizontal well production of segmentation, it is achieved reduce in fine and close oil scale development and invest Risk, improves and produces the target building benefit.

It is described in further details below in conjunction with accompanying drawing.

Accompanying drawing explanation

Fig. 1 is oil reservoir horizontal coefficient of heterogeneity HE_fRelation curve with the horizontal well individual well daily output；

Fig. 2 is mean permeability k and the relation curve of the horizontal well individual well daily output at individual well crack；

Fig. 3 is average well logging porosity at individual well crackRelation curve with the horizontal well individual well daily output；

Average well-log oil saturation S at Fig. 4 individual well crack_oRelation curve with the horizontal well individual well daily output；

Average rock brittleness index BI and the relation curve of the horizontal well individual well daily output at Fig. 5 individual well crack；

Average total hydrocarbon in gas logging value Q at Fig. 6 crack_TRelation curve with the horizontal well individual well daily output；

Fig. 7 is the relation curve of well logging permeability enveloping surface area K' and the horizontal well individual well daily output；

Fig. 8 is the relation curve of well logging porosity enveloping surface area Φ ' and the horizontal well individual well daily output；

Fig. 9 is well-log oil saturation enveloping surface area S_o' and the relation curve of the horizontal well individual well daily output；

Figure 10 is the relation curve of well logging brittleness index enveloping surface area B I' and the horizontal well individual well daily output；

Figure 11 is standard total hydrocarbon in gas logging value enveloping surface area Q_T' and the relation curve of the horizontal well individual well daily output；

Many bunches of volume fracturings of Figure 12 well segmentation enter ground liquid measure L_vRelation curve with the horizontal well individual well daily output；

Many bunches of volume fracturing seam hop count STA of Figure 13 well segmentation and the relation curve of the horizontal well individual well daily output；

Figure 14 is the relation curve of many bunches of average sand feeding amount SAN of volume fracturing of individual well segmentation and the horizontal well individual well daily output；

Figure 15 is the relation curve of many bunches of volume fracturing average discharge PR of individual well segmentation and the horizontal well individual well daily output；

Figure 16 is the many bunches of average sand of volume fracturing of the individual well segmentation relation curves than PC Yu the horizontal well individual well daily output；

Figure 17 is the relation curve of enveloping surface coefficient and the horizontal well individual well daily output；

Figure 18 is that GP47-64 borehole logging tool explains individual well card.

Detailed description of the invention

Embodiment 1:

Present embodiments provide a kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation, including following step Rapid:

Step 1) after horizontal well drilling terminates, obtain oil reservoir horizontal coefficient of heterogeneity HE_f；

Step 2) according to horizontal segment in well logging result along the permeability of journey oil reservoir, obtain horizontal segment along the well logging infiltration of journey oil reservoir Rate enveloping surface area K'；

Step 3) according to well logging result real-time during horizontal well drilling, it is eliminated when boring and entrance mud flow rate Standard total hydrocarbon in gas logging value Q_T, application definite integral obtains net horizontal section along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T'；

Step 4) with Barnett shale brittleness index formula, by compressional wave time difference, shear wave slowness and rock volume density value Calculate rock brittleness index；

Step 5) according to calculated rock brittleness index, apply definite integral principle, obtain horizontal segment crisp along journey oil reservoir Sex index enveloping surface area B I'；

Step 6) according to during many bunches of volume fracturings of segmentation every section enter ground liquid measure add up read group total individual well always enter ground liquid Amount L_v；

Step 7) obtain many bunches of volume fracturing horizontal well production Q predictor formulas of segmentation:

Q=1.628*Ln (HE_f×K'×Q_T'×BI'×L_V)-50.717。

The horizontal segment that the present invention obtains with well logging reflects oiliness along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Size, with the change of well log interpretation permeability reflection reservoir permeability, to log well brittleness index and man-made fracture enters ground liquid measure and combines Close reflection man-made fracture development degree, instruct and set up a kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation, The method has the advantage that parameter is easy to get and reliability is high.

Embodiment 2:

(1) fine and close oil level well production sensitive parameter screening

The influence factor of the fine and close profit many bunches of volume fracturing yield of horizontal well segmentation of impact has geologic(al) factor and engineering factor.Ground Considering in quality factor that reservoir oiliness, reservoir permeability and reservoir seam net form complexity, wherein reservoir oiliness can be used The horizontal coefficient of heterogeneity of oil reservoir, well logging total hydrocarbon in gas logging value, well-log oil saturation, 4 parameters of well logging porosity characterize, storage Layer permeability well logging permeability characterizes, and reservoir seam net forms complexity rock brittleness index and characterizes；Engineering Factor mainly affects reservoir seam net development degree after volume fracturing transformation, mainly includes into ground liquid measure, artificial fracture transformation section Number, sand feeding amount, discharge capacity and sand are than 5 parameters；

Step one: geologic(al) factor is screened

(1) oil reservoir horizontal coefficient of heterogeneity HE_f: according to well log interpretation conclusion, add up every mouthful of sample level well horizontal segment and bore Meet oil reservoir length L_oAnd horizontal segment oil layer section number N, obtain the oil reservoir horizontal coefficient of heterogeneity HE of every mouthful of sample level well_f；

${\mathrm{HE}}_f=\frac{L_o}{N}=\frac{\underset{i=1}{\overset{N}{\mathrm{\Σ}}}{L}_i}{N}---\left(1\right)$

(2) sample level well horizontal coefficient of heterogeneity HE is calculated according to formula (1)_f, do horizontal coefficient of heterogeneity with corresponding The scatterplot of horizontal well production, as it is shown in figure 1, the good relationship of regression formula and sample point, therefore by horizontal for oil reservoir non-all Matter coefficient is as the main geologic parameter of the fine and close oil level well production factor of impact；

Geologic parameter at conventional levels well Analysis On Production-affecting Fact-ors horizontal segment pressure-break represents whole level The parameter of well, for huge discharge volume fracturing densification profit horizontal well, applies same parameter statistic to fine and close profit horizontal well Influence factor is analyzed；

(3) averagely log well at individual well crack permeability k；

$k=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{k}_i}{n}---\left(2\right)$

(4) average well logging porosity at individual well crack

(5) average well-log oil saturation S at individual well crack_o；

$S_o=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{S}_{oi}}{n}---\left(4\right)$

(6) average total hydrocarbon in gas logging value Q at individual well crack_T；

$Q_T=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{Q}_{Ti}}{n}---\left(5\right)$

Total hydrocarbon in gas logging value q that well logging real-time during horizontal well drilling is obtained_t, application antitrigonometric function method is corrected Correction, is eliminated when boring and standard total hydrocarbon in gas logging value Q of entrance mud flow rate_T；

Standard total hydrocarbon in gas logging value

$Q_{Ti}=\frac{2}{\mathrm{\π}}\×100\×arctg\[\frac{{\mathrm{\ν}}_i{t}_i}{{\mathrm{\ν}}_o{t}_o}tg(\frac{\mathrm{\π}}{2}{q}_{ti}\×\frac{1}{100})\]$

(7) average rock brittleness index BI at individual well crack；

$BI=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{BI}}_i}{n}---\left(6\right)$

Rock brittleness index BI can not be directly obtained by log, by Barnett shale brittleness index computing formula Long 7 compact reservoir brittleness indexs are calculated according to compressional wave time difference, shear wave slowness and rock volume density value；

Brittleness index

Normalization Young's modulus

Normalization Poisson's ratio

Young's modulus

Poisson's ratio

The empirical equation that application is applicable to gas-free sandstone or argillaceous sandstone stratum is estimated for shear wave is calculated ground Layer shear wave slowness Δ t_s；

${\mathrm{\Δt}}_{si}=\frac{{\mathrm{\Δt}}_{pi}}{{\[1-1.15\frac{1/{\mathrm{\ρ}}_{bi}+{(1/{\mathrm{\ρ}}_{bi})}^3}{e^{1/{\mathrm{\ρ}}_{ib}}}\]}^{1.5}}$

(8) according to formula (1)～the average geologic parameter of (6) calculated level Jing Liefengchu, and each geologic parameter and water are done The scatterplot of the horizontal well individual well daily output, as shown in Fig. 2～Fig. 6, by table 1 regressive trend line relative coefficient it can be seen that various places Meansigma methods at matter parameter crack is poor with the dependency of yield, therefore analyzes and considers that whole horizontal segment is along journey reservoir heterogeneity The relation of geologic parameter enveloping surface area and yield；

(9), after terminating according to horizontal well drilling, horizontal segment reservoir permeability in well log interpretation conclusion, application definite integral is former Reason, obtain horizontal segment along journey oil reservoir log well permeability enveloping surface area K':

$K^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{k}_{ij}\×\mathrm{\Δ}l)---\left(7\right)$

(10), after terminating according to horizontal well drilling, horizontal segment oil reservoir porosity in well log interpretation conclusion, application definite integral is former Reason, obtains horizontal segment along journey oil reservoir well logging porosity enveloping surface area Φ ':

(11) after terminating according to horizontal well drilling, horizontal segment oil reservoir oil saturation in well log interpretation conclusion, apply constant volume Point principle, obtains horizontal segment along journey oil reservoir well-log oil saturation enveloping surface area S_o':

${S_o}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{S}_{oij}\×\mathrm{\Δ}l---\left(9\right)$

(12) total hydrocarbon in gas logging value q that well logging real-time during horizontal well drilling is obtained_t, application antitrigonometric function method is carried out Correction correction, is eliminated when boring and standard total hydrocarbon in gas logging value Q of entrance mud flow rate_T, apply definite integral principle, obtain water Flat section is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T':

Standard total hydrocarbon in gas logging value

$Q_{Tij}=\frac{2}{\mathrm{\π}}\×100\×arctg\[\frac{{\mathrm{\ν}}_{ij}{t}_{ij}}{{\mathrm{\ν}}_o{t}_o}tg(\frac{\mathrm{\π}}{2}{q}_{tij}\×\frac{1}{100})\]$

Horizontal segment is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area:

${Q_T}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{Q}_{Tij}\×\mathrm{\Δ}l)---\left(10\right)$

(13) rock brittleness index BI can not be directly obtained by log, calculates by Barnett shale brittleness index Formula calculates long 7 compact reservoir brittleness indexs according to compressional wave time difference, shear wave slowness and rock volume density value；

Rock brittleness index

Normalization Young's modulus

Normalization Poisson's ratio

Young's modulus

Poisson's ratio

The empirical equation that application is applicable to gas-free sandstone or argillaceous sandstone stratum is estimated for shear wave is calculated ground Layer shear wave slowness Δ t_s:

${\mathrm{\Δt}}_{sij}=\frac{{\mathrm{\Δt}}_{p_{ij}}}{{\[1-1.15\frac{1/{\mathrm{\ρ}}_{bij}+{(1/{\mathrm{\ρ}}_{bij})}^3}{e^{1/{\mathrm{\ρ}}_{b_{ij}}}}\]}^{1.5}}$

According to calculated rock brittleness index, apply definite integral principle, obtain horizontal segment along journey oil reservoir well logging index Enveloping surface area B I'；

${\mathrm{BI}}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{BI}}_{ij}\×\mathrm{\Δ}l)---\left(11\right)$

(14) the logging program sampling interval used by Changqing oilfields is 0.125m, i.e. Δ l=0.125m, therefore formula ～(11) can be to be simplified to following form (7):

$K^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{k}_{ij}\×\mathrm{\Δ}l)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{k}_{ij}\×0.125)---\left(12\right)$

${S_o}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{S}_{oij}\×\mathrm{\Δ}l)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{S}_{oij}\×0.125)---\left(14\right)$

${Q_T}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{Q}_{Tij}\×\mathrm{\Δ}l)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{Q}_{Tij}\×0.125)---\left(15\right)$

${\mathrm{BI}}^{\′}=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{BI}}_{ij}\×\mathrm{\Δ}l)=\underset{i=1}{\overset{N}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{BI}}_{ij}\×0.125)---\left(16\right)$

(15) each sample level well geologic parameter horizontal segment is calculated along journey oil reservoir enveloping surface according to formula (12)～(16) Area, does the scatterplot (Fig. 7～Figure 13) of these parameters and the individual well daily output.The dependency of Trendline is obtained according to regression formula Coefficients R²(being shown in Table 1), it can be seen that compared to the relation of meansigma methods at crack Yu the individual well daily output, each geologic parameter enveloping surface face Long-pending more preferable with the relation of the individual well daily output.Filter out oil reservoir well logging permeability enveloping surface area, the standard of good relationship simultaneously Total hydrocarbon in gas logging value enveloping surface area and brittleness index enveloping surface area are as the principal element affecting fine and close oil level well production.

Step 2: engineering factor screens

Engineering factor mainly affects reservoir seam net development degree after volume fracturing transformation, mainly includes into ground liquid measure, artificial Pressure-break hop count, sand feeding amount, discharge capacity and sand are than 5 parameters；

(16) many bunches of volume fracturings of individual well segmentation enter ground liquid measure L_v；

$L_v=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{L}_{vi}---\left(17\right)$

(17) many bunches of volume fracturing seam hop count STA of individual well segmentation；

STA=n (18)

(18) many bunches of average sand feeding amount SAN of volume fracturing of individual well segmentation；

$SAN=\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{SAN}}_i---\left(19\right)$

(19) many bunches of volume fracturing average discharge PR of individual well segmentation；

$PR=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{PR}}_i}{n}---\left(20\right)$

(20) many bunches of average sand of volume fracturing of individual well segmentation compare PC；

$PC=\frac{\underset{i=1}{\overset{n}{\mathrm{\Σ}}}{\mathrm{PC}}_i}{n}---\left(21\right)$

The each sample horizontal well Single well project parameter calculated according to formula (17)～(21), does each engineering parameter and individual well day The scatterplot (Figure 12～Figure 16) of yield, enters ground liquid measure according to the size (table 1) of relative coefficient, preferably individual well and causes as impact The principal element of close oil level well production；

(2) fine and close oil level well production formulary regression

Being analyzed the relation of each influence factor with the individual well daily output, as shown in figure 17,5 parametric regressions obtain densification Oil many bunches of volume fracturing horizontal well production predictor formulas of segmentation；

Q=1.628*Ln (HE_f×K'×Q_T'×BI'×L_V)-50.717 (22)

Q is volume fracturing horizontal well production, t/d.

The fine and close oil level well production influence factor of table 1 and the correlation coefficient of yield scatterplot

Embodiment 3:

Heshui length 7 is controlled by western, southwest and many things source, south, based on thing source, the west and south, grows delta and lake Two kinds of type of sedimentary facies；The turbidite Reservoir Body that major developmental delta front slump is formed.Turbidity channel sand is as this district Skeleton matching, in a shape, lumps distribution, sand thickness is big, and oil-layer distribution is the most stable, and thickness is relatively big, and oil reservoir has bigger Scale, resource potential is big.

GP47-64 is Heshui length 7 oil reservoir development horizontal well (Figure 18), horizontal section length 680m, bores and meets oil layer section 647m, bores and meets oil layer section number 7, uses water-jet annular space to add sand volume fracturing reforming mode, enters ground liquid measure 3550m³, oil Interval well logging permeability enveloping surface area is 143.8mD m, and total hydrocarbon in gas logging value enveloping surface area is 19137.1% m, fragility Exponential envelope face area is 34938.2% m, and the horizontal coefficient of heterogeneity of oil reservoir is 96.3m/ section.Produce according to fine and close profit horizontal well Amount predictor formula, calculates GP47-64 well individual well and has a daily output of 11.1t/d, and actual individual well has a daily output of 10.5t/d, it was predicted that error Being 5.9%, effect is preferable, and concrete calculating process is shown in GP47-64 production forecast process.

GP47-64 production forecast process:

Step one: calculate oil reservoir horizontal coefficient of heterogeneity HE_f:

Step 2: calculated level section is along journey oil reservoir well logging permeability enveloping surface area (being shown in Table 2):

$K^{\′}=\underset{i=1}{\overset{7}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{k}_{ij}\×0.125)=143.8mD\·m$

Table 2 horizontal segment is along journey oil reservoir well logging permeability enveloping surface area calculation step

Step 3: calculated level section is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area (being shown in Table 3):

${Q_T}^{\′}=\underset{i=1}{\overset{7}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{Q}_{Tij}\×0.125)=17933.3\%\·m$

Table 3 horizontal segment is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area calculation step

Step 4: calculated level section is along journey oil reservoir well logging brittleness index enveloping surface area (being shown in Table 4):

${\mathrm{BI}}^{\′}=\underset{i=1}{\overset{7}{\mathrm{\Σ}}}(\underset{j=1}{\overset{m}{\mathrm{\Σ}}}{\mathrm{BI}}_{ij}\×0.125)=34938.2\%\·m$

Table 4 is along journey oil reservoir brittleness index enveloping surface area calculation step

Step 5: always enter ground liquid measure L according to the ground liquid measure calculating individual well that enters of every section during many bunches of volume fracturings of segmentation_v:

$L_v=\underset{i=1}{\overset{8}{\mathrm{\Σ}}}{L}_{vi}=386+324+443+553+520+464+494+366=3550m^3$

Step 6: calculating GP47-64 yield:

Q=1.628*Ln (HE_f×K'×Q_T'×BI'×L_V)-50.717

=1.628 × Ln (96.3 × 143.8 × 19137.1 × 34938.2 × 3550)-50.717

=11.1t/d

GP47-64 actual production is 10.5t/d, and value of calculation and actual value error are 5.9%.

Exemplified as above is only the illustration to the present invention, is not intended that the restriction to protection scope of the present invention, all It is within design same or analogous with the present invention belongs to protection scope of the present invention.

Claims (9)

1. fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation, it is characterised in that comprise the following steps:

Step 1) after horizontal well drilling terminates, obtain oil reservoir horizontal coefficient of heterogeneity HE_f；

Step 2) according to horizontal segment in well logging result along the permeability of journey oil reservoir, obtain horizontal segment and log well permeability bag along journey oil reservoir Network face area K'；

Step 3) according to well logging result real-time during horizontal well drilling, it is eliminated when boring and the standard of entrance mud flow rate Total hydrocarbon in gas logging value Q_T, application definite integral obtains net horizontal section along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T'；

Step 4) with Barnett shale brittleness index formula, calculated by compressional wave time difference, shear wave slowness and rock volume density value Rock brittleness index；

Step 5) according to calculated rock brittleness index, apply definite integral principle, obtain horizontal segment and refer to along journey oil reservoir fragility Number enveloping surface area B I'；

Step 6) according to during many bunches of volume fracturings of segmentation every section enter ground liquid measure add up read group total individual well always enter ground liquid measure L_v；

Step 7) obtain many bunches of volume fracturing horizontal well production Q predictor formulas of segmentation:

Q=1.628*Ln (HE_f×K'×Q_T'×BI'×L_V)-50.717。

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature It is: step 1) described oil reservoir horizontal coefficient of heterogeneity HE_fHorizontal segment oil reservoir total length L is met for boring_oWith discontinuous oil layer section number The ratio of N:

In formula, HE_fFor the horizontal coefficient of heterogeneity of oil reservoir, m/ section；L_oHorizontal segment oil reservoir total length, m is met for boring；N meets horizontal segment for boring Middle oil layer section number, individual；L_iI-th oil layer section length in horizontal segment, m is met for boring.

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature Be: step 2) horizontal segment along journey oil reservoir well logging permeability enveloping surface area K' be calculated by following formula:

Wherein, K' is horizontal segment to log well permeability enveloping surface area along journey oil reservoir, mD m；

k_ijFor the permeability of jth sample point, mD in i-th oil layer section in brill chance horizontal segment；

Δ l is well logging sampling interval, m.

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature It being: step 3) net horizontal section is along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area Q_T' it is calculated by following formula:

Wherein, Q_T' for horizontal segment along journey oil reservoir standard total hydrocarbon in gas logging value enveloping surface area, % m；Q_TijBore and meet in horizontal segment i-th The standard total hydrocarbon in gas logging value that individual oil layer section jth sample point is corresponding, %, %.

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature It being: step 4) rock brittleness index is calculated by following formula:

Wherein,

In formula, BI_ijThe well logging brittleness index that in horizontal segment, i-th oil layer section jth sample point is corresponding, % is met for boring；ΔE_ij Bore and meet Young's modulus after the normalization that in horizontal segment, i-th oil layer section jth sample point is corresponding, dimensionless；Δμ_ijMeet for boring Poisson's ratio after the normalization that in horizontal segment, i-th oil layer section jth sample point is corresponding, dimensionless；E_ijBore and meet in horizontal segment The Young's modulus that i-th oil layer section jth sample point is corresponding, 10⁴MPa；μ_ijI-th oil layer section jth in horizontal segment is met for boring The Poisson's ratio that individual sample point is corresponding, dimensionless；Corresponding for boring i-th oil layer section jth sample point in chance horizontal segment Rock volume density, g/cm³；The shear wave of the rock that i-th oil layer section jth sample point is corresponding in horizontal segment is met for boring The time difference, μ s/m；The compressional wave time difference of the rock that i-th oil layer section jth sample point is corresponding, μ s/ in horizontal segment is met for boring m。

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature It being: step 5) horizontal segment obtained along journey oil reservoir brittleness index enveloping surface area B I' by following formula:

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 1, its feature Be: step 6) described in individual well always enter ground liquid measureWherein, L_vGround liquid measure, m is always entered for individual well³；L_viFor dividing Many bunches of volume fracturings of section i-th section enter ground liquid measure, m³。

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 4, its feature It is, described standard total hydrocarbon in gas logging value Q_TThe gas obtained well logging real-time during horizontal well drilling by antitrigonometric function method is surveyed Total hydrocarbon value q_tIt is corrected obtaining:

In formula, ν_ijI-th oil layer section jth sample point actual well drilled flow quantity in horizontal segment, m is met for boring³/min；t_ijFor boring When meeting the actual brill of i-th oil layer section jth sample point in horizontal segment, min/m；ν_oFor block standard well liquid entrance, m³/ min；t_oWhen boring for block standard, min/m；q_tijFor bore meet i-th oil layer section jth sample point in horizontal segment corresponding with Bore total hydrocarbon in gas logging value, %.

A kind of fine and close oil many bunches of volume fracturing horizontal well production Forecasting Methodologies of segmentation the most according to claim 5, its feature It is: the shear wave slowness of described rockCalculating formula when being used for logging well only to compressional wave time differenceTest, not test The situation of shear wave slowness.