Background
Beach dam sand is one of the very important oil and gas reservoirs in the broken basin, and beach dam sand quantitative description has been the focus of attention of geologist both at home and abroad. With the continuous advancement of exploration progress, deep beach-dam sand reservoirs gradually become important fields of oil and gas exploration of land-phase fracture basin in the eastern part of China, and by taking victory oil fields as an example, beach-dam sandstone exploration reserves are accumulated and reported only in western oil areas to be approximately 2 hundred million tons, so that development potential is huge. At present, researchers at home and abroad mostly consider that the deposition of the broken sink lake basin beach dam is mainly controlled by three palos, including paleo-terrain, paleo-hydrodynamic force and paleo-reference surface. The paleo-topography can provide evidence for representing the spatial relationship between a deposition area and a deposition center, and the carrying channel for reacting the deposition has important significance for predicting a deposition system. By searching domestic and foreign documents, domestic scholars have mostly qualitative predictions on beach and dam sand reservoir descriptions, lack of applicability technical research on quantitative prediction of sand thickness, and do not form an effective method for quantitatively describing the boundary and thickness of the sand.
The method comprises the steps of dividing sequence levels of different single wells according to historical data and logging data in a target area, comparing the sand bodies of different single wells of the same sequence level to obtain sand body comparison results, preprocessing and standardizing well logging curves of the single wells, screening inversion sample curves from the well logging curves according to lithology data of the well logging curves, carrying out band-pass filtering on seismic data according to the sand body comparison results to obtain thin reservoir seismic reflection data, carrying out waveform phase control random inversion according to the inversion sample curves, the acoustic wave curves and density curves in the seismic data and the thin reservoir seismic reflection data to obtain thin reservoir inversion results, and determining the position of the thin reservoir according to the thin reservoir inversion results. The method and the device can improve the precision of predicting the thin reservoir in the sand body of the lake-phase beach dam.
The application number of the Chinese patent application CN201710711093.2 relates to a reservoir sand prediction method and device based on a wind field, an object source and a basin system, wherein the method comprises the steps of obtaining geological data of an area to be predicted, inputting the geological data into a preset wind field, the object source and the basin system model to generate beach-dam sand formation process data of the area to be predicted, wherein the wind field, the object source and the basin system model at least comprise various of an ancient-object-source recovery tool, an ancient-wind-direction recovery tool, an ancient-landform recovery tool and an ancient-water-depth recovery tool, and predicting specific distribution positions of beach-dam sand in the area to be predicted by adopting a geological method and a geophysical method according to the beach-dam sand formation process data. The method can effectively identify and predict the distribution position and range of the shallow water thin layer beach dam sand body, and improves the feasibility and accuracy of a reservoir sand body prediction mode.
In the Chinese patent application with the application number of CN202010691793.1, a beach-dam sand reservoir prediction method and a beach-dam sand reservoir prediction device are related, wherein the method comprises the steps of dividing a reservoir corresponding to logging data into a plurality of sand development periods of beach-dam sand according to the logging data of a target well in a beach-dam sand research area, selecting the development period with the largest oil-gas content as a target layer section based on oil-gas exploration results corresponding to each development period, determining seismic horizon data of the target layer section, determining a beach-dam sand development area in the beach-dam sand research area according to the seismic horizon data of the target layer section, and determining beach-dam sand reservoir distribution results of the beach-dam sand research area based on the beach-dam sand development area. The method and the device can effectively improve the identification precision of acquiring the beach-dam sand reservoir, can effectively improve the efficiency, reliability and accuracy of predicting the beach-dam sand reservoir, and can further provide an effective and accurate data basis for beach-dam sand exploration and oil gas exploitation.
In the Chinese patent application with the application number of CN201210212603.9, a method for establishing a beach and dam sandstone microphase identification mode is related. The method comprises the steps of accurately identifying lithology according to a logging curve by combining the logging curve with geological logging data, dividing a small layer in a geological microphase layer section by combining the logging curve with lithology, extracting curve values and form values of the small layer, accurately dividing the lithology of the beach-dam sandstone sediment microphase according to the curve values, the form values and the lithology, and establishing a beach-dam sandstone sediment microphase modeling identification standard library and a beach-dam sandstone microphase identification mode. The beach bar sandstone microphase identification mode is adopted, based on a beach bar sandstone sediment microphase mode identification standard library, the beach bar sandstone sediment microphase can be accurately identified by combining the well logging curve value and the form value, the method can be suitable for all beach bar sandstone sediment microphase division, and a judgment basis can be provided for reservoir properties and oil gas content.
The prior art is greatly different from the method, the technical problem which is needed to be solved by the user cannot be solved, and the method for quantitatively describing the sand body of the beach and dam based on the ancient topography is invented.
Disclosure of Invention
The invention aims to provide a quantitative description method of beach and dam sand bodies based on ancient terrains, which can improve prediction accuracy, is simple and convenient to operate.
The invention aims to realize the quantitative description method of the beach and dam sand body based on the ancient topography, which comprises the following technical measures:
step 1, restoring paleotopography in a certain deposition period by synthesizing earthquake, geology and logging data;
Step 2, selecting a plurality of representative well point samples;
Step 3, obtaining a reduction stratum depth value of the position of the single sample well point and drilling sand thickness data;
Step 4, establishing a sand thickness prediction formula based on the depth of the reduced stratum;
And 5, quantitatively predicting the sand thickness of the beach and dam in the non-moving area according to an established sand thickness prediction formula based on the depth of the reduced stratum.
The aim of the invention can be achieved by the following technical measures:
in step 1, determining the top and bottom interfaces of a stratum with accurate representation in a certain geological period, and restoring the paleo-topography in the geological period by comprehensively utilizing the methods of residual thickness calculation, stratum filling and alignment and deposit compaction correction.
And 2, selecting a plurality of representative well point samples at the underwater highland, the underwater uplift gentle slope zone and the underwater low-lying place respectively based on the ancient topography data analysis.
In step 2, the paleo-topography of a certain geological period is restored, classification is carried out according to different paleo-topography characteristics, three categories are divided into a subsurface plateau, an underwater raised gentle slope zone and an underwater low-lying position, and then a plurality of representative wells are selected at different paleo-topography positions, so that the number diversity of samples is ensured, and the reliability of analysis results is improved.
And 3, based on the paleo-topography recovery result of a certain geological period, counting the depth value of the restored stratum at the position of the single sample well point, and based on the fine stratum comparison result, counting the thickness value of the sand body drilled by the single sample well point in the geological period.
And 4, calculating the correlation coefficient of the thickness of the sand body encountered by the single sample well point drill and the depth of the reduced stratum on the basis of classifying the positions of different paleo-topography, obtaining a sand body thickness prediction formula based on the depth of the reduced stratum, and drawing a fitting trend graph of the thickness of the sand body along with the depth change of the reduced stratum.
And 4, fitting stratum depth values of single well points at different paleotopography positions counted in the previous step with the thickness of the sand body when drilling through a rock physical simulation experiment to respectively obtain a sand body thickness-reduction stratum depth fitting trend equation at the underwater high-speed land, the underwater raised gentle slope zone and the underwater low-lying position, and drawing a fitting trend graph of the sand body thickness along with the reduction stratum depth.
The trend equation of the sand thickness-reduction stratum depth fitting of the underwater highland and the underwater uplift gentle slope belt is as follows:
Y1=a1e-b1x1
Wherein Y 1 is the thickness of the sand body actually drilled at the well point, the unit is m, x 1 is the depth of the reduced stratum at the position where the well point is located, and the unit is m. a 1 is a fitting coefficient between the thickness of the actual drilling sand body at the well point and the depth of the reduced stratum, and b 1 is a constant.
In step 4, from the trend equation of the underwater plateau, the underwater raised gentle slope zone sand thickness-reduction stratum depth fitting, it can be seen that the underwater plateau and the gentle slope zone sand thickness are exponentially inversely related to the reduction stratum depth, and the higher the paleo-terrain is, the less the sand thickness is deposited.
In step4, the trend equation of sand thickness-reduced stratum depth fitting at the underwater depression is:
Y2=a2eb2x2
wherein Y 2 is the thickness of the sand body actually drilled at the well point in the underwater depression, the unit is m, and x 2 is the depth of the reduced stratum at the well point in the underwater depression, and the unit is m. a 2 is a fitting coefficient between the thickness of the actual drilling sand body at the well point and the depth of the reduced stratum, and b 2 is a constant.
In step 4, from the trend equation of the sand thickness-reduction stratum depth fitting of the underwater low-lying place, it can be seen that the sand thickness of the underwater low-lying place is exponentially and positively correlated with the reduction stratum depth, and the higher the paleo-topography is, the larger the deposited sand thickness is.
The quantitative description method of the beach and dam sand body based on paleo-topography is based on earthquake, geology and logging data, and is corrected through sedimentation and compaction. The comprehensive application of the method reduces the paleo-topography in the sand four-deposition period, extracts the reduced stratum depth value and the drilling sand body thickness data of the position of a single sample well point respectively on the basis of classifying the positions of different paleo-topography, takes the reduced stratum depth as an independent variable, takes the drilling sand body thickness of the well point as a dependent variable, and establishes a sand body thickness prediction formula based on the reduced stratum depth by regression fitting, thereby quantitatively predicting the sand body thickness of the beach dam in the non-moving area. Compared with the prior art, the method has the advantages that aiming at the difficult problem of beach-dam sand description, under the guidance of a former three-paleo sand control deposition mode, the method is deepened, based on paleo-topography recovery results, the quantitative description of the beach-dam sand in an inactive area is guided by utilizing the nonlinear index relation between the thickness of a reduced stratum at the position of a well point of a dense well pattern area and the thickness of a sand body in which the well point is drilled, the accuracy of beach-dam sand description and the reliability of results can be effectively improved, no quantitative prediction of the beach-dam sand thickness is guided based on paleo-topography recovery results in the conventional beach-dam sand fine description method, and the method has stronger innovation.
Detailed Description
It should be noted that the following detailed description is exemplary and is intended to provide further explanation of the invention. Unless defined otherwise, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs.
It is noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to be limiting of exemplary embodiments according to the present invention. As used herein, the singular forms also are intended to include the plural forms unless the context clearly indicates otherwise, and furthermore, it should be understood that when the terms "comprises" and/or "comprising" are used in this specification, they specify the presence of stated features, steps, operations, and/or combinations thereof.
The quantitative description method of the beach-dam sand body based on the paleo-topography is characterized in that paleo-topography is restored by comprehensively applying residual thickness calculation, stratum filling and leveling and sedimentation compaction correction methods based on earthquake, geology and logging data, a plurality of representative well point samples are selected at underwater high-altitude, underwater raised gentle slope zones and underwater low-altitude positions based on paleo-topography data analysis to obtain reduction stratum depth values and drilling sand body thickness data of a single sample well point, and on the basis of classifying different paleo-topography positions, a sand body thickness-reduction stratum depth fitting trend equation is established according to correlation coefficients of single sample well point drilling sand body thickness and reduction stratum depth obtained through experimental calculation, and the beach-dam sand body thickness of an inactive area is quantitatively predicted according to the fitting trend equation. The method disclosed by the invention has high prediction precision, effectively overcomes the limitation that the conventional beach-dam sand prediction method cannot be quantitatively described, is simple and easy to operate, and can realize the characterization of the thickness of the beach-dam sand body.
The following are several embodiments of the invention
Example 1
In the embodiment 1 of the invention, as shown in fig. 1, fig. 1 is a flow chart of a quantitative description method of beach and dam sand based on ancient topography of the invention. The quantitative description method of the beach and dam sand body based on the ancient topography comprises the following steps:
step 101, restoring paleotopography in a certain deposition period by synthesizing earthquake, geology and logging data;
102, based on ancient topography data analysis, selecting a plurality of representative well point samples at the underwater highland, the underwater raised gentle slope zone and the underwater low-lying place respectively;
Step 103, obtaining a reduction stratum depth value of the position of the single sample well point and drilling sand thickness data;
104, on the basis of classifying different paleo-terrain positions, taking the reduced stratum depth as an independent variable, taking the well point drilling sand body thickness as a dependent variable, establishing a sand body thickness prediction formula based on the reduced stratum depth by regression fitting, and drawing a fitting trend chart of the sand body thickness along with the reduced stratum depth;
And step 105, quantitatively predicting the sand thickness of the beach and dam in the non-moving area according to an established sand thickness prediction formula based on the depth of the reduced stratum.
The method for restoring the paleo-topography of a certain sedimentary period by synthesizing earthquake, geology and logging data specifically comprises the steps of determining the top and bottom interfaces of a stratum with accurate representation of the certain geology period, and restoring the paleo-topography of the geology period by comprehensively utilizing residual thickness calculation, stratum filling and sedimentary compaction correction methods.
The method is characterized in that a plurality of representative well point samples are selected from underwater high lands, underwater raised gentle slope zones and underwater low-lying places based on ancient topography data analysis, and concretely comprises the steps of restoring ancient topography in a certain geological period by using professional software, classifying according to different ancient topography characteristics, and dividing the ancient topography into three categories of underwater high lands, underwater raised gentle slope zones and underwater low-lying places, so that a plurality of representative wells are selected from different ancient topography positions, the number diversity of the samples is guaranteed, and the reliability of analysis results is improved.
The method for obtaining the reduction stratum depth value and the drilling sand thickness data of the position of the single sample well point comprises the steps of utilizing professional software to count the reduction stratum depth value of the position of the single sample well point based on the paleo-topography recovery result of a certain geological period, and counting the sand thickness value of the single sample well point in the geological period based on the fine stratum comparison result.
On the basis of classifying different paleo-topography positions, calculating the correlation coefficient of the thickness of a single sample well point drilling sand body and the depth of a reduced stratum, obtaining a sand body thickness prediction formula based on the depth of the reduced stratum, and drawing a fitting trend graph of the thickness of the sand body along with the depth of the reduced stratum. The trend equation of the sand thickness-reduction stratum depth fitting of the underwater highland/underwater uplift gentle slope belt is as follows:
Y1=a1e-b1x1
Wherein Y 1 is the thickness of the sand body actually drilled at the well point, the unit is m, x 1 is the depth of the reduced stratum at the position where the well point is located, and the unit is m.
According to the equation of trend of the underwater plateau/underwater uplift gentle slope zone sand thickness-reduction stratum depth fitting, the underwater plateau and the gentle slope zone sand thickness are inversely related to the reduction stratum depth in an exponential manner, and the higher the paleo-terrain is, the lower the deposited sand thickness is.
The trend equation of sand thickness-reduction stratum depth fitting of the underwater depression is as follows:
Y2=a2eb2x2
wherein Y 3 is the thickness of the sand body actually drilled at the well point in the underwater depression, the unit is m, and x 3 is the depth of the reduced stratum at the well point in the underwater depression, and the unit is m.
According to the trend equation of the sand thickness of the underwater low-lying part and the depth fitting of the reduction stratum, the sand thickness of the underwater low-lying part and the depth of the reduction stratum are exponentially and positively correlated, and the higher the paleo-topography is, the larger the deposited sand thickness is.
Example 2
The quantitative description method of the beach and dam sand body based on the ancient topography comprises the following steps:
1. The seismic, geological and well logging data are synthesized, and the paleo-terrain of a certain block in a sand four-sediment period is restored by comprehensively utilizing the residual thickness calculation, stratum filling and sediment compaction correction method. As shown in fig. 2, the underwater raised gentle slope zone, the underwater high land, the structural turning zone and other favorable paleo-terrain can be accurately depicted, and the object is distributed from southwest to northeast.
2. And (3) establishing a sand thickness prediction formula based on the reduced stratum depth by taking the reduced stratum depth as an independent variable and taking the thickness of the sand encountered by the well point drill as an independent variable and performing regression fitting.
Based on ancient topography data analysis, respectively selecting a plurality of representative well point samples at the underwater highland, the underwater raised gentle slope zone and the underwater low-lying place to obtain a reduction stratum depth value and drilling sand body thickness data of the position of a single sample well point;
On the basis of classifying different paleo-topography positions, the reduced stratum depth is taken as an independent variable, the well point drilling sand body thickness is taken as an independent variable, regression fitting is carried out to establish a sand body thickness prediction formula based on the reduced stratum depth, and a fitting trend chart of the sand body thickness changing along with the reduced stratum depth is drawn, as shown in figure 3. The fitting trend equation of the thickness of the sand body of the underwater uplift gentle slope belt and the depth of the reduction stratum is as follows:
Y1=2.638e-0.005x1
R2=0.9568
Wherein Y 1 is the thickness of the actual drilling sand at the position of the underwater rising gentle slope zone of the well point, the unit is m, x 1 is the depth of the reduced stratum at the position of the underwater rising gentle slope zone of the well point, the unit is m, and R 2 is the correlation between the thickness of the actual drilling sand at the position of the underwater rising gentle slope zone of the well point and the depth of the reduced stratum and the constant.
The sand thickness-reduced formation depth fitting trend equation at the underwater elevation is:
Y2=1.996e-0.005x2
R2=0.9266
Wherein Y 2 is the thickness of the actual drilling sand body at the underwater high place, the unit is m, x 2 is the depth of the reduction stratum at the underwater high place, the unit is m, and R 2 is the correlation and constant between the thickness of the actual drilling sand body at the underwater high place and the depth of the reduction stratum.
The trend equation of sand thickness-reduction stratum depth fitting of the underwater depression is as follows:
Y3=1529.5e0.009x3
R2=0.9699
Wherein Y 3 is the thickness of the actual drilling sand body at the well point in the underwater low-lying place, the unit is m, x 3 is the depth of the reduced stratum at the well point in the underwater low-lying place, the unit is m, and R 2 is the correlation between the thickness of the actual drilling sand body at the well point in the underwater low-lying place and the depth of the reduced stratum and the constant.
3. Drawing a sand thickness plane distribution prediction graph
And drawing a sand thickness plane distribution prediction graph according to a sand thickness prediction formula, as shown in fig. 4. The fitting trend graph of the sand thickness along with the depth change of the reduction stratum shows that the sand thickness of the paleo-terrain at different positions is obviously different, and the deposited sandstone thickness is respectively from large to small, namely, an underwater raised slow slope zone > an underwater high land > an underwater low depression. In the research example, the sand thickness of the underwater low-lying zone is exponentially and positively correlated with the depth of the reduced stratum, the higher the paleo-topography is, the larger the thickness of the deposited sand body is, and the sand thickness of the underwater high-lying zone and the gentle slope zone is exponentially and negatively correlated with the depth of the reduced stratum, the higher the paleo-topography is, and the smaller the thickness of the deposited sand body is.
Example 3
In a specific embodiment 3 to which the present invention is applied, the method comprises the steps of:
1. the seismic, geological and well logging data are synthesized, and the paleo-topography of a certain block in the sand four-sedimentation period is restored by comprehensively utilizing the residual thickness calculation, stratum filling and sedimentation compaction correction method, and the paleo-topography comprises favorable paleo-topography such as underwater raised gentle slope zones, underwater high lands and the like.
2. And (3) establishing a sand thickness prediction formula based on the reduced stratum depth by taking the reduced stratum depth as an independent variable and taking the thickness of the sand encountered by the well point drill as an independent variable and performing regression fitting.
Based on ancient topography data analysis, respectively selecting a plurality of representative well point samples at the underwater highland, the underwater raised gentle slope zone and the underwater low-lying place to obtain a reduction stratum depth value and drilling sand body thickness data of the position of a single sample well point;
On the basis of classifying different paleo-topography positions, taking the reduced stratum depth as an independent variable, taking the thickness of a well point drilling sand body as an independent variable, establishing a sand body thickness prediction formula based on the reduced stratum depth by regression fitting, and drawing a fitting trend chart of the sand body thickness along with the reduced stratum depth, wherein the fitting trend equation of the underwater uplift gentle slope belt sand body thickness-reduced stratum depth is as follows:
Y1=3.726e-0.005x1
R2=0.9007
Wherein Y 1 is the thickness of the actual drilling sand at the position of the underwater rising gentle slope zone of the well point, the unit is m, x 1 is the depth of the reduced stratum at the position of the underwater rising gentle slope zone of the well point, the unit is m, and R 2 is the correlation between the thickness of the actual drilling sand at the position of the underwater rising gentle slope zone of the well point and the depth of the reduced stratum and the constant.
The sand thickness-reduced formation depth fitting trend equation at the underwater elevation is:
Y2=2.132e-0.005x2
R2=0.9308
Wherein Y 2 is the thickness of the actual drilling sand body at the underwater high place, the unit is m, x 2 is the depth of the reduction stratum at the underwater high place, the unit is m, and R 2 is the correlation and constant between the thickness of the actual drilling sand body at the underwater high place and the depth of the reduction stratum.
The trend equation of sand thickness-reduction stratum depth fitting of the underwater depression is as follows:
Y3=1426.3e0.009x3
R2=0.9661
Wherein Y 3 is the thickness of the actual drilling sand body at the well point in the underwater low-lying place, the unit is m, x 3 is the depth of the reduced stratum at the well point in the underwater low-lying place, the unit is m, and R 2 is the correlation between the thickness of the actual drilling sand body at the well point in the underwater low-lying place and the depth of the reduced stratum and the constant.
3. Drawing a sand thickness plane distribution prediction graph
And drawing a sand thickness plane distribution prediction graph according to the sand thickness prediction formula. The fitting trend graph of the sand thickness along with the depth change of the reduction stratum shows that the sand thickness of the paleo-terrain at different positions is obviously different, and the deposited sandstone thickness is respectively from large to small, namely, an underwater raised slow slope zone > an underwater high land > an underwater low depression. In the research example, the sand thickness of the underwater low-lying zone is exponentially and positively correlated with the depth of the reduced stratum, the higher the paleo-topography is, the larger the thickness of the deposited sand body is, and the sand thickness of the underwater high-lying zone and the gentle slope zone is exponentially and negatively correlated with the depth of the reduced stratum, the higher the paleo-topography is, and the smaller the thickness of the deposited sand body is.
It should be noted that the above description is only a preferred embodiment of the present invention, and not intended to limit the present invention, but although the present invention has been described in detail with reference to the above embodiment, it will be apparent to those skilled in the art that modifications may be made to the technical solutions described in the above embodiment, or equivalents may be substituted for some of the technical features thereof. Any modification, equivalent replacement, improvement, etc. made within the spirit and principle of the present invention should be included in the protection scope of the present invention.
Other than the technical features described in the specification, all are known to those skilled in the art.