CN110873896B - Method and system for intelligently analyzing seismic trace head information - Google Patents
Method and system for intelligently analyzing seismic trace head information Download PDFInfo
- Publication number
- CN110873896B CN110873896B CN201811017860.0A CN201811017860A CN110873896B CN 110873896 B CN110873896 B CN 110873896B CN 201811017860 A CN201811017860 A CN 201811017860A CN 110873896 B CN110873896 B CN 110873896B
- Authority
- CN
- China
- Prior art keywords
- track
- head information
- information
- data
- list
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000000007 visual effect Effects 0.000 claims description 56
- 238000010586 diagram Methods 0.000 claims description 34
- 238000012986 modification Methods 0.000 claims description 16
- 230000004048 modification Effects 0.000 claims description 16
- 230000001174 ascending effect Effects 0.000 claims description 7
- 238000012545 processing Methods 0.000 claims description 5
- 238000005070 sampling Methods 0.000 description 13
- 238000013523 data management Methods 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000007726 management method Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000003908 quality control method Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a method and a system for intelligently analyzing seismic trace head information, and belongs to the field of oil-gas geophysical exploration. The method comprises the following steps: the first step is as follows: scanning the channel head information of the seismic data to obtain the effective field of each channel head information; the second step is that: determining the track head information to be modified according to the interpretation data, determining the position of the track head information to be modified according to the effective field of each track head information obtained in the first step, and then modifying the track head information to be matched with the interpretation data; the third step: displaying a seismic data profile according to the modified trace head information, and determining a storage mode of seismic data according to the seismic data profile; the fourth step: and displaying the seismic data profile according to the storage mode determined in the third step. The method and the system can efficiently and intelligently analyze the seismic data, can acquire the position of the key parameter, not only improve the efficiency, but also ensure the correctness of the data.
Description
Technical Field
The invention belongs to the field of oil and gas geophysical exploration, and particularly relates to a method and a system for intelligently analyzing seismic trace head information.
Background
Seismic data are generally organized in units of seismic traces (trace) and stored in SEG-Y file format. The SEG-Y format is one of the standard tape data formats proposed by SEG (Society of Exploration geography scientists), which is one of the most common formats of seismic data in the oil Exploration industry. The standard SEG-Y file generally comprises three parts, the first part is an EBCDIC file header (3200 bytes) consisting of 40 cards (for example: 80 characters per line: 40 lines) for storing some information describing the seismic data volume; the second part is a binary file header (400 bytes) used for storing some key information describing the SEG-Y file, including some information such as the data format, the sampling point number, the sampling interval, the measurement unit and the like of the SEG-Y file, and the information is generally stored at a fixed position of the binary file header; the third part is the actual seismic traces, each containing 240 bytes of trace header information and seismic trace data. The trace header information generally stores information such as line number, trace number, sampling point number, geodetic coordinates and the like corresponding to the seismic trace, but some key parameter positions (such as the line number and the position of the trace number in the trace header) are not fixed, which brings inconvenience to subsequent explanation work. The seismic trace data is obtained by sampling the waveform of the seismic signal at a certain time interval delta t and recording the series of discrete amplitude values in a certain mode. The seismic data format can be IBM floating point type, IEEE floating point type, integer type, long integer type and the like, and the seismic data format processed by one three-dimensional seismic work area at the same time is unique. The number of seismic channel sampling points is determined by the number of sampling points in the seismic channel header information, all the seismic channel sampling points of most SEG-Y files are consistent, but the situation that the number of different seismic channel sampling points is different exists, the SEG-Y file is generally called as a SEG-Y file with a variable channel length format, and the file can cause the problem of warehousing.
The method is characterized in that the introduction of seismic data is the key point of explanation for seismic interpreters, and at present, TB-level data are frequently seen no matter pre-stack seismic data or post-stack seismic data, so that how to accurately analyze the trace information of the seismic data, reduce the data management cost and reduce decision deviation and loss caused by unreliable data is important.
Disclosure of Invention
The invention aims to solve the problems in the prior art and provides a method and a system for intelligently analyzing seismic trace header information, which are used for analyzing the storage mode of seismic data, analyzing the seismic data trace header information, and ensuring the correct storage of the seismic data through multi-mode quality monitoring such as graphic display, trace header modification, seismic data preview and the like.
The invention is realized by the following technical scheme:
a method for intelligently analyzing seismic trace header information comprises the following steps:
the first step is as follows: scanning the channel head information of the seismic data to obtain the effective field of each channel head information;
the second step is that: determining the track head information to be modified according to the interpretation data, determining the position of the track head information to be modified according to the effective field of each track head information obtained in the first step, and then modifying the track head information to be matched with the interpretation data;
the third step: displaying a seismic data profile according to the modified trace head information, and determining a storage mode of seismic data according to the seismic data profile;
the fourth step: and displaying the seismic data profile according to the storage mode determined in the third step.
The operation of the first step includes:
(11) reading a predetermined amount of track head information from the seismic data file in the SEGY format each time, and drawing an intuitive graph of each track head information and a list of the track head information;
(12) determining the byte position of each channel of header information by using the visual map and the list;
(13) and obtaining the effective field of each type of header information according to the byte position of each type of header information and the field byte number of the header information.
The operation of step (11) comprises:
reading a preset amount of track head information from the seismic data file in the SEGY format each time;
displaying the data of the channel head information on an intuitive graph by taking the channel number of a preset number as an abscissa and taking the value of the channel head information as an ordinate; one kind of heading information corresponds to one visual map;
displaying the value of each track corresponding to each track of each track head information in a list by taking the number of tracks of a preset number as a column and taking a byte number as a row; all track header information is displayed in a list.
The operation of step (12) comprises:
the header information includes: line number, track number, X coordinate, Y coordinate and coordinate proportion;
if the visual diagram of the track head information shows at least one horizontal straight line and the data of the track head information in the list keeps unchanged for a period of time, judging that the track head information is a line number, and determining the byte position corresponding to the track head information in the list as the byte position of the line number;
if a gradually ascending oblique line is displayed on the visual diagram of the track head information and the data of the track head information in the list shows a trend of increasing in sequence, judging that the track head information is a track number, and determining the byte position corresponding to the track head information in the list as the byte position of the track number;
if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 6 digits, judging that the track header information is an X coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the X coordinate;
if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 7 digits, judging that the track header information is Y coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the Y coordinate;
if the visual diagram of the track header information is always displayed as a straight line and the data of the track header information in the list is always kept unchanged, the track header information is judged to be in a coordinate proportion, and the byte position corresponding to the track header information in the list is determined to be the byte position of the coordinate proportion.
The operation of step (13) comprises:
and (4) obtaining the effective field of each kind of header information by utilizing the byte position of each kind of header information determined in the step (12) and according to the field byte number of each kind of header information specified by SEGY of the standard.
The second step of operation includes:
comparing the value of each channel head information obtained by scanning in the first step with the value corresponding to the channel head information in the interpretation data in sequence, and if the value of the channel head information is different from the value in the interpretation data, determining that the channel head information is the channel head information needing to be modified;
and determining the position of the track head information needing to be modified by using the effective field obtained in the first step, and modifying the value of the track head information at the position to be equal to, namely matched with, the value corresponding to the track head information in the interpretation data.
The second step further comprises:
processing the X coordinate and the Y coordinate according to the coordinate proportion: if the coordinate ratio is-1 or 0, the X and Y coordinates remain unchanged; if the coordinate ratio is 100, dividing the X and Y coordinates by 100, respectively; if the coordinate ratio is-100, then the X and Y coordinates are multiplied by 100, respectively, and so on.
The operation of the third step includes:
displaying a seismic data profile of the seismic data after the channel head information is modified according to the line direction, and if discontinuous lines appear in the seismic data profile, preliminarily determining that the seismic data are stored in the channel direction;
displaying a seismic data profile of the seismic data after the channel head information is modified according to the channel direction, and preliminarily determining that the seismic data are stored in the line direction if a vacant channel appears in the seismic data profile;
and acquiring the storage mode of the seismic data from the interpreted horizon data, comparing the storage mode with the preliminarily determined storage mode, judging that the preliminarily determined storage mode is correct if the storage mode is consistent with the preliminarily determined storage mode, and judging again if the storage mode is inconsistent with the preliminarily determined storage mode.
The invention also provides an intelligent analysis system for the seismic channel head information, which comprises:
seismic channel head information scanning unit: the method is used for scanning the channel head information of the seismic data and obtaining the effective field of each channel head information, and comprises the following steps:
the seismic channel head information scanning unit reads a preset amount of channel head information from the seismic data file in the SEGY format each time, and draws an intuitive graph and a list of the channel head information of each channel head information: reading a preset amount of track head information from the seismic data file in the SEGY format each time; displaying the data of the channel head information on an intuitive graph by taking the channel number of a preset number as an abscissa and taking the value of the channel head information as an ordinate; one kind of heading information corresponds to one visual map; displaying the value of each track corresponding to each track of each track head information in a list by taking the number of tracks of a preset number as a column and taking a byte number as a row; all track header information is displayed in a list;
the seismic channel head information scanning unit determines the byte position of each channel head information by using the visual chart and the list: if the visual diagram of the track head information shows at least one horizontal straight line and the data of the track head information in the list keeps unchanged for a period of time, judging that the track head information is a line number, and determining the byte position corresponding to the track head information in the list as the byte position of the line number; if a gradually ascending oblique line is displayed on the visual diagram of the track head information and the data of the track head information in the list shows a trend of increasing in sequence, judging that the track head information is a track number, and determining the byte position corresponding to the track head information in the list as the byte position of the track number; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 6 digits, judging that the track header information is an X coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the X coordinate; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 7 digits, judging that the track header information is Y coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the Y coordinate; if the visual graph of the road header information is always displayed as a straight line and the data of the road header information in the list is always kept unchanged, judging that the road header information is in a coordinate proportion, and determining the byte position corresponding to the road header information in the list as the byte position of the coordinate proportion;
the seismic channel head information scanning unit obtains the effective field of each channel head information according to the byte position of each channel head information and the field byte number of the channel head information: obtaining an effective field of each kind of header information by using the determined byte position of each kind of header information and according to the field byte number of each kind of header information specified by the SEGY standard;
a track head data modification unit: the device is connected with the seismic channel head information scanning unit and used for determining channel head information needing to be modified according to the interpretation data, determining the position of the channel head information needing to be modified according to the effective field of each channel head information obtained by the seismic channel head information scanning unit, and then modifying the channel head information to be matched with the interpretation data, wherein the specific steps are as follows: comparing the value of each channel head information scanned by the seismic channel head information scanning unit with the value corresponding to the channel head information in the interpretation data in sequence, and if the value of the channel head information is different from the value in the interpretation data, determining that the channel head information is the channel head information needing to be modified; determining the position of the channel head information needing to be modified by using the effective field obtained by the seismic channel head information scanning unit, and modifying the value of the channel head information at the position to enable the value to be equal to, namely matched with, the value corresponding to the channel head information in the interpretation data;
the track head data modification unit can also process the X coordinate and the Y coordinate according to the coordinate proportion: if the coordinate ratio is-1 or 0, the X and Y coordinates remain unchanged; if the coordinate ratio is 100, dividing the X and Y coordinates by 100, respectively; if the coordinate ratio is-100, then the X and Y coordinates are multiplied by 100, respectively, and so on.
Storage mode determination means: and the head data modification unit is connected with the head data modification unit and is used for displaying the seismic data profile according to the modified head information and determining the storage mode of the seismic data according to the seismic data profile, wherein the storage mode specifically comprises the following steps: displaying a seismic data profile of the seismic data after the channel head information is modified according to the line direction, and if discontinuous lines appear in the seismic data profile, preliminarily determining that the seismic data are stored in the channel direction; displaying a seismic data profile of the seismic data after the channel head information is modified according to the channel direction, and preliminarily determining that the seismic data are stored in the line direction if a vacant channel appears in the seismic data profile; obtaining the storage mode of the seismic data from the interpreted horizon data, comparing the storage mode with the preliminarily determined storage mode, if the two storage modes are consistent, judging that the preliminarily determined storage mode is correct, and if the two storage modes are not consistent, judging again;
seismic data profile display unit: and the storage mode determining unit is connected with the seismic data acquisition unit and is used for displaying the seismic data profile according to the storage mode determined by the storage mode determining unit.
The present invention also provides a computer-readable storage medium storing at least one program executable by a computer, the at least one program, when executed by the computer, causing the computer to perform the steps in the method of intelligent analysis of seismic header information of the present invention.
Compared with the prior art, the invention has the beneficial effects that:
the method and the system can efficiently and intelligently analyze the seismic data, can acquire the position of the key parameter, not only improve the efficiency, but also ensure the correctness of the data.
Drawings
FIG. 1 is a heading information display interface of the system of the present invention;
FIG. 2 is a seismic data profile display interface of the system of the present invention;
FIG. 3 is a diagram of a header data modification interface of the system of the present invention;
FIG. 4 is an interface for backup after header data modification in the system of the present invention;
FIG. 5 is a cross-sectional view of seismic data after modification of the trace header data as displayed by the system of the present invention;
FIG. 6 is a block diagram of the steps of the method of the present invention;
fig. 7 is a block diagram of the system of the present invention.
Detailed Description
The invention is described in further detail below with reference to the accompanying drawings:
as shown in fig. 6, the method of the present invention comprises:
the first step is as follows: scanning seismic channel head information to obtain an effective field of each channel head information;
analyzing a track header information file (namely track header information of a third part in an SEG-Y file of pre-stack seismic data or post-stack seismic data obtained by performing conventional processing on original ground seismic data acquired in the field), and determining the number of bytes of fields of each track header information such as a line number and a track number through data and graphic multimode (including data, an intuitive graph and a seismic data profile in the track header information file) quality monitoring, wherein the number of bytes is as follows:
the SEGY-formatted seismic data file is analyzed, a predetermined number of pieces of header information are read each time, and a visual map of the header information (as shown in a coordinate diagram below the "header information" in FIG. 1) is drawn, each visual map of the header information representing the value of the header information for each of a predetermined number of (2000 in this embodiment) tracks. The abscissa of the visual map is 2000 tracks, and the ordinate is the corresponding field value range, for example: if the line number is a visual map of the line number, the ordinate is the range of the line number in each line, if the line number is a visual map of the line number, the ordinate is the range of the line number, if the line number is a visual map of the X coordinate, the ordinate is the range of the value of the X coordinate, and so on.
For example, every time 2000 tracks are checked, the track number is used as a column, the track head information data of each track is arranged in a list mode by taking a byte number as a row, as shown in a table below a coordinate diagram in fig. 1, the track head information data of each track corresponds to a column of data in the table, and it can be found that the track number shows a trend of ascending on the data in sequence due to 2000 seismic tracks, while the line number does not change when being kept on the data for a period of time, a visual diagram of the track number is drawn, the track number of the 2000 track is shown on the visual diagram, and the track number shows a trend of ascending in sequence due to the track number, so that the track number is displayed as an oblique line in the visual diagram; drawing a line number visual graph, wherein the line number of 2000 traces is displayed on the visual graph, at least one horizontal straight line appears on the visual graph because the line number is kept unchanged for a period of time, if the line number is a plurality of straight lines, the straight lines are in a step shape, as shown in the visual graph in fig. 1, the first straight line in the visual graph in fig. 1 represents all seismic traces with the line number of 105, the second straight line represents all seismic traces with the line number of 106, the third straight line represents all seismic traces with the line number of 107, and it can be seen that the line number of 2000 seismic traces is from 105 to 107. Shown in the X-coordinate diagram are the X-coordinates of each of the 2000 tracks, which are discrete points, and shown in the Y-coordinate diagram are the Y-coordinates of each of the 2000 tracks, which are also discrete points. The X coordinate is typically 6 digits and the Y coordinate is typically 7 digits, and the X and Y coordinates can be found directly from the list based on the digits of the number. The coordinate scale map is shown as a straight line because the coordinate scale value in an SEGY file does not change (e.g., -1 or-100 or-1000 or 1000, etc.). In this way, each kind of the track head information is displayed by adopting the visual map, and the type of the track head information can be clearly judged by the difference of the graphics presented by each kind of the track head information in the visual map, for example, if a slash is presented in the visual map, the track head information is the track number, and if a plurality of parallel straight lines are presented in the visual map, the track head information is the line number. Through the intuitive multi-mode quality control, the byte positions of the line number, the track number, the X coordinate, the Y coordinate, the coordinate proportion and the like can be determined. Then, the effective field of each track head information is obtained according to the determined byte position and the field byte number. The number of field bytes is well specified in the standard SEGY.
The result of the first step auto-scan is the list of fig. 1 and a visual representation of each header information.
For non-standard SEGY data, after a track head information list is obtained by scanning, the seismic data section of the SEGY data shows incomplete data, an interpreter basically judges the positions of the line number, the track number, an X coordinate, a Y coordinate and the like by looking at a visual map and a section map, the number of bytes of a field and the type of the data can be modified at the moment, and the modified data can normally display the seismic data section.
For the coordinate scale, if the coordinate scale in the header writes a value of-1 or 0, then the X and Y coordinates remain unchanged, if 100, then both the X and Y coordinates need to be divided by 100, if-100, then the X and Y coordinates need to be multiplied by 100, and so on. The processing of the coordinate values can be performed before the input into the system of the invention, or can be performed in the system. The processed coordinates can be matched with the interpretation data.
The second step is that: modifying the channel head information according to the requirement of subsequent explanation, modifying the seismic data to be matched with the explained data:
the interpretation data refers to interpreted horizon data, fault data, destination layer data and the like, the data can be stored in an X coordinate, a Y coordinate and a Time value, the X coordinate and the Y coordinate are required to be consistent with SEGY data, if the X coordinate and the Y coordinate are inconsistent with the SEGY data, subsequent interpretation work cannot be carried out, and in order to keep the data consistent, appropriate modification needs to be carried out on trace header information of the seismic data according to the interpretation data. Such as: seismic data and interpretation horizon data are kept consistent, at least three points being consistent: XYT or XLT, i.e. X-coordinate, Y-coordinate, Time value or line number, track number, Time value. Only valid values are obtained from the SEGY file, and the position of each value is determined in order to obtain their valid values, and if the value does not match with the subsequent interpretation horizon data value, the corresponding value in the track header needs to be modified.
Examples are as follows: after seismic data are collected, sometimes the seismic data are area codes, but interpretation data are area codes, sometimes the seismic data are not area codes, but the interpretation data are area codes, and the data of the seismic data and the interpretation data are inconsistent, at this time, the seismic data need to be modified. And finding the field needing to be modified through the valid field determined in the first step, and then modifying the value of the field.
In the standard SEGY file, the X and Y coordinates are at 73 and 77 bytes, respectively, and the X coordinate is: 15234567.0, the Y coordinate is: 152345678.0, the area code 15 before the value is added or removed from the geodetic coordinates (i.e., X coordinates, Y coordinates) in the seismic data, i.e., the area code 15 before the coordinates is added or removed. For example, the X coordinate shown in fig. 3 is 18351364, the first 18 is the area code, and the area code is deleted by directly subtracting 18000000 from each X coordinate. In this embodiment, the area code is removed or added by way of example, and the value of any track header information, such as the line number, the track number, the coordinate ratio, and the like, may also be modified by using the present invention.
In addition, for the SEG-Y file with the variable track length format, according to the number of the sampling points in the track header information, under the condition that the number of a certain sampling point is insufficient, the number of the sampling points is complemented, and an invalid value (-9999.99) is complemented at the value corresponding to the complemented sampling points, so that the supplemented sampling points become seismic data with the consistent track length, and the storage management is facilitated.
The header modification function may be performed on the original file, or may generate a new file, that is, generate new seismic data in the SEGY format, as shown in fig. 4, where fig. 4 shows a selection interface for selecting "modify on original file" and "save as new seismic data", if "save as" is selected, a new SEGY file is generated, the old file is retained, and if "modify on original file" is selected, the new SEGY file is directly modified on original file.
The third step: determining the storage mode of the seismic data by combining the profile of the seismic data
A seismic data profile showing a predetermined number of traces, as shown in fig. 2, and fig. 2 is a seismic data profile corresponding to the 2000 traces of seismic data in fig. 1, from which it can be basically determined whether the seismic data is stored in the segy file in the line direction or the trace direction, as follows:
displaying a seismic data section diagram according to the line direction, and if discontinuous lines appear in the seismic data section, preliminarily judging that the seismic data are stored in the channel direction;
and displaying the seismic data section according to the channel direction, and if a vacant channel appears in the seismic data section, preliminarily judging that the seismic data is stored in the line direction.
And finally, sequentially reading each channel of seismic data according to a storage mode (if the seismic data are stored in the line direction, reading according to the line number, and if the seismic data are stored in the channel direction, reading according to the channel number).
After the storage mode of the seismic data is basically determined according to the seismic data profile, whether the storage mode of the seismic data is accurate can be verified through subsequent interpretation of horizon data. X, Y, a line number, a track number and a Time value are contained in the explained level data, and a storage mode can be obtained according to X, Y, the line number and the track number in the explained level data and is used for checking whether the judgment is accurate or not, if the judgment is consistent with the judgment, the primarily determined storage mode is judged to be accurate, and if the judgment is inconsistent with the judgment, the judgment is carried out again.
As shown in FIG. 2, there are many missing traces of the 2000 traces, indicating that the seismic data was stored in an abnormal line direction.
The fourth step: seismic data profile shows: and displaying the seismic data profile according to the storage mode determined in the third step.
FIG. 2 is an original profile, and an adjusted seismic data profile is shown in FIG. 5. comparing the two profiles, it can be seen that the invention, through intelligent analysis, ensures the integrity of data, reduces the data management cost, and reduces decision deviation and loss caused by unreliable data.
The invention provides an intelligent analysis system for seismic trace header information, as shown in fig. 7, comprising:
the seismic-channel-head information scanning unit 10: the method is used for scanning the channel head information of the seismic data and obtaining the effective field of each channel head information, and comprises the following steps:
the seismic channel head information scanning unit 10 reads a predetermined amount of channel head information from the seismic data file in the SEGY format each time, and draws a visual chart of each channel head information and a list of the channel head information: reading a preset amount of track head information from the seismic data file in the SEGY format each time; displaying the data of the channel head information on an intuitive graph by taking the channel number of a preset number as an abscissa and taking the value of the channel head information as an ordinate; one kind of heading information corresponds to one visual map; displaying the value of each track corresponding to each track of each track head information in a list by taking the number of tracks of a preset number as a column and taking a byte number as a row; all track header information is displayed in a list;
the seismic channel head information scanning unit 10 determines the byte position of each channel head information by using the visual chart and the list: if the visual diagram of the track head information shows at least one horizontal straight line and the data of the track head information in the list keeps unchanged for a period of time, judging that the track head information is a line number, and determining the byte position corresponding to the track head information in the list as the byte position of the line number; if a gradually ascending oblique line is displayed on the visual diagram of the track head information and the data of the track head information in the list shows a trend of increasing in sequence, judging that the track head information is a track number, and determining the byte position corresponding to the track head information in the list as the byte position of the track number; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 6 digits, judging that the track header information is an X coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the X coordinate; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 7 digits, judging that the track header information is Y coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the Y coordinate; if the visual graph of the road header information is always displayed as a straight line and the data of the road header information in the list is always kept unchanged, judging that the road header information is in a coordinate proportion, and determining the byte position corresponding to the road header information in the list as the byte position of the coordinate proportion;
the seismic channel head information scanning unit 10 obtains the effective field of each channel head information according to the byte position of each channel head information and the field byte number of the channel head information: obtaining an effective field of each kind of header information by using the determined byte position of each kind of header information and according to the field byte number of each kind of header information specified by the SEGY standard;
the header data modification unit 20: the seismic channel head information scanning unit 10 is connected to determine the channel head information to be modified according to the interpretation data, and the effective field of each channel head information obtained by the seismic channel head information scanning unit 10 determines the position of the channel head information to be modified, and then the channel head information is modified to match with the interpretation data, which is as follows: comparing the value of each channel head information scanned by the seismic channel head information scanning unit 10 with the value corresponding to the channel head information in the interpretation data in sequence, and if the value of the channel head information is different from the value in the interpretation data, determining that the channel head information is the channel head information needing to be modified; determining the position of the channel head information needing to be modified by using the effective field obtained by the seismic channel head information scanning unit 10, and modifying the value of the channel head information at the position to be equal to, namely matched with, the value corresponding to the channel head information in the interpretation data;
the track head data modification unit 20 is further capable of processing the X-coordinate and the Y-coordinate according to a coordinate ratio: if the coordinate ratio is-1 or 0, the X and Y coordinates remain unchanged; if the coordinate ratio is 100, dividing the X and Y coordinates by 100, respectively; if the coordinate ratio is-100, then the X and Y coordinates are multiplied by 100, respectively, and so on.
Storage manner determination unit 30: and the head data modification unit 20 is connected to display the seismic data profile according to the modified head information, and determine the storage mode of the seismic data according to the seismic data profile, which includes the following specific steps: displaying a seismic data profile of the seismic data after the channel head information is modified according to the line direction, and if discontinuous lines appear in the seismic data profile, preliminarily determining that the seismic data are stored in the channel direction; displaying a seismic data profile of the seismic data after the channel head information is modified according to the channel direction, and preliminarily determining that the seismic data are stored in the line direction if a vacant channel appears in the seismic data profile; obtaining the storage mode of the seismic data from the interpreted horizon data, comparing the storage mode with the preliminarily determined storage mode, if the two storage modes are consistent, judging that the preliminarily determined storage mode is correct, and if the two storage modes are not consistent, judging again;
seismic-data profile display unit 40: connected to the deposit mode determining unit 30 for displaying the seismic data profile according to the deposit mode determined by the deposit mode determining unit.
The present invention also provides a computer-readable storage medium storing at least one program executable by a computer, the at least one program, when executed by the computer, causing the computer to perform the steps of the method for intelligent analysis of seismic head information provided by any of the embodiments of the present invention.
The above-described embodiment is only one embodiment of the present invention, and it will be apparent to those skilled in the art that various modifications and variations can be easily made based on the application and principle of the present invention disclosed in the present application, and the present invention is not limited to the method described in the above-described embodiment of the present invention, so that the above-described embodiment is only preferred, and not restrictive.
Claims (8)
1. A method for intelligently analyzing seismic trace header information is characterized by comprising the following steps: the method comprises the following steps:
the first step is as follows: scanning the channel head information of the seismic data to obtain the effective field of each channel head information;
the second step is that: determining the track head information to be modified according to the interpretation data, determining the position of the track head information to be modified according to the effective field of each track head information obtained in the first step, and then modifying the track head information to be matched with the interpretation data;
the third step: displaying a seismic data profile according to the modified trace head information, and determining a storage mode of seismic data according to the seismic data profile;
the fourth step: displaying the seismic data profile according to the storage mode determined in the third step;
the operation of the first step includes:
(11) reading a predetermined amount of track head information from the seismic data file in the SEGY format each time, and drawing an intuitive graph of each track head information and a list of the track head information;
(12) determining the byte position of each channel of header information by using the visual map and the list;
(13) obtaining an effective field of each channel head information according to the byte position of each channel head information and the field byte number of the channel head information;
the operation of step (12) comprises:
the header information includes: line number, track number, X coordinate, Y coordinate and coordinate proportion;
if the visual diagram of the track head information shows at least one horizontal straight line and the data of the track head information in the list keeps unchanged for a period of time, judging that the track head information is a line number, and determining the byte position corresponding to the track head information in the list as the byte position of the line number;
if a gradually ascending oblique line is displayed on the visual diagram of the track head information and the data of the track head information in the list shows a trend of increasing in sequence, judging that the track head information is a track number, and determining the byte position corresponding to the track head information in the list as the byte position of the track number;
if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 6 digits, judging that the track header information is an X coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the X coordinate;
if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 7 digits, judging that the track header information is Y coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the Y coordinate;
if the visual diagram of the track header information is always displayed as a straight line and the data of the track header information in the list is always kept unchanged, the track header information is judged to be in a coordinate proportion, and the byte position corresponding to the track header information in the list is determined to be the byte position of the coordinate proportion.
2. The method for intelligent analysis of seismic header information according to claim 1, wherein: the operation of step (11) comprises:
reading a preset amount of track head information from the seismic data file in the SEGY format each time;
displaying the data of the channel head information on an intuitive graph by taking the channel number of a preset number as an abscissa and taking the value of the channel head information as an ordinate; one kind of heading information corresponds to one visual map;
displaying the value of each track corresponding to each track of each track head information in a list by taking the number of tracks of a preset number as a column and taking a byte number as a row; all track header information is displayed in a list.
3. The method for intelligent analysis of seismic header information according to claim 2, wherein: the operation of step (13) comprises:
and (4) obtaining the effective field of each kind of header information by utilizing the byte position of each kind of header information determined in the step (12) and according to the field byte number of each kind of header information specified by SEGY of the standard.
4. The method for intelligent analysis of seismic header information according to claim 3, wherein: the second step of operation includes:
comparing the value of each channel head information obtained by scanning in the first step with the value corresponding to the channel head information in the interpretation data in sequence, and if the value of the channel head information is different from the value in the interpretation data, determining that the channel head information is the channel head information needing to be modified;
and determining the position of the track head information needing to be modified by using the effective field obtained in the first step, and modifying the value of the track head information at the position to be equal to, namely matched with, the value corresponding to the track head information in the interpretation data.
5. The method for intelligent analysis of seismic header information according to claim 4, wherein: the second step further comprises:
processing the X coordinate and the Y coordinate according to the coordinate proportion: if the coordinate ratio is-1 or 0, the X and Y coordinates remain unchanged; if the coordinate ratio is 100, dividing the X and Y coordinates by 100, respectively; if the coordinate ratio is-100, then the X and Y coordinates are multiplied by 100, respectively, and so on.
6. The method for intelligent analysis of seismic header information according to claim 5, wherein: the operation of the third step includes:
displaying a seismic data profile of the seismic data after the channel head information is modified according to the line direction, and if discontinuous lines appear in the seismic data profile, preliminarily determining that the seismic data are stored in the channel direction;
displaying a seismic data profile of the seismic data after the channel head information is modified according to the channel direction, and preliminarily determining that the seismic data are stored in the line direction if a vacant channel appears in the seismic data profile;
and acquiring the storage mode of the seismic data from the interpreted horizon data, comparing the storage mode with the preliminarily determined storage mode, judging that the preliminarily determined storage mode is correct if the storage mode is consistent with the preliminarily determined storage mode, and judging again if the storage mode is inconsistent with the preliminarily determined storage mode.
7. The utility model provides a seismic channel head information intelligent analysis system which characterized in that: the system comprises:
seismic channel head information scanning unit: the method comprises the steps of scanning the channel head information of the seismic data, and obtaining the effective field of each channel head information;
a track head data modification unit: the earthquake channel head information scanning unit is connected with the earthquake channel head information scanning unit and is used for determining channel head information to be modified according to the interpretation data, determining the position of the channel head information to be modified according to the effective field of each channel head information obtained by the earthquake channel head information scanning unit, and then modifying the channel head information to be matched with the interpretation data;
storage mode determination means: the system is connected with the head data modification unit and used for displaying a seismic data profile according to the modified head information and determining a storage mode of the seismic data according to the seismic data profile;
seismic data profile display unit: the storage mode determining unit is connected with the seismic data acquisition unit and used for determining the storage mode of the seismic data acquisition unit;
the seismic channel head information scanning unit reads a preset amount of channel head information from the seismic data file in the SEGY format each time, and draws an intuitive graph and a list of the channel head information of each channel head information: reading a preset amount of track head information from the seismic data file in the SEGY format each time; displaying the data of the channel head information on an intuitive graph by taking the channel number of a preset number as an abscissa and taking the value of the channel head information as an ordinate; one kind of heading information corresponds to one visual map; displaying the value of each track corresponding to each track of each track head information in a list by taking the number of tracks of a preset number as a column and taking a byte number as a row; all track header information is displayed in a list;
the seismic channel head information scanning unit determines the byte position of each channel head information by using the visual chart and the list: if the visual diagram of the track head information shows at least one horizontal straight line and the data of the track head information in the list keeps unchanged for a period of time, judging that the track head information is a line number, and determining the byte position corresponding to the track head information in the list as the byte position of the line number; if a gradually ascending oblique line is displayed on the visual diagram of the track head information and the data of the track head information in the list shows a trend of increasing in sequence, judging that the track head information is a track number, and determining the byte position corresponding to the track head information in the list as the byte position of the track number; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 6 digits, judging that the track header information is an X coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the X coordinate; if a visual diagram of the track header information shows a plurality of discrete points and the data of the track header information in the list is 7 digits, judging that the track header information is Y coordinate, and determining the byte position corresponding to the track header information in the list as the byte position of the Y coordinate; if the visual diagram of the track header information is always displayed as a straight line and the data of the track header information in the list is always kept unchanged, the track header information is judged to be in a coordinate proportion, and the byte position corresponding to the track header information in the list is determined to be the byte position of the coordinate proportion.
8. A computer-readable storage medium characterized by: the computer-readable storage medium stores at least one program executable by a computer, the at least one program, when executed by the computer, causing the computer to perform the steps of the method for intelligent analysis of seismic header information of any of claims 1-6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811017860.0A CN110873896B (en) | 2018-09-03 | 2018-09-03 | Method and system for intelligently analyzing seismic trace head information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811017860.0A CN110873896B (en) | 2018-09-03 | 2018-09-03 | Method and system for intelligently analyzing seismic trace head information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110873896A CN110873896A (en) | 2020-03-10 |
| CN110873896B true CN110873896B (en) | 2021-12-31 |
Family
ID=69716505
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811017860.0A Active CN110873896B (en) | 2018-09-03 | 2018-09-03 | Method and system for intelligently analyzing seismic trace head information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110873896B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111983686A (en) * | 2020-09-01 | 2020-11-24 | 天津大学 | SEGY-based visual display method for shallow stratum profile original data |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5920828A (en) * | 1997-06-02 | 1999-07-06 | Baker Hughes Incorporated | Quality control seismic data processing system |
| CN105320513A (en) * | 2014-07-31 | 2016-02-10 | 中国石油化工股份有限公司 | I/O (input/output) method suitable for any SEG-Y seismic data |
| CN106291703A (en) * | 2016-08-31 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of binary system processing method of seismic velocity data |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8825409B2 (en) * | 2010-09-08 | 2014-09-02 | International Business Machines Corporation | Tracing seismic sections to convert to digital format |
-
2018
- 2018-09-03 CN CN201811017860.0A patent/CN110873896B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5920828A (en) * | 1997-06-02 | 1999-07-06 | Baker Hughes Incorporated | Quality control seismic data processing system |
| CN105320513A (en) * | 2014-07-31 | 2016-02-10 | 中国石油化工股份有限公司 | I/O (input/output) method suitable for any SEG-Y seismic data |
| CN106291703A (en) * | 2016-08-31 | 2017-01-04 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | A kind of binary system processing method of seismic velocity data |
Non-Patent Citations (3)
| Title |
|---|
| Converting scanned images of seismic reflection data into SEG-Y format;Daniel Sopher et al.;《Earth Sci Inform》;20171113;第241-255页 * |
| 基于VC++的SEGY数据格式地震剖面图绘制;刘继承等;《计算机应用》;20121231;第31卷(第3期);第25页 * |
| 磁盘SEG－D地震数据格式解编方法;王建锋等;《石油管材与仪器》;20160831;第2卷(第4期);第95-97页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110873896A (en) | 2020-03-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| WO2009018525A2 (en) | Method and system for dynamic, three-dimensional geological interpretation and modeling | |
| CN110874846A (en) | Color curve bitmap vectorization method, computer equipment and storage medium | |
| CN110873896B (en) | Method and system for intelligently analyzing seismic trace head information | |
| CN115131486B (en) | Engineering exploration data acquisition system and method | |
| CN113989414A (en) | Method and system for drawing cross-section diagram and computer storage medium | |
| US10422904B2 (en) | Data processing system for quantifying geologic growth history of subsurface oil field structures based on structural growth indications | |
| CN106294525A (en) | A kind of well logging columnar section information extracting method and system | |
| CN104297794A (en) | Seismic data mapping method and system | |
| CN115659700B (en) | Waveform automatic contrast analysis method, apparatus and medium based on signal dependency relationship | |
| CN114092538B (en) | Method and system for processing breakpoint area in single-well geological section | |
| CN109686077B (en) | Vehicle running condition monitoring method and device | |
| CN110346647B (en) | Three-dimensional map expression method for wide-area relative current ratio | |
| CN111751880B (en) | Well earthquake data combined display method and device | |
| CN117971931A (en) | Seismic data coordinate conversion method and device, electronic equipment and medium | |
| CN113608271A (en) | Geological exploration data processing method and system | |
| GB2588685A (en) | Hydrocarbon exploration method | |
| CN110990936A (en) | Method and system for determining rock mass quality index | |
| CN112925041A (en) | Method, device and equipment for determining development degree of reservoir fracture and storage medium | |
| CN116521777B (en) | Remote interface method for tunnel engineering monitoring data | |
| Guide | Profile Viewing and Analysis Software | |
| CN112014876B (en) | Reservoir prediction method and device based on pseudo-three-dimensional post-stack multi-attribute inversion | |
| US20220397690A1 (en) | Detection and Removal of Delayed Seismic Travel Times Produced by Velocity Inversions | |
| CN114398077A (en) | Processing method for adding known geological information into coal field seismic interpretation software | |
| CN104656135A (en) | Method and system for accurately distinguishing seismic shear waves | |
| CN113376690A (en) | Reservoir parameter prediction method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |