CN108614296B - Method and device for determining repeatability of observation system - Google Patents

Method and device for determining repeatability of observation system Download PDF

Info

Publication number
CN108614296B
CN108614296B CN201810573479.6A CN201810573479A CN108614296B CN 108614296 B CN108614296 B CN 108614296B CN 201810573479 A CN201810573479 A CN 201810573479A CN 108614296 B CN108614296 B CN 108614296B
Authority
CN
China
Prior art keywords
observation system
channel
repeatability
observation
cmp
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810573479.6A
Other languages
Chinese (zh)
Other versions
CN108614296A (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.)
Tianjin Hailong Petroleum Geophysical Exploration Co Ltd
BGP Inc
Original Assignee
Tianjin Hailong Petroleum Geophysical Exploration Co Ltd
BGP Inc
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 Tianjin Hailong Petroleum Geophysical Exploration Co Ltd, BGP Inc filed Critical Tianjin Hailong Petroleum Geophysical Exploration Co Ltd
Priority to CN201810573479.6A priority Critical patent/CN108614296B/en
Publication of CN108614296A publication Critical patent/CN108614296A/en
Application granted granted Critical
Publication of CN108614296B publication Critical patent/CN108614296B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/14—Signal detection
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/10—Aspects of acoustic signal generation or detection
    • G01V2210/16—Survey configurations

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 discloses a method and a device for determining repeatability of an observation system, and belongs to the technical field of time-lapse seismic. The method comprises the following steps: determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element; acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system; when the repeatability of the multi-channel observation system corresponding to each matching relation is calculated according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair, the repeatability of the single-channel observation system with the numerical value larger than the first numerical value is re-assigned to be the first numerical value; acquiring the minimum value in the repeatability of the multi-channel observation systems corresponding to the multiple matching relations as the repeatability of the multi-channel observation systems of each CMP surface element; and screening out a plurality of CMP surface elements with the repeatability of the observation system larger than a first value for displaying. The method not only improves the accuracy and reliability of the repeatability determination of the observation system, but also can intuitively display the surface element with the improved repeatability.

Description

