CN106569288B - Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis - Google Patents

Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis Download PDF

Info

Publication number
CN106569288B
CN106569288B CN201610896467.8A CN201610896467A CN106569288B CN 106569288 B CN106569288 B CN 106569288B CN 201610896467 A CN201610896467 A CN 201610896467A CN 106569288 B CN106569288 B CN 106569288B
Authority
CN
China
Prior art keywords
reservoir
fracture
effectiveness
well
research area
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.)
Active
Application number
CN201610896467.8A
Other languages
Chinese (zh)
Other versions
CN106569288A (en
Inventor
祖克威
曾大乾
陈建
程秀申
范凌霄
孙力
崔长鹏
万龙
刘卉
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
Original Assignee
China Petroleum and Chemical Corp
Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Exploration and Development Research Institute of Sinopec Zhongyuan Oilfield Co filed Critical China Petroleum and Chemical Corp
Priority to CN201610896467.8A priority Critical patent/CN106569288B/en
Publication of CN106569288A publication Critical patent/CN106569288A/en
Application granted granted Critical
Publication of CN106569288B publication Critical patent/CN106569288B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V9/00—Prospecting or detecting by methods not provided for in groups G01V1/00 - G01V8/00
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00—Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24—Earth materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Geology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention discloses a fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis, and belongs to the technical field of fractured oil and gas reservoir evaluation. The method creatively establishes a reservoir fracture effectiveness index model by utilizing various factors influencing the effectiveness of the reservoir fracture, determines the reservoir fracture effectiveness index, establishes the corresponding relation between the fracture reservoir quality and each corresponding influence factor for representing the reservoir fracture effectiveness under the constraint of the plurality of reservoir fracture effectiveness influence factor factors, quantifies the fractured reservoir quality, and provides scientific basis for the economic and effective development of the low-permeability oil-gas reservoir. The evaluation method has strong adaptability and high precision.

Description

Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis
Technical Field
The invention relates to a fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis, and belongs to the technical field of fractured oil and gas reservoir evaluation.
Background
Currently, the center of gravity of world oil and gas exploration has been shifted from conventional reservoirs to unconventional reservoirs, wherein exploration and development of fractured reservoirs are particularly important. According to statistics, the reserves of fractured oil and gas reservoirs in the world account for about half of the total reserves which are proved to be important components for increasing the oil and gas production and storage in the future. China's North China fractured buried hill hydrocarbon reservoirs, Sichuan fractured cave carbonate rock and fractured clastic rock hydrocarbon reservoirs, Tarim basin Ordovician fractured cave stratum lithologic hydrocarbon reservoirs, Ordos basin fractured hydrocarbon reservoirs and the like are important hydrocarbon production bases in China. The natural fracture system is a main control factor for controlling and forming the low-permeability hydrocarbon reservoir, and the low-permeability hydrocarbon reservoir cannot be developed economically without natural fractures.
Although natural fractures are widely developed in the underground, the fractures are filled or closed through the transformation of later tectonic movement and diagenesis, and only the unfilled fractures can be used as a reservoir space and a migration channel in the underground. In the process of developing the oil and gas field, the quality of the fractured reservoir depends on the development condition of the reservoir fracture, and the foothold of reservoir fracture research is the evaluation on the effectiveness of the fracture. The effectiveness of reservoir fractures refers to that fractures are in an open state in the ground and can be used as channels for effective flow of fluid, the effectiveness not only refers to the opening degree of the fractures and the development degree of the fractures, and in the development process, the connectivity of the fractures and the opening capacity of the fractures are also main factors influencing the effectiveness of the fractures.
Although the research on the effectiveness of the fracture starts in the middle of the last century, the evaluation on the effectiveness of the fracture is still in an exploration stage, and the research results of predecessors are combined, wherein the influence factors of the effectiveness of the fracture comprise the relationship of the fracture trend to the same-stress-level maximum principal stress direction, the fracture inclination angle, the fracture opening degree, the fracture development strength, the fracture filling property, the fracture connectivity, the later-stage corrosion modification, the formation lifting, the formation corrosion, the abnormal high pressure and the like. Currently, the evaluation of fracture effectiveness is mainly performed from the following aspects: (1) considering the size of an included angle between the ground stress direction and the fracture strike, the imaging logging, the microseism and other production actual measurement data are generally utilized to analyze the horizontal maximum main stress direction of a single well, the effective fracture azimuth identified by the imaging logging is determined, and the smaller the included angle between the imaging logging and the microseism is, the higher the fracture effectiveness is. The evaluation method is a main evaluation method at present, but the cracks are three-dimensionally spread in space, so that not only is the relation between the crack trend and the principal stress azimuth considered, but also the relation between the dip angle of the crack and the principal stress azimuth cannot be ignored; (2) considering the influence of tectonic movement on the effectiveness of the fracture, and reflecting the seepage and communication capacity of the fracture by analyzing the density of the fracture development, namely the greater the fracture density is, the better the effectiveness is; (3) analyzing the filling property of the fracture, and reflecting the effectiveness of the fracture by analyzing and counting the proportion of unfilled fractures under a single well scope and a core fracture; (4) the relationship between the erosion effect and the effectiveness of the fractures is mainly applied to carbonate rocks, and the constructive effect of the erosion holes on reservoir fractures is qualitatively analyzed and identified.
In summary, the existing effectiveness evaluation method for reservoir fractures mainly performs single qualitative analysis on each influence factor, and cannot comprehensively evaluate the effectiveness of the reservoir fractures and the influence degrees of different influence factors. Therefore, the traditional effectiveness evaluation method cannot accurately indicate the quality of the fractured reservoir.
Disclosure of Invention
The invention aims to provide a fractured reservoir quality evaluation method based on reservoir fracture effectiveness clustering analysis, which overcomes the defects that the existing evaluation method is limited by a single factor, can only qualitatively depict fractured reservoir quality, and has low evaluation precision and poor applicability.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
the fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis determines factors influencing reservoir fracture effectiveness by utilizing the representing reservoir fracture parameters and the production measured data of the analysis and calculation ground stress state, establishes an index model for evaluating the reservoir fracture effectiveness, and quantitatively evaluates the fractured reservoir quality according to the corresponding relation between each influencing factor in the index model influencing the reservoir fracture effectiveness and the reservoir quality. The method specifically comprises the following steps:
1) determining fracture parameters and ground stress states of a single-well reservoir which characterize factors influencing the effectiveness of reservoir fractures according to basic data of a research area;
2) establishing a model for evaluating the effectiveness index of the reservoir fractures by using factors influencing the effectiveness of the reservoir fractures, and determining the effectiveness indexes of the reservoir fractures corresponding to different well sections of a single well in a research area;
3) repeating the step 1) and the step 2), and determining the effectiveness index of the reservoir fracture influenced by the reservoir fracture corresponding to each single well in the research area;
4) and grading the fractured reservoir quality of the research area according to the yield of each single well corresponding to each single well reservoir fracture of the research area, establishing the corresponding relation between each influence factor in the effective index model of the influenced reservoir fracture and the reservoir quality by utilizing a cluster analysis method, and quantitatively evaluating the fractured reservoir quality.
The method for determining the fracture parameters and the ground stress state of the single-well reservoir which characterize the factors influencing the effectiveness of the reservoir fractures in the step 1) comprises the following steps:
determining the parameters of the single-well reservoir fracture in the characterization research area: determining reservoir fracture parameters of the research area by using single-well imaging logging or directional coring data of the research area, wherein the reservoir fracture parameters comprise at least 2 parameters in fracture occurrence, fracture opening or fracture density of a characteristic reservoir corresponding to the depth of a single-well reservoir fracture development section of the research area;
determining the ground stress state of a single-well reservoir in a research area: the method comprises the steps of obtaining the current crustal stress direction of a reservoir by utilizing the well diameter collapse information of a single well in a research area and the single well drilling induced joint information, and determining the maximum and minimum principal stress of the single well in the research area by utilizing a single well hydraulic fracturing construction curve or a single well acoustic logging curve.
In the step 2), an index model for evaluating the effectiveness of the reservoir fractures, which is shown as the formula 1, is established according to factors influencing the effectiveness of the reservoir fractures, including the characteristic reservoir fracture parameters and the ground stress state:
formula 1:
Figure BDA0001130325400000031
in the formula: e is reservoir fracture effectiveness index, FiFactors for influencing reservoir fracture effectiveness, including fracture opening factor FOpening degreeFracture density factor FDensity ofGround stress factor FStress of earthFracture morphogenetic factor FThe state of birth，AiA weighting factor for each influencing factor.
Determining the effectiveness indexes of reservoir fractures corresponding to different well sections of a single well in a research area in the step 2), specifically: determining by using the fracture parameters and the ground stress state of the single-well reservoir determined in the step 1) and combining the established index model for evaluating the effectiveness of the reservoir fracture. The effectiveness index corresponding to the single-well reservoir fracture comprises FOpening degree、FDensity of、FStress of earthAnd FThe state of birthThe calculation is shown in formulas 2 to 5.
Formula 2: fOpening degree＝(b-bmin)/(bmax-bmin)；
In the formula: b is the single well reservoir fracture opening, bmaxAnd bminThe maximum value and the minimum value of the reservoir fracture opening of the research area are respectively.
Formula 3: fDensity of＝(ρ-ρmin)/(ρmax-ρmin)；
In the formula: rho is the fracture density of the single well reservoir, rhomaxAnd ρminThe maximum value and the minimum value of reservoir fracture density in the research area are respectively.
Formula 4: fStress of earth＝(Pmax-Pi)/(Pmax-Pmin)；
In the formula: piSingle well reservoir fracture opening pressure, PmaxAnd PminThe maximum value and the minimum value of the reservoir fracture opening pressure of the research area are respectively.
Formula 5: fThe state of birth＝(90-θ)/90；
In the formula: theta is the angle between the normal direction of the crack and the direction of the minimum principal stress.
The corresponding relation between each influence factor and the reservoir quality in the step 4) is shown as the following formula:
formula 6:
Figure BDA0001130325400000032
in the formula: p is a characteristic parameter of reservoir quality determined by fracture effectiveness index clustering analysis calculation, BiAnd B is a constant term for each influence factor Fi coefficient of each single well influencing the effectiveness index of the reservoir fracture.
The invention has the beneficial effects that:
the method creatively establishes a reservoir fracture effectiveness index model by utilizing various factors influencing the effectiveness of the reservoir fracture, determines the reservoir fracture effectiveness index, establishes the corresponding relation between the fracture reservoir quality and each corresponding influence factor for representing the reservoir fracture effectiveness under the constraint of the plurality of reservoir fracture effectiveness influence factor factors, quantifies the fractured reservoir quality, and provides scientific basis for the economic and effective development of the low-permeability oil-gas reservoir. The evaluation method has strong adaptability and high precision.
Drawings
FIG. 1 is a schematic flow chart of the technical solution of the present invention;
FIG. 2 is the current horizontal maximum principal stress pattern of each individual well in a certain layer of a study area in an example;
FIG. 3 is a cross-sectional view of the horizontal maximum and minimum principal stress calculated for a single well acoustic log of a study area;
FIG. 4 is a graph showing the relationship between the normal direction of the fracture surface and the principal stress;
FIG. 5 is a graph of correlation between fracture validity index cluster analysis calculation results and reservoir quality.
Detailed Description
The following examples are intended to illustrate the invention in further detail, but are not to be construed as limiting the invention in any way.
Example 1
Taking a certain crack type gas reservoir in the eastern part of China as an example, the effectiveness of the reservoir crack is evaluated, the flow schematic diagram is shown in figure 1, and the specific steps are as follows:
1. and (4) according to the actual data of the research area, performing statistical representation and determination of the ground stress state on all well fracture parameters, and preparing for calculating the effectiveness index of the subsequent fracture.
1.1 characterization and statistics of fracture parameters
And determining reservoir fracture parameters of the research area by using single-well imaging logging or directional coring well data of the research area, wherein the reservoir fracture parameters comprise fracture tendency, fracture inclination angle, fracture opening and fracture density of a characteristic reservoir corresponding to the depth of a single-well reservoir fracture development section of the research area. Combining the statistical information of the reservoir fracture parameters of different depths explained by the imaging logging of 69 wells in the research area, the maximum opening of the underground fractures in the research area is 1.228mm, the minimum opening is 0.001mm, the maximum fracture density after the fracture dip angle correction is 6.92 pieces/m, and the minimum fracture density is 0.108 pieces/m.
1.2 determination of the ground stress State
1.2.1 analysis of the ground stress Direction
The well wall collapse direction is vertical to the current horizontal maximum main stress direction, the well drilling induced joint direction is consistent to the current horizontal maximum main stress direction, and the single well ground stress direction is analyzed by using the data collected in the research area. As can be seen from the present horizontal maximum principal stress pattern of each single well in a certain layer of the research area shown in fig. 2, the overall crustal stress direction of the research area is the east-west direction, and the local crustal stress direction is deflected due to the fault and lithology difference.
1.2.2 calculation of the magnitude of the ground stress
And (3) calculating the maximum and minimum main stress of the single well level in the research area by using the acoustic logging, and selecting a Newberry model for the geostress calculation model according to the geological characteristics of the research area for calculation. And acquiring a horizontal maximum and minimum principal stress profile calculated by a single well acoustic logging curve of the research area as shown in the attached figure 3. The maximum value of the opening pressure of the deep cracks in the research area is 141.57MPa, and the minimum value is 31.23 MPa.
2. Determination of fracture effectiveness index
Analyzing and selecting factors influencing the effectiveness of the cracks in the research area, determining the magnitude of each influence factor of different crack development sections of a single well by considering the relationship between the crack opening, the crack density, the ground stress state and the current maximum level main stress and the crack occurrence, and then multiplying the magnitudes of the influence factor by weighting coefficients to obtain a crack effectiveness index, namely:
Figure BDA0001130325400000051
wherein E is a crack effectiveness index, Fi is each influence factor, and Ai is a weighting coefficient of each influence factor. When only four factors of fracture opening, fracture density, ground stress state and the relation between the maximum level principal stress and the fracture occurrence are considered, the above formula can be specifically written as follows:
E＝A1×Fopening degree+A2×FDensity of+A3×FStress of earth+A4×FThe state of birth。
Taking a single well 3277-3283 depth section of a research area as an example, the method for solving the factors of four influencing factors in the well depth section comprises the following steps:
(1) factor of crack opening
The fracture opening is an important parameter for representing the effectiveness of the fracture, the larger the fracture opening is, the better the connectivity of the fracture is, and the stronger the migration capacity of the fracture to the fluid is. The fracture opening index is defined to reflect the magnitude of a certain fracture opening value relative to all fracture opening values in the study area. The evaluation of the crack opening degree is mainly to count the interpretation results of the underground crack opening degrees of different depth sections of each well in the research area according to the logging information, and the maximum crack opening degree b of the underground crack in the research area is calculatedmaxIs 1.228mm, minimum opening bminThe diameter is 0.001mm, the imaging logging identification result of the opening of the 3277-3283m internal fracture is 0.251, and the opening factor is calculated as follows:
Fopening degree＝(b-bmin)/(bmax-bmin)；
Thus, FOpening degreeThe calculation result was 0.204.
(2) Fracture density factor
The linear density of the cracks can represent the development degree of the cracks, and the higher the development degree of the cracks is, the better the connectivity among crack group systems isThe greater its percolation capacity and, therefore, the greater its effectiveness. Fracture density factors are defined to reflect the relative magnitude of the density value of a certain section of fracture relative to the density values of all fractures in different depth sections in different regions of the study area. Taking the density of well interpretation fracture lines of different depth sections of each well in a research area as a maximum value rhomaxAnd minimum value of linear density ρminFor reference, the fracture density ρ of any interval in the deep underground layer has a density factor of:
Fdensity of＝(ρ-ρmin)/(ρmax-ρmin)；
The statistical analysis shows the maximum value rho of the fracture density in different depth sections of each single well in the research areamaxIs 6.92 bars/m, minimum value rhominThe number of the cracks is 0.108, and the density of the cracks in the single well 3277-3283m depth section is 0.312. Substituting into a formula to calculate to obtain F in the sectionDensity ofIs 0.03.
(3) Factor of earth stress
The seepage capability of the cracks is obviously influenced by the ground stress state, particularly in the oil and gas development process, the magnitude and the direction of the ground stress directly determine the strength of the opening capability of the underground deep cracks, and further influence the communication among crack groups. The ground stress factor is defined to reflect the relative ease of fracture opening within a certain depth segment. According to a calculation formula of fracture opening pressure in a tight reservoir:
Figure RE-GDA0001204678450000061
in the formula, PiFor fracture opening pressure, μ is the rock Poisson ratio, H is the depth of burial, ρsIs the volume weight of the rock, theta is the angle of inclination of the fracture, rhowIs the specific gravity of water; sigmaHAnd σhThe maximum and minimum horizontal principal stress are respectively, and beta is the included angle between the maximum horizontal principal stress and the fracture trend at present.
According to the calculation result of the mechanical parameters of the rock calculated by logging, the Poisson ratio mu of the rock is 0.23, the burial depth H is 3280m according to the average burial depth, and the average volume weight of the rock in the stratum is uniformly taken as 2.5g/cm3The crack angle theta is37 deg. of water specific gravity ρwIs 1g/cm3Horizontal maximum and minimum principal stress σHAnd σhRespectively 72MPa and 55.8MPa, the fracture trend is 20 degrees to the current maximum main stress included angle, the fracture opening pressure calculation formula is substituted, and the fracture opening pressure P of the single well 3277-3283m well section can be obtainediIs 92 MPa.
And calculating the fracture opening pressures of different groups of wells at different depths in the research area by combining the quantitative description result of the fracture occurrence of the well logging explanation in the research area. Counting the maximum opening pressure P of the target intervalmaxAnd a minimum opening pressure PminThe smaller the opening pressure in the underground depth is, the easier the fracture is to open, and further the fracture network is easier to form, and the seepage capacity of the reservoir is increased. The geostress factor for fracture effectiveness is:
Fstress of earth＝(Pmax-Pi)/(Pmax-Pmin)。
Calculating and counting the fracture opening pressure of different regions and depth sections of the research area to obtain the maximum value P of the fracture opening pressure of the research areamax141.57MPa, minimum value Pmin31.23MPa, the crack opening pressure of the single well 3277-3283m depth is 92MPa, and the calculation is carried out by substituting the formulaStress of earthIs 0.45.
(4) Fracture morphogenetic factor
When the dip angle of the crack is fixed, the effectiveness of the crack is mainly controlled by the trend of the crack, and when the included angle between the trend of the crack and the horizontal maximum principal stress direction is smaller, namely when the included angle between the normal direction of the crack surface and the horizontal minimum principal stress direction is smaller, the effectiveness of the crack is better; when the crack trends to be fixed, the effectiveness of the crack is mainly controlled by the inclination angle of the crack, taking a research area as an example, when the vertical principal stress of the crack is the maximum principal stress, and the included angle between the high-angle crack and the maximum principal stress is smaller, namely, when the included angle between the normal direction of the crack surface and the horizontal plane is smaller, the effectiveness of the crack is better. The maximum principal stress of the research area is in the vertical direction, the middle principal stress and the minimum principal stress level, therefore, in a three-dimensional space, the effectiveness of the crack depends on the included angle between the normal direction of the crack surface and the horizontal minimum principal stress, and the effectiveness of the crack is better when the included angle is smaller. Therefore, we define the fracture occurrence index to reflect the relative magnitude between the fracture face normal direction and the present minimum principal stress.
And calculating the included angle between the normal direction of the fracture of the single well and the minimum principal stress by taking the fracture azimuth data of different depths explained by logging as the basis and combining the geostress azimuth data simulated by the present geostress in the research area. FIG. 4 shows a graph of the relationship between the normal direction of the crack surface and the principal stress, wherein an arrow n indicates the normal direction of the crack surface, the normal direction of the crack surface and the horizontal minimum principal stress direction are theta, and the inclination of the crack surface and the horizontal minimum principal stress angle are theta1The angle between the normal line of the crack and the projection of the normal line on the horizontal plane is theta2According to the spatial functional relationship, the three angles satisfy the formula:
cosθ＝cosθ1×cosθ2。
the crack surface of the single well 3277-3283m section is inclined to be the same as the horizontal minimum principal stress included angle theta1Is 20 degrees, and the normal line of the crack surface forms an included angle theta with the horizontal plane projection of the normal line252.7 deg., so that the angle theta between the normal direction of the fracture surface and the horizontal direction of the minimum principal stress is 55.3 deg.. The greater the included angle, the less effective the fracture, and therefore, the shape factor of the fracture effectiveness is:
Fthe state of birth＝(90-θ)/90；
Substituting theta into the formula to obtain FThe state of birthIs 0.614.
When the stress state, the fracture occurrence, the opening and the density of a certain interval at the deep part of the underground part of the research area are known, the factors of four influencing factors can be calculated.
3. And (3) repeating the step 1 and the step 2 on the fracture and the ground stress parameters of the 69 wells in the research area, calculating all effectiveness indexes of different well sections of each single well in the research area, and compiling the corresponding relation between the effectiveness indexes of the reservoir fractures of each single well and the yield as shown in the table 1.
4. And (3) grading the reservoirs according to the gas production of the fractured reservoirs in the research area and the table 1, and performing cluster analysis on each influence factor of the fracture effectiveness indexes of different intervals and the fractured reservoir grades. And obtaining the corresponding relation between each influence factor of the reservoir fracture effectiveness index and the reservoir quality as shown in the formula 3, and evaluating the reservoir quality of the fractured reservoirs in other areas of the research area.
TABLE 1 relationship between fracture effectiveness index and gas production for different intervals
Figure BDA0001130325400000081
Counting the gas production of different intervals, grading the quality of the reservoir according to the production requirement, and making the gas production less than 1 x 10 in the research area4m3The gas production layer of the/d is defined as a 1-grade reservoir layer, and the gas production rate is 1-10 multiplied by 104m3The/d is defined as 2-grade reservoir, and the gas production is larger than 10 multiplied by 104m3 Stage 3 reservoirs were assigned/d (see Table 1). And matching the reservoir quality with the fracture effectiveness index influence factors of the corresponding depth sections, and performing cluster analysis by using SPSS software to obtain an optimal solution. The analysis result is as follows:
P＝2.90175×Fstress of earth+2.37337×FThe state of birth+0.657891×FOpening degree+2.72257×FDensity of-1.15316；
In the formula, P is a characteristic parameter of reservoir quality determined by fracture effectiveness index cluster analysis calculation, a graph 5 is drawn according to the analysis result, a correlation graph of the fracture effectiveness index cluster analysis calculation result P and the reservoir quality is drawn, as shown in FIG 5, the reservoir quality determined by calculation and the actual reservoir quality have a linear correlation relationship, and the correlation coefficient R is larger than 0.9. The value of each coefficient represents the weight of each influence factor and represents the relative strength of the influence degree on the effectiveness of the crack.
In conclusion, analysis by combining with the embodiment shows that the factors influencing the effectiveness evaluation of the reservoir fracture have mutual relevance, on the basis of establishing the reservoir fracture effectiveness evaluation model, the quality of the fractured reservoir is quantitatively evaluated under the common constraint of the influence factors of the evaluation model, and the consistency of the evaluation result and the grading result of the actual fractured reservoir is good. Can meet the requirement of economic and effective development of low-permeability oil-gas reservoirs.

Claims (3)

1. The fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis is characterized by comprising the following steps of: the method comprises the steps of determining factors influencing the effectiveness of reservoir fractures by using production measured data representing reservoir fracture parameters and analyzing and calculating the ground stress state, establishing an index model for evaluating the effectiveness of the reservoir fractures, and quantitatively evaluating the quality of fractured reservoirs according to the corresponding relation between each influencing factor in the index model influencing the effectiveness of the reservoir fractures and the reservoir quality;
the method specifically comprises the following steps:
1) determining fracture parameters and ground stress states of a single-well reservoir which characterize factors influencing the effectiveness of reservoir fractures according to basic data of a research area;
2) establishing a model for evaluating the effectiveness index of the reservoir fractures by using factors influencing the effectiveness of the reservoir fractures, and determining the effectiveness indexes of the reservoir fractures corresponding to different well sections of a single well in a research area;
3) repeating the step 1) and the step 2), and determining the effectiveness index of the reservoir fracture influenced by the reservoir fracture corresponding to each single well in the research area;
4) grading the fractured reservoir quality of the research area according to the yield of each single well corresponding to each single well reservoir fracture of the research area, establishing a corresponding relation between each influence factor in an influence reservoir fracture effectiveness index model and the reservoir quality by utilizing a cluster analysis method, and quantitatively evaluating the fractured reservoir quality;
determining the effectiveness indexes of reservoir fractures corresponding to different well sections of a single well in a research area in the step 2), specifically: determining by using the fracture parameters and the ground stress state of the single-well reservoir determined in the step 1) and combining the established index model for evaluating the effectiveness of the reservoir fractures;
the effectiveness index includes FOpening degree、FDensity of、FStress of earthAnd FThe state of birthSee formulas 2 to 5:
formula 2: fOpening degree＝(b-bmin)/(bmax-bmin)，
In the formula: b is the single well reservoir fracture opening, bmaxAnd bminRespectively the maximum value and the minimum value of the reservoir fracture opening of the research area;
formula 3: fDensity of＝(ρ-ρmin)/(ρmax-ρmin)，
In the formula: rho is the fracture density of the single well reservoir, rhomaxAnd ρminRespectively the maximum value and the minimum value of reservoir fracture density in a research area;
formula 4: fStress of earth＝(Pmax-Pi)/(Pmax-Pmin)，
In the formula: piSingle well reservoir fracture opening pressure, PmaxAnd PminRespectively the maximum value and the minimum value of reservoir fracture opening pressure in a research area;
formula 5: fThe state of birth＝(90-θ)/90，
In the formula: theta is an included angle between the normal direction of the crack and the direction of the minimum principal stress;
the corresponding relation between each influence factor and the reservoir quality in the step 4) is shown as the following formula:
formula 6:
Figure FDA0002742856490000021
in the formula: p is a characteristic parameter of reservoir quality determined by fracture effectiveness index clustering analysis calculation, BiAnd B is a constant term for each influence factor Fi coefficient of each single well influencing the effectiveness index of the reservoir fracture.
2. The evaluation method according to claim 1, characterized in that: the method for determining the fracture parameters and the ground stress state of the single-well reservoir which characterize the factors influencing the effectiveness of the reservoir fractures in the step 1) comprises the following steps:
determining the parameters of the single-well reservoir fracture in the characterization research area: determining reservoir fracture parameters of the research area by using single-well imaging logging or directional coring data of the research area, wherein the reservoir fracture parameters comprise at least 2 parameters in fracture occurrence, fracture opening or fracture density of a characteristic reservoir corresponding to the depth of a single-well reservoir fracture development section of the research area;
determining the ground stress state of a single-well reservoir in a research area: the method comprises the steps of obtaining the current crustal stress direction of a reservoir by utilizing the well diameter collapse information of a single well in a research area and the single well drilling induced joint information, and determining the maximum and minimum principal stress of the single well in the research area by utilizing a single well hydraulic fracturing construction curve or a single well acoustic logging curve.
3. The evaluation method according to claim 1, characterized in that: in the step 2), an index model for evaluating the effectiveness of the reservoir fractures, which is shown as the formula 1, is established according to factors influencing the effectiveness of the reservoir fractures, including the characteristic reservoir fracture parameters and the ground stress state:
formula 1:
Figure FDA0002742856490000022
in the formula: e is reservoir fracture effectiveness index, FiFactors for influencing reservoir fracture effectiveness, including fracture opening factor FOpening degreeFracture density factor FDensity ofGround stress factor FStress of earthFracture morphogenetic factor FThe state of birth，AiA weighting factor for each influencing factor.
CN201610896467.8A 2016-10-13 2016-10-13 Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis Active CN106569288B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610896467.8A CN106569288B (en) 2016-10-13 2016-10-13 Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610896467.8A CN106569288B (en) 2016-10-13 2016-10-13 Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis

Publications (2)

Publication Number Publication Date
CN106569288A CN106569288A (en) 2017-04-19
CN106569288B true CN106569288B (en) 2021-02-05

Family

ID=58532093

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610896467.8A Active CN106569288B (en) 2016-10-13 2016-10-13 Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis

Country Status (1)

Country Link
CN (1) CN106569288B (en)

Families Citing this family (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107315203A (en) * 2017-08-04 2017-11-03 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 A kind of tight sandstone reservoir dessert identification and integrated evaluating method
CN107908919B (en) * 2017-11-01 2020-03-31 中国矿业大学 Method for evaluating effectiveness of tight oil reservoir bedding joint
CN109339776B (en) * 2018-10-30 2021-03-09 中国石油大学(北京) Method for measuring anisotropic formation ground stress azimuth
CN112780253B (en) * 2020-01-20 2022-05-10 中国石油天然气集团有限公司 Method for predicting and evaluating fractured reservoir
CN114183121B (en) * 2020-09-14 2024-02-27 中国石油化工股份有限公司 Quantitative evaluation method and device for effectiveness of crack, electronic equipment and storage medium
CN112859171B (en) * 2020-11-19 2024-01-30 中国石油天然气股份有限公司 Method and apparatus for reservoir connectivity analysis using effective fractures
CN112987106B (en) * 2021-02-22 2023-05-16 河南理工大学 Coalbed methane well productivity potential evaluation method based on microseism static monitoring
CN112835124B (en) * 2021-03-10 2024-01-26 长江大学 Crack effectiveness evaluation method based on imaging logging and array acoustic logging data

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102200008A (en) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 Reservoir effectiveness identification method based on electrical imaging logging
WO2014205162A1 (en) * 2013-06-21 2014-12-24 Schlumberger Canada Limited Determining change in permeability caused by a hydraulic fracture in reservoirs
CN105134189A (en) * 2015-08-24 2015-12-09 西南石油大学 Logging GeoMechanics Identify Reservoir (LogGMIR) method
CN105334536A (en) * 2015-12-01 2016-02-17 中国石油大学(华东) Effectiveness evaluation method for compact sandstone reservoir map cracking system
CN105842754A (en) * 2016-03-24 2016-08-10 中国石油天然气股份有限公司 Method for quantitatively determining effectiveness of ultra-deep-layer compact sandstone crack

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102200008A (en) * 2010-03-26 2011-09-28 中国石油天然气股份有限公司 Reservoir effectiveness identification method based on electrical imaging logging
WO2014205162A1 (en) * 2013-06-21 2014-12-24 Schlumberger Canada Limited Determining change in permeability caused by a hydraulic fracture in reservoirs
CN105134189A (en) * 2015-08-24 2015-12-09 西南石油大学 Logging GeoMechanics Identify Reservoir (LogGMIR) method
CN105334536A (en) * 2015-12-01 2016-02-17 中国石油大学(华东) Effectiveness evaluation method for compact sandstone reservoir map cracking system
CN105842754A (en) * 2016-03-24 2016-08-10 中国石油天然气股份有限公司 Method for quantitatively determining effectiveness of ultra-deep-layer compact sandstone crack

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
缝洞型碳酸盐岩油藏开发描述及评价;鲁新便;《中国博士学位论文全文数据库 基础科学辑》;20050715(第7期);第4章 *

Also Published As

Publication number Publication date
CN106569288A (en) 2017-04-19

Similar Documents

Publication Publication Date Title
CN106569288B (en) Fractured reservoir quality evaluation method based on reservoir fracture effectiveness cluster analysis
CN105334536B (en) Tight sandstone reservoir chicken-wire cracking system effectiveness evaluation method
CN104992468B (en) Fracture-cavity type carbonate oil-gas reservoir three-dimensional geological modeling method
CN105651676B (en) A kind of reservoir heterogeneity characterizing method under horizontal well rule well pattern
CN110644980B (en) Comprehensive classification evaluation method for ultra-low permeability oil reservoir
CN109388817B (en) Three-dimensional modeling method for reservoir fractures
CN113901681B (en) Three-dimensional compressibility evaluation method for dual desserts of shale gas reservoir in whole life cycle
CN106869911B (en) Evaluation method for describing compressibility of shale reservoir
CN106569267A (en) Multi-scale crack model of compact low-penetration reservoir and modeling method of model
CN107038313B (en) Layered crustal stress fine description method based on numerical value core
CN103857876A (en) System and method for performing wellbore fracture operations
CN108825223A (en) A kind of texture feature extracting method of shale formation
CN112502701A (en) Classification evaluation method for comprehensive geology-engineering of low-permeability reservoir
CN110532507B (en) Method for improving drilling rate of compact oil reservoir of horizontal well
CN110390154A (en) A method of improving Complex reservoir reservoir numerical simulation efficiency
CN107704646B (en) Modeling method after compact reservoir volume modification
RU2761935C1 (en) Method for localization of promising zones in oil-producing strata
CN110989034B (en) Method for inverting logging transverse wave time difference by regression-fractal interpolation method
CN113971334A (en) Method and device for predicting fracture effectiveness based on three-dimensional stress field numerical simulation
Xiao et al. Perforation location optimization considering microscopic structure for multi-cluster fracturing technology
CN110244357A (en) A kind of structure destruction type Foregone pool sentences knowledge method
Lamberghini et al. Performance Based Reservoir Characterization in a Tight Gas Reservoir System-Case Study from Lajas and Punta Rosada Formations in the Neuquen Basin, Argentina
Lozano et al. Integrated Geocellular Static Model for Geomechanical and Dynamic Simulations in the Vaca Muerta Formation
Li et al. Three-dimensional reservoir architecture modeling by geostatistical techniques in BD block, Jinhu depression, northern Jiangsu Basin, China
CN109541683B (en) Construction method and system of underground point dam inclined stratum model grid

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