CN111337974B - Method and device for determining fault closure - Google Patents

Method and device for determining fault closure Download PDF

Info

Publication number
CN111337974B
CN111337974B CN202010096204.5A CN202010096204A CN111337974B CN 111337974 B CN111337974 B CN 111337974B CN 202010096204 A CN202010096204 A CN 202010096204A CN 111337974 B CN111337974 B CN 111337974B
Authority
CN
China
Prior art keywords
fault
dip angle
effective
mud ratio
target interval
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202010096204.5A
Other languages
Chinese (zh)
Other versions
CN111337974A (en
Inventor
苏玉平
卫平生
陈广坡
石兰亭
景紫岩
邓松涛
李娟�
谢明贤
马丽娜
张斌
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202010096204.5A priority Critical patent/CN111337974B/en
Publication of CN111337974A publication Critical patent/CN111337974A/en
Application granted granted Critical
Publication of CN111337974B publication Critical patent/CN111337974B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/30Analysis
    • G01V1/306Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

Abstract

The invention provides a fault closure determination method and a fault closure determination device, wherein the method comprises the following steps: acquiring seismic data of a target interval; determining the position of a spatial projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the spatial projection window; obtaining a fault dip angle value of each calculation point according to the seismic data; obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point; and determining the fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR. According to the method, after the effective fault mud ratio ESGR is obtained, influence factors brought by the change of the fault dip angle are introduced, so that the subsequent fault rock mud content is changed by the fault dip angle, the method is more in line with the actual situation, and the accuracy of the fault closure determination result is improved.

Description

Method and device for determining fault closure
Technical Field
The invention relates to the technical field of oil and gas field exploration, in particular to a fault closure determination method and device.
Background
Fault fractures are filled with fault rocks with uneven distribution such as thickness and lithology, the lateral sealing capacity of the fault is not determined by the butt joint condition of two lithologies, but is determined by the displacement pressure difference between the fault rocks and a target reservoir, and for a clastic rock stratum, the displacement pressure of the fault rocks is mainly determined by the mud content of the fault rocks, so that the larger the mud content of the fault rocks in the clastic rock stratum is, the stronger the lateral sealing capacity of the fault is; the weaker the opposite. Therefore, researchers have attempted to evaluate the ability of a fault to laterally trap hydrocarbons by studying the mudcontent of the fault rock.
The muddiness in the fault rocks are mainly mudstone smear and sand-mud mixed filler of shale layers. The two methods for evaluating the mudstone smearing development degree in the prior art mainly comprise a mudstone smearing factor SSF and a mudstone smearing potential CSP, but the two methods only consider the fault distance and the thickness of the mudstone layer, so that the evaluation on the mudstone smearing development degree is not accurate, and the final fault closure determination result is not accurate.
Disclosure of Invention
The embodiment of the invention provides a fault closure determination method, which is used for improving the accuracy of a fault closure determination result and comprises the following steps:
acquiring seismic data of a target interval;
determining the position of a space projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the space projection window;
obtaining a fault dip angle value of each calculation point according to the seismic data;
obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point;
determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
obtaining a fault dip correction factor of each calculation point according to the fault dip value of each calculation point, wherein the fault dip correction factor comprises the following steps: and obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point according to the following formula:
Figure GDA0003481963550000021
wherein f (θ) represents a fault dip correction factor; theta represents a fault dip angle value;
determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR, wherein the fault closure of the target interval comprises the following steps:
obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value to determine fault closure of the target interval; and the lower threshold of the effective fault mud ratio correction value is obtained according to historical fault logging data statistics.
The embodiment of the invention also provides a fault sealing performance determining device, which is used for improving the accuracy of the fault sealing performance determining result, and comprises the following components:
the data acquisition module is used for acquiring seismic data of a target interval;
the calculation domain selection module is used for determining the position of a space projection window of the target interval on the fault plane according to the seismic data and selecting a calculation point on the fault plane on the position of the space projection window;
the fault dip angle determining module is used for obtaining a fault dip angle value of each calculation point according to the seismic data;
the correction factor calculation module is used for obtaining a fault dip correction factor of each calculation point according to the fault dip value of each calculation point;
the sealing performance determining module is used for determining the sealing performance of the fault of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
wherein the correction factor calculation module is specifically configured to:
and obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point according to the following formula:
Figure GDA0003481963550000022
wherein f (θ) represents a fault dip correction factor; theta represents a fault dip angle value;
wherein the closure determination module is specifically configured to:
obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value to determine fault closure of the target interval; and the lower threshold of the effective fault mud ratio correction value is obtained according to historical fault logging data statistics.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the fault closure determination method when executing the computer program.
Embodiments of the present invention also provide a computer-readable storage medium storing a computer program for executing the fault seal determination method described above.
In the embodiment of the invention, the seismic data of the target interval is obtained; determining the position of a spatial projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the spatial projection window; obtaining a fault dip angle value of each calculation point according to the seismic data; obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point; determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR; after the effective fault mud ratio ESGR is obtained, influence factors caused by changes of fault dip angles are introduced, so that subsequent fault rock mud content is changed by the fault dip angles, actual conditions are met, and the accuracy of a fault closure determination result is improved.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a schematic diagram of a fault seal determination method in an embodiment of the present invention.
Fig. 2 is a schematic diagram of a method for implementing step 105 in an embodiment of the present invention.
FIG. 3 is a schematic diagram illustrating a variation of a fault dip correction factor according to an embodiment of the present invention.
FIG. 4 is a schematic diagram illustrating the relationship between dip angle variation and fault seal in an embodiment of the present invention.
Fig. 5 is a block diagram showing a structure of a fault closure determining apparatus according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
The applicant finds that neither of the two methods for evaluating the mudstone smearing development degree in the prior art takes into account the objective phenomenon that the dip angles of different positions of the fault are changed. In the prior art, for the research of sand-mud mixed fillings in faults, an effective fault mud ratio ESGR method is mainly adopted, namely, the ratio of the accumulated thickness of a mud-rock layer at a certain point, which slides through the fault in the process of fault movement, to the fault distance, the contribution of mud (silicate) in sandstone to the ratio of the mud in the fault zone in the actual fracture action process is considered, and the ESGR can more accurately predict the mud content in the fault zone and the smearing degree of the mud-rock, but the method does not consider the condition that the dip angles of different positions of the same fault are different, particularly the fault with obvious change of dip angle occurrence of a plow-type fault.
The determination of the shale content of the fault rock mainly depends on the shale smearing of the shale layer and the determination of the sand-mud mixed filling, but because the method does not consider the change of the fault dip angle, the influence factors brought to the shale content calculation of the fault rock are not considered, and the method has limitation, so that the final fault closure determination result is inaccurate.
In order to solve the problem that the fault blocking performance determination result is inaccurate, an embodiment of the present invention provides a fault blocking performance determination method, which is used for improving the accuracy of the fault blocking performance determination result, as shown in fig. 1, and the method includes:
step 101: acquiring seismic data of a target interval;
step 102: determining the position of a spatial projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the spatial projection window;
step 103: obtaining a fault dip angle value of each calculation point according to the seismic data;
step 104: obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point;
step 105: and determining the fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR.
As can be known from the flow shown in FIG. 1, in the embodiment of the present invention, the seismic data of the target interval is obtained; determining the position of a spatial projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the spatial projection window; obtaining a fault dip angle value of each calculation point according to the seismic data; obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point; determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR; after the effective fault mud ratio ESGR is obtained, influence factors caused by the change of the fault dip angle are introduced, so that the subsequent fault rock mud content is changed by the fault dip angle, the actual situation is better met, and the accuracy of the fault closure determination result is improved.
When the method is specifically implemented, firstly, the seismic data of the target interval is obtained, and the seismic data of the target interval is obtained from the data obtained by petroleum seismic exploration.
After seismic data of the target interval are obtained, the position of a space projection window of the target interval on the fault plane is determined, and a calculation point is selected on the fault plane on the position of the space projection window. During specific implementation, the fault plane on the position of the space projection window is subjected to gridding division according to the preset grid size, and a plurality of grids obtained through division are used as selected calculation points. For example, the size of the selected spatial projection window is 100m × 100m, the spatial projection window may be divided into 10000 grids of 1m × 1m, and the 10000 grids obtained by the division are used as the selected computation points.
And after the selection of the calculation points is finished, obtaining the fault dip angle value of each calculation point according to the seismic data. During specific implementation, an underground fault model is established according to seismic data, and a fault dip angle value of each calculation point on a section and a vertical fault distance of the fault are calculated.
Obtaining a fault dip correction factor of each calculation point according to the fault dip value of each calculation point, and specifically, calculating according to the following formula:
Figure GDA0003481963550000051
wherein f (θ) represents a fault dip correction factor; θ represents a fault dip value.
And after the fault dip angle correction factor of each calculation point is obtained, determining the fault closure of the target layer section according to the fault dip angle correction factor and the effective fault mud ratio ESGR. The specific process is shown in fig. 2, and includes:
step 201: obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
step 202: and comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value to determine the fault closure of the target interval.
Wherein, the effective fault mud ratio ESGR in step 201 is calculated according to the following formula:
Figure GDA0003481963550000052
wherein i represents the number of formations that have slipped past the calculation point; t isiRepresents the thickness, m, of the ith layer of rock sliding over the computation point; vsh(i)Denotes the shale content,%, of the i-th layer sliding over the calculation point; d represents the vertical fault distance of the fault, m.
The above parameter TiAnd Vsh(i)Can be obtained according to the analysis of seismic data and well drilling data.
After the effective fault mud ratio ESGR is obtained, an effective fault mud ratio correction value DSGR is obtained as f (θ) × ESGR.
And the lower threshold of the effective fault mud ratio correction value is obtained according to the statistics of historical fault logging data. The discovered reservoir drilling information associated with faults in the interval of interest over the area of interest is counted, and then a lower threshold of valid fault mud ratio corrections for faults in the counted interval of interest that have been identified as closed is determined.
When the step 202 is implemented, the method includes:
if the effective fault mud ratio correction value DSGR is larger than the lower limit threshold, the fault of the target interval is laterally closed, and the closing capacity is determined according to the difference between the effective fault mud ratio correction value DSGR and the lower limit threshold; wherein, the larger the difference between the DSGR and the lower threshold value is, the stronger the fault sealing capability is.
And if the effective fault mud ratio correction value DSGR is smaller than the lower limit threshold, the fault of the target interval is laterally opened.
A specific example is given below to illustrate how embodiments of the present invention perform fault seal determination.
(1) Establishing an underground fault model according to data obtained by petroleum seismic exploration;
(2) obtaining a stratum model according to petroleum seismic exploration and well drilling data, coupling the stratum model with a fault model, then determining the position of a spatial projection window of an evaluated target interval on a section, and selecting a calculation point to be subjected to closure evaluation on the section;
(3) and calculating the fault dip angle value and the vertical fault distance of each calculation point according to the seismic data, and calculating the fault dip angle correction factor of each calculation point on the section. As shown in fig. 3, the fault inclination correction factor f (θ) gradually decreases as the fault inclination increases. When the fault dip angle is changed between 0-90 degrees, the range corresponding to f (theta) is 1-0.841.
(4) The shale content ESGR value of the effective fault rock of each calculation point is obtained, and the T required in the formula is determined according to petroleum seismic exploration and well drilling data in the stepiAnd Vsh(i)The required fault vertical offset D in the formula has been obtained in step (3).
(5) And applying the fault dip angle correction factor to ESGR correction and solving a DSGR value. The DSGR value considering the fault dip angle factor is smaller than the original ESGR value on the whole; meanwhile, the smaller the fault dip angle is, the larger the corresponding fault dip angle correction factor f (theta) is, the larger the corresponding DSGR is. This is because the flatter the fracture surface (the smaller the fault dip), the greater the positive pressure of the fracture surface due to the gravity of the upper tray, the stronger the crushing and grinding action on the filler in the fracture zone, the more the filler component of the produced fine particles, the higher the shale content of the fault rock, the stronger the fault blocking capability reflected by the DSGR, and vice versa. As shown in FIG. 4, the method can quantitatively explain the relationship between the inclination angle change and the sealing property of the plow type normal fault with the steep upper part (the fault inclination angle is theta 1) and the gentle lower part (the fault inclination angle is theta 2) well, so that the effect brought by the fault sealing property determining method provided by the invention is more practical and meets the geological knowledge.
(6) And (4) carrying out quantitative judgment on the lower threshold and the closure of the statistical effective fault mud ratio correction value. In this step, the discovered reservoir well drilling information associated with faults in the interval of interest in the study area is first counted, then a lower threshold of effective fault mud ratio correction values for faults in the counted interval of interest that have been identified as closed is determined, and the computed value of DSGR in (5) is compared with the lower threshold: if the DSGR calculated value is larger than the closed lower limit threshold value, the fault is laterally closed, otherwise, the fault is laterally opened; the size of the sealing capacity of the fault.
After the fault closure determination results of all the calculation points are obtained, the general situation of the target interval can be analyzed according to the fault closure determination results of all the calculation points, and the subsequent industrial production is guided.
The implementation of the above specific application is only an example, and the rest of the embodiments are not described in detail.
Based on the same inventive concept, embodiments of the present invention further provide a fault seal determination apparatus, and because the principle of the problem solved by the fault seal determination apparatus is similar to that of the fault seal determination method, the implementation of the fault seal determination apparatus may refer to the implementation of the fault seal determination method, and repeated parts are not repeated, and the specific structure is as shown in fig. 5:
the data acquisition module 501 is used for acquiring seismic data of a target interval;
a calculation domain selection module 502, configured to determine, according to the seismic data, a spatial projection window position of the target interval on the fault plane, and select a calculation point on the fault plane at the spatial projection window position;
a fault dip angle determination module 503, configured to obtain a fault dip angle value of each calculation point according to the seismic data;
a correction factor calculation module 504, configured to obtain a fault dip correction factor for each calculation point according to the fault dip value of each calculation point;
and the sealing performance determining module 505 is used for determining the fault sealing performance of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR.
In a specific embodiment, the calculation domain selecting module 502 is specifically configured to:
and carrying out gridding division on the fault plane at the position of the space projection window according to the preset grid size, wherein a plurality of grids obtained by division are used as selected calculation points.
In the embodiment of the present invention, the correction factor calculation module 504 is specifically configured to:
and obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point according to the following formula:
Figure GDA0003481963550000081
wherein f (θ) represents a fault dip correction factor; θ represents a fault dip value.
In this embodiment of the present invention, the closure determining module 505 is specifically configured to:
obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value, and determining the fault closure of the target interval; and the lower threshold of the effective fault mud ratio correction value is obtained according to historical fault logging data statistics.
Specifically, the closure determination module 505 is specifically configured to:
if the effective fault mud ratio correction value is larger than the lower limit threshold, the fault of the target interval is laterally closed, and the closing capacity is determined according to the difference between the effective fault mud ratio correction value and the lower limit threshold;
and if the effective fault mud ratio correction value is smaller than the lower limit threshold, the fault of the target interval is laterally opened.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the fault closure determination method when executing the computer program.
An embodiment of the present invention further provides a computer-readable storage medium storing a computer program for executing the fault seal determination method.
In summary, the fault sealing determination method and device provided by the embodiment of the invention have the following advantages:
acquiring seismic data of a target interval; determining the position of a spatial projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the spatial projection window; obtaining a fault dip angle value of each calculation point according to the seismic data; obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point; determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR; after the effective fault mud ratio ESGR is obtained, influence factors caused by the change of the fault dip angle are introduced, so that the subsequent fault rock mud content is changed by the fault dip angle, the actual situation is better met, and the accuracy of the fault closure determination result is improved.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, apparatus, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above description is only a preferred embodiment of the present invention, and is not intended to limit the present invention, and various modifications and changes may be made to the embodiment of the present invention by those skilled in the art. Any modification, equivalent replacement, or improvement made within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (8)

1. A fault seal determination method, comprising:
acquiring seismic data of a target interval;
determining the position of a space projection window of the target interval on the fault plane according to the seismic data, and selecting a calculation point on the fault plane on the position of the space projection window;
obtaining a fault dip angle value of each calculation point according to the seismic data;
obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point;
determining the fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
obtaining a fault dip correction factor of each calculation point according to the fault dip value of each calculation point, wherein the fault dip correction factor comprises the following steps: and obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point according to the following formula:
Figure FDA0003481963540000011
wherein f (θ) represents a fault dip correction factor; theta represents a fault dip angle value;
determining fault closure of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR, wherein the fault closure of the target interval comprises the following steps:
obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value to determine fault closure of the target interval; and the lower threshold of the effective fault mud ratio correction value is obtained according to historical fault logging data statistics.
2. The method of claim 1, wherein selecting a computed point on a fault plane at the location of the spatial projection window comprises:
and carrying out gridding division on the fault plane on the position of the space projection window according to a preset grid size, wherein a plurality of grids obtained by division are used as selected calculation points.
3. The method of claim 1, wherein comparing the effective fault mud ratio correction value to a lower threshold of effective fault mud ratio correction values to determine fault seal for the interval of interest comprises:
if the effective fault mud ratio correction value is larger than the lower threshold, the fault of the target interval is laterally closed, and the closing capacity is determined according to the difference between the effective fault mud ratio correction value and the lower threshold;
and if the effective fault mud ratio correction value is smaller than the lower threshold, the fault of the target interval is laterally opened.
4. A fault seal determination apparatus, comprising:
the data acquisition module is used for acquiring seismic data of a target interval;
the calculation domain selection module is used for determining the position of a space projection window of the target interval on the fault plane according to the seismic data and selecting a calculation point on the fault plane on the position of the space projection window;
the fault dip angle determining module is used for obtaining a fault dip angle value of each calculation point according to the seismic data;
the correction factor calculation module is used for obtaining a fault dip correction factor of each calculation point according to the fault dip value of each calculation point;
the sealing performance determining module is used for determining the sealing performance of the fault of the target interval according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
wherein the correction factor calculation module is specifically configured to:
and obtaining a fault dip angle correction factor of each calculation point according to the fault dip angle value of each calculation point according to the following formula:
Figure FDA0003481963540000021
wherein f (θ) represents a fault dip correction factor; theta represents a fault dip angle value;
wherein the closure determination module is specifically configured to:
obtaining an effective fault mud ratio correction value according to the fault dip angle correction factor and the effective fault mud ratio ESGR;
comparing the effective fault mud ratio correction value with a lower limit threshold of the effective fault mud ratio correction value to determine fault closure of the target interval; and the lower threshold of the effective fault mud ratio correction value is obtained according to historical fault logging data statistics.
5. The apparatus of claim 4, wherein the calculation domain selection module is specifically configured to:
and carrying out gridding division on the fault plane on the position of the space projection window according to a preset grid size, wherein a plurality of grids obtained by division are used as selected calculation points.
6. The apparatus of claim 4, wherein the seal determination module is specifically configured to:
if the effective fault mud ratio correction value is larger than the lower threshold, the fault of the target interval is laterally closed, and the closing capacity is determined according to the difference between the effective fault mud ratio correction value and the lower threshold;
and if the effective fault mud ratio correction value is smaller than the lower threshold, the fault of the target interval is laterally opened.
7. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 3 when executing the computer program.
8. A computer-readable storage medium, characterized in that it stores a computer program for executing the method of any one of claims 1 to 3.
CN202010096204.5A 2020-02-17 2020-02-17 Method and device for determining fault closure Active CN111337974B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010096204.5A CN111337974B (en) 2020-02-17 2020-02-17 Method and device for determining fault closure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010096204.5A CN111337974B (en) 2020-02-17 2020-02-17 Method and device for determining fault closure

Publications (2)

Publication Number Publication Date
CN111337974A CN111337974A (en) 2020-06-26
CN111337974B true CN111337974B (en) 2022-05-06

Family

ID=71181565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010096204.5A Active CN111337974B (en) 2020-02-17 2020-02-17 Method and device for determining fault closure

Country Status (1)

Country Link
CN (1) CN111337974B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112083144B (en) * 2020-09-01 2021-10-15 中国科学院地质与地球物理研究所 Fault on-off prediction method and device, computer equipment and storage medium
CN113466931B (en) * 2021-05-31 2023-08-22 中国石油天然气股份有限公司 Geological period fault closure evaluation method and system

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
MY125603A (en) * 2000-02-25 2006-08-30 Shell Int Research Processing seismic data
GB0900687D0 (en) * 2009-01-16 2009-02-25 Rock Deformation Res Ltd Fault analysis system
CN106324676B (en) * 2016-10-21 2018-09-04 中国石油天然气股份有限公司 A kind of method and device of determining fault sealing property
CN106873030B (en) * 2017-02-07 2019-05-07 中国石油天然气股份有限公司 The determination method and apparatus of Fault Seal
CN108343430A (en) * 2017-07-21 2018-07-31 中国石油大学胜利学院 Fault sealing property fuzzy mathematics comprehensive quantitative evaluation method
CN108241181B (en) * 2018-01-30 2019-07-09 中国石油化工股份有限公司 A kind of evaluation method of fault sealing property
CN110632656B (en) * 2019-08-16 2021-08-03 中国石油天然气股份有限公司 Method and device for analyzing transverse sealing property of fault
CN110632655B (en) * 2019-08-16 2021-08-03 中国石油天然气股份有限公司 Method and device for analyzing sealing property of heterogeneous fault
CN110632651B (en) * 2019-08-16 2021-09-28 中国石油天然气股份有限公司 Method and device for quantitatively analyzing fault closure

Also Published As

Publication number Publication date
CN111337974A (en) 2020-06-26

Similar Documents

Publication Publication Date Title
CN106526693B (en) Crack identification method and apparatus
CN110632651B (en) Method and device for quantitatively analyzing fault closure
CN103744109B (en) Method for identifying a weathering crust structure of clastic rock in a area of covering no well
CN111337974B (en) Method and device for determining fault closure
US20220291418A1 (en) An integrated geomechanics model for predicting hydrocarbon and migration pathways
CN110632652B (en) Method and device for determining oil-gas migration mode of fault block oil-gas reservoir
CN110632656B (en) Method and device for analyzing transverse sealing property of fault
CN110632657B (en) Mudstone smearing type fault sealing analysis method and device
CN110633557B (en) Identification method for favorable area of coal bed gas structure
CN110632654B (en) Method and device for determining oil-containing boundary of broken block trap
CN112394401A (en) Clastic rock trap hydrocarbon-containing evaluation method based on risk probability
Pontiggia et al. New integrated approach for diagenesis characterization and simulation
CN116930023A (en) Fine interpretation method and device for dense sandstone phase-control classified porosity logging
Rotimi et al. Reservoir characterization and modeling of lateral heterogeneity using multivariate analysis
CN114779333A (en) Method and system for determining location of hydrocarbon reservoir
CN113534248A (en) Method, device and system for quantitatively analyzing closure of trap fault
CN113031056B (en) Fault closure analysis method and device under construction constraint
CN113466931B (en) Geological period fault closure evaluation method and system
Faraji et al. Estimation and mapping of vitrinite reflectance from seismic data in South Pars gas field, Persian Gulf
CN112147677A (en) Oil and gas reservoir parameter tag data generation method and device
CN106405638A (en) Method and apparatus for determining gas saturation
US10460051B2 (en) Computationally-efficient modeling of viscous fingering effect for enhanced oil recovery (EOR) agent injected at multiple injection concentrations
CN108318936A (en) A kind of STRATIGRAPHIC DIVISION treating method and apparatus
CN114139328B (en) Prediction method for favorable ore-forming zone of sandstone-type uranium ore in interlayer oxidation zone
Vu et al. 3D dynamic fault sealing capacity modelling to improve history matching: a case study in Oligocene reservoir, Tay Ho Field, Blocks A, Cuu Long basin, Offshore Vietnam

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