CN110688728A - Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment - Google Patents
Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment Download PDFInfo
- Publication number
- CN110688728A CN110688728A CN201910769556.XA CN201910769556A CN110688728A CN 110688728 A CN110688728 A CN 110688728A CN 201910769556 A CN201910769556 A CN 201910769556A CN 110688728 A CN110688728 A CN 110688728A
- Authority
- CN
- China
- Prior art keywords
- deposition
- information
- dimensional
- data
- sediment
- 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
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Abstract
The invention relates to a method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments in a time domain and a water body environment, which aims at solving the problems that the existing sedimentary strata is weak to research in the time domain and lacks an effective depth domain-to-time domain conversion method, and provides a novel sedimentary strata one-dimensional depth domain-to-time domain conversion method based on a sedimentary forward modeling result for the first time, so that the sedimentary strata sedimentation, erosion, deposition history and duration time can be quantitatively researched, and the sedimentation characteristics, evolution processes and the relation between different types of sediments and the water body environment can be quantitatively analyzed.
Description
Technical Field
The invention belongs to the technical field of petroleum and natural gas exploration and development, and particularly relates to a quantitative analysis method for one-dimensional sediment deposition characteristics and a water body environment in a time domain.
Background
The discontinuity characteristics of sedimentary formations in the time domain have been widely recognized (Barrell, 1917; Sadler, 1981; Burgess and Wright, 2003; Straub and Foreman, 2018). Currently, the field of sedimentology mainly uses sedimentary formation records to research in a depth domain, and an effective method for quantitatively converting a formation from the depth domain to a time domain is lacked. The radioisotope dating method can determine the absolute age of the stratum and further convert a certain section of the stratum from a depth domain to a time domain, but on one hand, the method has a large error in dating the stratum older in the age and on the other hand, cannot accurately determine the deposition time of the stratum with smaller spacing (Sadler, 1981; Weedon, 2003). The physical simulation experiment of the water tank can recover the whole deposition process in a time domain in a frequent process recording mode, and the integrity of the stratum is researched (Strauband Foreman,2018), but the actual geology cannot be recovered through the physical simulation experiment.
Barrell,1917：Barrell,J.,1917,Rhythms and the measurements of geologictime:Geological Society of America Bulletin,v.28,p.745-904.
Sadler,1981：Sadler,P.M.,1981,Sediment accumulation rates and thecompleteness of stratigraphic sections:The Journal of Geology,v.89,p.569–584.
Burgess and Wright,2003：Burgess,P.M.,Wright,V.P.,2003,Numericalforward modelling of carbonate platform dynamics:an evaluation of complexityand completeness in carbonate strata:Journal of Sedimentary Research,v.73,no.5,p.637-652.
Straub and Foreman,2018：Straub,K.M.,Foreman,B.Z.,2018,Geomorphicstasis and spatiotemporal scales of stratigraphic completeness:Geology,v.46,no.4,p.311-314.
Disclosure of Invention
Aiming at the technical problems, the invention provides a deposition forward modeling method through process constraint, which comprises the steps of firstly establishing a reasonable three-dimensional deposition forward modeling numerical simulation model in a research area, selecting a target position on the basis, converting a deposition stratum at the position from a depth domain into a time domain by utilizing an autonomous research and development program, then quantitatively analyzing the deposition, deposition-free, ablation history and duration of the stratum in the time domain, analyzing the characteristics of deposition thickness, water depth, offshore distance and the like of different types of sediments, and establishing a novel method for quantitatively analyzing the deposition evolution process of the sediments from the time domain.
Specifically, the technical scheme adopted by the invention is as follows:
the method for quantitatively analyzing the sedimentation characteristics of the one-dimensional sediments and the water body environment in the time domain comprises the following steps:
(1) selecting a research area range and a target layer, establishing a target layer deposition evolution conceptual model of the research area by means of basic geological data analysis, literature research, current analogy and the like, and quantitatively determining the initial deposition bottom shape form, the structural settlement amount, the sea/lake plane lifting curve, the sediment source and supply rate, the growth period, conditions and rate of the carbonate rock and organic matters.
The basic geological data in the step (1) comprise seismic data, well drilling data, well logging data, core analysis data, various test data, sedimentary facies data, stratum thickness, sand shale content, previous research result reports and the like.
(2) On the basis of the step (1), performing three-dimensional deposition forward numerical simulation on a target layer in a research area by using software, comparing a simulation result with actual geological data, if the error is large, readjusting input parameters, repeating the simulation and the correction until the error of the simulation result reaches a reasonable range (the error is less than 10%), and finally establishing a reasonable three-dimensional deposition simulation model;
the software is deposition field professional forward modeling software which takes a simplified hydrodynamic momentum equation and a continuous equation as a core algorithm and considers various geological action processes.
The simplified hydrodynamic momentum equation expression is as follows:
wherein q is a fluid velocity vector, m/s;
t is time, s;
Φ is the ratio of pressure to fluid density, i.e.Φ ═ p-Rho, wherein the unit of p is Pa, and the unit of rho is kg/m3；
v is the kinematic viscosity, i.e. ν ═ μ/ρ, where μ is the fluid viscosity in pa.s and ρ is the fluid density in kg/m3；
g is the acceleration of gravity, m/s2。
The simplified equation of continuity expression is as follows:
wherein the content of the first and second substances,
is the gradient operator;
q is the fluid velocity vector, m/s.
The multiple geological processes include: the method comprises the following steps of sediment erosion, carrying and sedimentation processes, wave, tide and storm wave action, isopycnic flow action, aeolian sedimentation action, slumping gravity flow action, carbonate rock and organic matter growth process, carbonate rock diagenesis action, tectonic lifting action, sea level change, compaction action after sedimentation and crustal equilibrium sedimentation action.
The actual geological data refer to two-dimensional seismic profiles, field outcrop data, the sedimentary thickness of the current target layer, the sandstone percentage content, sedimentary facies, the thickness of the single-well target layer, lithology and logging data.
The software is deposition forward modeling software Sedsim;
(3) selecting a certain grid point to be analyzed, and expressing the grid point by a row number and a column number of the grid point;
(4) for a file with a GRAPH suffix output by a deposition simulation result, sequentially searching the grid point deposition information in each time interval by taking fixed time as a step length from the initial deposition stage of a target layer, wherein the following three conditions exist: a) if the grid point information is retrieved from the GRAPH file in the first time interval and only one line exists, writing the line information into a new text document and starting with the sequence number of '1', wherein the sequence number represents the deposit information in the 1 st time interval; b) if the grid point information is retrieved from the GRAPH file in the first time interval but two lines of information are corresponding, the two lines of information are written into two different text documents respectively, the text corresponding to the line with the positive numerical value is named as 'deposition history', the text corresponding to the line with the negative numerical value is named as 'ablation history', and the texts begin of the text with the sequence number of '1' respectively represents the deposition or ablation information in the 1 st time interval; c) if the grid point information is not retrieved from the GRAPH file in the first time interval, assigning the line of sediment information to be 0 in the newly-built text document, and writing the line of sediment information into the newly-built text document at the beginning of the sequence number 1; after the retrieval of the sediment information in the first time interval is finished, the steps are repeated by taking fixed time as a step length, the sediment information in the 2 nd, 3 rd and 4 … th time intervals is sequentially retrieved, the corresponding time sequence number is taken as the beginning, the information is sequentially written into a newly-built text file, and finally, a text file of the sediment information in different time intervals is generated;
the fixed time step in the step (4) corresponds to the display time interval data in the deposition simulation input parameters, and the user can customize the data.
The sediment information in the step (4) comprises the respective deposition thickness information of 8 types of sediments including 4 types of particle-size clastic rocks (coarse sandstone, medium sandstone, fine sandstone and mudstone), 2 types of carbonate rocks (carbonate rock 1 and carbonate rock 2) and 2 types of organic matters (organic matters 1 and organic matters 2), and the total thickness information of the sediments.
The data No. 4 to No. 11 in the GRAPH file in the step (4) respectively correspond to: and the sediment thickness data of the coarse sandstone, the medium sand debris rock, the fine sandstone, the mudstone, the carbonate rock 1, the carbonate rock 2, the organic matter 1 and the organic matter 2.
In the step (4), n in the deposition information in the 2 nd, 3 rd and 4 … th time intervals is sequentially searched, wherein the n represents the last time interval obtained by dividing the total deposition time of the target layer according to a fixed time step.
(5) Sequentially searching the ancient water depth and offshore distance data of the grid point in different deposition periods by taking fixed time as a step length from the initial deposition stage of a target layer for a file with a GRAPH suffix output by a deposition simulation result, writing the search result information into a newly-built text document at the beginning of a time interval sequence number if the grid information is searched in the GRAPH file in a certain time interval, assigning the ancient water depth and the offshore distance to be 0 if the grid information is not searched, writing the ancient water depth and the offshore distance into the newly-built text document at the beginning of the time interval sequence number, and finally establishing the text file of the ancient water depth and the offshore distance of the deposit deposition in different time intervals;
(6) and establishing a one-dimensional sediment deposition characteristic and water body environment comprehensive analysis graph in a time domain by respectively taking the sediment information, the ancient water depth, the offshore distance as a horizontal axis and the deposition time as a vertical axis, comparing the one-dimensional sediment deposition characteristic and the water body environment comprehensive analysis graph with the depth domain one-dimensional lithologic stratum simulation result, and quantitatively analyzing the deposition evolution process of the lithologic sediment at different grid point positions and the relation between the deposition evolution process and the water body environment.
Has the advantages that:
aiming at the problems that the existing sedimentary strata is weak in research in a time domain and lacks of an effective depth domain to time domain conversion method, the invention firstly provides a novel sedimentary strata one-dimensional depth domain to time domain conversion method based on a sedimentary forward modeling result, can quantitatively research the sedimentary strata deposition, the denudation, the non-deposition history and the duration time in the time domain, and quantitatively analyze the sedimentary characteristics, the evolution process and the relation between different types of sediments and the water environment.
The analysis method is simple and easy to implement, low in cost and high in operability, and provides a new technical means for extracting time information in the stratum.
Description of the drawings:
FIG. 1 is a three-dimensional sedimentary forward modeling model of debris rocks in Shanxi group in Ordos basin according to an embodiment of the invention;
FIG. 2 is a graph of the total thickness of sediments and the thickness of different types of sediments within a one-dimensional time domain of a clastic rock sedimentary formation according to an embodiment of the present invention;
FIG. 3 is a diagram illustrating the relationship between the sedimentary evolution process of the clastic rock sedimentary formation in one-dimensional time domain and the water body environment according to the embodiment of the invention;
FIG. 4 is a forward modeling model of the three-dimensional deposition of carbonate rock in the middle seismic-denier system lamp shadow system of the Sichuan basin in accordance with the embodiment of the present invention;
FIG. 5 is a total thickness of sediments and different types of sediments within a one-dimensional time domain of a carbonate sedimentary formation according to an embodiment of the present invention;
FIG. 6 shows the relationship between the sedimentary evolution process of carbonate sedimentary formations in one-dimensional time domain and the water environment in the embodiment of the present invention.
Detailed Description
The invention is described in detail below by way of exemplary embodiments.
The embodiment of the invention provides a method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water environment in a time domain, which comprises the following steps:
s1: selecting a research area range and a target layer, establishing a target layer deposition evolution conceptual model of the research area by means of basic geological data analysis, literature research, current analogy and the like, and quantitatively determining the initial deposition bottom shape form, the structural settlement amount, the sea/lake plane lifting curve, the sediment source and supply rate, the growth period, conditions and rate of the carbonate rock and organic matters.
In the step, the basic geological data comprises seismic data, well drilling data, well logging data, core analysis data, various test data, sedimentary facies data, stratum thickness, sand shale content, previous research result reports and the like.
S2: on the basis of S1, performing three-dimensional deposition forward numerical simulation on a target layer of a research area by using software, comparing a simulation result with actual geological data, if the error is large, readjusting input parameters, repeating the simulation and the correction until the error of the simulation result reaches a reasonable range (the error is less than 10%), and finally establishing a reasonable three-dimensional deposition simulation model;
in the step, the software is deposition field professional forward modeling software which takes a simplified hydrodynamic momentum equation and a continuous equation as core algorithms and considers various geological action processes. Wherein the simplified hydrodynamic momentum equation expression is as follows:
wherein the content of the first and second substances,
comprises the following steps: differentiation of the parameter t;
represents the rate of change of the function q due to a single change in t, i.e. the partial differential of q over t;
q is the fluid velocity vector, m/s;
t is time, s;
Φ is the ratio of pressure and fluid density, i.e., Φ is p/ρ, where p is in Pa and ρ is in kg/m3；
V is kinematic viscosity, that is, v is μ/ρ, where μ is fluid viscosity in pa.s and ρ is fluid density in kg/m3；
g is the acceleration of gravity, m/s2。
The simplified equation of continuity expression is as follows:
wherein the content of the first and second substances,is the gradient operator;
q is the fluid velocity vector, m/s.
The various geologic processes considered by the software include: the method comprises the following steps of sediment erosion, carrying and sedimentation processes, wave, tide and storm wave action, isopycnic flow action, aeolian sedimentation action, slumping gravity flow action, carbonate rock and organic matter growth process, carbonate rock diagenesis action, tectonic lifting action, sea level change, compaction action after sedimentation and crustal equilibrium sedimentation action.
The general steps of establishing a reasonable three-dimensional deposition forward numerical simulation model in a research area by utilizing Sedsim software are as follows: firstly, establishing a conceptual model of a reasonable deposition process in a research area on the basis of literature investigation and basic geological analysis; secondly, acquiring main input parameters required by deposition simulation, wherein the main input parameters comprise a work area range, a grid number, a simulation time range, a deposition initial bottom shape, a constructed settlement amount, a sea/lake plane change curve, a source direction, a supply rate and a deposition component; and then, performing three-dimensional deposition forward numerical simulation by using software, comparing the simulation result with actual geological data, if the error is large, readjusting the input parameters, repeating the simulation and the correction until the error of the simulation result reaches a reasonable range, and establishing a final reasonable three-dimensional deposition simulation model.
In the step, the actual geological data mainly refers to two-dimensional seismic profiles, field outcrop data, sedimentary thickness of the current target layer, sandstone percentage content, sedimentary facies, thickness of the single-well target layer, lithology and logging data.
S3: selecting a certain grid point to be analyzed, and expressing the grid point by a row number and a column number of the grid point;
s4: for a file with a GRAPH suffix output by a deposition simulation result, sequentially searching the grid point deposition information in each time interval by taking fixed time as a step length from the initial deposition stage of a target layer, wherein the following three conditions exist: a) if the grid point information is retrieved from the GRAPH file in the first time interval and only one line exists, writing the line information into a new text document and starting with the sequence number of '1', wherein the sequence number represents the deposit information in the 1 st time interval; b) if the grid point information is searched in the GRAPH file in the first time interval but two lines of information are corresponding, adding corresponding No. 4 to No. 11 data in the two lines of information to form a new line of data, writing the new line of data into a new text document, starting with a sequence number of '1', and representing the sediment information in the 1 st time interval; c) if the grid point information is not retrieved from the GRAPH file in the first time interval, assigning the line of sediment information to be 0 in the newly-built text document, and writing the line of sediment information into the newly-built text document at the beginning of the sequence number 1; after the retrieval of the sediment information in the first time interval is finished, the steps are repeated by taking fixed time as a step length, the sediment information in the 2 nd, 3 rd and 4 … th time intervals is sequentially retrieved, the corresponding time sequence number is taken as the beginning, the information is sequentially written into a newly-built text file, and finally, a text file of the sediment information in different time intervals is generated;
in this step, the fixed time step corresponds to the "display time interval" data in the deposition simulation input parameters, and the user can customize the data.
The sediment information comprises sediment thickness information of 4 types of particle size clastic rock, 2 types of carbonate rock and 2 types of organic matters, and 8 types of sediment, and total sediment thickness information. Outputting a suffix to the deposition simulation result as a GRAPH file, wherein the data No. 4 to 11 correspond to: and the sediment thickness data of the coarse sandstone, the medium sand debris rock, the fine sandstone, the mudstone, the carbonate rock 1, the carbonate rock 2, the organic matter 1 and the organic matter 2.
And sequentially searching the 'n' in the deposition information in the 2 nd, 3 rd and 4 … th time intervals to obtain the last time interval, wherein the 'n' in the deposition information indicates that the total deposition time of the target layer is divided according to the fixed time step.
S5: sequentially searching the ancient water depth and offshore distance data of the grid point in different deposition periods by taking fixed time as a step length from the initial deposition stage of a target layer for a file with a GRAPH suffix output by a deposition simulation result, writing the search result information into a newly-built text document at the beginning of a time interval sequence number if the grid information is searched in the GRAPH file in a certain time interval, assigning the ancient water depth and the offshore distance to be 0 if the grid information is not searched, writing the ancient water depth and the offshore distance into the newly-built text document at the beginning of the time interval sequence number, and finally establishing the text file of the ancient water depth and the offshore distance of the deposit deposition in different time intervals;
s6: and establishing a one-dimensional sediment deposition characteristic and water body environment comprehensive analysis graph in a time domain by respectively taking the sediment information, the ancient water depth, the offshore distance as a horizontal axis and the deposition time as a vertical axis, comparing the one-dimensional sediment deposition characteristic and the water body environment comprehensive analysis graph with the depth domain one-dimensional lithologic stratum simulation result, and quantitatively analyzing the deposition evolution process of the lithologic sediment at different grid point positions and the relation between the deposition evolution process and the water body environment.
In order to more clearly and specifically describe the method for quantitatively analyzing the one-dimensional sediment deposition characteristics and the water environment in the time domain according to the embodiments of the present invention, the following description will be made with reference to specific embodiments.
The analysis was performed by taking the clastic rock stratum in the Shanxi group in the Ordos basin and the carbonate rock stratum in the Shandan series lamp shadow group in the Sichuan basin as an example.
(1) Shanxi group clastic rock in middle of Ordos basin
S1: selecting Shanxi group stratums in the middle of an Ordos basin as a research target layer, establishing a target layer deposition evolution conceptual model of a research area by means of basic geological data (the basic geological data comprises seismic data, well drilling data, well logging data, core analysis data, various test data, sedimentary facies data, stratum thickness, sandstone content, previous research result report) analysis, literature investigation, current analogy and the like, and quantitatively determining the bottom form, the structural settlement amount, a sea/lake plane lifting curve, a sediment source and supply rate, and the growth period, conditions and rate of carbonate rocks and organic matters at the early deposition stage of the target layer.
S2: on the basis of S1, carrying out three-dimensional deposition forward numerical simulation and result correction on the Shanxi group stratum in the middle of the Ordos basin by using deposition forward simulation software Sedsim, and establishing a reasonable three-dimensional deposition simulation model as shown in the attached drawing 1;
s3: selecting a virtual well PW-1 of the delta deposition environment in the work area, wherein the grid point of the well is the 60 th row and the 40 th column;
s4: sequentially searching the grid point sediment information in each time interval by taking 40000 years as step length from the initial stage of deposition of the target layer to the file with GRAPH as the suffix output by the deposition simulation result to obtain text files of the sediment deposition information in 285-280 Ma in different time intervals;
s5: sequentially searching the ancient water depth and offshore distance data of the grid point at different deposition periods by taking 40000 years as step length from the initial deposition stage of the target layer according to a file with a GRAPH suffix output by the deposition simulation result, and establishing text files of the ancient water depth and the offshore distance of the deposit in different time intervals;
s6: and respectively making a total thickness of the sediments in the one-dimensional time domain of the position point and a thickness map (figure 2) of different types of sediments and a relation map (figure 3) between the deposition evolution process in the one-dimensional time domain and the water environment by using the sediment information, the ancient water depth, the offshore distance as a horizontal axis and the deposition time as a vertical axis and using a mapping software graph, comparing the relation map with the depth domain lithology column, and quantitatively analyzing the deposition evolution process of the sediments in the time domain.
As shown in fig. 2: the pseudo-well PW-1 stratum of the Shanxi group (285-280 Ma) in the middle of the Ordos basin is mainly formed by clastic rock deposition and contains a small amount of organic matters, the thickness of the deposited stratum reaches more than 150 meters, and the thickness of the medium clastic rock deposition is the largest. The two sections (285-282.4 Ma) of the mountain deposit a large amount of medium-grained clastic rock and a small amount of fine-grained clastic rock, clay and organic matters, and the one section (282.4-280 Ma) of the mountain still mainly deposits the medium-grained clastic rock and simultaneously starts to deposit a large amount of coarse-grained clastic rock. As shown in fig. 3: when the stratum is deposited, the stratum is degraded, and the stratum degradation of the virtual well PW-1 of the Shanxi group in the middle of the Ordos basin mainly occurs in the second section of the mountain (285-282.4 Ma) and is small in thickness; the thickness of the stratum deposited at a certain moment is related to the rise and fall of the ancient water depth and the sea level at the moment, and the stratum deposition of the clastic rock is not facilitated when the water depth is too deep or too shallow; the distance of the virtual well PW-1 offshore changes always in accordance with the change trend of the ancient water depth.
(2) Middle seismic denier lamp shadow group carbonate rock of Sichuan basin
S1: selecting a carbonate rock stratum of a seismic-denier lamp shadow group in the middle of the Sichuan basin as a research target layer, establishing a target layer deposition evolution conceptual model of a research area by means of basic geological data analysis, literature investigation, current analogy and the like, and quantitatively determining the bottom shape form, the structural sedimentation amount, a sea/lake plane lifting curve, a sediment source and supply rate, and the growth period, conditions and rate of carbonate rock and organic matters at the initial deposition stage of the target layer.
S2: on the basis of S1, performing three-dimensional sedimentation forward modeling numerical simulation and result correction on the carbonate rock stratum of the camerawork group of the lamp shadow system in the middle of the Sichuan basin by using sedimentation forward modeling software Sedsim, and establishing a reasonable three-dimensional sedimentation simulation model as shown in an attached figure 4;
s3: selecting a virtual well PW-2 of a delta deposition environment in the work area, wherein the grid point of the well is the 128 th row and 64 th column;
s4: sequentially searching the grid point sediment information in each time interval from the initial deposition stage of the target layer by taking 20000 years as step length for a file with a GRAPH suffix output by the deposition simulation result to obtain text files of the sediment deposition information in different time intervals in 551.1-539 Ma;
s5: sequentially searching the ancient water depth and offshore distance data of the grid point at different deposition periods by using 20000 years as step length from the initial deposition stage of the target layer for the file with the GRAPH suffix output by the deposition simulation result, and establishing text files of the ancient water depth and the offshore distance of the deposit in different time intervals;
s6: and respectively making a total thickness of the sediments in the one-dimensional time domain of the position point and a thickness map (figure 5) of different types of sediments and a relation map (figure 6) between the deposition evolution process in the one-dimensional time domain and the water environment by using the sediment information, the ancient water depth, the offshore distance as a horizontal axis and the deposition time as a vertical axis and using a mapping software graph, comparing the relation map with the depth domain lithology column, and quantitatively analyzing the deposition evolution process of the sediments in the time domain.
As shown in fig. 5: the PW-2 stratum of the seismic-denier lamp shadow group (551.1-539 Ma) in the middle of the Sichuan basin is mainly carbonate rock sediment and contains a small amount of clastic rock and organic matter sediment, and the thickness of the sediment stratum reaches 1120 meters. The lamp shadow component is 4 layer sections (a first lamp section, a second lamp section, a third lamp section and a fourth lamp section), wherein the deposition thickness of the first lamp section (551.1-549.1 Ma) and the third lamp section (554.5-542.5 Ma) is smaller, and the thick carbonate rock stratum is deposited by the second lamp section (549.1-545.5 Ma) and the fourth lamp section (542.5-540.5 Ma). As shown in fig. 6: during stratum deposition, a stratum ablation effect exists, the middle seismic denier lamp shadow group virtual well PW-2 stratum ablation of the Sichuan basin mainly occurs in the movement period of the Tung Bay, a small amount of stratum ablation exists in the second lamp section and the fourth lamp section, wherein the stratum is not deposited and a large amount of stratum ablation exists in the I curtain period (545.5-544.5 Ma) of the Tung Bay, the stratum is deposited and ablated at the same time in the II curtain period (540.5-539 Ma) of the Tung Bay, and the ablation thickness is larger than the deposition thickness. The development of the carbonate rock stratum is closely related to the water depth, the water depth of the area is-30 meters, the development of the carbonate rock stratum is facilitated, and the development of the stratum is not facilitated if the water depth exceeds 30 meters. The elevation of the sea level also has important influence on the development of the stratum, the first lamp section and the third lamp section of the lamp shadow group in the area are a sea invasion area, the second lamp section and the fourth lamp section are a high-level area, and the high-level area is a main development layer section of the reservoir, so that the lithofacies ancient geographic characteristics of the lamp shadow group can be clarified.
Claims (9)
1. A method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in a time domain is characterized by comprising the following steps of: the method comprises the following steps:
(1) selecting a research area range and a target layer, establishing a target layer deposition evolution conceptual model of the research area by means of basic geological data analysis, literature research, current analogy and the like, and quantitatively determining the initial deposition bottom shape form, the structural settlement amount, the sea/lake plane lifting curve, the sediment source and supply rate, the growth period, conditions and rate of the carbonate rock and organic matters.
(2) On the basis of the step (1), carrying out three-dimensional deposition forward modeling numerical simulation on a target layer in a research area by using deposition forward modeling software, comparing a simulation result with actual geological data, if the error is large, readjusting input parameters, repeatedly simulating and correcting until the error of the simulation result reaches a reasonable range, and finally establishing a reasonable three-dimensional deposition modeling model;
(3) selecting a certain grid point to be analyzed, and expressing the grid point by a row number and a column number of the grid point;
(4) for a file with a GRAPH suffix output by a deposition simulation result, sequentially searching the grid point deposition information in each time interval by taking fixed time as a step length from the initial deposition stage of a target layer, wherein the following three conditions exist: a) if the grid point information is retrieved from the GRAPH file in the first time interval and only one line exists, writing the line information into a new text document and starting with the sequence number of '1', wherein the sequence number represents the deposit information in the 1 st time interval; b) if the grid point information is retrieved from the GRAPH file in the first time interval but two lines of information are corresponding, the two lines of information are written into two different text documents respectively, the text corresponding to the line with the positive numerical value is named as 'deposition history', the text corresponding to the line with the negative numerical value is named as 'ablation history', and the texts begin of the text with the sequence number of '1' respectively represents the deposition or ablation information in the 1 st time interval; c) if the grid point information is not retrieved from the GRAPH file in the first time interval, assigning the line of sediment information to be 0 in the newly-built text document, and writing the line of sediment information into the newly-built text document at the beginning of the sequence number 1; after the retrieval of the sediment information in the first time interval is finished, the steps are repeated by taking fixed time as a step length, the sediment information in the 2 nd, 3 rd and 4 … th time intervals is sequentially retrieved, the corresponding time sequence number is taken as the beginning, the information is sequentially written into a newly-built text file, and finally, a text file of the sediment information in different time intervals is generated;
(5) sequentially searching the ancient water depth and offshore distance data of the grid point in different deposition periods by taking fixed time as a step length from the initial deposition stage of a target layer for a file with a GRAPH suffix output by a deposition simulation result, writing the search result information into a newly-built text document at the beginning of a time interval sequence number if the grid information is searched in the GRAPH file in a certain time interval, assigning the ancient water depth and the offshore distance to be 0 if the grid information is not searched, writing the ancient water depth and the offshore distance into the newly-built text document at the beginning of the time interval sequence number, and finally establishing the text file of the ancient water depth and the offshore distance of the deposit deposition in different time intervals;
(6) and establishing a one-dimensional sediment deposition characteristic and water body environment comprehensive analysis graph in a time domain by respectively taking the sediment information, the ancient water depth, the offshore distance as a horizontal axis and the deposition time as a vertical axis, comparing the one-dimensional sediment deposition characteristic and the water body environment comprehensive analysis graph with the depth domain one-dimensional lithologic stratum simulation result, and quantitatively analyzing the deposition evolution process of the lithologic sediment at different grid point positions and the relation between the deposition evolution process and the water body environment.
2. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein: the basic geological data in the step (1) comprise seismic data, well drilling data, well logging data, core analysis data, various test data, sedimentary facies data, stratum thickness, sand shale content, previous research result reports and the like.
3. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the deposition forward modeling software is deposition field professional forward modeling software which takes a simplified hydrodynamic momentum equation and a continuous equation as core algorithms and considers various geological action processes;
the simplified hydrodynamic momentum equation expression is as follows:
wherein the content of the first and second substances,
comprises the following steps: differentiation of the parameter q;
q is the fluid velocity vector, m/s;
t is time, s;
Φ is the ratio of pressure and fluid density, i.e., Φ is p/ρ, where p is in Pa and ρ is in kg/m3；
v is the kinematic viscosity, i.e. ν ═ μ/ρ, where μ is the fluid viscosity in pa.s and ρ is the fluid density in kg/m3；
g is the acceleration of gravity, m/s2。
The simplified equation of continuity expression is as follows:
wherein the content of the first and second substances,
is the gradient operator;
q is the fluid velocity vector, m/s.
4. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the multiple geological processes include: the method comprises the following steps of sediment erosion, carrying and sedimentation processes, wave, tide and storm wave action, isopycnic flow action, aeolian sedimentation action, slumping gravity flow action, carbonate rock and organic matter growth process, carbonate rock diagenesis action, tectonic lifting action, sea level change, compaction action after sedimentation and crustal equilibrium sedimentation action.
The actual geological data refer to two-dimensional seismic profiles, field outcrop data, the sedimentary thickness of the current target layer, the sandstone percentage content, sedimentary facies, the thickness of the single-well target layer, lithology and logging data.
5. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the fixed time step in the step (4) corresponds to the display time interval data in the deposition simulation input parameters, and the user can define the data.
6. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the sediment information in the step (4) comprises the respective deposition thickness information of 8 types of sediments including 4 types of particle-size clastic rocks (coarse sandstone, medium sandstone, fine sandstone and mudstone), 2 types of carbonate rocks (carbonate rock 1 and carbonate rock 2) and 2 types of organic matters (organic matters 1 and organic matters 2), and the total thickness information of the sediments.
7. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the data No. 4 to No. 11 in the GRAPH file in the step (4) respectively correspond to: and the sediment thickness data of the coarse sandstone, the medium sand debris rock, the fine sandstone, the mudstone, the carbonate rock 1, the carbonate rock 2, the organic matter 1 and the organic matter 2.
8. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
and (4) sequentially searching for 'n' in the deposition information in the 2 nd, 3 rd and 4 … th time intervals in the step (4), wherein the 'n' represents the final time interval obtained by dividing the total deposition time of the target layer according to the fixed time step.
9. The method for quantitatively analyzing the sedimentation characteristics of one-dimensional sediments and the water body environment in the time domain as claimed in claim 1, wherein:
the deposition forward modeling software is software Sedsim.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910769556.XA CN110688728B (en) | 2019-08-20 | 2019-08-20 | Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910769556.XA CN110688728B (en) | 2019-08-20 | 2019-08-20 | Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110688728A true CN110688728A (en) | 2020-01-14 |
| CN110688728B CN110688728B (en) | 2022-08-26 |
Family
ID=69108578
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910769556.XA Active CN110688728B (en) | 2019-08-20 | 2019-08-20 | Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110688728B (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112100878A (en) * | 2020-08-18 | 2020-12-18 | 重庆地质矿产研究院 | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation |
| CN117195511A (en) * | 2023-08-23 | 2023-12-08 | 中国科学院南海海洋研究所 | Quantitative calculation method for initial crust thickness and expansion coefficient |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104932031A (en) * | 2014-03-19 | 2015-09-23 | 中国石油化工股份有限公司 | Paleo-water-depth quantitative calculation method aiming at lake facies deposition |
-
2019
- 2019-08-20 CN CN201910769556.XA patent/CN110688728B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104932031A (en) * | 2014-03-19 | 2015-09-23 | 中国石油化工股份有限公司 | Paleo-water-depth quantitative calculation method aiming at lake facies deposition |
Non-Patent Citations (2)
| Title |
|---|
| 倪斌等: "塔里木盆地玉北地区埋藏史及热史分析", 《现代地质》 * |
| 张文彪等: "综合沉积正演与多点地质统计模拟碳酸盐岩台地――以巴西Jupiter油田为例", 《石油学报》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112100878A (en) * | 2020-08-18 | 2020-12-18 | 重庆地质矿产研究院 | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation |
| CN117195511A (en) * | 2023-08-23 | 2023-12-08 | 中国科学院南海海洋研究所 | Quantitative calculation method for initial crust thickness and expansion coefficient |
| CN117195511B (en) * | 2023-08-23 | 2024-04-30 | 中国科学院南海海洋研究所 | Quantitative calculation method for initial crust thickness and expansion coefficient |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110688728B (en) | 2022-08-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bhattacharya et al. | Estimation of source area, river paleo-discharge, paleoslope, and sediment budgets of linked deep-time depositional systems and implications for hydrocarbon potential | |
| Steel et al. | Clinoforms, clinoform trajectories and deepwater sands | |
| US7742875B2 (en) | Method for geologic modeling through hydrodynamics-based gridding (Hydro-Grids) | |
| CA2692425C (en) | Method, program and computer system for scaling hydrocarbon reservoir model data | |
| CN110609319B (en) | Method for establishing time domain sequence stratigraphic section by using forward modeling of sedimentation | |
| CN102243678A (en) | Method for analyzing sand bodies in reservoirs based on inversion technique of sedimentary dynamics | |
| Burton et al. | Geologically-based permeability anisotropy estimates for tidally-influenced reservoirs using quantitative shale data | |
| CN105467464A (en) | Dividing and contrasting method of near-shore underwater fan glutenite effective connected bodies | |
| CN110688728B (en) | Method for quantitatively analyzing sedimentation characteristics of one-dimensional sediments in time domain and water environment | |
| Steventon et al. | Evolution of a sand-rich submarine channel–lobe system, and the impact of mass-transport and transitional-flow deposits on reservoir heterogeneity: Magnus Field, Northern North Sea | |
| Abdelwahhab et al. | 3D-static reservoir and basin modeling of a lacustrine fan-deltaic system in the Gulf of Suez, Egypt | |
| Abdelwahhab et al. | Machine learning-supported seismic stratigraphy of the Paleozoic Nubia Formation (SW Gulf of Suez-rift): implications for paleoenvironment− petroleum geology of a lacustrine-fan delta | |
| Shabani et al. | 3D forward stratigraphic modeling of the Albian succession in a part of the northeastern margin of the Arabian Plate and its implications for exploration of subtle traps | |
| Trippetta et al. | Carbonate-ramp reservoirs modelling best solutions: Insights from a dense shallow well database in Central Italy | |
| Hampson et al. | Quantitative analysis of net‐transgressive shoreline trajectories and stratigraphic architectures: mid‐to‐late Jurassic of the North Sea Rift Basin | |
| Whitaker et al. | Simulating flow through forward sediment model stratigraphies: insights into climatic control of reservoir quality in isolated carbonate platforms | |
| CUI et al. | A volumetric model for evaluating tight sandstone gas reserves in the Permian Sulige gas field, Ordos Basin, Central China | |
| Chudi et al. | Reservoir quality prediction of deep-water Oligocene sandstones from the west Niger Delta by integrated petrological, petrophysical and basin modelling | |
| CN114482995A (en) | Fine determination method for argillaceous content of fine-grain sediment | |
| CN109917489B (en) | Novel method for determining underground pressure-bearing water level | |
| Andresen | The Norwegian–Danish Basin: a dynamic setting for subsurface sand remobilization–established concepts on distribution and controlling factors | |
| Song et al. | A workflow for turbidite reservoir characterization—a case study of the Macedon member, Northern Carnarvon Basin, NW Australia | |
| CN113341480B (en) | Frequency-division RGB slice and multi-attribute fusion sandy hydrate reservoir prediction method | |
| Ayodele | An integrated study of the early cretaceous (Valanginian) reservoir from the Gamtoos Basin, offshore South Africa with special reference to seismic cacies, formation evaluation and static reservoir modeling | |
| Sagen | Spatial reservoir characterisation with focus on distribution of porosity and permeability properties:-A study of the Tarbert and Ness Formations at the Gullfaks Field |
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 |