CN111155980A - Water flow dominant channel identification method and device - Google Patents
Water flow dominant channel identification method and device Download PDFInfo
- Publication number
- CN111155980A CN111155980A CN201811223943.5A CN201811223943A CN111155980A CN 111155980 A CN111155980 A CN 111155980A CN 201811223943 A CN201811223943 A CN 201811223943A CN 111155980 A CN111155980 A CN 111155980A
- Authority
- CN
- China
- Prior art keywords
- water flow
- water
- model
- change rate
- channel
- 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
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 305
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000008859 change Effects 0.000 claims abstract description 110
- 230000035699 permeability Effects 0.000 claims abstract description 72
- 238000010521 absorption reaction Methods 0.000 claims abstract description 61
- 230000008901 benefit Effects 0.000 claims abstract description 37
- 238000011161 development Methods 0.000 claims abstract description 25
- 238000004088 simulation Methods 0.000 claims abstract description 17
- 230000008569 process Effects 0.000 claims abstract description 14
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- IYLGZMTXKJYONK-ACLXAEORSA-N (12s,15r)-15-hydroxy-11,16-dioxo-15,20-dihydrosenecionan-12-yl acetate Chemical compound O1C(=O)[C@](CC)(O)C[C@@H](C)[C@](C)(OC(C)=O)C(=O)OCC2=CCN3[C@H]2[C@H]1CC3 IYLGZMTXKJYONK-ACLXAEORSA-N 0.000 claims description 21
- IYLGZMTXKJYONK-UHFFFAOYSA-N ruwenine Natural products O1C(=O)C(CC)(O)CC(C)C(C)(OC(C)=O)C(=O)OCC2=CCN3C2C1CC3 IYLGZMTXKJYONK-UHFFFAOYSA-N 0.000 claims description 21
- 239000011435 rock Substances 0.000 claims description 8
- 238000004590 computer program Methods 0.000 claims description 4
- 239000011800 void material Substances 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 6
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 4
- 239000003345 natural gas Substances 0.000 abstract description 2
- 239000003209 petroleum derivative Substances 0.000 abstract description 2
- 239000003921 oil Substances 0.000 description 30
- 239000004576 sand Substances 0.000 description 12
- 230000014509 gene expression Effects 0.000 description 6
- 230000007774 longterm Effects 0.000 description 6
- 239000011148 porous material Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004891 communication Methods 0.000 description 3
- 238000010276 construction Methods 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000000700 radioactive tracer Substances 0.000 description 3
- 238000005070 sampling Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000000704 physical effect Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000005540 biological transmission Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 238000011010 flushing procedure Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 208000019585 progressive encephalomyelitis with rigidity and myoclonus Diseases 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000012764 semi-quantitative analysis Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Geology (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Fluid Mechanics (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Geochemistry & Mineralogy (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Geophysics (AREA)
- Agronomy & Crop Science (AREA)
- Animal Husbandry (AREA)
- Marine Sciences & Fisheries (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention relates to a method and a device for identifying a water flow dominant channel, and belongs to the technical field of petroleum and natural gas exploration and development. The method utilizes an indoor core water drive experiment to establish the relation between the water passing multiple and the water absorption strength and the relation between the water passing multiple change rate and the permeability change rate, establishes a final water passing multiple model through numerical reservoir simulation on the basis, determines a water absorption strength change rate model and a permeability change rate model, and utilizes the two models to obtain a comprehensive judgment value so as to realize the identification of a water flow advantage channel. In the whole process, data are obtained through objective acquisition and calculation, the artificial influence is small, the water flow advantage channel can be identified quantitatively and objectively, and the accuracy of current water flow advantage identification is improved.
Description
Technical Field
The invention relates to a method and a device for identifying a water flow dominant channel, and belongs to the technical field of petroleum and natural gas exploration and development.
Background
The sandstone oil is hidden in the long-term water drive development process, the reservoir characteristics are changed continuously, particularly after entering a high water-cut period, under the comprehensive influence of multiple factors such as stratum heterogeneity, fluid geological acting force, gravity difference and the like, a large number of water flow dominant seepage channels are formed in sand bodies with high permeability, thin thickness and strong heterogeneity, so that the production well has the disadvantages of sudden water logging, high water content and low extraction degree, and meanwhile, invalid circulation formed by injected water along the water flow dominant seepage channels wastes a large amount of hydraulic and electric power resources, is high in cost and poor in benefit, and cannot realize benefit development under low oil price.
The key of effective water control of the sandstone oil reservoir in the high water cut period is the treatment of water flow advantage seepage, and how to accurately depict a water flow advantage channel is a difficult point. At present, a plurality of methods are used for describing the water flow dominant channel, including a tracer method, a well logging data method, a grey correlation method and the like. The tracer method mainly identifies the water flow dominant channel by monitoring parameters such as the concentration of the tracer, breakthrough time, peak time and the like, is greatly influenced by factors such as a well measuring instrument, field construction and the like, has long period and high cost, can only identify qualitatively and cannot determine a developed horizon. The logging data method utilizes logging data of wells in different periods to establish an identification template to identify the water flow dominant seepage channel, is limited by differences of logging data acquisition standards in different periods, has low representation precision and cannot represent all parameters of the logging data. The grey correlation method is used for normalizing and judging each index step by screening main factors influencing and marking the formation of the water flow advantage seepage channel, finally obtaining a comprehensive judgment index value of a well pair and establishing a qualitative identification expert system model.
In general, the existing identification method for the water flow advantage seepage channel mainly takes qualitative analysis or semi-quantitative analysis as a main part, the influence of artificial subjective factors is large, meanwhile, the development degree, the scale size and the grade of the water flow advantage seepage channel are not analyzed, and a method for quantitatively, objectively and accurately representing the water flow advantage channels of different grades is lacked.
Disclosure of Invention
The invention aims to provide a water flow advantage channel identification method, which aims to solve the problems that the water flow advantage channel cannot be accurately identified due to strong subjectivity and great artificial influence in the current water flow advantage channel identification process; based on the problems, the invention also provides a water flow dominant channel recognition device.
The invention provides a method for identifying a water flow dominant channel to solve the technical problems, which comprises the following steps:
1) collecting a rock core in a target area, and determining the relation between the water passing multiple and the change rate of the water absorption strength and the relation between the water passing multiple and the change rate of the permeability by performing a water flooding experiment on the rock core;
2) performing numerical reservoir simulation on a three-dimensional geological model in an initial state of a reservoir in a target area to obtain a water passing multiple model, and respectively determining a water absorption strength change rate model and a permeability change rate model according to the relationship between the water passing multiple and the water absorption strength change rate as well as the relationship between the water passing multiple and the permeability change rate;
3) determining the water absorption strength change rate and the permeability change rate of each layer of the target area according to the determined water absorption strength change rate model and the determined permeability change rate model;
4) and determining a comprehensive discrimination value of the corresponding layer according to the water absorption intensity variation rate and the permeability variation rate of each layer, determining a water flow advantage channel according to the comprehensive discrimination value, wherein the larger the comprehensive discrimination value is, the stronger the corresponding water flow advantage channel development area is.
The method utilizes an indoor core water drive experiment to establish the relation between the water passing multiple and the water absorption strength and the relation between the water passing multiple change rate and the permeability change rate, establishes a final water passing multiple model through numerical reservoir simulation on the basis, determines a water absorption strength change rate model and a permeability change rate model, and utilizes the two models to obtain a comprehensive judgment value so as to realize the identification of a water flow advantage channel. In the whole process, data are obtained through objective acquisition and calculation, the artificial influence is small, the water flow advantage channel can be identified quantitatively and objectively, and the accuracy of current water flow advantage identification is improved.
Further, in order to characterize the water flow dominant channels of different scales, the method also comprises the step of grading the water flow dominant channels according to the comprehensive discrimination values.
Furthermore, in order to realize accurate classification, the invention also provides a specific calculation mode of boundary values of each grade, wherein the classification of the water flow dominant channel is to divide the water flow dominant channel into a strong level, a medium level and a weak level, and the boundary value T of each grade1、T2And T3Respectively as follows:
T1＝1.4×Fmin
T2＝(Fmax-Fmin)/3
T3＝2×(Fmax-Fmin)/3
wherein FminIs the minimum value, F, of the comprehensive discrimination value of the water flow dominant channel in the target areamaxThe maximum value of the comprehensive judgment value of the water flow dominant channel of the target area is obtained, and if the comprehensive judgment value F of the water flow dominant channel of the target area is more than or equal to T1And is less than T2If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a weak water flow dominant channel2And is less than T3Then, the water flow is describedThe dominant channel is a medium water flow dominant channel development area, and if the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T3And then, the water flow dominant channel is indicated to be a strong water flow dominant channel.
Further, in order to obtain the water passing multiple model, the invention also provides a specific model establishing process, and the determination process of the water passing multiple model in the step 2) is as follows:
A. establishing a three-dimensional geological model in an initial state of an oil deposit in a target area, wherein the three-dimensional geological model comprises a structural model, a porosity model and a permeability model;
B. carrying out numerical reservoir simulation on the established geological model to obtain an accumulated water passing model;
C. and dividing the water passing model by the effective void volume model to obtain a water passing multiple model.
Furthermore, in order to identify the water flow dominant channel by fusing the water absorption strength change rate and the permeability change, the invention also provides how to obtain a comprehensive judgment value, wherein the comprehensive judgment value is equal to the sum of the water absorption strength change rate and the permeability change rate, and the step 4) is to superpose the water absorption strength change rate isoline plan and the permeability change rate isoline plan to obtain a comprehensive judgment value isoline plan so as to determine the comprehensive judgment value.
The invention also provides a water flow dominant channel identification device, which comprises a processor, a memory and a computer program stored on the memory and running on the processor, wherein the processor executes the computer program to realize the following steps:
1) determining the relation between the water passing multiple and the change rate of the water absorption strength and the relation between the water passing multiple and the change rate of the permeability according to the water drive experimental data of the rock core in the target area;
2) performing numerical reservoir simulation on a three-dimensional geological model in an initial state of a reservoir in a target area to obtain a water passing multiple model, and respectively determining a water absorption strength change rate model and a permeability change rate model according to the relationship between the water passing multiple and the water absorption strength change rate as well as the relationship between the water passing multiple and the permeability change rate;
3) determining the water absorption strength change rate and the permeability change rate of each layer of the target area according to the determined water absorption strength change rate model and the determined permeability change rate model;
4) and determining a comprehensive discrimination value of the corresponding layer according to the water absorption intensity variation rate and the permeability variation rate of each layer, determining a water flow advantage channel according to the comprehensive discrimination value, wherein the larger the comprehensive discrimination value is, the stronger the corresponding water flow advantage channel development area is.
In the identification process, data are obtained through objective acquisition and calculation, the artificial influence is small, the water flow advantage channel can be identified quantitatively and objectively, and the accuracy of the current water flow advantage identification is improved.
Further, in order to represent the water flow dominant channels with different scales, the device also comprises a step of grading the water flow dominant channels according to the comprehensive discrimination value.
Furthermore, in order to realize accurate classification, the invention also provides a specific calculation mode of boundary values of each grade, wherein the classification of the water flow dominant channel is to divide the water flow dominant channel into a strong level, a medium level and a weak level, and the boundary value T of each grade1、T2And T3Respectively as follows:
T1＝1.4×Fmin
T2＝(Fmax-Fmin)/3
T3＝2×(Fmax-Fmin)/3
wherein FminIs the minimum value, F, of the comprehensive discrimination value of the water flow dominant channel in the target areamaxThe maximum value of the comprehensive judgment value of the water flow dominant channel of the target area is obtained, and if the comprehensive judgment value F of the water flow dominant channel of the target area is more than or equal to T1And is less than T2If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a weak water flow dominant channel2And is less than T3If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a medium water flow dominant channel3And then, the water flow dominant channel is indicated to be a strong water flow dominant channel.
Further, in order to obtain the water passing multiple model, the invention also provides a specific model establishing process, and the determination process of the water passing multiple model in the step 2) is as follows:
A. establishing a three-dimensional geological model in an initial state of an oil deposit in a target area, wherein the three-dimensional geological model comprises a structural model, a porosity model and a permeability model;
B. carrying out numerical reservoir simulation on the established geological model to obtain an accumulated water passing model;
C. and dividing the water passing model by the effective void volume model to obtain a water passing multiple model.
Furthermore, in order to comprehensively consider the change rate of the water absorption strength and the change of the permeability to identify the water flow dominant channel, the invention also provides how to obtain a comprehensive judgment value, wherein the comprehensive judgment value is equal to the sum of the change rate of the water absorption strength and the change rate of the permeability, and the step 4) is to superpose the isoline planar graph of the water absorption strength change rate and the isoline planar graph of the permeability change rate to obtain an isoline planar graph of the comprehensive judgment value so as to determine the comprehensive judgment value.
Drawings
FIG. 1 is a flow chart of a water flow dominant channel identification method of the present invention;
FIG. 2 is a graph showing the relationship between the water passing times and the change rate of the water absorption strength in the example of the method of the present invention;
FIG. 3 is a graph of water passage multiple versus permeability change rate for a process embodiment of the present invention;
FIG. 4 is a three-dimensional view of a Puchenxi region 2+3 reservoir formation model in a method embodiment of the present invention;
FIG. 5 is a three-dimensional plot of a Puchenxi region 2+3 reservoir porosity model in a method embodiment of the present invention;
FIG. 6 is a three-dimensional plot of a Puchenxi region 2+3 reservoir permeability model in a method embodiment of the present invention;
FIG. 7 is a three-dimensional plot of the water-crossing multiple model of Puchenxi region 2+3 reservoir in an example of the method of the present invention;
FIG. 8 is a three-dimensional plot of a Puchenxi region 2+3 reservoir water absorption strength change rate model in a method embodiment of the present invention;
FIG. 9 is a three-dimensional plot of a Puchenxi region 2+3 reservoir permeability change rate model in a method embodiment of the present invention;
FIG. 10 is a plan view of a Puchenxi region 2+3 reservoir S25-1 small water absorption strength change rate contour in accordance with an embodiment of the method of the present invention;
FIG. 11 is a plan view of a Puchenxi region 2+3 reservoir S25-1 small zone permeability change rate contour according to an embodiment of the method of the present invention;
FIG. 12 is a plan view of integrated discriminatory values for the S25-1 small layer water flow dominant channel for Puchenxi region 2+3 reservoir in accordance with an embodiment of the method of the present invention;
fig. 13 is a graded plan view of a minor water flow dominant channel of the Puchenxi 2+3 reservoir S25-1 in accordance with an embodiment of the method of the present invention.
Detailed Description
Embodiment of water flow dominant channel identification method
The sandstone oil is hidden in the long-term water drive development process, the permeability is a main control factor for forming the water flow dominant channel, and the water absorption strength is a main parameter for reflecting the forming and development degree of the water flow dominant channel. In the embodiment, an indoor core water flooding experiment is utilized to establish a relational expression of water passing times and water absorption strength and a relational expression of water passing time change rate and permeability change rate, on the basis, a final water passing time model is established through numerical reservoir simulation, a water absorption strength change rate model and a permeability change rate model are established, the water absorption strength change rate and the permeability change rate of each layer of a research area are determined, a water flow advantage channel grading standard is established through comprehensive analysis, and a water flow advantage channel is graded and carved. The flow of the method is shown in figure 1, and the sand of a 2+3 oil reservoir in a certain sunken western region is 85-1The small reservoir structure is illustrated as an example, and the specific implementation steps of the method are as follows.
1. Based on indoor experiments, the relationship between the water passing multiple and the water absorption strength and the relationship between the water passing multiple and the permeability are determined.
1.1 collecting and organizing the data of the target oil reservoir core well, selecting the wells within five years after the start of industrial development to ensure that the samples can represent the initial state of the oil reservoir, and sampling in a target layer. For this embodiment, the data of the 2+3 core well of the oil reservoir in a depressed western region is obtained, and in order to ensure that the core well sample can represent the initial state of the oil reservoir, the new Pu 44 well of the core well at the early stage of the oil reservoir exploitation is selected for sampling.
1.2 carrying out long-term water flushing experiments on rock sample saturated formation water by using injected water, carrying out constant-speed water flooding, wherein the injected water is not lower than 1000PV, and acquiring water absorption strength change rate and permeability change rate data under different water passing multiples. In the initial stage of water flooding development, the physical properties of the reservoir change greatly, and the experimental data needs to be recorded intensively, so that in the initial stage of experiment, when the injection multiple does not reach 150PV, the experimental data are collected at 0.1PV, 0.2PV, 0.3PV, 0.4PV, 0.5PV, 0.6PV, 0.7PV, 0.8PV, 0.9PV, 1.0PV, 2.0PV, 5.0PV, 10.0PV, 20.0PV, 30.0PV, 50.0PV, 70.0PV, 100PV and 150.0PV respectively; when the water filling multiple is more than 150PV, experimental data are collected at every 100PV interval. The water absorption strength change and permeability change data obtained in the manner described above for this example at different water breakthrough multiples after sampling the fresh pu 44 well are shown in table 1.
TABLE 1
1.3 fitting a relational expression of the water passing multiple and the change rate of the water absorption strength and a relational expression of the water passing multiple and the change rate of the permeability respectively.
For this example, the water passage multiple versus permeability change was fitted to the data in Table 1 as follows:
ΔPERM＝4.1514ln(PV)+5.9617
wherein, Δ PERM is the permeability change rate, PV is the water passing multiple, and the relationship is shown in FIG. 2.
The fitted relation between the water passing multiple and the change rate of the water absorption strength is as follows:
ΔPEQ＝2.2226ln(PV)+3.2259
wherein, Δ PEQ is the change rate of the water absorption strength, and PV is the water passing multiple, and the relationship is shown in FIG. 3.
2. And (4) carrying out numerical simulation on the oil reservoir, and establishing a water absorption strength change rate model and a permeability change rate model of the current state of the oil reservoir.
2.1, establishing a three-dimensional geological model in the initial state of the oil deposit in the target area, in the embodiment, establishing the three-dimensional geological model in the initial state of the oil deposit in the target area by using PETREL software, wherein the three-dimensional geological model comprises a structural model, a porosity model and a permeability model, and data and steps required for establishing the model are as follows:
a. the method comprises the steps of establishing a construction model needing fault data, stratum interpretation data and well point layering data, firstly adjusting the spatial relationship among faults, establishing the fault model, generating a stratum surface by a Krigin method by utilizing the stratum interpretation data, correcting the stratum surface according to the well point layering data, adjusting the contact relationship and the fault distance between the faults and the stratum surface, and completing the establishment of the construction model.
b. Because the well pattern of the oil reservoir after long-term development is perfect, the well distance is close, the precision of the porosity plane graph and the permeability plane graph drawn by well point data is high, and the physical property characteristics of the oil reservoir can be accurately reflected, the porosity model and the permeability model are established by a deterministic modeling method.
For this example, Puchenxi region 2+3 reservoir Sand top 8 was established5-1The model of the formation in the small layer undeveloped state is shown in fig. 4, the porosity model is shown in fig. 5, and the permeability model is shown in fig. 6.
And 2.2, carrying out numerical reservoir simulation based on the geological model, deriving an accumulated excess water model in the last time step, and dividing the excess water model by the effective pore volume model to obtain a target reservoir excess water multiple model.
Specifically, the established geological model is led into Tnavigator software for numerical reservoir simulation to obtain a water passing multiple model, and the steps are as follows:
a. and (3) importing the oil-water interface depth, temperature and relative permeability data of the oil deposit, the high-pressure physical parameters of the rock, the density, viscosity and pressure coefficient parameters of water and crude oil into the model, fitting the reserves of the oil deposit, and finishing the initialization of the model.
b. And importing the production history data and the perforation data actually recorded by the target oil reservoir into the model, simulating the development and production history of the target oil reservoir, and completing history fitting.
c. And (5) deriving an oil reservoir numerical simulation accumulated water passing field of the last time step to obtain an accumulated water passing model.
d. And multiplying the model grid volume model by the porosity model to obtain an effective pore volume model, and dividing the water passing capacity model by the effective pore volume model to obtain a water passing multiple model.
For this embodiment, a numerical reservoir simulation is performed on the 85-1 small layer geological model on the second top of the 2+3 reservoir sand in the Puncheng region, and after the simulation is completed, there are 477 time steps, wherein the model at the 477 time step can represent the characteristics of the current state of the reservoir, and the derived accumulated water-passing capacity model at the 477 time step is derived and divided by the effective pore volume model, so as to obtain the current water-passing multiple model of the 85-1 small layer reservoir on the second top of the 2+3 reservoir sand in the Puncheng region, and the result is shown in FIG. 7.
And 2.3, establishing a current water absorption strength change rate model and a current permeability change rate model of the target area by using the target oil reservoir accumulated water passing multiple model obtained by calculation in 2.2 and the relational expression of the water passing multiple and the water absorption strength and the water passing multiple and the permeability change rate fitted in 1.3.
For this embodiment, based on the water passing multiple model of the 2+3 oil reservoir in the Puchenxi region, a current water absorption strength change rate model (as shown in FIG. 8) and a current permeability change rate model (as shown in FIG. 9) of the target region are established through fitted relational expressions of the water passing multiple and the water absorption strength and the water passing multiple and the permeability change rate.
3. And drawing a water absorption intensity change rate isoline plan view and a permeability change rate isoline plan view of each layer of the research area based on the water absorption intensity model and the permeability change rate model established in 2.3. For this example, Puchenxi region 2+3 reservoir Sand top 85-1The contour plane of the change rate of water absorption strength of the small layer is shown in FIG. 10, and the contour plane of the change rate of permeability is shown in FIG. 11.
4. The water absorption strength change rate belongs to a dynamic change index of the oil reservoir after long-term water drive, the permeability change rate reflects a reservoir change index after long-term water drive, a comprehensive judgment value F is obtained by utilizing the two parameters, the comprehensive judgment value is equal to the sum of the water absorption strength change rate and the permeability change rate, the water flow advantage channel is identified according to the judgment value, in order to accurately represent the water flow advantage channels with different scales, the grade of the identified water flow advantage channel is also divided, and the water flow advantage channel is graded and carved.
4.1 superposing the water absorption strength change rate isoline plan drawing in the step 3 with the permeability change rate isoline plan drawing an integrated discrimination value isoline plan drawing, and counting the minimum value Fmin and the maximum value Fmax of the integrated discrimination value of the water flow advantage channel of the target block; for this example, Puchenxi region 2+3 reservoir sand two was increased by 85-1Superposing the isoline plan of the change rate of the small layer water absorption intensity and the isoline plan of the change rate of the permeability, wherein the superposed water flow dominant channel comprehensive discrimination value plan is shown in FIG. 12, and the minimum value F in the statistical chart min2, maximum value Fmax＝55.98。
4.2 dividing the water flow dominant channel into a strong level, a medium level and a weak level, and dividing values T of each level1、T2And T3Respectively as follows:
T1＝1.4×Fmin
T2＝(Fmax-Fmin)/3
T3＝2×(Fmax-Fmin)/3
wherein FminIs the minimum value, F, of the comprehensive discrimination value of the water flow dominant channel in the target areamaxThe maximum value of the comprehensive discrimination value of the water flow dominant channel in the target area is obtained. For the embodiment, the boundary value T of the dominant channel of the three levels of water flows1＝2.8、T2＝17.99、T3＝35.98。
4.3 carrying out water flow advantage channel grading plan view based on the comprehensive discrimination value isoline plan view drawn in the step 4.1, if the comprehensive discrimination value F of the water flow advantage channel in the target area is more than or equal to T1And is less than T2Then, the water flow dominant channel is a development area of the weak water flow dominant channelIf the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T2And is less than T3If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is a medium water flow dominant channel development area3And then, the water flow dominant channel is a strong water flow dominant channel development area.
2+3 oil reservoir sand second-upper 8 based on Puchenxi region5-1Analyzing a small layer comprehensive discriminant value isoline plan, wherein a region larger than 2.8 and smaller than 17.99 in the plan is a weak water flow dominant channel development region, a region larger than 17.99 and smaller than 35.98 is a medium water flow dominant channel development region, and a region larger than 35.98 is a strong water flow dominant channel development region, and completing the steps to obtain the oil pool sand 8 in the 2+3 oil pool sand II in the Pupusilu region5-1And drawing a grading plan view of the laminar flow dominant channel, as shown in fig. 13.
Embodiment of water flow dominant channel recognition device
The invention discloses a water flow dominant channel recognition device, which is a computer or other equipment with data processing capability, and comprises: memory, processor and network module. The memory, processor, and network module are electrically connected to each other, directly or indirectly, to enable transmission or interaction of data. For example, the components may be electrically connected to each other via one or more communication buses or signal lines. The memory comprises at least one software functional module which can be stored in the memory in the form of software or firmware (firmware), and the processor executes various functional applications and data processing by running the software programs and modules stored in the memory, namely, the water flow dominant channel identification method in the method embodiment of the invention is realized.
The Memory may be, but is not limited to, a Random Access Memory (RAM), a Read Only Memory (ROM), a Programmable Read-Only Memory (PROM), an Erasable Read-Only Memory (EPROM), an electrically Erasable Read-Only Memory (EEPROM), and the like. The memory is used for storing programs, and the processor executes the programs after receiving the execution instructions.
The processor may be an integrated circuit chip having signal processing capabilities. The processor may be a general-purpose processor including a Central Processing Unit (CPU), a Network Processor (NP), and the like. But may also be a Digital Signal Processor (DSP)), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other programmable logic device, discrete gate or transistor logic device, discrete hardware components. The various methods, steps and logic blocks disclosed in embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like.
The network module is used for establishing communication connection between the electronic equipment and an external communication terminal through a network, and realizing the transceiving operation of network signals and data. The network signal may include a wireless signal or a wired signal.
Claims (10)
1. A water flow dominant channel identification method is characterized by comprising the following steps:
1) collecting a rock core in a target area, and determining the relation between the water passing multiple and the change rate of the water absorption strength and the relation between the water passing multiple and the change rate of the permeability by performing a water flooding experiment on the rock core;
2) performing numerical reservoir simulation on a three-dimensional geological model in an initial state of a reservoir in a target area to obtain a water passing multiple model, and respectively determining a water absorption strength change rate model and a permeability change rate model according to the relationship between the water passing multiple and the water absorption strength change rate as well as the relationship between the water passing multiple and the permeability change rate;
3) determining the water absorption strength change rate and the permeability change rate of each layer of the target area according to the determined water absorption strength change rate model and the determined permeability change rate model;
4) and determining a comprehensive discrimination value of the corresponding layer according to the water absorption intensity variation rate and the permeability variation rate of each layer, determining a water flow advantage channel according to the comprehensive discrimination value, wherein the larger the comprehensive discrimination value is, the stronger the corresponding water flow advantage channel development area is.
2. The method as claimed in claim 1, further comprising the step of ranking the dominant water flow paths according to the integrated discrimination values.
3. The method for identifying a dominant water flow channel as claimed in claim 2, wherein the classification of the dominant water flow channel is performed by dividing the channel into three classes of strong, medium and weak, and the boundary value T of each class1、T2And T3Respectively as follows:
T1＝1.4×Fmin
T2＝(Fmax-Fmin)/3
T3＝2×(Fmax-Fmin)/3
wherein FminIs the minimum value, F, of the comprehensive discrimination value of the water flow dominant channel in the target areamaxThe maximum value of the comprehensive judgment value of the water flow dominant channel of the target area is obtained, and if the comprehensive judgment value F of the water flow dominant channel of the target area is more than or equal to T1And is less than T2If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a weak water flow dominant channel2And is less than T3If the target area is the medium water flow dominant channel development area, the water flow dominant channel is indicated as the medium water flow dominant channel development areaThe comprehensive judgment value F of the water flow advantage channel is more than or equal to T3And then, the water flow dominant channel is indicated to be a strong water flow dominant channel.
4. The method for identifying the water flow dominant channel according to claim 1 or 3, wherein the determination process of the water flow multiple model in the step 2) is as follows:
A. establishing a three-dimensional geological model in an initial state of an oil deposit in a target area, wherein the three-dimensional geological model comprises a structural model, a porosity model and a permeability model;
B. carrying out numerical reservoir simulation on the established geological model to obtain an accumulated water passing model;
C. and dividing the water passing model by the effective void volume model to obtain a water passing multiple model.
5. The method for identifying the water flow dominant channel according to claim 1 or 3, wherein the step 4) is to superpose the water absorption intensity change rate isoline plan and the permeability change rate isoline plan to obtain a comprehensive discrimination value isoline plan so as to determine the comprehensive discrimination value.
6. A water flow dominant channel identification apparatus, comprising a processor and a memory, and a computer program stored on the memory and executable on the processor, wherein the processor executes the computer program to perform the steps of:
1) determining the relation between the water passing multiple and the change rate of the water absorption strength and the relation between the water passing multiple and the change rate of the permeability according to the water drive experimental data of the rock core in the target area;
2) performing numerical reservoir simulation on a three-dimensional geological model in an initial state of a reservoir in a target area to obtain a water passing multiple model, and respectively determining a water absorption strength change rate model and a permeability change rate model according to the relationship between the water passing multiple and the water absorption strength change rate as well as the relationship between the water passing multiple and the permeability change rate;
3) determining the water absorption strength change rate and the permeability change rate of each layer of the target area according to the determined water absorption strength change rate model and the determined permeability change rate model;
4) and determining a comprehensive discrimination value of the corresponding layer according to the water absorption intensity variation rate and the permeability variation rate of each layer, determining a water flow advantage channel according to the comprehensive discrimination value, wherein the larger the comprehensive discrimination value is, the stronger the corresponding water flow advantage channel development area is.
7. The apparatus as claimed in claim 6, further comprising a step of classifying the dominant water flow channel according to the integrated discrimination value.
8. The apparatus as claimed in claim 7, wherein the classification of the dominant water flow channel is divided into three stages of strong, medium and weak, and the boundary value T of each stage is1、T2And T3Respectively as follows:
T1＝1.4×Fmin
T2＝(Fmax-Fmin)/3
T3＝2×(Fmax-Fmin)/3
wherein FminIs the minimum value, F, of the comprehensive discrimination value of the water flow dominant channel in the target areamaxThe maximum value of the comprehensive judgment value of the water flow dominant channel of the target area is obtained, and if the comprehensive judgment value F of the water flow dominant channel of the target area is more than or equal to T1And is less than T2If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a weak water flow dominant channel2And is less than T3If the comprehensive judgment value F of the water flow dominant channel in the target area is more than or equal to T, the water flow dominant channel is indicated to be a medium water flow dominant channel3And then, the water flow dominant channel is indicated to be a strong water flow dominant channel.
9. The water flow dominant channel identification apparatus as claimed in claim 6 or 8, wherein the determination process of the water flow multiple model in step 2) is as follows:
A. establishing a three-dimensional geological model in an initial state of an oil deposit in a target area, wherein the three-dimensional geological model comprises a structural model, a porosity model and a permeability model;
B. carrying out numerical reservoir simulation on the established geological model to obtain an accumulated water passing model;
C. and dividing the water passing model by the effective void volume model to obtain a water passing multiple model.
10. The device for identifying the water flow dominant channel according to claim 6 or 8, wherein the step 4) is to superpose the water absorption intensity change rate isoline plan and the permeability change rate isoline plan to obtain a comprehensive discrimination value isoline plan so as to determine the comprehensive discrimination value.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811223943.5A CN111155980B (en) | 2018-10-19 | 2018-10-19 | Water flow dominant channel identification method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811223943.5A CN111155980B (en) | 2018-10-19 | 2018-10-19 | Water flow dominant channel identification method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111155980A true CN111155980A (en) | 2020-05-15 |
| CN111155980B CN111155980B (en) | 2023-09-01 |
Family
ID=70554462
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811223943.5A Active CN111155980B (en) | 2018-10-19 | 2018-10-19 | Water flow dominant channel identification method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111155980B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111749688A (en) * | 2020-08-10 | 2020-10-09 | 西南石油大学 | Method for predicting development position and direction of dominant seepage channel |
| CN112081582A (en) * | 2020-09-21 | 2020-12-15 | 中国石油大学(北京) | Prediction method, system and device for dominant channel in water-drive oil reservoir development |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102392636A (en) * | 2011-11-03 | 2012-03-28 | 中国石油大学(华东) | Underground fluid channel parameter description method for water drive development of sandstone reservoir |
| CN102419787A (en) * | 2011-10-19 | 2012-04-18 | 中国石油天然气股份有限公司 | Method and system for identifying water flow dominant channel of high-water-content oil field |
| CN103742130A (en) * | 2013-01-11 | 2014-04-23 | 中国石油化工股份有限公司 | High permeability band characterization method based on time varying |
| CN104727789A (en) * | 2013-12-23 | 2015-06-24 | 中国石油化工股份有限公司 | Dynamic description method of medium-high-permeability sandstone reservoir water flooding sweep efficiency and water flooding pore volumes |
| CN104806212A (en) * | 2015-02-16 | 2015-07-29 | 中国石油大学(北京) | Water drive oilfield advantage channel forming mechanism and development situation analysis method |
| CN105041307A (en) * | 2015-07-23 | 2015-11-11 | 中国石油大学(华东) | Industrial procedure for identifying preferential seepage channels of oil and gas reservoirs of clastic rock |
| CN105740563A (en) * | 2016-02-05 | 2016-07-06 | 北京华夏晓能石油技术有限公司 | Preferential channel identification technology for secondary development of mature oil field |
| CN107092719A (en) * | 2017-03-17 | 2017-08-25 | 中国石油天然气股份有限公司 | Method and device for identifying water injection dominant channel and selecting microsphere plugging particle size |
-
2018
- 2018-10-19 CN CN201811223943.5A patent/CN111155980B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102419787A (en) * | 2011-10-19 | 2012-04-18 | 中国石油天然气股份有限公司 | Method and system for identifying water flow dominant channel of high-water-content oil field |
| CN102392636A (en) * | 2011-11-03 | 2012-03-28 | 中国石油大学(华东) | Underground fluid channel parameter description method for water drive development of sandstone reservoir |
| CN103742130A (en) * | 2013-01-11 | 2014-04-23 | 中国石油化工股份有限公司 | High permeability band characterization method based on time varying |
| CN104727789A (en) * | 2013-12-23 | 2015-06-24 | 中国石油化工股份有限公司 | Dynamic description method of medium-high-permeability sandstone reservoir water flooding sweep efficiency and water flooding pore volumes |
| CN104806212A (en) * | 2015-02-16 | 2015-07-29 | 中国石油大学(北京) | Water drive oilfield advantage channel forming mechanism and development situation analysis method |
| CN105041307A (en) * | 2015-07-23 | 2015-11-11 | 中国石油大学(华东) | Industrial procedure for identifying preferential seepage channels of oil and gas reservoirs of clastic rock |
| CN105740563A (en) * | 2016-02-05 | 2016-07-06 | 北京华夏晓能石油技术有限公司 | Preferential channel identification technology for secondary development of mature oil field |
| CN107092719A (en) * | 2017-03-17 | 2017-08-25 | 中国石油天然气股份有限公司 | Method and device for identifying water injection dominant channel and selecting microsphere plugging particle size |
Non-Patent Citations (2)
| Title |
|---|
| 丁帅伟等: "水驱砂岩油藏优势通道识别综述", 《石油地质与工程》 * |
| 杜庆龙: "长期注水开发砂岩油田储层渗透率变化规律及微观机理", 《石油学报》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111749688A (en) * | 2020-08-10 | 2020-10-09 | 西南石油大学 | Method for predicting development position and direction of dominant seepage channel |
| CN112081582A (en) * | 2020-09-21 | 2020-12-15 | 中国石油大学(北京) | Prediction method, system and device for dominant channel in water-drive oil reservoir development |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111155980B (en) | 2023-09-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104533400B (en) | Method for reconstructing logging curve | |
| CN111425193B (en) | Reservoir compressibility evaluation method based on clustering analysis logging rock physical facies division | |
| CN112343587A (en) | Ultra-low permeability reservoir dominant seepage channel identification and characterization method | |
| CN104899411B (en) | A kind of reservoir productivity prediction model method for building up and system | |
| CN104564006A (en) | Hypotonic gas well fracturing water-producing capacity judgment method | |
| CN106353820A (en) | Method for interpreting well logging on basis of knowledge models | |
| US20230160304A1 (en) | Method and system for predicting relative permeability curve based on machine learning | |
| CN113033648A (en) | Method for realizing logging interpretation by using machine learning algorithm | |
| CN116427915A (en) | Conventional logging curve crack density prediction method and system based on random forest | |
| CN111155980B (en) | Water flow dominant channel identification method and device | |
| CN110988997A (en) | Hydrocarbon source rock three-dimensional space distribution quantitative prediction technology based on machine learning | |
| CN117251802B (en) | Heterogeneous reservoir parameter prediction method and system based on transfer learning | |
| CN110566196B (en) | Reservoir connectivity analysis method | |
| CN115629188B (en) | Core productivity simulation experiment system | |
| CN104880737A (en) | Multivariate Logistic method using logging information to identify type of underground fluid | |
| CN117131971A (en) | Reservoir dominant seepage channel prediction method based on Xgboost algorithm | |
| CN111236934B (en) | Method and device for determining flooding level | |
| CN114202036A (en) | Static sounding automatic soil layer dividing technology based on pattern recognition technology | |
| CN115705452A (en) | Novel recovery ratio prediction method for middle and later stages of integrated sandstone reservoir development | |
| CN111706323A (en) | Water flooded layer fine interpretation and evaluation method based on GWO-LSSVM algorithm | |
| CN118135103B (en) | Quantitative assessment method for water saturation of strong heterogeneous clastic rock reservoir | |
| CN115961939B (en) | Multistage fracturing crack identification method based on geological engineering integration | |
| CN117449846B (en) | Comprehensive diagnosis method for water production source of tight sandstone gas reservoir | |
| CN117634301A (en) | Method and equipment for predicting clastic rock reservoir flow unit type | |
| CN107944749B (en) | Shale gas block development potential evaluation method based on relative preference index |
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 |