Method and device for determining repeatability of observation system
Technical Field
The invention relates to the technical field of time-lapse seismic, in particular to a method and a device for determining repeatability of an observation system.
Background
The time-lapse earthquake refers to an earthquake which is repeatedly carried out in different time in the same exploration area, the difference between earthquake data in different time can reflect the dynamic change of the characteristic of an oil-gas reservoir, and the change can guide related personnel to carry out the development of oil and gas, thereby improving the exploitation benefit of the oil and gas. The repeatability of the seismic observation system usually has direct influence on the repeatability of seismic data, and the repeatability of the seismic data can be fundamentally improved by controlling the repeatability of the observation system, so that how to determine the repeatability of the observation system is very important.
The method for determining the repeatability of the observation system in the related art comprises the following steps: and determining the repeatability of the observation system according to shot point deviation and feather angle deviation between the single-channel monitoring observation system and the baseline observation system. Or, for a single shot-geophone pair in the observation system, determining the repeatability of the observation system according to shot point deviation and geophone point deviation between the monitoring observation system and the baseline observation system.
In the process of implementing the invention, the inventor finds that the related art has at least the following problems:
the method for determining the repeatability of the observation system through the single-channel observation system only utilizes partial data of the observation system, and the determined repeatability of the observation system is poor in accuracy and reliability.
Disclosure of Invention
The embodiment of the invention provides a method and a device for determining repeatability of an observation system, which can solve the problem of poor accuracy and reliability of the repeatability of the observation system in the related art. The technical scheme is as follows:
in a first aspect, a method for determining repeatability of an observation system is provided, including:
determining a multi-channel baseline observation system and a multi-channel monitoring observation system for each CMP (Common Mid-Point) bin of a plurality of CMP bins of a specified area;
acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system;
performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in the first numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relation;
performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval;
repeatedly acquiring a single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value;
for each matching relationship, when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair in each matching relationship, reassigning the repeatability of the single-channel observation system with the value larger than the first value as the first value;
acquiring the minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relations of each CMP surface element as the repeatability of the multiple observation systems of each CMP surface element;
and screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements for displaying.
In one possible implementation, the determining a multi-channel baseline observation system and a multi-channel monitor observation system for each of the plurality of CMP bins comprises:
and determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element in the plurality of CMP surface elements according to the plurality of CMP surface elements and the positions of the shot points and the demodulator probes.
In one possible implementation, the method further includes:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, distributing virtual monitoring observation systems for mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element, wherein the mismatched baseline observation systems are residual baseline observation systems after matching;
at the execution calculation multichannel observation system repeatability that every kind of matching relation corresponds and obtain when the step of multichannel observation system repeatability of every CMP surface element, will mismatch baseline observation system with virtual monitoring observation system corresponds single-channel observation system repeatability assignment is the second numerical value, obtains the first multichannel observation system repeatability of arbitrary CMP surface element.
In one possible implementation, the method further includes:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, rejecting mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element;
and after the step of calculating the repeatability of the multi-channel observation system corresponding to each matching relation and acquiring the repeatability of the multi-channel observation system of each CMP surface element is executed, the repeatability of the second multi-channel observation system of any CMP surface element is obtained.
In one possible implementation, the method further includes:
obtaining the difference value of the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP surface element;
according to the difference value corresponding to the CMP surface element, wherein the number of the multi-channel baseline observation systems is larger than that of the multi-channel monitoring observation systems, the target CMP surface element in the plurality of CMP surface elements is obtained and displayed, and the difference value corresponding to the target CMP surface element is larger than a preset threshold value.
In a second aspect, an observation system repeatability determination apparatus is provided, including:
the determining module is used for determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element in a plurality of common center point CMP surface elements of the designated area;
the acquisition module is used for acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system;
the acquisition module is further used for performing linear fitting on relation data between the repeatability of the single-channel observation system with the numerical value in the first numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relation; performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval; repeatedly acquiring a single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value;
the calculation module is used for reassigning the repeatability of the single-channel observation system with the value larger than the first value to the first value when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair in each matching relationship;
the obtaining module is further configured to obtain a minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relationships of each CMP bin as the repeatability of the multiple observation systems of each CMP bin;
and the display module is used for screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements and displaying the CMP surface elements.
In one possible implementation, the determining module is configured to determine a multi-channel baseline observation system and a multi-channel monitoring observation system for each of the plurality of CMP bins according to the plurality of CMP bins and the locations of the shot points and the detector points.
In one possible implementation, the apparatus further includes:
the system comprises a distribution module, a monitoring module and a monitoring module, wherein the distribution module is used for distributing a virtual monitoring observation system for a mismatched baseline observation system in the multi-baseline observation system of any CMP surface element when the number of the multi-baseline observation systems of the any CMP surface element is larger than that of the multi-monitoring observation systems, and the mismatched baseline observation system is a residual baseline observation system after matching;
the calculation module is also used for calculating the repeatability of the multichannel observation system corresponding to each matching relationship and acquiring the repeatability of the multichannel observation system of each CMP surface element, and the repeatability of the mismatch baseline observation system and the repeatability of the single-channel observation system corresponding to the virtual monitoring observation system are assigned as a second numerical value to obtain the repeatability of the first multichannel observation system of any CMP surface element.
In one possible implementation, the apparatus further includes:
a rejecting module, configured to reject a mismatched baseline observation system in the multi-channel baseline observation system of any CMP bin when the number of the multi-channel baseline observation systems of any CMP bin is greater than the number of the multi-channel monitoring observation systems;
the calculation module is further used for calculating the repeatability of the multiple observation systems corresponding to each matching relationship and obtaining the repeatability of the multiple observation systems of each CMP surface element, and then the repeatability of the second multiple observation systems of any CMP surface element is obtained.
In a possible implementation manner, the obtaining module is further configured to obtain a difference between the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP bin;
the display module is further used for acquiring a difference value corresponding to a CMP surface element, wherein the number of the multi-channel baseline observation systems is larger than that of the multi-channel monitoring observation systems, in the plurality of CMP surface elements, the target CMP surface element in the plurality of CMP surface elements is displayed, and the difference value corresponding to the target CMP surface element is larger than a preset threshold value.
In a third aspect, a computer device is provided, comprising a processor and a memory; the memory is used for storing a computer program; the processor is configured to execute the computer program stored in the memory to implement the method steps of any one of the implementation manners of the first aspect.
In a fourth aspect, a computer-readable storage medium is provided, in which a computer program is stored which, when being executed by a processor, carries out the method steps of any of the implementations of the first aspect.
The technical scheme provided by the embodiment of the invention has the beneficial effects that at least:
after the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element are matched, the repeatability of the multi-channel observation system is calculated for each matching relation, when the repeatability of the multi-channel observation system is calculated, the repeatability of the single-channel observation system which is larger than a first value is re-assigned, and finally, the repeatability of the multi-channel observation system of the CMP surface element is screened and displayed in a difference mode according to the repeatability of the multi-channel observation system of each CMP surface element. The method determines the repeatability of the observation system through a plurality of observation systems, comprehensively utilizes the data of the observation systems, improves the accuracy and reliability of determining the repeatability of the observation systems, and can visually display the surface elements with the improved repeatability.
Drawings
In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings needed to be used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on these drawings without creative efforts.
Fig. 1 is a flowchart of a method for determining repeatability of an observation system according to an embodiment of the present invention;
FIG. 2 is a flowchart of a method for determining repeatability of an observation system according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a relationship between observation system repeatability and seismic data repeatability provided by an embodiment of the invention;
FIG. 4 is a schematic diagram of data display during an X-relational schema processing process according to an embodiment of the present invention;
FIG. 5 is a schematic diagram illustrating a repetitive calculation result of a multi-channel observation system according to an embodiment of the present invention;
FIG. 6 is a schematic diagram of data display of repeatability of a multi-channel observation system without X node convergence according to an embodiment of the present invention;
FIG. 7 is a schematic diagram of a data display of repeatability of a multi-channel observation system without mismatch extrapolation processing according to an embodiment of the present invention;
fig. 8 is a schematic structural diagram of an observation system repeatability determination apparatus according to an embodiment of the present invention;
fig. 9 is a schematic structural diagram of an observation system repeatability determination apparatus according to an embodiment of the present invention;
fig. 10 is a schematic structural diagram of an observation system repeatability determination apparatus according to an embodiment of the present invention;
fig. 11 is a schematic structural diagram of a computer device 1100 according to an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, embodiments of the present invention will be described in detail with reference to the accompanying drawings.
Fig. 1 is a flowchart of a method for determining repeatability of an observation system according to an embodiment of the present invention. Referring to fig. 1, the method includes:
101. a multi-channel baseline observation system and a multi-channel monitor observation system for each CMP bin of a plurality of CMP bins of a specified area are determined.
102. And acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system.
103. And performing linear fitting on the relation data between the repeatability of the single-channel observation system with the numerical value in the first numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relation.
104. And performing linear fitting on the relation data between the repeatability of the single-channel observation system with the numerical value in the second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval.
105. And repeatedly acquiring the single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value.
106. For each matching relationship, when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the calculation formula of the repeatability of the single-channel observation system and the repeatability of the multi-channel observation system corresponding to each matching pair in each matching relationship, reassigning the repeatability of the single-channel observation system with the value larger than the first value as the first value.
107. And acquiring the minimum value in the repeatability of the multiple observation systems corresponding to the multiple matching relations of each CMP surface element as the repeatability of the multiple observation systems of each CMP surface element.
108. And screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements for displaying.
According to the method provided by the embodiment of the invention, after the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP (chemical mechanical polishing) surface element are matched, the repeatability of the multi-channel observation system is respectively calculated for each matching relation, when the repeatability of the multi-channel observation system is calculated, the repeatability of the single-channel observation system which is greater than a first value is re-assigned, and finally, the repeatability of the multi-channel observation system of each CMP surface element is screened and displayed in a difference mode according to the repeatability of the multi-channel observation system of each CMP surface element. The method determines the repeatability of the observation system through a plurality of observation systems, comprehensively utilizes the data of the observation systems, improves the accuracy and reliability of determining the repeatability of the observation systems, and can visually display the surface elements with the improved repeatability.
In one possible implementation, the determining a multi-channel baseline observation system and a multi-channel monitor observation system for each of the plurality of CMP bins includes:
and determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP bin in the plurality of CMP bins according to the plurality of CMP bins and the positions of the shot points and the demodulator probes.
In one possible implementation, the method further comprises:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, distributing virtual monitoring observation systems for mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element, wherein the mismatched baseline observation systems are residual baseline observation systems after matching;
and when the steps of calculating the repeatability of the multi-channel observation system corresponding to each matching relation and acquiring the repeatability of the multi-channel observation system of each CMP surface element are executed, assigning the repeatability of the single-channel observation system corresponding to the mismatch baseline observation system and the virtual monitoring observation system as a second numerical value to obtain the repeatability of the first multi-channel observation system of any CMP surface element.
In one possible implementation, the method further comprises:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, rejecting mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element;
and after the steps of calculating the repeatability of the multiple observation systems corresponding to each matching relation and acquiring the repeatability of the multiple observation systems of each CMP surface element are executed, obtaining the repeatability of a second multiple observation systems of any CMP surface element.
In one possible implementation, the method further comprises:
obtaining the difference value of the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP surface element;
and acquiring a target CMP surface element in the plurality of CMP surface elements for displaying according to the difference corresponding to the CMP surface element of which the number of the multi-channel baseline observation systems is greater than that of the multi-channel monitoring observation systems in the plurality of CMP surface elements, wherein the difference corresponding to the target CMP surface element is greater than a preset threshold value.
All the above-mentioned optional technical solutions can be combined arbitrarily to form the optional embodiments of the present invention, and are not described herein again.
Fig. 2 is a flowchart of a method for determining repeatability of an observation system according to an embodiment of the present invention. The execution subject of the method may be a computer device, see fig. 2, the method comprising:
201. and determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP bin in the plurality of CMP bins according to the plurality of CMPs of the designated area and the positions of the shot points and the demodulator probes.
In the embodiment of the invention, the designated area can be an area needing oil reservoir exploration and exploitation. For example, a given area may be provided with multiple sources (e.g., air guns) and multiple detectors, each source and each detector being located differently, the sources and different detectors may define different observation systems. The location of the source is called the shot and the location of the receivers is called the geophone. The seismic detector can convert the vibration signal into an electric signal and send the electric signal to the computer equipment for recording, and the computer equipment can obtain seismic waveform data according to the electric signal. Due to the plurality of detectors, the vibration signal emitted by the seismic source can be reflected back through a plurality of reflection points in the ground or underwater and received by the plurality of detectors.
By taking the example of one way without loss of generality every time seismic data is acquired, by moving the positions of the seismic source and the detectors multiple times, for example, moving the seismic source and each detector 50m to the left every time, a reflection point can be covered multiple times, that is, a vibration signal emitted by the seismic source after each movement can be reflected by the reflection point and received by one of the detectors. Taking the reflection point 1 as an example, when the moving object is not moving, the vibration signal emitted by the seismic source is reflected by the reflection point 1 and may be received by the detector 1, and after the first moving object, the vibration signal emitted by the seismic source is reflected by the reflection point 1 and may be received by the detector 3. Accordingly, for each reflection point, the computer device may determine one observation system according to the shot point and the detector point when the reflection point is covered each time, so that a plurality of observation systems of the CMP bin where the reflection point is located may be determined by the shot point and the detector point when the reflection point is covered for a plurality of times.
According to a simplified condition, the vertical projection of the reflection point on the shot point plane and the demodulator probe plane is the central point of each pair of the shot point and the demodulator probe, theoretically, the central point corresponding to one reflection point falls on the same point, which is called a CMP point, a rectangular (including square) area with the CMP point intervals as the side lengths in the corresponding directions is called a CMP bin, and represents the actual range of the reflection point corresponding to the corresponding relationship, and therefore, the multi-channel observation system of each CMP bin is the multi-channel observation system in which the central point corresponding to each CMP bin falls.
The baseline observation system is an observation system determined when the seismic data are acquired for the first time, and the monitoring observation system is an observation system determined when the seismic data are acquired again after a period of time (such as 2 years), wherein in the period of time, the specified area may have undergone oil and gas exploitation, so that the fluid condition in the reservoir of the specified area is changed.
It should be noted that this step 201 is one possible implementation of determining a multi-channel baseline observation system and a multi-channel monitoring observation system for each of a plurality of CMP bins of a designated area, and repeatability of seismic imaging (seismic data) can be ensured by performing repeatability analysis and evaluation on the observation systems of the bins.
202. And for each CMP surface element, acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system.
In the embodiment of the invention, for each CMP bin, the computer device may match the multi-channel baseline observation system and the multi-channel monitoring observation system of the CMP bin according to the number of the multi-channel baseline observation system and the multi-channel monitoring observation system of the CMP bin, so as to obtain a plurality of matching relationships.
In one possible implementation, when the number of the multi-channel baseline observation systems is equal to the number of the multi-channel monitoring observation systems, the multi-channel baseline observation systems are matched with the multi-channel monitoring observation systems to obtain a set of matching relationships. For example, if the number of single-channel observation systems (hereinafter, briefly described as: the number of multi-channel baseline observation systems) of the multi-channel baseline observation systems is 4 and the number of multi-channel monitoring observation systems is 4, the computing device may match the 4-channel baseline observation systems with the 4-channel monitoring observation systems one by one to obtain a set of matching relationships, where each matching relationship includes 4 matching pairs.
And when the number of the multi-channel baseline observation systems is smaller than that of the multi-channel monitoring observation systems, removing the monitoring observation systems left after matching. For example, if the number of the multi-channel baseline observation systems is 4 and the number of the multi-channel monitoring observation systems is 5, the computing device may select 4 monitoring observation systems from the 5 monitoring observation systems to match with the 4 baseline observation systems, so that 120 matching relationships may be obtained, where each matching relationship includes 4 matching pairs. The computer equipment can remove the rest 1 monitoring observation system after obtaining each matching relation, and the removed monitoring observation systems are different along with the difference of the matching relations.
In a possible implementation manner, when the number of the multi-channel baseline observation systems is greater than the number of the multi-channel monitoring observation systems, a virtual monitoring observation system is allocated to a mismatched baseline observation system in the multi-channel baseline observation systems, where the mismatched baseline observation system is a remaining baseline observation system after matching. For example, if the number of the multi-channel baseline observation systems is 4 and the number of the multi-channel monitoring observation systems is 3, the computing device may select 3 baseline observation systems from the 4 baseline observation systems to match with the 3 monitoring observation systems, so that 24 matching relationships may be obtained, each matching relationship includes 3 matching pairs, and the computer device may allocate a virtual monitoring observation system to the remaining 1 baseline observation system after obtaining each matching relationship. Further, the computer device may assign a constant value for repeatability between the mismatched baseline observation system and the assigned virtual surveillance observation system.
The matching process comprises rejection processing of the mismatch monitoring observation system, distribution processing of the mismatch baseline observation system and repeated assignment.
In one possible implementation, when the number of the multi-channel baseline observation systems is greater than the number of the multi-channel monitoring observation systems, the computer device may further reject a mismatched baseline observation system in the multi-channel baseline observation systems, and the visually mismatched baseline observation system does not exist.
It should be noted that, since each observation system includes a shot-geophone pair, that is, a shot point and a geophone point, the process of matching the multi-channel baseline observation system and the multi-channel monitoring observation system by the computer device is also a process of matching the shot-geophone pair.
203. And for each matching relationship, calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair in each matching relationship.
In the embodiment of the invention, for each matching relationship, the computer equipment can calculate the repeatability of the single-channel observation system corresponding to each matching pair. Each matching pair includes a single pass of the baseline observation system and the monitoring observation system. For example, the computer device may calculate the repeatability (denoted by d) of the single-channel observation system corresponding to each matching pair using the following formula:
d＝|ΔS|+|ΔR| (1)
wherein, | Δ S | is a model of a distance between a shot point in the monitoring observation system and a shot point in the baseline observation system, | Δ R | is a model of a distance between a geophone point in the monitoring observation system and a geophone point in the baseline observation system.
For a repetitive calculation formula of a multi-channel observation system, in a possible implementation manner, the repetitive calculation formula of the multi-channel observation system is obtained through the following process:
a measure of the repeatability of single-pass seismic data (imaging) is NRMS (Normalized rms difference), which is expressed as follows:
wherein Monitor represents the seismic data monitored by the monitoring observation system, Baseline represents the seismic data observed by the Baseline observation system, and the rms operator is defined as:
where N represents the number of sampling points of the seismic data, xiThe ith value representing the array acted upon,i.e., the difference between the ith sample of the monitored seismic data and the ith sample of the baseline seismic data.
According to the formula (2) and the formula (3) and the ideal energy assumption, a multi-channel seismic data repeatability calculation formula can be obtained as follows:
wherein f represents a matching relationship, Di(f) Representing the repeatability, p, of the single-track seismic data corresponding to the ith matching pairiThe weighting coefficient is a single-channel weighting coefficient, and has different relations with the shot-checking distance Offset according to the difference of the upper limit beta of the dynamic correction stretch coefficient; m represents the number of matching pairs. The computer equipment can calculate p by adopting specific seismic wave velocity function and stretching coefficient beta being 0.2iAnd Offset from the detected shot distance.
The computer device may establish a relationship between the seismic data repeatability of a single trace and the observation system repeatability based on experimental data of the relationship between the seismic data repeatability and the observation system repeatability, which may be a linear relationship as shown in the following equation:
Di＝kdi (5)
wherein D isiRepresenting seismic data repeatability of the ith matched pair, k being a constant, diIndicating the observation system repeatability of the ith matched pair.
Combining the above equation (4) and equation (5), the following equation can be obtained:
wherein, the calculation formula of d (f) is as follows:
wherein f represents a matching relationship; d (f) reflecting the repeatability of the multi-channel observation system of the NRMS under the matching relation f; di(f) Represents the single corresponding to the ith matching pairObserving the repeatability of the system; p is a radical ofiThe weighting coefficient is a single-channel weighting coefficient, the weighting coefficient can be determined by the data quantity of each channel participating in the operation in the seismic data repeatability calculation, and p is determined according to the difference of the upper limit beta of the dynamic correction stretch coefficientiHave different relationships with offset; m represents the number of matching pairs.
In the embodiment of the present invention, the computer device may calculate the repeatability of the multi-channel observation system corresponding to each matching relationship by using the formula (7). Wherein, the repeatability of the multi-channel observation system is the weighted root mean square of the repeatability of each single-channel observation system.
It should be noted that the expressions (6) and (7) are only reasonably present if the expression (5) is satisfied. However, the inventor finds that the relation fitting and application of the formula (5) are accurate only when the repeatability value of the observation system is within a certain range, and the formula (5) can not well accord with experimental data within an unlimited range. When the linear relationship is obtained through fitting, if the adopted data range is too large, although the overall fitting effect on a large range is better, the k value in the linear relationship is reduced, the deviation from the actual data is larger in a small range of repeatability of an observation system, and the numerical range of repeatability of the observation system serving as fitting data is increased, so that the k value in the linear relationship is directly reduced monotonously, and the selection of the k value is not objective. Therefore, the linear relationship between the seismic data repeatability and the observation system repeatability represented by the formula (5) is only suitable for the limited range of single-channel observation system repeatability.
In this regard, in the embodiment of the present invention, the computer device may perform piecewise fitting on the relation data between the single-channel seismic data repeatability and the observation system repeatability, and specifically, the computer device may perform linear fitting on the relation data between the single-channel observation system repeatability and the corresponding single-channel seismic data repeatability, the numerical value of which is in the first numerical value interval, to obtain the first linear relation. For example, the computer device may have the log value at [0, dx]Linear fitting is carried out on relation data of repeatability of the single-channel observation system and repeatability of the single-channel seismic data in the interval to obtain an inclined straight lineThe tilted straight line or the proportional function represented by the tilted straight line is denoted as a first linear relationship. In the fitting process, when the intercept of the linear relation (the numerical value of the repeatability of the single-channel seismic data when the repeatability of the single-channel observation system is 0) is not zero, fitting is carried out by forcibly restricting the intercept to be zero.
The computer equipment can also perform linear fitting on relation data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval. For example, the computer device may log a value at (d)xAnd infinity) interval, fitting a horizontal straight line, and recording a constant function represented by the horizontal straight line or the horizontal straight line as a second linear relation. Referring to FIG. 3, a graphical representation of the relationship between observation system repeatability and seismic data repeatability is provided, as shown in FIG. 3, at 0 to some offset dxA straight line with a direct proportional relation is fitted between the two, based on dxTo the value range between theoretical infinity (∞), the characteristic of small overall change of seismic data repeatability is given to dxFitting a horizontal straight line to infinity or data in the range between the maximum values of the experimental and application data, intersecting the straight line in the direct proportional relationship at the X node, and the abscissa corresponding to the X node is dx. The repeatability of the single-channel observation system is within the interval [0, dx]In the previous process, the relation between the repeatability of the single-channel seismic data and the repeatability of the single-channel observation system is a fitted first linear relation; repeatability of single channel observation system in interval (d)xInfinity), the relationship between the repeatability of the single-channel seismic data and the repeatability of the single-channel observation system is a second linear relationship, and the repeatability of the seismic data keeps a constant value and is not related to the repeatability of the observation system. The first linear relationship and the second linear relationship constitute an X-relationship model reflecting a relationship model between the single-channel seismic data repeatability of the single-channel observation system in the interval [0, ∞ ] and the single-channel observation system repeatability.
Particularly at the backIn the process of calculating the repeatability of the multi-channel observation system, the computer equipment can control (d)xInfinity) to d at the X nodexValue, the definition domain convergence of the single-pass observation system repeatability is [0, dx]So that the repeatability of the seismic data of a single channel and the repeatability of an observation system are kept to be [0, dx]The above direct proportional linear relationship, thus making equations (5), (6) and (7) generally true. That is, when the repetitive original value of the single-channel observation system is in the interval [0, dx]When the method is used, the method is used according to the original value; when the original value of the repeatability of the single-channel observation system is in the interval (dx, infinity), the repeatability of the single-channel observation system is reassigned to dx。
In the embodiment of the present invention, when the computer device calculates the repeatability of the multi-channel observation system corresponding to each matching relationship according to the multi-channel observation system repeatability calculation formula (7)), the repeatability of the single-channel observation system is evaluated according to the above characteristics, specifically, the computer device may re-assign the repeatability of the single-channel observation system having a value greater than a first value to the first value, where the first value may be a value of the repeatability of the single-channel observation system at an intersection (i.e., an X node) of a first linear relationship and a second linear relationship, such as dx. That is, when the repeatability value of the single-channel observation system corresponding to any matching pair is greater than the first value, the computer device may assign the repeatability value of the single-channel observation system corresponding to any matching pair to the first value again (referred to as convergence to the first value or convergence to the X node), and then substitute the first value into the above formula (7) to perform calculation.
In addition, for the case that the computer device allocates a virtual monitoring observation system to the mismatched baseline observation system in the multi-channel baseline observation system when the number of the multi-channel baseline observation systems is greater than the number of the multi-channel monitoring observation systems in step 202, when the computer device calculates the repeatability of the multi-channel observation system corresponding to each matching relationship according to the above formula (7), the repeatability of the single-channel observation system corresponding to the mismatched baseline observation system and the virtual monitoring observation system can be assigned to a second numerical value (e.g., the constant mentioned in step 202) (this assignment is referred to as mismatch extrapolation), and then the second numerical value is substituted into the above formula (7) for calculation, so as to obtain the repeatability of the first multi-channel observation system of the CMP bin. The first multi-channel observation system repeatability is the multi-channel observation system repeatability of the CMP surface element obtained through calculation under the condition that a virtual monitoring observation system is distributed for the mismatch baseline observation system.
In a possible implementation, the second value may be obtained by a mismatch extrapolation process, and specifically, the obtaining of the second value includes: the computer equipment repeatedly assigns a value to a third numerical value for the single-channel seismic data between the mismatch baseline observation system and the virtual monitoring observation system; and obtaining the second value according to the third value and the first linear relation. For example, a computer device may reduce noise repeatabilityAs a third value, assign to the single-channel seismic data repeatability, extrapolate the single-channel observation system repeatability to the second value according to the first linear relationship, e.g.Wherein k' is the slope of the line represented by the first linear relationship. As shown in fig. 3, the abscissa d of the extrapolated point0I.e. the second value, and the ordinate is the third value.
The repeatability of the single-channel observation system under the condition of overlarge value (larger than the first numerical value) is repeatedly assigned to be the first numerical value, and the repeatability of the single-channel observation system under the condition of mismatch is assigned to be the second numerical value, so that when the repeatability of the multi-channel observation system is calculated, the repeatability under the condition of mismatch can be distinguished from the repeatability under the condition of overlarge value. When the second numerical value is larger than the first numerical value, that is, the numerical value of the repeatability of the single-channel observation system under the mismatch condition is larger, it can be known from formula (1) that the larger the numerical value of the repeatability of the single-channel observation system is, the larger the deviation of the shot point and the deviation of the wave detection point are, that is, the lower the repeatability of the single-channel monitoring observation system and the baseline observation system is, and the repeatability of the single-channel observation system under the mismatch condition is quantitatively reflected.
In step 202, when the number of the multi-channel baseline observation systems is greater than the number of the multi-channel monitoring observation systems, the condition of mismatched baseline observation systems in the multi-channel baseline observation systems is eliminated, and the computer device calculates the repeatability of the multi-channel observation systems corresponding to each matching relationship according to the formula (2), so that the repeatability of the second multi-channel observation systems of the CMP surface element can be obtained. The second multi-channel observation system repeatability is the multi-channel observation system repeatability of the CMP surface element obtained through calculation under the condition that the mismatched baseline observation system is removed.
204. And acquiring the minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relations of each CMP surface element as the repeatability of the multiple observation systems of each CMP surface element.
In the embodiment of the present invention, for each CMP bin, after the computer device calculates the multiple matching relationships respectively and obtains the multiple observation system repeatability corresponding to each matching relationship, the optimal matching relationship in the matching relationships can be determined, the multiple observation system repeatability corresponding to the optimal matching relationship is the minimum value among the multiple observation system repeatability corresponding to all matching relationships, and the minimum value is used as the multiple observation system repeatability of the CMP bin.
For example, the minimum value is expressed as:
d＝min(d(f)) (8)
wherein d is referred to as observation system mismatch extrapolation optimal orientation matching weighted root mean square repeatability, and d (f) is multi-channel observation system repeatability corresponding to various matching relations.
The X relation mode, the first linear relation, the second linear relation and the multi-channel observation system repeatability measurement provided by the embodiment of the invention can well distinguish the influence of single-channel matching with large repeatability values and mismatching in multi-channel repeatability, and is beneficial to selecting the optimal matching.
It should be noted that the computer device performs the above steps 202 to 204 for each CMP bin, so that a multi-view system repeatability can be obtained for all CMP bins. Referring to fig. 4, a schematic diagram of data display during processing of an X-relationship model is provided, where both graphs (a) and (b) in fig. 4 are only X-node convergence processing (when calculating the repeatability of a multi-channel observation system, only a single graph with a larger value is usedThe repeatability of the road observation system is converged to d at the X nodexValue) is the repetitive data display condition of the single-channel observation system, except that (a) mismatch is marked by pure black and (b) mismatch is marked by pure white; (c) the figure is the data display situation of the single-channel observation system repeatability only in the mismatch extrapolation processing (when the multi-channel observation system repeatability is calculated, the single-channel seismic data repeatability between the mismatch baseline observation system and the virtual monitoring observation system is assigned through the mismatch extrapolation processing); (d) the figure shows the repetitive data display condition of the single-channel observation system during the processing of the complete X-relation mode.
In the embodiment of the invention, the computer equipment performs repetitive calculation on the surface element according to the relation mode (X relation mode) between the seismic data repeatability and the observation system repeatability, and researches the performance characteristics of the observation system repeatability. The application of this method is limited to multi-channel observation systems in which the center point corresponding to each bin falls.
205. And screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements for displaying.
In the embodiment of the present invention, it can be known from the above calculation method that the bin in which the repeatability of the multi-channel observation system is greater than the X node value (first value) is determined to have the baseline mismatch. The repeatability and the mismatch of the large-value single-channel observation system can increase the repeatability of the multi-channel observation system, but under the condition that the repeatability of the multi-channel observation system does not exceed the X node value, the situation that whether the mismatch exists or the repeatability of the single-channel observation system reaches the X node value cannot be judged from the repeatability of the multi-channel observation system of the bin group. When the repeatability of the multichannel observation system of the surface element exceeds the X node value, the fact that mismatch exists can be shown, and when the repeatability of the single channel observation system of all the matched pairs reaches the X node value, the repeatability of the multichannel observation system can also be the X node value. Thus, for a multi-track observation system repeatability, the X node value is a particular critical point, reaching or exceeding this value indicates that the repeatability falls within a particular range, which can be considered an unacceptable range. On the other hand, as the meaning of the repeatability value of the multi-channel observation system is the same as that of the single-channel observation system, the point that the X node value represents the worst repeatability value is reached, and the repeatability exceeding the value represents that the repeatability passes through the worst repeatability of the observation system and the mismatching is dominant. Thus, the X node value may be considered a quality level demarcation point, and the bins for which the observation system repeatability meets or exceeds this value are bins for which the observation system repeatability needs to be improved.
Accordingly, the computer device may use the first value as a critical point, and the CMP bins whose repeatability of the multi-vision system exceeds the critical point are the bins that need improvement, and the computer device may screen the CMP bins for display. And screening out the surface element with the repeatability greater than or equal to a first value from the calculated repeatability of the multi-channel observation system for displaying, thereby visually reflecting the surface element with the total influence of the deviation and the mismatch of the observation system exceeding a certain degree.
Evaluation and analysis by the method of the invention, with established relatively objective criteria, shows a small number of local areas where observation system repeatability remains to be improved. The region shown in fig. 5, in which the reproducibility is greater than the X node value, is objectively defined in display control, and is significantly distinguished from the subjectively determined and perceived region having poor reproducibility, according to which different reproducibility improvement targets can be defined.
Referring to fig. 5, a schematic diagram for displaying the repetitive calculation results of the multichannel observation system is provided, in which in fig. 5, (a) the surface elements with the repeatability greater than the X node value are marked with pure black, (b) the surface elements with the repeatability greater than the X node value are marked with pure white, and (c) the surface elements with the repeatability greater than the X node value are screened out and displayed separately.
The computer device performs a specific filtering on the calculated data, the repeatability greater than the X node value is filtered, and the corresponding positions in (a) and (b) of fig. 5 are marked with pure black and pure white. Further, the computer device may filter out bins with repeatability less than the X node value, and display bins with repeatability greater than the X node value and a distinction of their degree, as shown in fig. 5 (c).
The processing result of the embodiment of the invention can reasonably apply the relation between the single-channel seismic data repeatability and the observation system repeatability, and form a relation mode with wider application range, so that the mismatched single-channel observation system repeatability and the large-value single-channel observation system repeatability are reasonably measured and distinguished. By clearly showing the repeatability of the bin observation system where the mismatch is dominant, the basis is given for bins where the repeatability clearly needs improvement. The computer device may also display all bins affected by the mismatch and their degree of effect.
In a possible implementation manner, for the computer device in step 204, obtaining the first multi-channel observation system repeatability and the second multi-channel observation system repeatability of the CMP bin under different conditions, further, the computer device may obtain a difference between the first multi-channel observation system repeatability and the second multi-channel observation system repeatability of the CMP bin; and acquiring a target CMP surface element in the plurality of CMP surface elements for displaying according to the difference corresponding to the CMP surface element, wherein the number of the multi-channel baseline observation systems is greater than that of the multi-channel monitoring observation systems, and the difference corresponding to the target CMP surface element is greater than a preset threshold value. In this way all bins affected by the baseline mismatch and their degree of effect can be given.
Referring to fig. 6, a schematic data display diagram of a multi-channel observation system repeatability without X node convergence is provided, and the diagram (a) of fig. 6 is an observation system repeatability without X node convergence based on an X mode; (b) the figure is the repeatability of the observation system without X node convergence marked in black for the area above the extrapolated point value (the value of the repeatability of the extrapolated point observation system); (c) the graph is the repeatability of the observation system without X node convergence marked by white in the area above the extrapolation value; (d) the figure is a separate display of the repeatability of the observation system above the extrapolated point values screened; (e) the figure is a single display of the repeatability of the observation system above the screened X node value; (f) the figure is a multi-channel observation system repeatability value of the surface element between the X node and the extrapolation point and the distribution of the repeatability value.
It can be seen by comparison that if the repeatability of the X node value is not converged, the repeatability value increases greatly from the X node value, and there may be observation system repeatability values exceeding the extrapolation point, so that the mismatched repeatability value is submerged in the observation system repeatability value, as shown in (a) of fig. 6. Thus, the computer device may mark the display of the areas of repetitiveness values outside the extrapolated points, as shown in fig. 6 (b) and (c) where the display of the areas of repetitiveness values outside the extrapolated points are marked with solid black and solid white colors, respectively, and the filtered out data is shown in fig. 6 (d). Fig. 6 (e) and (f) show the multi-channel observation system repeatability values and distributions of filtered X node exterior bins and the multi-channel observation system repeatability values and distributions of bins between the X node and the exterior bin, respectively. Of course, the computer device may also display the display result of the repetitive data excluding the X node.
In addition, the computer equipment can also process the data in a mismatch-free extrapolation mode. Referring to FIG. 7, a schematic data display of repeatability of a multi-channel observation system without mismatch extrapolation is provided, and FIG. 7 (a) shows repeatability of a multi-channel observation system without mismatch extrapolation; (b) the figure shows the repeatability of a multi-channel observation system without mismatch extrapolation (mismatch is marked in pure black); (c) the figure shows the repeatability of a multi-channel observation system without mismatch extrapolation (mismatch is marked in pure white); (d) the figure shows the difference between the full X-mode and its mismatch-free extrapolated multi-channel observation system repeatability. Fig. 7 (d) shows a difference plot between the full X-mode and the repeatability of the multi-channel observation system without mismatch extrapolation, where all areas of mismatch, along with the degree of mismatch, can be shown.
The embodiment of the invention further researches the relation between the seismic data repeatability and the observation system repeatability, and provides an X relation mode for embodying the relation, so that the research on the repeatability of a plurality of observation systems is more reasonably applied, the application range of the research is more generalized, and the overall accuracy is higher. In addition, the repeatability of the multi-channel observation system calculated according to the X relation mode provided by the embodiment of the invention not only keeps the repeatability of the optimal orientation matching weighted root mean square based on the mismatching extrapolation of the observation system in the surface element, but also reflects the characteristics of the effective repeatability and the coverage degree of all the observation systems collected by the monitoring observation system, and also enables the self to contain richer information, so that the mismatching of a certain degree can be revealed, and the standard of the repeatability evaluation of the multi-channel observation system can be naturally provided. The original display of the data obtained by the method solves the general technical problem of numerical value display limitation in the data display, and has engineering technical value and certain scientific significance.
According to the method provided by the embodiment of the invention, after the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP (chemical mechanical polishing) surface element are matched, the repeatability of the multi-channel observation system is respectively calculated for each matching relation, when the repeatability of the multi-channel observation system is calculated, the repeatability of the single-channel observation system which is greater than a first value is re-assigned, and finally, the repeatability of the multi-channel observation system of each CMP surface element is screened and displayed in a difference mode according to the repeatability of the multi-channel observation system of each CMP surface element. The method determines the repeatability of the observation system through a plurality of observation systems, comprehensively utilizes the data of the observation systems, improves the accuracy and reliability of determining the repeatability of the observation systems, and can visually display the surface elements with the improved repeatability.
Fig. 8 is a schematic structural diagram of an observation system repeatability determination apparatus according to an embodiment of the present invention. Referring to fig. 8, the apparatus includes:
a determining module 801, configured to determine a multi-channel baseline observation system and a multi-channel monitoring observation system for each CMP bin of a plurality of common-center CMP bins of a designated area;
an obtaining module 802, configured to obtain multiple matching relationships between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP bin, where each matching relationship includes multiple matching pairs, and each matching pair includes a single-channel baseline observation system and a single-channel monitoring observation system;
the obtaining module 802 is further configured to perform linear fitting on relationship data between the repeatability of the single-channel observation system with a value in the first value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relationship; performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval; repeatedly acquiring a single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value;
a calculating module 803, configured to, for each matching relationship, reassign the repeatability of the single-channel observation system having a value greater than the first value to the first value when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the single-channel observation system repeatability and the multi-channel observation system repeatability calculation formula corresponding to each matching pair in each matching relationship;
the obtaining module 802 is further configured to obtain a minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relationships of each CMP bin as the repeatability of the multiple observation systems of each CMP bin;
a display module 804, configured to screen a plurality of CMP bins from the plurality of CMP bins, where the repeatability of the observation system is greater than the first value, and display the CMP bins.
In one possible implementation, the determining module 801 is configured to determine a multi-channel baseline observation system and a multi-channel monitor observation system for each CMP bin of the plurality of CMP bins according to the locations of the plurality of CMP bins and shot and detector points.
In one possible implementation, referring to fig. 9, the apparatus further includes:
the allocating module 805 is configured to, when the number of the multi-channel baseline observation systems of any CMP bin is greater than the number of the multi-channel monitoring observation systems, allocate a virtual monitoring observation system to a mismatched baseline observation system in the multi-channel baseline observation systems of the any CMP bin, where the mismatched baseline observation system is a remaining baseline observation system after matching;
the calculating module 803 is further configured to assign the repeatability of the mismatch baseline observation system and the repeatability of the single-channel observation system corresponding to the virtual monitoring observation system to a second value to obtain the repeatability of the first multi-channel observation system of any CMP bin when the steps of calculating the repeatability of the multi-channel observation system corresponding to each matching relationship and obtaining the repeatability of the multi-channel observation system of each CMP bin are performed.
In one possible implementation, referring to fig. 10, the apparatus further includes:
a rejecting module 806, configured to reject a mismatched baseline observation system in the multi-baseline observation systems of any CMP bin when the number of the multi-baseline observation systems of any CMP bin is greater than the number of the multi-monitoring observation systems;
the calculating module 803 is further configured to obtain a second multi-channel observation system repeatability of any CMP bin after the steps of calculating the multi-channel observation system repeatability corresponding to each matching relationship and obtaining the multi-channel observation system repeatability of each CMP bin are performed.
In a possible implementation manner, the obtaining module 802 is further configured to obtain a difference between the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP bin;
the display module 804 is further configured to obtain a target CMP bin of the plurality of CMP bins for displaying according to a difference corresponding to the CMP bin, where the number of the plurality of baseline observation systems is greater than the number of the plurality of monitoring observation systems, and the difference corresponding to the target CMP bin is greater than a preset threshold.
In the embodiment of the invention, after the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element are matched, the repeatability of the multi-channel observation system is respectively calculated for each matching relation, when the repeatability of the multi-channel observation system is calculated, the repeatability of a single-channel observation system with the value larger than the first value is re-assigned, and finally, the repeatability of the multi-channel observation system of the CMP surface element is screened and displayed in a difference mode according to the repeatability of the multi-channel observation system of each CMP surface element. The method determines the repeatability of the observation system through a plurality of observation systems, comprehensively utilizes the data of the observation systems, improves the accuracy and reliability of determining the repeatability of the observation systems, and can visually display the surface elements with the improved repeatability.
It should be noted that: the observation system repeatability determining apparatus provided in the above embodiment is only illustrated by dividing the functional modules when determining the repeatability of the observation system, and in practical applications, the function distribution may be completed by different functional modules according to needs, that is, the internal structure of the device is divided into different functional modules to complete all or part of the functions described above. In addition, the observation system repeatability determination device and the observation system repeatability determination method provided by the above embodiments belong to the same concept, and specific implementation processes thereof are described in the method embodiments in detail and are not described herein again.
Fig. 11 is a schematic structural diagram of a computer device 1100 according to an embodiment of the present invention, where the computer device 1100 may have a relatively large difference due to different configurations or performances, and may include one or more processors (CPUs) 1101 and one or more memories 1102, where the memory 1102 stores at least one instruction, and the at least one instruction is loaded and executed by the processors 1101 to implement the methods provided by the above method embodiments. Certainly, the computer device may further have components such as a wired or wireless network interface, a keyboard, and an input/output interface, so as to perform input and output, and the computer device may further include other components for implementing the functions of the device, which is not described herein again.
In an exemplary embodiment, a computer-readable storage medium, such as a memory, storing a computer program is also provided, which when executed by a processor implements the observation system repeatability determination method in the above embodiments. For example, the computer-readable storage medium may be a Read-Only Memory (ROM), a Random Access Memory (RAM), a Compact Disc Read-Only Memory (CD-ROM), a magnetic tape, a floppy disk, an optical data storage device, and the like.
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, and 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 present invention is not limited to the above preferred embodiments, and any modifications, equivalent replacements, improvements, etc. within the spirit and principle of the present invention should be included in the protection scope of the present invention.

Claims (10)

1. An observation system repeatability determination method, characterized in that the method comprises:
determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element in a plurality of CMP surface elements of a designated area;
acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system;
performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in the first numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relation;
performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval;
repeatedly acquiring a single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value;
for each matching relationship, when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair in each matching relationship, reassigning the repeatability of the single-channel observation system with the value larger than the first value as the first value;
acquiring the minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relations of each CMP surface element as the repeatability of the multiple observation systems of each CMP surface element;
and screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements for displaying.
2. The method of claim 1, wherein determining a multi-channel baseline observation system and a multi-channel monitor observation system for each CMP bin of the plurality of CMP bins comprises:
and determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element in the plurality of CMP surface elements according to the plurality of CMP surface elements and the positions of the shot points and the demodulator probes.
3. The method of claim 1, further comprising:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, distributing virtual monitoring observation systems for mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element, wherein the mismatched baseline observation systems are residual baseline observation systems after matching;
at the execution calculation multichannel observation system repeatability that every kind of matching relation corresponds and obtain when the step of multichannel observation system repeatability of every CMP surface element, will mismatch baseline observation system with virtual monitoring observation system corresponds single-channel observation system repeatability assignment is the second numerical value, obtains the first multichannel observation system repeatability of arbitrary CMP surface element.
4. The method of claim 3, further comprising:
when the number of the multi-channel baseline observation systems of any CMP surface element is larger than that of the multi-channel monitoring observation systems, rejecting mismatched baseline observation systems in the multi-channel baseline observation systems of any CMP surface element;
and after the step of calculating the repeatability of the multi-channel observation system corresponding to each matching relation and acquiring the repeatability of the multi-channel observation system of each CMP surface element is executed, the repeatability of the second multi-channel observation system of any CMP surface element is obtained.
5. The method of claim 4, further comprising:
obtaining the difference value of the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP surface element;
according to the difference value corresponding to the CMP surface element, wherein the number of the multi-channel baseline observation systems is larger than that of the multi-channel monitoring observation systems, the target CMP surface element in the plurality of CMP surface elements is obtained and displayed, and the difference value corresponding to the target CMP surface element is larger than a preset threshold value.
6. An observation system repeatability determination apparatus, the apparatus comprising:
the determining module is used for determining a multi-channel baseline observation system and a multi-channel monitoring observation system of each CMP surface element in a plurality of CMP surface elements of the designated area;
the acquisition module is used for acquiring a plurality of matching relations between the multi-channel baseline observation system and the multi-channel monitoring observation system of each CMP surface element, wherein each matching relation comprises a plurality of matching pairs, and each matching pair comprises a single-channel baseline observation system and a single-channel monitoring observation system;
the acquisition module is further used for performing linear fitting on relation data between the repeatability of the single-channel observation system with the numerical value in the first numerical value interval and the repeatability of the corresponding single-channel seismic data to obtain a first linear relation; performing linear fitting on relational data between the repeatability of the single-channel observation system with the numerical value in a second numerical value interval and the corresponding repeatability of the single-channel seismic data to obtain a second linear relation, wherein the left boundary of the second numerical value interval is the right boundary of the first numerical value interval; repeatedly acquiring a single-channel observation system corresponding to the intersection point of the two straight lines represented by the first linear relation and the second linear relation as a first numerical value;
the calculation module is used for reassigning the repeatability of the single-channel observation system with the value larger than the first value to the first value when calculating the repeatability of the multi-channel observation system corresponding to each matching relationship according to the repeatability calculation formula of the single-channel observation system and the multi-channel observation system corresponding to each matching pair in each matching relationship;
the obtaining module is further configured to obtain a minimum value of the repeatability of the multiple observation systems corresponding to the multiple matching relationships of each CMP bin as the repeatability of the multiple observation systems of each CMP bin;
and the display module is used for screening a plurality of CMP surface elements with the repeatability of the observation system larger than the first value from the plurality of CMP surface elements and displaying the CMP surface elements.
7. The apparatus of claim 6, wherein the determining module is configured to determine a multi-channel baseline observation system and a multi-channel monitor observation system for each of the plurality of CMP bins based on the locations of the plurality of CMP bins and shot and geophone points.
8. The apparatus of claim 6, further comprising:
the system comprises a distribution module, a monitoring module and a monitoring module, wherein the distribution module is used for distributing a virtual monitoring observation system for a mismatched baseline observation system in the multi-baseline observation system of any CMP surface element when the number of the multi-baseline observation systems of the any CMP surface element is larger than that of the multi-monitoring observation systems, and the mismatched baseline observation system is a residual baseline observation system after matching;
the calculation module is also used for calculating the repeatability of the multichannel observation system corresponding to each matching relationship and acquiring the repeatability of the multichannel observation system of each CMP surface element, and the repeatability of the mismatch baseline observation system and the repeatability of the single-channel observation system corresponding to the virtual monitoring observation system are assigned as a second numerical value to obtain the repeatability of the first multichannel observation system of any CMP surface element.
9. The apparatus of claim 8, further comprising:
a rejecting module, configured to reject a mismatched baseline observation system in the multi-channel baseline observation system of any CMP bin when the number of the multi-channel baseline observation systems of any CMP bin is greater than the number of the multi-channel monitoring observation systems;
the calculation module is further used for calculating the repeatability of the multiple observation systems corresponding to each matching relationship and obtaining the repeatability of the multiple observation systems of each CMP surface element, and then the repeatability of the second multiple observation systems of any CMP surface element is obtained.
10. The apparatus of claim 9,
the acquisition module is further used for acquiring a difference value between the repeatability of the first multi-channel observation system and the repeatability of the second multi-channel observation system of any CMP surface element;
the display module is further used for acquiring a difference value corresponding to a CMP surface element, wherein the number of the multi-channel baseline observation systems is larger than that of the multi-channel monitoring observation systems, in the plurality of CMP surface elements, the target CMP surface element in the plurality of CMP surface elements is displayed, and the difference value corresponding to the target CMP surface element is larger than a preset threshold value.
CN201810573479.6A 2018-06-06 2018-06-06 Method and device for determining repeatability of observation system Active CN108614296B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810573479.6A CN108614296B (en) 2018-06-06 2018-06-06 Method and device for determining repeatability of observation system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810573479.6A CN108614296B (en) 2018-06-06 2018-06-06 Method and device for determining repeatability of observation system

Publications (2)

Publication Number Publication Date
CN108614296A CN108614296A (en) 2018-10-02
CN108614296B true CN108614296B (en) 2019-12-31

Family

ID=63664677

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810573479.6A Active CN108614296B (en) 2018-06-06 2018-06-06 Method and device for determining repeatability of observation system

Country Status (1)

Country Link
CN (1) CN108614296B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112198550B (en) * 2019-07-08 2023-09-26 中国石油天然气集团有限公司 Time-lapse seismic data repeatability measurement method and device
CN110824561B (en) * 2019-11-14 2021-11-23 中国石油集团东方地球物理勘探有限责任公司 Method for measuring repeatability of time-lapse seismic data
CN113568042B (en) * 2021-07-21 2024-01-26 中海石油(中国)有限公司 Streamer acquisition method and system for improving repeatability of time-lapse seismic data

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004264B (en) * 2010-10-18 2015-09-23 中国石油化工股份有限公司 A kind of earthquake-capturing data quality quantitative test and evaluation method
EP2815254A4 (en) * 2012-02-14 2016-06-15 Halliburton Energy Services Inc Shear wave source for vsp and surface seismic exploration
US20140003191A1 (en) * 2012-06-28 2014-01-02 Chevron U.S.A. Inc System and method for evaluating a time-lapse seismic signal recording using shifted normalized root mean square metric
EP3014304A2 (en) * 2013-06-28 2016-05-04 CGG Services SA System and method for estimating repeatability using base data
CN105301650A (en) * 2015-10-09 2016-02-03 中国石油天然气集团公司 Quality monitoring method and quality monitoring device for time-lapse seismic data acquired by marine towline
CN105549083B (en) * 2015-12-25 2018-03-13 中国石油天然气集团公司 A kind of method and apparatus for determining seismic observation system repeatability and integrally measuring
CN105652344A (en) * 2016-01-08 2016-06-08 中国神华能源股份有限公司 Coal mining-based quantitative evaluation method for four-dimensional seismic observation system

Also Published As

Publication number Publication date
CN108614296A (en) 2018-10-02

Similar Documents

Publication Publication Date Title
CN108614296B (en) Method and device for determining repeatability of observation system
Marshall et al. Modelling the Galactic interstellar extinction distribution in three dimensions
MacGillivray et al. New results with the COSMOS machine
US5787201A (en) High order fractal feature extraction for classification of objects in images
CA2515345C (en) Estimating the time of arrival of a seismic wave
Buyer et al. Calculation the spacing of discontinuities from 3D point clouds
JP3581149B2 (en) Method and apparatus for identifying an object using a regular sequence of boundary pixel parameters
US4868883A (en) Analysis of thin section images
EP3588061B1 (en) Optical inspection system and method including accounting for variations of optical path length within a sample
US7119909B2 (en) Film thickness and boundary characterization by interferometric profilometry
WO2018151680A1 (en) Methods and devices for identifying population clusters in data
JP2020514762A (en) Particle characterization
JP2020169994A5 (en)
WO2017107862A1 (en) Method and device for determining overall measurement of seismic observation system repeatability
JPH0521180B2 (en)
Lin et al. BEHAVIOR OF APPARENT DISPERSION CURVE AND ITS IMPLICATION TO MASW TESTING.
US8360979B2 (en) Apparatus and method for ultrasonic color imaging
JP4398819B2 (en) Method and apparatus for identifying surface properties
CN109272491A (en) The recognition methods of crack tip, device and equipment under experimental enviroment
Barry et al. Observer reproducibility during computer-assisted planimetric measurements of nuclear features
CN116359185A (en) Stream type lattice data analysis method and device and computer equipment
US20140324358A1 (en) Surface multiple prediction
CN111626979B (en) Pipe diameter measuring method and device
TWI677323B (en) Method for evaluating x-ray imaging quality
CN114235347A (en) Lens quality evaluation method and device

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