CN111650637A - Seismic horizon interpretation method and device - Google Patents

Seismic horizon interpretation method and device Download PDF

Info

Publication number
CN111650637A
CN111650637A CN201910161351.3A CN201910161351A CN111650637A CN 111650637 A CN111650637 A CN 111650637A CN 201910161351 A CN201910161351 A CN 201910161351A CN 111650637 A CN111650637 A CN 111650637A
Authority
CN
China
Prior art keywords
seismic
point
target
horizon
determining
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
Application number
CN201910161351.3A
Other languages
Chinese (zh)
Other versions
CN111650637B (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
Daqing Oilfield Co Ltd
Original Assignee
Petrochina Co Ltd
Daqing Oilfield 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, Daqing Oilfield Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201910161351.3A priority Critical patent/CN111650637B/en
Publication of CN111650637A publication Critical patent/CN111650637A/en
Application granted granted Critical
Publication of CN111650637B publication Critical patent/CN111650637B/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

Abstract

The application discloses a seismic horizon interpretation method and device, and belongs to the technical field of seismic data processing. The method comprises the following steps: selecting a set of seismic data bodies from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, wherein a plurality of seismic sections in each set of seismic data bodies are in one-to-one correspondence; tracking a first homophasic axis and a second homophasic axis of a first seismic section of a reference data volume to respectively obtain a target layer horizon and a standard layer horizon of the first seismic section; and determining the target horizon of a second seismic profile of the target data volume according to the target horizon and the standard horizon of the first seismic profile, wherein the target data volume is any one of the seismic data volumes except the reference data volume in the plurality of sets of seismic data volumes, and the second seismic profile is the seismic profile corresponding to the first seismic profile in the plurality of seismic profiles of the target data volume. The method and the device can improve the seismic horizon interpretation efficiency and the interpretation precision.

Description

Seismic horizon interpretation method and device
Technical Field
The application relates to the technical field of seismic data processing, in particular to a seismic horizon interpretation method and device.
Background
Seismic wave signals reflected back to the ground in a region to be researched of a stratum can be collected through an OVT (Offset Vector Tile) domain pre-stack seismic reservoir prediction technology, and a plurality of sets of seismic data bodies can be obtained after the collected seismic wave signals are processed. In order to clearly understand the reservoir physical property and lithology changes of the layer to be researched, namely the target layer in the area to be researched, the target layer of each seismic data body in the multiple seismic data bodies needs to be sequentially subjected to horizon interpretation.
Currently, the horizons of interest for these sets of seismic data volumes are often interpreted in two ways. The first mode is as follows: and performing target layer horizon interpretation on each seismic data body in the multiple sets of seismic data bodies one by one according to a conventional seismic horizon interpretation method to obtain target layer horizon interpretation results of the multiple sets of seismic data bodies. The second mode is as follows: and performing target layer horizon interpretation on one set of seismic data bodies in the multiple sets of seismic data bodies according to a conventional seismic horizon interpretation method, and taking the target layer horizon interpretation result of the set of seismic data bodies as the target layer horizon interpretation result of the other seismic data bodies in the multiple sets of seismic data bodies.
However, the first method described above is very labor intensive and time consuming to perform horizon interpretation on multiple sets of seismic data volumes. When the second mode is adopted, because the plurality of sets of seismic data volumes have differences, when the layer position interpretation result of one set of seismic data volume is directly used as the layer position interpretation result of the other seismic data volumes, the interpretation precision is low for the other seismic data volumes, and even a local layer crossing situation can occur.
Disclosure of Invention
The embodiment of the application provides a seismic horizon interpretation method and device, and the problems of low seismic horizon interpretation efficiency and low interpretation precision in the related art can be solved. The technical scheme is as follows:
in a first aspect, a seismic horizon interpretation method is provided, the method comprising:
selecting a set of seismic data bodies from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, wherein each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data bodies are in one-to-one correspondence;
tracking a first homodyne axis of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, and tracking a second homodyne axis of the first seismic section to obtain a standard horizon of the first seismic section, wherein the first seismic section is any one of a plurality of seismic sections of the reference data volume;
and determining the target horizon of a second seismic profile of the target data volume according to the target horizon and the standard horizon of the first seismic profile, wherein the target data volume is any one of the plurality of sets of seismic data volumes except the reference data volume, and the second seismic profile is a seismic profile corresponding to the first seismic profile in the plurality of seismic profiles of the target data volume.
Optionally, each seismic section of the plurality of seismic sections comprises a plurality of seismic traces;
the tracking a first co-axis of the first seismic section to obtain a target horizon of the first seismic section comprises:
obtaining a plurality of first reference points on a first co-axial of the first seismic section;
determining amplitude maximum points corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
Optionally, the tracking a second in-phase axis of the first seismic section to obtain a standard layer horizon of the first seismic section includes:
obtaining a plurality of second reference points on a second event axis of the first seismic section;
determining amplitude maximum points corresponding to each second reference point in a second time window on the plurality of seismic traces according to the plurality of second reference points;
and determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
Optionally, the determining a target layer horizon of a second seismic section of the target data volume according to the target layer horizon and the standard layer horizon of the first seismic section includes:
determining a standard layer position of the second seismic section according to the standard layer position of the first seismic section;
determining a difference value between each point on the standard layer level of the first seismic profile and a corresponding point on the standard layer level of the second seismic profile according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile, and correspondingly storing the coordinate of each point on the standard layer level of the second seismic profile and the determined corresponding difference value to obtain a corresponding relation between the coordinate and the difference value;
and determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference value.
Optionally, the determining the standard layer horizon of the second seismic section according to the standard layer horizon of the first seismic section includes:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining an amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
Optionally, the determining, according to the standard layer level of the first seismic section and the standard layer level of the second seismic section, a difference between each point on the standard layer of the first seismic section and a corresponding point on the standard layer level of the second seismic section includes:
acquiring coordinates of a first target point on a standard layer horizon of the first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on a standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
determining a time difference between a first time of the first target point and a second time of the second target point, and determining the time difference as a difference between the first target point and the second target point.
Optionally, the determining a target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relationship between the coordinate and the difference includes:
acquiring coordinates of each point on a target horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
In a second aspect, there is provided a seismic horizon interpretation apparatus comprising:
the system comprises a selection module, a data analysis module and a data analysis module, wherein the selection module is used for selecting a set of seismic data body from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data body correspond to one another;
a tracking module, configured to track a first event of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, and track a second event of the first seismic section to obtain a standard horizon of the first seismic section, where the first seismic section is any one of a plurality of seismic sections of the reference data volume;
and the determining module is used for determining the target horizon of a second seismic section of the target data volume according to the target horizon and the standard horizon of the first seismic section, wherein the target data volume is any one of the plurality of sets of seismic data volumes except the reference data volume, and the second seismic section is a seismic section corresponding to the first seismic section in the plurality of seismic sections of the target data volume.
Optionally, each seismic section of the plurality of seismic sections comprises a plurality of seismic traces;
the tracking module includes:
a first acquisition submodule for acquiring a plurality of first reference points on a first event axis of the first seismic section;
the first determining submodule is used for determining an amplitude maximum value point corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and the second determining submodule is used for determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
Optionally, the tracking module further comprises:
a second acquisition submodule for acquiring a plurality of second reference points on a second in-phase axis of the first seismic section;
a third determining submodule, configured to determine, according to the plurality of second reference points, an amplitude maximum point corresponding to each second reference point in a second time window on the plurality of seismic traces;
and the fourth determining submodule is used for determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
Optionally, the determining module includes:
the fifth determining submodule is used for determining the standard layer position of the second seismic section according to the standard layer position of the first seismic section;
a sixth determining submodule, configured to determine, according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile, a difference between each point on the standard layer level of the first seismic profile and a corresponding point on the standard layer level of the second seismic profile, and store coordinates of each point on the standard layer level of the second seismic profile and the corresponding difference determined for the corresponding point, so as to obtain a corresponding relationship between the coordinates and the difference;
and the seventh determining submodule is used for determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference values.
Optionally, the fifth determining submodule is configured to:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining an amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
Optionally, the sixth determining submodule is configured to:
acquiring coordinates of a first target point on a standard layer horizon of the first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on a standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
determining a time difference between a first time of the first target point and a second time of the second target point, and determining the time difference as a difference between the first target point and the second target point.
Optionally, the seventh determining sub-module is configured to:
acquiring coordinates of each point on a target horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
In a third aspect, there is provided a seismic horizon interpreting apparatus, the apparatus comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to perform the steps of the seismic horizon interpretation method of the first aspect described above.
In a fourth aspect, a computer readable storage medium is provided, having instructions stored thereon, which when executed by a processor, implement the steps of the seismic horizon interpretation method according to the first aspect described above.
The technical scheme provided by the embodiment of the application can at least bring the following beneficial effects:
according to the method and the device for determining the seismic profile of the target data volume, one set of seismic data volume can be selected from multiple sets of seismic data volumes to serve as a reference data volume, a first homophase axis and a second homophase axis of a first seismic profile of the reference data volume are tracked, a target layer position and a standard layer position of the first seismic profile are obtained respectively, and then the target layer position of a second seismic profile of the target data volume is determined according to the target layer position and the standard layer position of the first seismic profile of the reference data volume. Therefore, only the reference data bodies in the multiple sets of seismic data bodies need to be subjected to horizon interpretation in a conventional mode of in-phase axis tracking, the target layer horizon of any seismic data body except the reference data body in the multiple sets of seismic data bodies can be obtained through the target layer horizon and the standard layer horizon of the reference data body, and the consumption of a large amount of time caused by the conventional in-phase axis tracking of other seismic data bodies except the reference data body for multiple times is avoided. In addition, the method and the device can determine the target layers of other seismic data volumes in the multiple sets of seismic data volumes according to the standard layer position and the target layer position of the reference data volume to perform layer position interpretation, and compared with the method and the device which directly take the layer position interpretation result of one set of seismic data volume as the layer position interpretation result of other seismic data volumes in the related technology, the method and the device improve interpretation precision.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present application, 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 application, 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 representation of horizons of a first layer under study and a second layer under study of a first seismic section of a reference data volume in the related art;
FIG. 2 is a schematic view of a horizon of a first layer under study of a second seismic profile of a target data volume of the related art;
FIG. 3 is a flow chart of a seismic horizon interpreting method according to an embodiment of the present application;
FIG. 4 is a flow chart of another seismic horizon interpretation method provided by embodiments of the present application;
FIG. 5 is a diagram illustrating a comparison of two layers of horizons on a second seismic section of a target data volume according to the method of the present application and two layers of horizons on a second seismic section of a target data volume according to the method of the related art;
FIG. 6 is a schematic structural diagram of a seismic horizon interpreting apparatus according to an embodiment of the present application;
FIG. 7 is a schematic structural diagram of another seismic horizon interpreting apparatus according to an embodiment of the present application.
Detailed Description
To make the objects, technical solutions and advantages of the present application more clear, embodiments of the present application will be described in further detail below with reference to the accompanying drawings.
Before explaining the embodiments of the present application in detail, the application scenarios related to the embodiments of the present application are explained first.
Currently, the horizons of interest for these sets of seismic data volumes are often interpreted in two ways. The first mode is as follows: and performing target layer horizon interpretation on each seismic data body in the multiple sets of seismic data bodies one by one according to a conventional seismic horizon interpretation method to obtain target layer horizon interpretation results of the multiple sets of seismic data bodies. The second mode is as follows: and performing target layer horizon interpretation on one set of seismic data bodies in the multiple sets of seismic data bodies according to a conventional seismic horizon interpretation method, and taking the target layer horizon interpretation result of the set of seismic data bodies as the target layer horizon interpretation result of the other seismic data bodies in the multiple sets of seismic data bodies. However, the first method described above is very labor intensive and time consuming to perform horizon interpretation on multiple sets of seismic data volumes. When the second mode is adopted, because the plurality of sets of seismic data volumes have differences, when the layer position interpretation result of one set of seismic data volume is directly used as the layer position interpretation result of the other seismic data volumes, the interpretation precision is low for the other seismic data volumes, and even a local layer crossing situation can occur.
For example, a first seismic section of a reference data volume is shown in FIG. 1, L1 represents horizons of a first layer to be investigated of the first seismic section of the reference data volume interpreted according to a conventional seismic horizon interpretation method, and L2 represents horizons of a second layer to be investigated of the first seismic section of the reference data volume interpreted according to a conventional seismic horizon interpretation method. In fig. 1, neither L1 nor L2 has a cross-layer phenomenon, which indicates that the level of the first layer to be studied and the level of the second layer to be studied of the first seismic section of the reference data volume interpreted according to the conventional seismic level interpretation method are accurate. Fig. 2 shows the second seismic section corresponding to the first seismic section in the target data volume, and L3 is the horizon of the first layer to be studied of the second seismic section of the target data volume obtained by directly taking the horizon (L1) of the first layer to be studied of the first seismic section of the reference data volume as the horizon of the first layer to be studied of the second seismic section of the target data volume. As can be seen from fig. 2, L3 has a local cross-layer phenomenon on the second seismic section of the target data volume, which indicates that when the horizon interpretation result of one set of seismic data volume is directly used as the horizon interpretation result of other seismic data volumes, the interpretation accuracy is low for other seismic data volumes. Wherein, the channeling refers to the channeling from one in-phase axis to other in-phase axes on the seismic section.
Based on the above, the embodiment of the application provides a seismic horizon interpretation method, which is used for solving the problems of low seismic horizon interpretation efficiency and low interpretation precision in the related art.
FIG. 3 is a flowchart of a seismic horizon interpreting method according to an embodiment of the present application. Referring to fig. 3, the method includes:
step 301: selecting a set of seismic data bodies from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, wherein each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data bodies correspond to one another.
Step 302: and tracking a first homodyne axis of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, and tracking a second homodyne axis of the first seismic section to obtain a standard horizon of the first seismic section, wherein the first seismic section is any one of the plurality of seismic sections of the reference data volume.
Step 303: and determining the target horizon of a second seismic section of the target data volume according to the target horizon and the standard horizon of the first seismic section, wherein the target data volume is any one of the seismic data volumes except the reference data volume in the plurality of sets of seismic data volumes, and the second seismic section is the seismic section corresponding to the first seismic section in the plurality of seismic sections of the target data volume.
In the embodiment of the application, a set of seismic data volume is selected from a plurality of sets of seismic data volumes to serve as a reference data volume, a first homophase axis and a second homophase axis of a first seismic section of the reference data volume are tracked to obtain a target layer position and a standard layer position of the first seismic section respectively, and then the target layer position of a second seismic section of a target data volume is determined according to the target layer position and the standard layer position of the first seismic section of the reference data volume. Therefore, only the reference data bodies in the multiple sets of seismic data bodies need to be subjected to horizon interpretation in a conventional mode of in-phase axis tracking, the target layer horizon of any seismic data body except the reference data body in the multiple sets of seismic data bodies can be obtained through the target layer horizon and the standard layer horizon of the reference data body, and the consumption of a large amount of time caused by the conventional in-phase axis tracking of other seismic data bodies except the reference data body for multiple times is avoided. In addition, the method and the device can determine the target layers of other seismic data volumes in the multiple sets of seismic data volumes according to the standard layer position and the target layer position of the reference data volume to perform layer position interpretation, and compared with the method and the device which directly take the layer position interpretation result of one set of seismic data volume as the layer position interpretation result of other seismic data volumes in the related technology, the method and the device improve interpretation precision.
Optionally, each seismic section of the plurality of seismic sections comprises a plurality of seismic traces;
tracking a first co-axial line of the first seismic section to obtain a target horizon of the first seismic section, comprising:
obtaining a plurality of first reference points on a first event axis of a first seismic section;
determining an amplitude maximum value point corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
Optionally, tracking a second in-phase axis of the first seismic section to obtain a standard horizon of the first seismic section, including:
obtaining a plurality of second reference points on a second event axis of the first seismic section;
determining the amplitude maximum value point corresponding to each second reference point in a second time window on the plurality of seismic traces according to the plurality of second reference points;
and determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
Optionally, determining the target horizon of the second seismic section of the target data volume from the target horizon of the first seismic section and the standard horizon, comprising:
determining a standard layer position of a second seismic section according to the standard layer position of the first seismic section;
determining a difference value between each point on the standard layer position of the first seismic profile and a corresponding point on the standard layer position of the second seismic profile according to the standard layer position of the first seismic profile and the standard layer position of the second seismic profile, and correspondingly storing the coordinate of each point on the standard layer position of the second seismic profile and the determined corresponding difference value to obtain a corresponding relation between the coordinate and the difference value;
and determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference.
Optionally, determining the standard layer horizon of the second seismic section from the standard layer horizon of the first seismic section comprises:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining the amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
Optionally, determining a difference between each point on the standard layer of the first seismic profile and a corresponding point on the standard layer of the second seismic profile according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile comprises:
acquiring coordinates of a first target point on a standard layer horizon of a first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on the standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
a time difference between a first time of the first target point and a second time of the second target point is determined, and the time difference is determined as a difference between the first target point and the second target point.
Optionally, determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relationship between the coordinates and the difference, including:
acquiring coordinates of each point on a target layer horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
All the above optional technical solutions can be combined arbitrarily to form an optional embodiment of the present application, and the present application embodiment is not described in detail again.
For ease of understanding, the seismic horizon interpretation method provided by the embodiment of FIG. 3 is described below in conjunction with FIG. 4. Fig. 4 is a flowchart of a seismic horizon interpreting method according to an embodiment of the present application, where the method may be applied to a terminal, and the terminal may be a mobile phone, a tablet, a computer, or the like. Referring to fig. 4, the method includes:
step 401: selecting a set of seismic data bodies from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, wherein each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data bodies correspond to one another.
It should be noted that the area to be investigated may be any area of the formation that needs to be investigated. The multiple sets of seismic data bodies can be seismic data bodies obtained by processing seismic wave signals which are acquired by OVT domain pre-stack seismic reservoir prediction technology and are reflected to the ground in a region to be researched. After obtaining the sets of seismic data volumes, the user may number the sets of seismic data volumes to distinguish the sets of seismic data volumes. For example, when 4 sets of seismic data volumes are obtained after processing the seismic signals reflected back to the ground from the acquired region to be studied, the 4 sets of seismic data volumes may be numbered to obtain seismic data volume 1, seismic data volume 2, seismic data volume 3, and seismic data volume 4.
In addition, before the plurality of sets of seismic data volumes are interpreted, each set of seismic data volumes may be sectioned with a plurality of transverse sections or a plurality of longitudinal sections, so as to obtain a plurality of transverse seismic sections or a plurality of longitudinal seismic sections of each set of seismic data volumes, that is, each set of seismic data volumes may include a plurality of transverse seismic sections or a plurality of longitudinal sections. To facilitate seismic horizon interpretation, each of the sets of seismic data volumes may be sectioned with the same number of identically located cross sections or longitudinal sections to obtain the same number and the same location of seismic sections from the sets of seismic data volumes, and the plurality of seismic sections in each set of seismic data volumes may be numbered. For example, 4 sets of seismic data volumes are obtained by processing the seismic signals reflected back to the ground from the acquired region to be studied, and if each set of seismic data volume is sectioned through a transverse profile, the 4 sets of seismic data volumes may include 3 transverse seismic profiles. In this case, the 3 transverse seismic sections of the 4 sets of seismic data volumes may be identically numbered, resulting in the 3 transverse seismic sections LINE1, LINE2, and LINE 3. If each set of seismic data volumes is sectioned by longitudinal sections, the 4 sets of seismic data volumes may include 3 longitudinal seismic sections. In this case, 3 longitudinal seismic sections of the 4 sets of seismic data volumes may be identically numbered, resulting in 3 longitudinal seismic sections of TRACE1, TRACE2, and TRACE 3.
It is noted that for ease of seismic horizon interpretation, each of the sets of seismic data volumes in the embodiments of the present application includes the same number of transverse seismic sections and the same number of longitudinal seismic sections. And the transverse seismic sections of the multiple sets of seismic data bodies are in one-to-one correspondence, and the longitudinal seismic sections of the multiple sets of seismic data bodies are in one-to-one correspondence. By one-to-one correspondence of transverse seismic profiles, it is meant that one seismic profile in any one of the sets of seismic data volumes has a corresponding seismic profile in the other seismic data volumes, and the seismic profiles are all intended to be at the same location. In the embodiments of the present application, the corresponding seismic sections in each set of seismic data volumes may be identified by the same reference numerals. For example, the seismic data volume 1 includes 2 transverse seismic sections, LINE1 and LINE2, and 2 longitudinal seismic sections, LINE1 and LINE2, the seismic data volume 2 also includes 2 transverse seismic sections, LINE1 and LINE2, and 2 longitudinal seismic sections, LINE1 and LINE2, then the LINE1 transverse seismic section of the seismic data volume 1 corresponds to the LINE1 transverse seismic section of the seismic data volume 2, the LINE2 transverse seismic section of the seismic data volume 1 corresponds to the LINE2 transverse seismic section of the seismic data volume 2, the LINE1 longitudinal seismic section of the seismic data volume 1 corresponds to the LINE1 longitudinal section of the seismic data volume 2, and the LINE2 longitudinal seismic section of the seismic data volume 1 corresponds to the LINE2 longitudinal seismic section of the seismic data volume 2.
Finally, when one set of seismic data bodies is selected from the multiple sets of seismic data bodies as the reference data body, a user can lead all the multiple sets of seismic data bodies into the terminal, the number of any one seismic data body in the multiple sets of seismic data bodies is input or selected in the terminal, and correspondingly, the terminal can obtain the number of the seismic data body input or selected by the user, so that the seismic data body corresponding to the number can be used as the reference data body.
Step 402: and tracking a first co-axial of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, wherein the first seismic section is any one of the plurality of seismic sections of the reference data volume.
The first seismic section of the reference data volume may be any one of the plurality of transverse seismic sections of the reference data volume, or may be any one of the plurality of longitudinal seismic sections of the reference data volume. In addition, each seismic section, whether a transverse seismic section or a longitudinal seismic section, includes a plurality of seismic traces. The event axis refers to a connection line of extreme values (commonly called wave crest or wave trough) of the same vibration phase of each seismic channel on the seismic section, and each seismic section can be considered to be composed of a plurality of event axes. The first event is the event that the target zone to be investigated corresponds to on the first seismic section.
Additionally, the destination layer horizon of the first seismic section is used to indicate a location of the destination layer in the first seismic section. The target layer horizon of the reference data volume can be obtained by combining the target layers of all seismic sections of the reference data volume, and the target layer horizon of the reference data volume can reflect the fluctuation state of the target layer in the stratum.
It is noted that before step 402 is performed, the user may manually predict the continuity and stability of the first event and the amplitude intensity of the first reference points on the first event on each seismic trace of the first seismic section. If the continuity and stability of the first event are good, and the amplitude intensity of the first datum points on the first event on each seismic trace of the first seismic section meets the requirement of the amplitude intensity, the user can click the first datum points on the first event. Accordingly, the terminal may acquire a plurality of first reference points input by the user; determining an amplitude maximum value point corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points; and determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
It should be noted that, since the plurality of first reference points acquired by the terminal are determined according to the click operation of the user, there may be an error caused by human factors, and therefore, the user may input a time interval, and accordingly, the terminal may acquire the time interval input by the user and determine the first time window according to the time interval. After the terminal acquires the plurality of first reference points and the first time window, the amplitude maximum point corresponding to each first reference point in the first time window may be determined on the plurality of seismic traces by taking the plurality of first reference points as the center. In this way, the accuracy of the acquired amplitude maximum point can be improved. And finally, connecting the corresponding maximum amplitude points of each first reference point in the first time window by the terminal to obtain the target layer position of the first seismic section.
In practical application, the first time window may be a time range obtained by setting a time interval up and down with a certain time as a center. For example, when the time corresponding to one of the first reference points is 853ms (millisecond), and the time intervals set up above and below are 3ms and 7ms, respectively, the terminal determines that the first time window is 850ms to 860ms with 853ms as the center. The terminal may then determine the maximum amplitude point of the first reference point within 850ms to 860ms (the first time window) of the seismic trace on which the first reference point is located. Each point on each seismic trace can correspond to a coordinate, each coordinate can include 4 parameters, and the 4 parameters are X (east-west coordinate or abscissa), Y (north-south coordinate or ordinate), t (time), and a (amplitude), so that for all points located in the first time window on the seismic trace corresponding to the first reference point, the terminal can determine the amplitude maximum point according to the amplitude included in the coordinates of the points, and use the determined amplitude maximum point as the amplitude maximum point corresponding to the first reference point in the first time window.
Optionally, when the user determines that the continuity and stability of the first event are not good after analyzing the continuity and stability of the first event, or the amplitude intensity of the first reference points on the first event on each seismic trace of the first seismic profile does not meet the amplitude intensity requirement, the user may click the first reference points on the first event, and accordingly, the terminal may obtain the first reference points input by the user, and directly connect the obtained first reference points, thereby obtaining the target horizon of the first seismic profile.
Further, after the target horizon of the first seismic section is obtained in step 402, in order to facilitate subsequent use of the target horizon of the first seismic section, the target horizon of the first seismic section can be conveniently obtained, and the terminal can number and store the target horizon of the first seismic section.
The measurement areas corresponding to the sets of seismic data volumes are consistent, and only the relative relationship between the positions of the sets of seismic data volumes for receiving seismic waves and the positions for exciting seismic waves is changed, so that the sets of seismic data volumes are approximately similar. In this case, when it is determined that the continuity and stability of the first homodyne axis corresponding to the target interval on the first seismic section of the reference data volume are good and the amplitude intensities of the first reference points on the first homodyne axis on the seismic traces of the first seismic section satisfy the amplitude intensity requirement, it is considered that the continuity and stability of the homodyne axis corresponding to the target interval on each seismic section of the seismic data volumes other than the reference data volume among the plurality of sets of seismic data volumes are good and the amplitude intensities of the reference points on the homodyne axis on each seismic trace of each seismic section satisfy the amplitude intensity requirement. In this case, after the destination layer horizons of the first seismic section are determined by step 402, the destination layer horizons of the second seismic section of the target data volume may be determined directly from the destination layer horizons of the first seismic section without further performing steps 403 and 404.
When the target horizon of the second seismic profile of the target data volume is determined according to the target horizon of the first seismic profile, a plurality of points on the target horizon of the first seismic profile can be used as a plurality of reference points; determining the amplitude maximum point corresponding to each reference point in a preset time window on a plurality of seismic traces of the second seismic section according to the plurality of reference points; and determining the target layer horizon of the second seismic profile according to the corresponding amplitude maximum point of each reference point in the preset time window.
It is noted that the first seismic section includes a plurality of points at the horizon of the target layer. The terminal numbers and stores the target layer position of the first seismic section, so that a user can input or check the number of the target layer position of the first seismic section in the terminal, correspondingly, the terminal can obtain the number of the target layer position of the first seismic section input or checked by the user, and then a plurality of points on the target layer position of the first seismic section corresponding to the number can be used as a plurality of reference points. After the terminal acquires the plurality of reference points, the plurality of reference points can be mapped onto the second seismic profile to obtain a plurality of target points on the second seismic profile, wherein the target points correspond to the plurality of reference points one by one. Then, in order to improve the seismic horizon interpretation precision, the terminal may obtain a time interval preset by a user, determine a time window according to the time interval, and then determine the maximum amplitude point of each target point in the time window on the plurality of seismic traces of the second seismic profile. And finally, the terminal connects the corresponding maximum amplitude points of each target point in the time window, and the target layer horizon of the second seismic section can be obtained.
Alternatively, when it is determined that the continuity and stability of the first homodyne axis corresponding to the target interval on the first seismic section of the reference data volume is not good, and/or the amplitude intensities of the first reference points on the first homodyne axis on the seismic traces of the first seismic section do not meet the amplitude intensity requirement, it may be considered that the continuity and stability of the corresponding homodyne axes of the target interval on the seismic traces of the seismic data volumes other than the reference data volume in the plurality of sets of seismic data volumes are not good, and/or the amplitude intensities of the reference points on the homodyne axis on the seismic traces of the seismic sections do not meet the amplitude intensity requirement. In this case, after the target horizon for the first seismic section is determined by step 402, the target horizon for the second seismic section of the target data volume may be determined by steps 403 through 404 as follows.
Step 403: and tracking the second in-phase axis of the first seismic section to obtain a standard layer horizon of the first seismic section.
It should be noted that the terminal may obtain a second event determined by the user according to the first event, and then track the second event, so as to obtain the standard horizon of the first seismic section. The second homophasic axis of the first seismic section is a homophasic axis, the distance between the first homophasic axis and the second homophasic axis on the first seismic section is smaller than or equal to a preset distance threshold, the second homophasic axis is better in continuity and stability compared with the first homophasic axis, and the amplitude intensity of a plurality of datum points on the second homophasic axis on each seismic channel on the first seismic section meets the requirement of the amplitude intensity.
Additionally, the standard layer horizons of the first seismic section are locations of the standard layers of the first seismic section in the first seismic section.
Specifically, the implementation process of step 403 may be: obtaining a plurality of second reference points on a second event axis of the first seismic section; determining the amplitude maximum value point corresponding to each second reference point in a second time window on the plurality of seismic traces according to the plurality of second reference points; and determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
It should be noted that the terminal may acquire a plurality of second reference points on the second in-phase axis input by the user. In order to improve the accuracy of the obtained amplitude maximum point, the terminal may obtain a time interval input by the user and determine the second time window according to the time interval. After the terminal acquires the plurality of second reference points and the second time window, the maximum amplitude point of each second reference point in the second time window can be determined on the plurality of seismic traces by taking the plurality of second reference points as the center. And finally, connecting the corresponding amplitude maximum value points of each second reference point in the second time window by the terminal to obtain the standard layer horizon of the first seismic section.
In addition, the description of the second time window is similar to that of the first time window in step 402, and is not repeated here.
Further, after the standard horizon of the first seismic section is obtained in step 403, in order to facilitate the subsequent use of the standard horizon of the first seismic section, the standard horizon of the first seismic section can be conveniently obtained, and the terminal can number and store the standard horizon of the first seismic section.
After the target layer level and the standard layer level of the first seismic section of the reference data volume are obtained through steps 402 and 403, the target layer level of the second seismic section of the target data volume may be determined through step 404.
Step 404: and determining the target layer position of the second seismic section of the target data volume according to the target layer position and the standard layer position of the first seismic section.
The target data volume is any one of the seismic data volumes except the reference data volume, and the second seismic section is a seismic section corresponding to the first seismic section in the seismic sections of the target data volume.
It should be noted that the destination layer level of the second seismic section of the target data volume is the position of the destination layer in the second seismic section. And combining the target layers of all seismic sections of the target data volume to obtain the target layer horizon of the target data volume, wherein the target layer horizon of the target data volume can reflect the fluctuation state of the target layer in the stratum.
Specifically, step 404 may be implemented by steps 4041-4043 as follows:
step 4041: and determining the standard layer position of the second seismic section according to the standard layer position of the first seismic section.
In particular, the terminal may take a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points; determining the amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points; and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
It should be noted that, in the specific implementation process of this step, reference may be made to a related implementation manner of determining the target layer level of the second seismic profile of the target data volume according to the target layer level of the first seismic profile in step 402, and details of the embodiment of this application are not described herein again.
In addition, the description of the third time window is similar to that of the first time window in step 402, and is not repeated here.
Further, after the standard layer position of the second seismic profile is obtained in step 4041, in order to facilitate subsequent use of the standard layer position of the second seismic profile, the standard layer position of the second seismic profile can be conveniently obtained, and the terminal can number and store the standard layer position of the second seismic profile.
Step 4042: and determining the difference value between each point on the standard layer position of the first seismic profile and the corresponding point on the standard layer position of the second seismic profile according to the standard layer position of the first seismic profile and the standard layer position of the second seismic profile, and correspondingly storing the coordinate of each point on the standard layer position of the second seismic profile and the determined corresponding difference value to obtain the corresponding relation between the coordinate and the difference value.
And similarly, a plurality of intersection points also exist between the standard layer position of the second seismic section and the plurality of seismic channels included in the second seismic section. In an embodiment of the present application, the plurality of intersections of the standard layer horizon with the plurality of seismic traces may be taken as a plurality of points on the standard layer horizon. Based on this, after determining the standard layer level of the first seismic section and the standard layer level of the second seismic section, the terminal may determine a difference between each point on the standard layer level of the first seismic section and a corresponding point on the standard layer level of the second seismic section according to coordinates of each point on the standard layer level of the first seismic section and coordinates of each point on the standard layer level of the second seismic section.
It should be noted that, since the terminal numbers and stores both the standard layer position of the first seismic profile and the standard layer position of the second seismic profile, the user can input or select the number of the standard layer position of the first seismic profile and the number of the standard layer position of the second seismic profile in the terminal, and accordingly, the terminal can obtain the number of the standard layer position of the first seismic profile and the number of the standard layer position of the second seismic profile, which are input or selected by the user, so as to obtain the coordinates of each point on the standard layer position of the first seismic profile and the standard layer position of the second seismic profile, which correspond to the two numbers.
After the coordinates of each point on the standard layer level of the first seismic profile and the coordinates of each point on the standard layer level of the second seismic profile are obtained, for any point on the standard layer level of the first seismic profile, the terminal can determine and obtain the time difference between the two points according to the coordinates of the point and the coordinates of the point corresponding to the point on the standard layer level of the second seismic profile. The process of determining the difference between a first target point and the corresponding point of the first target point on the second seismic section will now be described, taking as an example the first target point on the standard horizon of the first seismic section. Wherein the first target point may be any point on a standard layer horizon of the first seismic section.
The terminal can acquire coordinates of a first target point on a standard layer horizon of a first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time; determining a second target point on the standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate; a time difference between a first time of the first target point and a second time of the second target point is determined, and the time difference is determined as a difference between the first target point and the second target point.
It should be noted that, when the second target point on the standard layer level of the second seismic profile is determined according to the first abscissa and the first ordinate, it may be determined that, among a plurality of points included on the standard layer level of the second seismic profile, an abscissa of the plurality of points is equal to the first abscissa, and a point whose ordinate is equal to the second ordinate is the second target point, according to the first abscissa and the first ordinate. After determining the second target point, the terminal may calculate a time difference between a first time included in the coordinates of the first target point and a second time included in the coordinates of the second target point, thereby determining the time difference as a difference between the first target point and the second target point.
For example, the first seismic section includes A (x) at the standard horizon1,y1,t1) And B (x)2,y2,t2) The standard layer horizon of the second seismic profile comprises C (x)1,y1,t3) And D (x)2,y2,t4) These two points. When the point A on the standard layer horizon of the first seismic section is the first target point, the first abscissa of the point A is x1The first ordinate of the point A is y1Since the abscissa of the point C is equal to the first abscissa of the point a and the ordinate of the point C is equal to the first ordinate of the point a, the point C may be determined as the second target point on the standard layer horizon of the second seismic profile. Thereafter, t may be adjusted1And t3Making a difference to obtain a time difference t between the first time of the first target point and the second time of the second target point1-t3This time difference t1-t3I.e. the difference between the first target point and the second target point.
After obtaining the time difference between each point on the first seismic section and the corresponding point of each point on the second seismic section, the terminal may store each time difference in correspondence with the coordinates of the corresponding point, thereby obtaining the correspondence between the coordinates and the difference.
It should be noted that, for each time difference, since the coordinates of the two points for calculating the time difference include the same abscissa and ordinate, the terminal may store the time difference in correspondence with the abscissa and ordinate of the corresponding two points.
Still taking the above example as an example, when the first target point is point a and the second target point is point C, the difference between the first target point and the second target point is t1-t3. After obtaining the time difference, the terminal may compare the time difference t because the points a and C have the same abscissa and the same ordinate1-t3Storing the abscissa and ordinate of the points A and C in correspondence, i.e. storing t1-t3And (x)1,y1) And correspondingly storing.
For each point on the standard layer horizon of the first seismic profile and each point on the standard layer horizon of the second seismic profile, the terminal can refer to the processing mode for the first target point and the second target point for processing, and correspondingly store each determined difference value and the coordinate of the point corresponding to each difference value, so as to obtain the corresponding relation between the coordinate and the difference value.
Step 4043: and determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference.
After the corresponding relation between the coordinates and the difference values is obtained, the terminal can obtain the coordinates of each point on the target layer position of the first seismic section, wherein the coordinates of each point comprise a third horizontal coordinate, a third vertical coordinate and a third time; determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section; determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point; determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point; and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
It should be noted that since the target horizon of the first seismic section also includes a plurality of intersections with the plurality of seismic traces, the target horizon of the first seismic section may be considered to include a plurality of points. The terminal numbers and stores the target layer position of the first seismic section, so that a user can input or select the number of the target layer position of the first seismic section in the terminal, and correspondingly, the terminal can obtain the number of the target layer position of the first seismic section input or selected by the user, and thus the coordinate of each point on the target layer position of the first seismic section corresponding to the number can be obtained.
After the coordinates of each point on the target layer level of the first seismic profile are obtained, for any point on the target layer level of the first seismic profile, the terminal can obtain a corresponding difference value from the corresponding relation between the coordinates and the difference value according to the coordinates of the point, and further determine the corresponding point of the point on the second seismic profile according to the coordinates of the point and the obtained difference value. The process of determining the difference values for the corresponding points will be described below using a target point on the target horizon of the first seismic section as an example. For the sake of convenience of explanation, this target point is referred to as a fourth target point.
The terminal may obtain, from the correspondence between the coordinates and the difference values, a difference value corresponding to the coordinate equal to the third abscissa and equal to the third ordinate according to the third abscissa and the third ordinate of the fourth target point, and determine the difference value as the difference value corresponding to the fourth target point.
After determining that the difference corresponding to the fourth target point is obtained, the terminal may calculate a difference between a third time included in the coordinates of the fourth target point and the obtained difference, so that the calculated difference is used as a time included in the coordinates of the third target point corresponding to the fourth target point on the second seismic profile. In addition, the terminal may use a third abscissa and a third ordinate included in the coordinates of the fourth target point as the abscissa and the ordinate of the third target point, so that the coordinates of the third target point may be obtained, and the third target point may be determined and obtained on the second seismic profile according to the coordinates of the third target point.
For example, when the fourth target point E has coordinates of (x)1,y1,t5) The third abscissa of the point E is x1The third ordinate of the point E is y1. If the difference t is1-t3The corresponding coordinate is (x)1,y1) Then according to the third abscissa x of the fourth target point E1And a third ordinate y1The coordinates (x) can be expressed1,y1) Corresponding difference t1-t3And determining the difference corresponding to the fourth target point E. Then, the third time t of the fourth target point E may be5The difference t corresponding to the third target point E1-t3Making a difference to obtain a difference value t between the third time of the fourth target point E and the difference value corresponding to the obtained target point E5-(t1-t3) That is, the time of obtaining the coordinates of the target point E in the third target point corresponding to the second seismic profile is: t is t5-(t1-t3). Thereafter, the abscissa x of the fourth target point E is plotted1As the abscissa of the third target point, the ordinate y of the fourth target point E1As the ordinate of the third target point, the coordinate of the third target point is obtained as (x)1,y1,t5-(t1-t3) And determining a third target point on the second seismic profile according to the coordinates of the third target point.
For each point on the target layer horizon of the first seismic profile, the terminal can process the fourth target point by referring to the processing of the fourth target point, so as to obtain a point corresponding to each point on the second seismic profile, and connect the points determined on the second seismic profile, so that the target layer horizon of the second seismic profile can be obtained.
FIG. 5 is a comparison of two-layer horizons on a second seismic profile of a target data volume interpreted by the methods provided in the embodiments of the present application with two-layer horizons of a second seismic profile of a target data volume interpreted according to methods in the related art. Wherein FIG. 5 shows the second seismic section of the target data volume, and L3 is the horizon of the first layer to be studied of the first seismic section of the reference data volume shown in FIG. 1, taken directly after the horizon of the first layer to be studied of the second seismic section of the target data volume (L1); l4 is the horizon of the second layer to be investigated of the second seismic section of the target data volume obtained directly after the horizon of the second layer to be investigated of the first seismic section of the reference data volume shown in fig. 1 (L2). L5 is the horizon of the first layer to be studied of the second seismic profile of the target data volume obtained according to the method of the embodiment of the present application, and L6 is the horizon of the second layer to be studied of the second seismic profile of the target data volume obtained according to the method of the embodiment of the present application. As can be seen from fig. 5, both L3 and L4 have a local layer crossing phenomenon on the second seismic section of the target data volume, and both L5 and L6 have no layer crossing phenomenon on the second seismic section of the target data volume, which indicates that when the layer interpretation result of one set of seismic data volume is directly used as the layer interpretation result of other seismic data volumes, the interpretation precision is low for other seismic data volumes, and the layer interpretation result obtained according to the embodiment of the present application has higher precision, and compared with the related art, the seismic layer interpretation precision is improved.
According to the method and the device for determining the seismic profile of the target data volume, one set of seismic data volume can be selected from multiple sets of seismic data volumes to serve as a reference data volume, a first homophase axis and a second homophase axis of a first seismic profile of the reference data volume are tracked, a target layer position and a standard layer position of the first seismic profile are obtained respectively, and then the target layer position of a second seismic profile of the target data volume is determined according to the target layer position and the standard layer position of the first seismic profile of the reference data volume. Therefore, only the reference data bodies in the multiple sets of seismic data bodies need to be subjected to horizon interpretation in a conventional mode of in-phase axis tracking, the target layer horizon of any seismic data body except the reference data body in the multiple sets of seismic data bodies can be obtained through the target layer horizon and the standard layer horizon of the reference data body, and the consumption of a large amount of time caused by the conventional in-phase axis tracking of other seismic data bodies except the reference data body for multiple times is avoided. In addition, the method and the device can determine the target layers of other seismic data volumes in the multiple sets of seismic data volumes according to the standard layer position and the target layer position of the reference data volume to perform layer position interpretation, and compared with the method and the device which directly take the layer position interpretation result of one set of seismic data volume as the layer position interpretation result of other seismic data volumes in the related technology, the method and the device improve interpretation precision.
Fig. 6 is a schematic structural diagram of a seismic horizon interpreting apparatus according to an embodiment of the present application. Referring to fig. 6, the apparatus includes: a selection module 601, a tracking module 602, and a determination module 603.
A selecting module 601, configured to select a set of seismic data volumes from multiple sets of seismic data volumes of an area to be studied as a reference data volume, where each set of seismic data volume in the multiple sets of seismic data volumes includes multiple seismic sections, and the multiple seismic sections in each set of seismic data volume correspond to each other one by one;
a tracking module 602, configured to track a first event of a first seismic section of a reference data volume to obtain a target horizon of the first seismic section, and track a second event of the first seismic section to obtain a standard horizon of the first seismic section, where the first seismic section is any one of a plurality of seismic sections of the reference data volume;
a determining module 603, configured to determine a target horizon of a second seismic profile of the target data volume according to the target horizon and the standard horizon of the first seismic profile, where the target data volume is any one of the plurality of sets of seismic data volumes except the reference data volume, and the second seismic profile is a seismic profile corresponding to the first seismic profile in the plurality of seismic profiles of the target data volume.
Optionally, each seismic section of the plurality of seismic sections comprises a plurality of seismic traces;
the tracking module 602 includes:
a first acquisition submodule for acquiring a plurality of first reference points on a first event axis of a first seismic section;
the first determining submodule is used for determining an amplitude maximum value point corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and the second determining submodule is used for determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
Optionally, the tracking module 602 further comprises:
the second acquisition submodule is used for acquiring a plurality of second datum points on a second in-phase axis of the first seismic section;
a third determining submodule, configured to determine, according to the plurality of second reference points, an amplitude maximum point corresponding to each second reference point in a second time window on the plurality of seismic traces;
and the fourth determining submodule is used for determining the standard layer horizon of the first seismic section according to the amplitude maximum point corresponding to each second reference point in the second time window.
Optionally, the determining module 603 includes:
the fifth determining submodule is used for determining the standard layer position of the second seismic section according to the standard layer position of the first seismic section;
a sixth determining submodule, configured to determine, according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile, a difference between each point on the standard layer level of the first seismic profile and a corresponding point on the standard layer level of the second seismic profile, and store coordinates of each point on the standard layer level of the second seismic profile and the determined corresponding difference, to obtain a correspondence between the coordinates and the differences;
and the seventh determining submodule is used for determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinate and the difference value.
Optionally, the fifth determining sub-module is configured to:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining the amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
Optionally, the sixth determining sub-module is configured to:
acquiring coordinates of a first target point on a standard layer horizon of a first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on the standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
a time difference between a first time of the first target point and a second time of the second target point is determined, and the time difference is determined as a difference between the first target point and the second target point.
Optionally, the seventh determining sub-module is configured to:
acquiring coordinates of each point on a target layer horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
According to the method and the device for determining the seismic profile of the target data volume, one set of seismic data volume can be selected from multiple sets of seismic data volumes to serve as a reference data volume, a first homophase axis and a second homophase axis of a first seismic profile of the reference data volume are tracked, a target layer position and a standard layer position of the first seismic profile are obtained respectively, and then the target layer position of a second seismic profile of the target data volume is determined according to the target layer position and the standard layer position of the first seismic profile of the reference data volume. Therefore, only the reference data bodies in the multiple sets of seismic data bodies need to be subjected to horizon interpretation in a conventional mode of in-phase axis tracking, the target layer horizon of any seismic data body except the reference data body in the multiple sets of seismic data bodies can be obtained through the target layer horizon and the standard layer horizon of the reference data body, and the consumption of a large amount of time caused by the conventional in-phase axis tracking of other seismic data bodies except the reference data body for multiple times is avoided. In addition, the method and the device can determine the target layers of other seismic data volumes in the multiple sets of seismic data volumes according to the standard layer position and the target layer position of the reference data volume to perform layer position interpretation, and compared with the method and the device which directly take the layer position interpretation result of one set of seismic data volume as the layer position interpretation result of other seismic data volumes in the related technology, the method and the device improve interpretation precision.
It should be noted that: in the seismic horizon interpretation apparatus provided in the above embodiment, when interpreting a seismic horizon, only the division of each functional module is exemplified, and in practical application, the function distribution may be completed by different functional modules according to needs, that is, the internal structure of the apparatus is divided into different functional modules, so as to complete all or part of the functions described above. In addition, the seismic horizon interpreting apparatus and the seismic horizon interpreting method provided by the above embodiments belong to the same concept, and the specific implementation process thereof is described in the method embodiments, and is not described herein again.
Fig. 7 is a schematic structural diagram of a seismic horizon interpreting apparatus according to an embodiment of the present application. Referring to fig. 7, the apparatus may be a terminal 700, and the terminal 700 may be: a smart phone, a tablet computer, an MP3 player (Moving picture Experts Group Audio Layer III, motion picture Experts compression standard Audio Layer 3), an MP4 player (Moving picture Experts Group Audio Layer IV, motion picture Experts compression standard Audio Layer 4), a notebook computer or a desktop computer. Terminal 700 may also be referred to by other names such as user equipment, portable terminal, laptop terminal, desktop terminal, and so on.
In general, terminal 700 includes: a processor 701 and a memory 702.
The processor 701 may include one or more processing cores, such as a 4-core processor, an 8-core processor, and so on. The processor 701 may be implemented in at least one hardware form of a DSP (Digital Signal Processing), an FPGA (Field-Programmable Gate Array), and a PLA (Programmable Logic Array). The processor 701 may also include a main processor and a coprocessor, where the main processor is a processor for processing data in an awake state, and is also called a Central Processing Unit (CPU); a coprocessor is a low power processor for processing data in a standby state. In some embodiments, the processor 701 may be integrated with a GPU (Graphics Processing Unit), which is responsible for rendering and drawing the content required to be displayed on the display screen. In some embodiments, the processor 701 may further include an AI (Artificial Intelligence) processor for processing computing operations related to machine learning.
Memory 702 may include one or more computer-readable storage media, which may be non-transitory. Memory 702 may also include high-speed random access memory, as well as non-volatile memory, such as one or more magnetic disk storage devices, flash memory storage devices. In some embodiments, a non-transitory computer readable storage medium in memory 702 is used to store at least one instruction for execution by processor 701 to implement a seismic horizon interpretation method as provided by method embodiments herein.
In some embodiments, the terminal 700 may further optionally include: a peripheral interface 703 and at least one peripheral. The processor 701, the memory 702, and the peripheral interface 703 may be connected by buses or signal lines. Various peripheral devices may be connected to peripheral interface 703 via a bus, signal line, or circuit board. Specifically, the peripheral device includes: at least one of radio frequency circuitry 704, touch screen display 705, camera 706, audio circuitry 707, positioning components 708, and power source 709.
The peripheral interface 703 may be used to connect at least one peripheral related to I/O (Input/Output) to the processor 701 and the memory 702. In some embodiments, processor 701, memory 702, and peripheral interface 703 are integrated on the same chip or circuit board; in some other embodiments, any one or two of the processor 701, the memory 702, and the peripheral interface 703 may be implemented on a separate chip or circuit board, which is not limited in this embodiment.
The Radio Frequency circuit 704 is used for receiving and transmitting RF (Radio Frequency) signals, also called electromagnetic signals. The radio frequency circuitry 704 communicates with communication networks and other communication devices via electromagnetic signals. The rf circuit 704 converts an electrical signal into an electromagnetic signal to transmit, or converts a received electromagnetic signal into an electrical signal. Optionally, the radio frequency circuit 704 includes: an antenna system, an RF transceiver, one or more amplifiers, a tuner, an oscillator, a digital signal processor, a codec chipset, a subscriber identity module card, and so forth. The radio frequency circuitry 704 may communicate with other terminals via at least one wireless communication protocol. The wireless communication protocols include, but are not limited to: metropolitan area networks, various generation mobile communication networks (2G, 3G, 4G, and 5G), Wireless local area networks, and/or WiFi (Wireless Fidelity) networks. In some embodiments, the radio frequency circuit 704 may also include NFC (Near Field Communication) related circuits, which are not limited in this application.
The display screen 705 is used to display a UI (User Interface). The UI may include graphics, text, icons, video, and any combination thereof. When the display screen 705 is a touch display screen, the display screen 705 also has the ability to capture touch signals on or over the surface of the display screen 705. The touch signal may be input to the processor 701 as a control signal for processing. At this point, the display 705 may also be used to provide virtual buttons and/or a virtual keyboard, also referred to as soft buttons and/or a soft keyboard. In some embodiments, the display 705 may be one, providing the front panel of the terminal 700; in other embodiments, the display 705 can be at least two, respectively disposed on different surfaces of the terminal 700 or in a folded design; in still other embodiments, the display 705 may be a flexible display disposed on a curved surface or on a folded surface of the terminal 700. Even more, the display 705 may be arranged in a non-rectangular irregular pattern, i.e. a shaped screen. The Display 705 may be made of LCD (liquid crystal Display), OLED (Organic Light-Emitting Diode), or the like.
The camera assembly 706 is used to capture images or video. Optionally, camera assembly 706 includes a front camera and a rear camera. Generally, a front camera is disposed at a front panel of the terminal, and a rear camera is disposed at a rear surface of the terminal. In some embodiments, the number of the rear cameras is at least two, and each rear camera is any one of a main camera, a depth-of-field camera, a wide-angle camera and a telephoto camera, so that the main camera and the depth-of-field camera are fused to realize a background blurring function, and the main camera and the wide-angle camera are fused to realize panoramic shooting and VR (Virtual Reality) shooting functions or other fusion shooting functions. In some embodiments, camera assembly 706 may also include a flash. The flash lamp can be a monochrome temperature flash lamp or a bicolor temperature flash lamp. The double-color-temperature flash lamp is a combination of a warm-light flash lamp and a cold-light flash lamp, and can be used for light compensation at different color temperatures.
The audio circuitry 707 may include a microphone and a speaker. The microphone is used for collecting sound waves of a user and the environment, converting the sound waves into electric signals, and inputting the electric signals to the processor 701 for processing or inputting the electric signals to the radio frequency circuit 704 to realize voice communication. For the purpose of stereo sound collection or noise reduction, a plurality of microphones may be provided at different portions of the terminal 700. The microphone may also be an array microphone or an omni-directional pick-up microphone. The speaker is used to convert electrical signals from the processor 701 or the radio frequency circuit 704 into sound waves. The loudspeaker can be a traditional film loudspeaker or a piezoelectric ceramic loudspeaker. When the speaker is a piezoelectric ceramic speaker, the speaker can be used for purposes such as converting an electric signal into a sound wave audible to a human being, or converting an electric signal into a sound wave inaudible to a human being to measure a distance. In some embodiments, the audio circuitry 707 may also include a headphone jack.
The positioning component 708 is used to locate the current geographic position of the terminal 700 to implement navigation or LBS (location based Service). The positioning component 708 may be a positioning component based on the GPS (global positioning System) in the united states, the beidou System in china, the graves System in russia, or the galileo System in the european union.
Power supply 709 is provided to supply power to various components of terminal 700. The power source 709 may be alternating current, direct current, disposable batteries, or rechargeable batteries. When power source 709 includes a rechargeable battery, the rechargeable battery may support wired or wireless charging. The rechargeable battery may also be used to support fast charge technology.
In some embodiments, terminal 700 also includes one or more sensors 710. The one or more sensors 710 include, but are not limited to: acceleration sensor 711, gyro sensor 712, pressure sensor 713, fingerprint sensor 714, optical sensor 715, and proximity sensor 716.
The acceleration sensor 711 can detect the magnitude of acceleration in three coordinate axes of a coordinate system established with the terminal 700. For example, the acceleration sensor 711 may be used to detect components of the gravitational acceleration in three coordinate axes. The processor 701 may control the touch screen 705 to display the user interface in a landscape view or a portrait view according to the gravitational acceleration signal collected by the acceleration sensor 711. The acceleration sensor 711 may also be used for acquisition of motion data of a game or a user.
The gyro sensor 712 may detect a body direction and a rotation angle of the terminal 700, and the gyro sensor 712 may cooperate with the acceleration sensor 711 to acquire a 3D motion of the terminal 700 by the user. From the data collected by the gyro sensor 712, the processor 701 may implement the following functions: motion sensing (such as changing the UI according to a user's tilting operation), image stabilization at the time of photographing, game control, and inertial navigation.
Pressure sensors 713 may be disposed on a side bezel of terminal 700 and/or an underlying layer of touch display 705. When the pressure sensor 713 is disposed on a side frame of the terminal 700, a user's grip signal on the terminal 700 may be detected, and the processor 701 performs right-left hand recognition or shortcut operation according to the grip signal collected by the pressure sensor 713. When the pressure sensor 713 is disposed at a lower layer of the touch display 705, the processor 701 controls the operability control on the UI interface according to the pressure operation of the user on the touch display 705. The operability control comprises at least one of a button control, a scroll bar control, an icon control and a menu control.
The fingerprint sensor 714 is used for collecting a fingerprint of a user, and the processor 701 identifies the identity of the user according to the fingerprint collected by the fingerprint sensor 714, or the fingerprint sensor 714 identifies the identity of the user according to the collected fingerprint. When the user identity is identified as a trusted identity, the processor 701 authorizes the user to perform relevant sensitive operations, including unlocking a screen, viewing encrypted information, downloading software, paying, changing settings, and the like. The fingerprint sensor 714 may be disposed on the front, back, or side of the terminal 700. When a physical button or a vendor Logo is provided on the terminal 700, the fingerprint sensor 714 may be integrated with the physical button or the vendor Logo.
The optical sensor 715 is used to collect the ambient light intensity. In one embodiment, the processor 701 may control the display brightness of the touch display 705 based on the ambient light intensity collected by the optical sensor 715. Specifically, when the ambient light intensity is high, the display brightness of the touch display screen 705 is increased; when the ambient light intensity is low, the display brightness of the touch display 705 is turned down. In another embodiment, processor 701 may also dynamically adjust the shooting parameters of camera assembly 706 based on the ambient light intensity collected by optical sensor 715.
A proximity sensor 716, also referred to as a distance sensor, is typically disposed on a front panel of the terminal 700. The proximity sensor 716 is used to collect the distance between the user and the front surface of the terminal 700. In one embodiment, when the proximity sensor 716 detects that the distance between the user and the front surface of the terminal 700 gradually decreases, the processor 701 controls the touch display 705 to switch from the bright screen state to the dark screen state; when the proximity sensor 716 detects that the distance between the user and the front surface of the terminal 700 gradually becomes larger, the processor 701 controls the touch display 705 to switch from the breath screen state to the bright screen state.
Those skilled in the art will appreciate that the configuration shown in fig. 7 is not intended to be limiting of terminal 700 and may include more or fewer components than those shown, or some components may be combined, or a different arrangement of components may be used.
It will be understood by those skilled in the art that all or part of the steps for implementing the above embodiments may be implemented by hardware, or may be implemented by a program instructing relevant hardware, where the program may be stored in a computer-readable storage medium, and the above-mentioned storage medium may be a read-only memory, a magnetic disk or an optical disk, etc.
The above description is only exemplary of the present application and should not be taken as limiting the present application, as any modification, equivalent replacement, or improvement made within the spirit and principle of the present application should be included in the protection scope of the present application.

Claims (16)

1. A seismic horizon interpretation method, the method comprising:
selecting a set of seismic data bodies from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, wherein each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data bodies are in one-to-one correspondence;
tracking a first homodyne axis of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, and tracking a second homodyne axis of the first seismic section to obtain a standard horizon of the first seismic section, wherein the first seismic section is any one of a plurality of seismic sections of the reference data volume;
and determining the target horizon of a second seismic profile of the target data volume according to the target horizon and the standard horizon of the first seismic profile, wherein the target data volume is any one of the plurality of sets of seismic data volumes except the reference data volume, and the second seismic profile is a seismic profile corresponding to the first seismic profile in the plurality of seismic profiles of the target data volume.
2. The method of claim 1, wherein each seismic section of the plurality of seismic sections includes a plurality of seismic traces;
the tracking a first co-axis of the first seismic section to obtain a target horizon of the first seismic section comprises:
obtaining a plurality of first reference points on a first co-axial of the first seismic section;
determining amplitude maximum points corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
3. The method of claim 2, wherein said tracking the second in-phase axis of the first seismic section to obtain a standard layer horizon for the first seismic section comprises:
obtaining a plurality of second reference points on a second event axis of the first seismic section;
determining amplitude maximum points corresponding to each second reference point in a second time window on the plurality of seismic traces according to the plurality of second reference points;
and determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
4. The method of claim 1, wherein determining the destination layer horizon for the second seismic section of the target data volume from the destination layer horizon and the standard layer horizon for the first seismic section comprises:
determining a standard layer position of the second seismic section according to the standard layer position of the first seismic section;
determining a difference value between each point on the standard layer level of the first seismic profile and a corresponding point on the standard layer level of the second seismic profile according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile, and correspondingly storing the coordinate of each point on the standard layer level of the second seismic profile and the determined corresponding difference value to obtain a corresponding relation between the coordinate and the difference value;
and determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference value.
5. The method of claim 4, wherein determining the standard layer horizons for the second seismic section from the standard layer horizons for the first seismic section comprises:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining an amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
6. The method of claim 4, wherein determining a difference between each point on the standard layer of the first seismic profile and a corresponding point on the standard layer horizon of the second seismic profile based on the standard layer horizons of the first seismic profile and the standard layer horizons of the second seismic profile comprises:
acquiring coordinates of a first target point on a standard layer horizon of the first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on a standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
determining a time difference between a first time of the first target point and a second time of the second target point, and determining the time difference as a difference between the first target point and the second target point.
7. The method of claim 6, wherein determining the destination layer level of the second seismic section from the destination layer level of the first seismic section and the correspondence of the coordinates and the difference comprises:
acquiring coordinates of each point on a target horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
8. A seismic horizon interpretation apparatus, the apparatus comprising:
the system comprises a selection module, a data analysis module and a data analysis module, wherein the selection module is used for selecting a set of seismic data body from a plurality of sets of seismic data bodies of an area to be researched as a reference data body, each set of seismic data body in the plurality of sets of seismic data bodies comprises a plurality of seismic sections, and the plurality of seismic sections in each set of seismic data body correspond to one another;
a tracking module, configured to track a first event of a first seismic section of the reference data volume to obtain a target horizon of the first seismic section, and track a second event of the first seismic section to obtain a standard horizon of the first seismic section, where the first seismic section is any one of a plurality of seismic sections of the reference data volume;
and the determining module is used for determining the target horizon of a second seismic section of the target data volume according to the target horizon and the standard horizon of the first seismic section, wherein the target data volume is any one of the plurality of sets of seismic data volumes except the reference data volume, and the second seismic section is a seismic section corresponding to the first seismic section in the plurality of seismic sections of the target data volume.
9. The apparatus of claim 8, wherein each seismic section of the plurality of seismic sections comprises a plurality of seismic traces;
the tracking module includes:
a first acquisition submodule for acquiring a plurality of first reference points on a first event axis of the first seismic section;
the first determining submodule is used for determining an amplitude maximum value point corresponding to each first reference point in a first time window on the plurality of seismic traces according to the plurality of first reference points;
and the second determining submodule is used for determining the target layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each first reference point in the first time window.
10. The apparatus of claim 9, wherein the tracking module further comprises:
a second acquisition submodule for acquiring a plurality of second reference points on a second in-phase axis of the first seismic section;
a third determining submodule, configured to determine, according to the plurality of second reference points, an amplitude maximum point corresponding to each second reference point in a second time window on the plurality of seismic traces;
and the fourth determining submodule is used for determining the standard layer horizon of the first seismic section according to the amplitude maximum value point corresponding to each second reference point in the second time window.
11. The apparatus of claim 8, wherein the determining module comprises:
the fifth determining submodule is used for determining the standard layer position of the second seismic section according to the standard layer position of the first seismic section;
a sixth determining submodule, configured to determine, according to the standard layer level of the first seismic profile and the standard layer level of the second seismic profile, a difference between each point on the standard layer level of the first seismic profile and a corresponding point on the standard layer level of the second seismic profile, and store coordinates of each point on the standard layer level of the second seismic profile and the corresponding difference determined for the corresponding point, so as to obtain a corresponding relationship between the coordinates and the difference;
and the seventh determining submodule is used for determining the target layer position of the second seismic section according to the target layer position of the first seismic section and the corresponding relation between the coordinates and the difference values.
12. The apparatus of claim 11, wherein the fifth determination submodule is to:
using a plurality of points on a standard layer horizon of the first seismic section as a plurality of third reference points;
determining an amplitude maximum point corresponding to each third reference point in a third time window on a plurality of seismic traces of the second seismic section according to the plurality of third reference points;
and determining the standard layer horizon of the second seismic profile according to the amplitude maximum point corresponding to each third reference point in the third time window.
13. The apparatus of claim 11, wherein the sixth determination submodule is to:
acquiring coordinates of a first target point on a standard layer horizon of the first seismic section, wherein the coordinates of the first target point comprise a first abscissa, a first ordinate and a first time, and the first target point is any point on the standard layer horizon of the first seismic section;
determining a second target point on a standard layer level of the second seismic profile according to the first abscissa and the first ordinate, wherein the second abscissa of the second target point is equal to the first abscissa, and the second ordinate of the second target point is equal to the first ordinate;
determining a time difference between a first time of the first target point and a second time of the second target point, and determining the time difference as a difference between the first target point and the second target point.
14. The apparatus of claim 13, wherein the seventh determination submodule is to:
acquiring coordinates of each point on a target horizon of the first seismic section, wherein the coordinates of each point comprise a third abscissa, a third ordinate and a third time;
determining a corresponding difference value from the corresponding relation between the coordinates and the difference value according to a third abscissa and a third ordinate which are included in the coordinates of each point on the target horizon of the first seismic section;
determining a difference value between the third time included in the third coordinate of each point and the obtained difference value corresponding to the corresponding point, and obtaining a fourth time corresponding to the corresponding point;
determining a plurality of third target points on the second seismic profile according to a third abscissa and a third ordinate included in the coordinates of each point and a fourth time corresponding to each point;
and determining the target layer horizon of the second seismic profile according to the plurality of third target points.
15. A seismic horizon interpretation apparatus, the apparatus comprising:
a processor;
a memory for storing processor-executable instructions;
wherein the processor is configured to perform the steps of any of the methods of claims 1-7.
16. A computer-readable storage medium having instructions stored thereon, wherein the instructions, when executed by a processor, implement the steps of any of the methods of claims 1-7.
CN201910161351.3A 2019-03-04 2019-03-04 Seismic horizon interpretation method and device Active CN111650637B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910161351.3A CN111650637B (en) 2019-03-04 2019-03-04 Seismic horizon interpretation method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910161351.3A CN111650637B (en) 2019-03-04 2019-03-04 Seismic horizon interpretation method and device

Publications (2)

Publication Number Publication Date
CN111650637A true CN111650637A (en) 2020-09-11
CN111650637B CN111650637B (en) 2022-12-02

Family

ID=72346169

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910161351.3A Active CN111650637B (en) 2019-03-04 2019-03-04 Seismic horizon interpretation method and device

Country Status (1)

Country Link
CN (1) CN111650637B (en)

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5056066A (en) * 1990-06-25 1991-10-08 Landmark Graphics Corporation Method for attribute tracking in seismic data
US20020042678A1 (en) * 2000-08-07 2002-04-11 Dimitri Bevc Typing picks to horizons in migration velocity analysis
US20020049540A1 (en) * 2000-08-07 2002-04-25 Dimitri Bevc Velocity analysis using angle-domain common image gathers
CN101359056A (en) * 2007-07-31 2009-02-04 中国石油天然气集团公司 Method for generating longitudinal wave time domain high precision converted wave sectional plane
CN104181596A (en) * 2014-08-27 2014-12-03 中国石油集团东方地球物理勘探有限责任公司 Geologic horizon automatic tracking method and device
CN105370269A (en) * 2015-05-13 2016-03-02 长江大学 Compaction-correction-based isochronous stratigraphic interface tracking and contrasting method
EP3067718A1 (en) * 2015-03-12 2016-09-14 CGG Services SA Boundary layer tomography method and device
CN106291681A (en) * 2015-05-30 2017-01-04 中国石油化工股份有限公司 Seismic horizon method for automatic tracking and device
CN107678060A (en) * 2017-08-29 2018-02-09 中国石油天然气股份有限公司 A kind of strata slicing method for optimizing
CN107783187A (en) * 2016-08-24 2018-03-09 中国石油化工股份有限公司 It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5056066A (en) * 1990-06-25 1991-10-08 Landmark Graphics Corporation Method for attribute tracking in seismic data
US20020042678A1 (en) * 2000-08-07 2002-04-11 Dimitri Bevc Typing picks to horizons in migration velocity analysis
US20020049540A1 (en) * 2000-08-07 2002-04-25 Dimitri Bevc Velocity analysis using angle-domain common image gathers
CN101359056A (en) * 2007-07-31 2009-02-04 中国石油天然气集团公司 Method for generating longitudinal wave time domain high precision converted wave sectional plane
CN104181596A (en) * 2014-08-27 2014-12-03 中国石油集团东方地球物理勘探有限责任公司 Geologic horizon automatic tracking method and device
EP3067718A1 (en) * 2015-03-12 2016-09-14 CGG Services SA Boundary layer tomography method and device
CN105370269A (en) * 2015-05-13 2016-03-02 长江大学 Compaction-correction-based isochronous stratigraphic interface tracking and contrasting method
CN106291681A (en) * 2015-05-30 2017-01-04 中国石油化工股份有限公司 Seismic horizon method for automatic tracking and device
CN107783187A (en) * 2016-08-24 2018-03-09 中国石油化工股份有限公司 It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field
CN107678060A (en) * 2017-08-29 2018-02-09 中国石油天然气股份有限公司 A kind of strata slicing method for optimizing

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
LILI LI ET AL.: "New Method of Horizon Recognition in Seismic Data", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 *
QI JINCHEGNG ET AL.: "An application of 3D seismic technology in tapping potential of remaining oil of narrow channel sand body in Saertu oil field", 《SPG/SEG BEIJING 2016 INTERNATIONAL GEOPHYSICAL CONFERENCE》 *
刘丽峰等: "三维精细构造解释的方法流程和关键技术", 《地球物理学进展》 *
曹来圣等: "吐哈地区三维地震资料连片处理方法研究", 《石油地球物理勘探》 *
郝兰英等: "地震沉积学在大庆长垣密井网条件下储层精细描述", 《地学前缘( 中国地质大学( 北京);北京大学)》 *
郝永庆等: "南四—八区井震联合构造精细建模及应用研究", 《中国优秀博硕士学位论文全文数据库(硕士) 基础科学辑》 *
韩宏伟等: "三维速度场建立方法", 《油气地质与采收率》 *

Also Published As

Publication number Publication date
CN111650637B (en) 2022-12-02

Similar Documents

Publication Publication Date Title
CN109166150B (en) Pose acquisition method and device storage medium
CN110174699B (en) Method, device and storage medium for determining fault closure
CN109696706B (en) Salt structure detection method, device and storage medium
CN110986930A (en) Equipment positioning method and device, electronic equipment and storage medium
CN111005715A (en) Method and device for determining gas well yield and storage medium
CN110570465A (en) real-time positioning and map construction method and device and computer readable storage medium
CN109281648B (en) Method and apparatus for determining a reasonable well pattern density of an oil reservoir
CN110633336B (en) Method and device for determining laser data search range and storage medium
CN111127541A (en) Vehicle size determination method and device and storage medium
CN112922589A (en) Pinch-out line determining method, device, terminal and storage medium
CN109281664B (en) Method and device for predicting oil production of horizontal well
CN109032492B (en) Song cutting method and device
CN110673214B (en) Method and device for predicting depths of inlet target point and end point of horizontal well
CN111650637B (en) Seismic horizon interpretation method and device
CN115545592A (en) Display positioning method, device, equipment and storage medium
CN109343113B (en) Method and device for predicting oil and gas reservoir position
CN110058306B (en) Method, apparatus and computer storage medium for determining three-dimensional velocity volume
CN112329355B (en) Method and device for determining single-well control area, computer equipment and storage medium
CN109116415B (en) Seismic wave data separation method, device and storage medium
CN113009573B (en) Seismic data acquisition method, device, computer equipment and storage medium
CN112525554B (en) Method and device for determining collision angle of automobile and computer storage medium
CN115263244B (en) Method and device for controlling drainage of air-water layer and computer storage medium
CN111597285B (en) Road network splicing method and device, electronic equipment and storage medium
CN115496828A (en) Method and device for determining stratum inclination angle parameters, computer equipment and storage medium
CN114254687A (en) Method, device and equipment for determining matching degree of drilling track and storage medium

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
CB03 Change of inventor or designer information

Inventor after: Jiang Yan

Inventor after: Li Hongxing

Inventor after: Song Baoquan

Inventor after: Qi Jincheng

Inventor after: Yang Huidong

Inventor after: Zhou Huajian

Inventor after: Pang Chunhong

Inventor after: Jiao Yanli

Inventor before: Li Hongxing

Inventor before: Song Baoquan

Inventor before: Qi Jincheng

Inventor before: Jiang Yan

Inventor before: Yang Huidong

Inventor before: Zhou Huajian

Inventor before: Pang Chunhong

Inventor before: Jiao Yanli

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant