WO2015100544A1 - 基于零偏垂直地震剖面数据估计品质因子的方法和装置 - Google Patents

基于零偏垂直地震剖面数据估计品质因子的方法和装置 Download PDF

Info

Publication number
WO2015100544A1
WO2015100544A1 PCT/CN2013/090914 CN2013090914W WO2015100544A1 WO 2015100544 A1 WO2015100544 A1 WO 2015100544A1 CN 2013090914 W CN2013090914 W CN 2013090914W WO 2015100544 A1 WO2015100544 A1 WO 2015100544A1
Authority
WO
WIPO (PCT)
Prior art keywords
vsp
quality factor
transmission coefficient
coefficient
objective function
Prior art date
Application number
PCT/CN2013/090914
Other languages
English (en)
French (fr)
Inventor
于永才
孙夕平
张明
崔兴福
侯连华
Original Assignee
中国石油天然气股份有限公司
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 中国石油天然气股份有限公司 filed Critical 中国石油天然气股份有限公司
Priority to PCT/CN2013/090914 priority Critical patent/WO2015100544A1/zh
Publication of WO2015100544A1 publication Critical patent/WO2015100544A1/zh
Priority to US14/919,053 priority patent/US20160091623A1/en

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
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48—Processing data
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2200/00—Details of seismic or acoustic prospecting or detecting in general
    • G01V2200/10—Miscellaneous details
    • G01V2200/14—Quality control
    • 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
    • G01V2210/161—Vertical seismic profiling [VSP]

Definitions

  • the present invention relates to the field of physical exploration technology, and in particular, to a method and apparatus for estimating a quality factor based on a zero-seismic vertical seismic profile (VSP) data.
  • VSP vertical seismic profile
  • the earth medium is incompletely elastic.
  • the absorption caused by the viscosity of the medium causes the energy attenuation and velocity dispersion of the seismic wave.
  • the inherent attenuation characteristics of the medium are usually described by the quality factor ⁇ . Due to the existence of the quality factor, the attenuation of the high-frequency energy of the seismic wave is more serious than the attenuation of the low-frequency energy.
  • the propagation speed of the high-frequency component is faster than the propagation velocity of the low-frequency component, and the deep energy of the seismic profile is weakened and the frequency band is narrowed. This leads to a reduction in the resolution of the data, which increases the difficulty of fine interpretation of seismic data.
  • inverse ⁇ filtering is an effective means for seismic data absorption attenuation compensation, which can compensate and correct energy attenuation and velocity dispersion during seismic wave propagation, thereby improving the resolution of data.
  • VSP Vertical Seismic Profiling
  • the former used VSP data to compare the time domain method such as amplitude attenuation method, rise time method, wavelet simulation method and analytical signal method, and the frequency domain method such as matching method, spectrum simulation method and spectral ratio method.
  • One method is applicable to any situation, and the effectiveness of each method depends on the quality of the data.
  • centroid frequency shift method which mainly uses the change of centroid frequency in the process of seismic wave propagation to obtain the quality factor.
  • the spectral ratio method is a commonly used method for estimating the quality factor.
  • the quality factor is estimated by using the linear relationship between the logarithm of the spectrum and the frequency, but the slope fitting is susceptible to the spectral logarithm error.
  • the influence of the quality factor estimation is affected.
  • the exponential method effectively avoids the disadvantages of the spectral logarithm.
  • the index method uses the forward matching method to estimate the quality factor, which can effectively improve the stability of the quality factor estimation.
  • the quality factor estimated by the exponential method often has a large error, and the accuracy of the quality factor estimation is difficult to guarantee. Summary of the invention
  • Embodiments of the present invention provide a method for estimating a quality factor based on zero-bias vertical seismic profile data, so as to achieve an objective of improving the estimation stability and estimation accuracy of a quality factor, and the method includes: Determining the transmission coefficient between adjacent two VSP channels based on the layer velocity of the seismic wave in the VSP data of the zero-bias vertical seismic section;
  • the transmission coefficient is determined as an indefinite coefficient in an exponential objective function, and a quality factor is estimated based on the exponential objective function.
  • the transmission velocity between adjacent two VSP tracks is determined based on the layer velocity of the seismic wave of the zero-bias VSP data, including:
  • the transmission coefficient is determined according to the following formula:
  • p k ⁇ + Vk ⁇ Vk ⁇ l
  • A represents the transmission coefficient
  • ⁇ represents the layer velocity between the i-th and the A-th VSP
  • ⁇ represents the layer velocity between the k-2th and the first VSP.
  • the transmission coefficient is determined as an indefinite coefficient in an exponential objective function
  • the quality factor is estimated according to the exponential objective function, including:
  • the quality factor corresponding to the minimum matching error is taken as the estimated quality factor.
  • the objective function is:
  • ; represents the lower frequency limit of the dominant frequency band
  • / 2 represents the upper frequency limit of the dominant frequency band
  • t k represents the A-th VSP channel
  • a represents the quality factor between the first and the fourth VSP tracks
  • G represents the matching error.
  • the transmission coefficients between adjacent two VSPs are determined based on the layer velocity of the seismic waves in the VSP data of the zero-bias vertical seismic profile, including:
  • a transmission coefficient between adjacent two VSP channels is determined according to a layer velocity of the seismic wave in the downward wave field.
  • Embodiments of the present invention provide a device for estimating a quality factor based on zero-bias vertical seismic profile data, so as to achieve an objective of improving the estimation stability and estimation accuracy of a quality factor, the device includes: a determining module, configured to determine a transmission coefficient between adjacent two VSP channels based on a layer velocity of the seismic wave in the VSP data of the zero-bias vertical seismic section;
  • an estimation module configured to determine the transmission coefficient as an indefinite coefficient in an exponential objective function, and estimate a quality factor according to the index method objective function.
  • the estimating module comprises:
  • a determining unit configured to determine a value range of the quality factor
  • a scanning unit configured to obtain a quality factor that causes a minimum matching error of the target function by means of quality factor scanning within the range of values
  • the estimating unit is configured to use the quality factor corresponding to the minimum matching error as the estimated quality factor.
  • the objective function is:
  • ; represents the lower frequency limit of the dominant frequency band
  • / 2 represents the upper frequency limit of the dominant frequency band
  • t k represents the first VSP channel
  • the propagation time of the seismic wave between the first VSP lanes represents an indefinite coefficient
  • ft represents the quality factor between the i-th and the A-th VSP
  • G represents the matching error.
  • the determining module comprises:
  • a wave field separating unit configured to perform geometric diffusion compensation on the zero-bias VSP data, and a wave field separation to obtain a downlink wave field
  • a transmission coefficient determining unit configured to determine two adjacent ones according to a layer velocity of the seismic wave in the downward wave field
  • the transmission coefficient between the two VSP channels is determined by the layer velocity of the seismic wave, and then the transmission coefficient is taken as an indefinite coefficient to the objective function of the exponential method to estimate the quality factor, thereby solving the existing
  • the estimated quality is caused by the existence of the indefinite coefficient.
  • FIG. 1 is a flow chart of a method for estimating a quality factor based on zero-bias VSP data according to an embodiment of the present invention
  • FIG. 2 is a structural block diagram of an apparatus for estimating a quality factor based on a zero-bias VSP data according to an embodiment of the present invention
  • FIG. 3 is a schematic diagram of a zero offset VSP observation system according to an embodiment of the present invention
  • FIG. 4 is a schematic diagram of a first-time recording of an attenuated VSP direct wave according to an embodiment of the present invention
  • FIG. 5 is a schematic diagram of a quality factor estimation result of a noiseless data index method according to an embodiment of the present invention.
  • FIG. 6 is a schematic diagram of a quality factor estimation result of a noisy data index method according to an embodiment of the present invention.
  • FIG. 7 is a schematic diagram showing matching error of an index method corresponding to different quality factors of noise-free data according to an embodiment of the present invention.
  • FIG. 9 is a schematic diagram showing matching errors of an index method corresponding to different quality factors of noisy data according to an embodiment of the present invention.
  • FIG. 10 is a schematic diagram showing coefficients of an index method corresponding to different quality factors of noise-free data according to an embodiment of the present invention.
  • Figure 12 is a schematic diagram showing the stability analysis of the exponential method when the coefficients of the noise-free data are known in the embodiment of the present invention.
  • Figure 13 is a schematic diagram showing the stability analysis of the exponential method when the coefficients of the noisy data are known in the embodiment of the present invention.
  • FIG. 14 is a schematic diagram of a velocity model and a quality factor model according to an embodiment of the present invention.
  • Figure 15 is a diagram showing the results of the quality factor estimation of the spectral ratio method, the index method and the improved index method for noise-free data according to an embodiment of the present invention
  • 16 is a schematic diagram of a spectral factor quality factor estimation result for noisy data according to an embodiment of the present invention
  • 17 is a schematic diagram showing an index quality factor estimation result for noisy data according to an embodiment of the present invention
  • FIG. 18 is a schematic diagram showing an improved index quality factor estimation result for noisy data according to an embodiment of the present invention.
  • the inventor's research on the exponential method found that the accuracy of the quality factor estimated by the exponential method is biased because there is an unknown coefficient in the exponential method, and the existence of the unknown coefficient affects the performance of the exponential method. Therefore, the inventor found that the unknown coefficient in the exponential method can be obtained first, and then the quality factor is estimated, thereby eliminating the influence of the unknown coefficient on the stability and accuracy of the method. By using the quality factor estimation in this way, it can be further improved. Stability of quality factor estimates.
  • the method includes the following steps:
  • Step 101 Determine a transmission coefficient between two adjacent VSP channels based on a layer velocity of a seismic wave of the VSP data
  • Step 102 Determine the transmission coefficient as an indefinite coefficient in an exponential objective function, and estimate a quality factor according to the exponential objective function.
  • the transmission coefficient between the two VSP channels is determined by the layer velocity of the seismic wave, and then the transmission coefficient is determined as an indefinite coefficient in the objective function of the exponential method, and the quality factor is estimated according to the objective function of the exponential method. Therefore, in the prior art, in the process of estimating the quality factor by the exponential method, the technical problem of the stability and the accuracy of the estimated quality factor caused by the existence of the indefinite coefficient is achieved, and the stability of the estimated quality factor is improved. The technical effect of precision.
  • the reason why the transmission coefficient is used as the unknown coefficient in the exponential method is mainly because:
  • the unknown coefficient generally has two aspects: the geometric diffusion effect in the propagation process and the transmission loss at the wave impedance interface. Among them, the amplitude attenuation caused by the geometric diffusion effect can be compensated by geometric diffusion correction, so the problem of finding the coefficient is transformed into how to solve the transmission loss caused by the wave impedance interface.
  • the transmission coefficient can be substituted into the objective function of the exponential method to obtain the quality factor, that is, the indeterminate coefficient of the transmission coefficient as the objective function is substituted into the objective function;
  • the objective function obtains a quality factor corresponding to the minimum matching error; and the quality factor corresponding to the minimum matching error is used as the obtained quality factor.
  • ; represents the lower frequency limit of the dominant frequency band
  • / 2 represents the upper frequency limit of the dominant frequency band
  • t k represents the first VSP channel
  • the propagation time of the seismic wave between the first VSP lanes represents an indefinite coefficient
  • ft represents the quality factor between the i-th and the A-th VSP
  • G represents the matching error.
  • determining a transmission coefficient between two adjacent VSP channels based on a layer velocity of the seismic wave in the VSP data of the zero-bias vertical seismic profile may include:: the zero-bias VSP data A geometrical diffusion compensation, and a wave field separation are performed to obtain a downward wave field; and a transmission coefficient between adjacent two VSP channels is determined according to a layer velocity of the seismic wave in the downward wave field. That is, the data input of the method for estimating the quality factor based on the zero-bias VSP data is the downward wave field of the geometrically diffused compensated zero-bias VSP data.
  • an embodiment of the present invention further provides an apparatus for estimating a quality factor based on zero-bias VSP data, as described in the following embodiments. Since the principle of the device solving problem based on the zero-bias VSP data estimation quality factor is similar to the method based on the zero-bias VSP data estimation quality factor, the implementation of the device for estimating the quality factor based on the zero-bias VSP data can be referred to the estimation quality based on the zero-bias VSP data. The implementation of the method of the factor, the repetition will not be repeated.
  • the term "unit” or "module" may implement a combination of software and/or hardware for a predetermined function.
  • FIG. 2 is a structural block diagram of an apparatus for estimating a quality factor based on zero-bias VSP data according to an embodiment of the present invention. As shown in FIG. 2, the method includes: a determining module 201 and an estimating module 202, which are described below.
  • a determining module 201 configured to determine a transmission coefficient between adjacent two VSP channels based on a layer velocity of the seismic wave of the zero-bias vertical seismic profile VSP data;
  • the estimating module 202 is configured to determine the transmission coefficient as an indefinite coefficient in an exponential objective function, and estimate a quality factor according to the exponential method objective function.
  • the estimation module 202 includes:
  • a determining unit configured to determine a value range of the quality factor
  • a scanning unit configured to obtain a quality factor that causes a minimum matching error of the target function by means of quality factor scanning within the range of values
  • the estimating unit is configured to use the quality factor corresponding to the minimum matching error as the estimated quality factor.
  • the objective function is:
  • ; represents the lower frequency limit of the dominant frequency band
  • / 2 represents the upper frequency limit of the dominant frequency band
  • t k represents the first VSP
  • Q k represents the quality factor between the ⁇ 1 and A A VSP tracks
  • G represents the matching error.
  • the determining module 201 includes: a wave field separating unit, configured to perform geometric diffusion compensation on the zero-bias VSP data, and a wave field separation to obtain a downlink wave field; and a transmission coefficient determining unit, configured to perform, according to the downlink The layer velocity of the seismic wave in the wave field determines the transmission coefficient between two adjacent VSP channels.
  • the present invention also provides a specific embodiment to describe the above method for estimating the quality factor based on the zero-bias VSP data, however, it is worth noting that the specific embodiment is only for better illustrating the present invention, and It does not constitute an undue limitation of the invention.
  • the amplitude spectrum of the seismic wavelet during propagation can be expressed as:
  • Equation 1 Q exp f (Equation 1)
  • / represents the frequency, and respectively represents the amplitude spectrum of the seismic wave at the depth of the detection point as shown in Figure 3, and ⁇ represents the seismic wave between the two detection points.
  • the spectral ratio method is the most commonly used quality factor estimation method.
  • the method uses a linear relationship between the spectral logarithm of the seismic wavelet at two different moments and the frequency to estimate the quality factor.
  • the spectral logarithm is easily affected by the amplitude spectrum notch and noise, which causes the linear relationship with frequency to oscillate, which affects the stability of the quality factor estimation.
  • the predecessors proposed an exponential method, which uses the forward matching method to obtain the inter-layer quality factor by solving the minimum value of the equation below:
  • the formula 2 contains only one unknown quantity, so that the ft of the minimum value of the formula 2 can be obtained by the quality factor scanning.
  • the attenuated VSP record as shown in Fig. 4 can be obtained by forward modeling.
  • only the first arrival wave of the downlink direct wave is simulated in Fig. 4. Field, and the wavelet waveform changes of adjacent channels are only caused by the inter-layer quality factor.
  • the scanning range of the quality factor is set from 1 to 300, and the matching error of the formula 2 in the coefficient formula 3 is calculated separately.
  • the quality factor corresponding to the minimum matching error is 150 as shown in Fig. 7.
  • the coefficient corresponding to the quality factor at this time is 1 as shown in Fig. 8, and completely kisses the true value.
  • the quality factor corresponding to the minimum matching error is 206 as shown in Fig. 9.
  • the coefficient corresponding to the quality factor at this time ( ⁇ is 0.9948 as shown in Fig. 10, both deviate from the true value.
  • the scan range is set from 1 to 300, and the scan range of the coefficient is
  • the coefficients generally include two aspects: the geometric diffusion effect during propagation and the transmission loss at the wave impedance interface.
  • the amplitude attenuation caused by the geometric diffusion effect can be compensated by geometric diffusion correction, so the problem of finding the coefficient is transformed into how to solve the transmission loss caused by the wave impedance interface.
  • the zero-bias VSP data according to the arrival time of each direct wave and the depth of the detection point, the average velocity of the surface to the detection point can be obtained, and then converted into the layer velocity, so that the adjacent two paths can be obtained.
  • Equation 5 The above formula 5 is the objective function of the improved exponential method, and the ⁇ value corresponding to the minimum matching error obtained by the ⁇ -scan is taken as the inter-layer quality factor sought.
  • the source is located at the surface, the receiving point depth is 100 ⁇ 800m, and the receiving point interval is 10m, using the propagation matrix method including absorption attenuation effect.
  • the zero-bias VSP full-wave field, the down-going wave field and the up-going wave field are obtained respectively.
  • the spectral ratio method, the index method, and the estimation methods in the present embodiment all agree well with the true values.
  • the spectral ratio method is shown in Figure 16. The deviation is very large.
  • the estimation result is difficult to accept.
  • the index method is shown in Figure 17, and the estimation result is relative to the spectrum.
  • the comparison method has been greatly improved.
  • the estimation method in this example is as shown in Fig. 18, and the estimation result of the quality factor is further improved.
  • the thick line represents the true value
  • the thin line represents the Estimated results.
  • Estimating the quality factor by the above method can effectively improve the stability and accuracy of the quality factor estimation.
  • a storage medium is also provided, the software being stored, including but not limited to: an optical disk, a floppy disk, a hard disk, a rewritable memory, and the like.
  • the embodiment of the present invention achieves the following technical effects: determining the transmission coefficient between two VSP channels by the layer velocity of the seismic wave, and then determining the transmission coefficient as the index function objective function Indefinite coefficient, the quality factor is estimated according to the objective function of the exponential method, thereby solving the prior art in the process of estimating the quality factor by using the exponential method, and the stability and accuracy of the estimated quality factor are low due to the existence of the indefinite coefficient.
  • the technical problem has reached the technical effect of improving the stability and accuracy of the estimated quality factor.
  • modules or steps of the embodiments of the present invention can be implemented by a general computing device, which can be concentrated on a single computing device or distributed in multiple computing devices. Alternatively, they may be implemented by program code executable by the computing device, such that they may be stored in the storage device by the computing device and, in some cases, may be different from The steps shown or described are performed sequentially, or they are separately fabricated into individual integrated circuit modules, or they are Multiple modules or steps in the fabrication are implemented as a single integrated circuit module. Thus, embodiments of the invention are not limited to any specific combination of hardware and software.

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Remote Sensing (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

提供了一种基于零偏垂直地震剖面数据估计品质因子的方法和装置，其中，该方法包括以下步骤：基于零偏垂直地震剖面VSP数据中的地震波的层速度，确定相邻两个VSP道之间的透射系数（101）；将所述透射系数确定为指数法目标函数中的不定系数，根据所述指数法目标函数估计品质因子（102）。本方法和装置解决了现有技术中在采用指数法估计品质因子的过程中，因为不定系数的存在而导致的估计品质因子的稳定性和精度较低的技术问题，达到了提高估计品质因子的稳定性和精度的技术效果。

Description

基于零偏垂直地震剖面数据估计品质因子的方法和装置 技术领域
本发明涉及物理勘探技术领域， 特别涉及一种基于零偏垂直地震剖面 （Vertical Seismic Profiling, VSP) 数据估计品质因子的方法和装置。 背景技术
地球介质是非完全弹性的， 地震波在地球介质中传播时， 介质粘滞性引起的吸收作 用会造成地震波的能量衰减和速度频散， 这种介质所固有的衰减特性通常用品质因子 ρ 来描述。 由于品质因子的存在， 使得地震波高频能量的衰减严重于低频能量的衰减， 高 频成分的传播速度快于低频成分的传播速度， 同时也会造成地震剖面深层能量变弱、 频 带变窄， 从而导致资料的分辨率降低， 增加了地震资料精细解释的困难。
针对上述问题， 反 ρ滤波是进行地震资料吸收衰减补偿的有效手段， 能够对地震波 传播过程中的能量衰减和速度频散进行补偿和校正， 从而提高资料的分辨率。
精确的品质因子估计是进行反 ρ滤波的前提条件， 相比于地面地震数据， 垂直地震 剖面 （Vertical Seismic Profiling, VSP) 数据由于受到的干扰较少， 被广泛应用于品质因 子估计。 前人用 VSP资料对振幅衰减法、 上升时间法、 子波模拟法、 解析信号法等时间 域方法和匹配法、 谱模拟法、 谱比法等频率域方法进行了对比研究后发现， 没有哪一种 方法适用于任何情况， 每种方法所产生的效果的好坏依赖于资料的质量。 目前有一种质 心频移法， 主要是利用地震波传播过程中质心频率的变化来获得品质因子。 在上述这些 方法中， 谱比法是一种进行品质因子估计的常用方法， 通过利用谱比对数与频率间的线 性关系来估计品质因子， 然而其斜率拟合容易受到谱比对数误差的影响， 导致品质因子 估计的稳定性受到影响， 指数法有效避免了谱比对数存在的弊端， 指数法是采用正演匹 配的方式来估计品质因子， 能够有效改善品质因子估计的稳定性。 然而， 在实际的操作 中发现指数法所估计的品质因子也经常会出现较大的误差， 品质因子估计的精度难以保 证。 发明内容
本发明实施例提供了一种基于零偏垂直地震剖面数据估计品质因子的方法， 以达到 提高品质因子的估计稳定性和估计精度的目的， 该方法包括： 基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻两个 VSP道之间 的透射系数；
将所述透射系数确定为指数法目标函数中的不定系数， 根据所述指数法目标函数估 计品质因子。
在一个实施例中， 基于零偏 VSP资料的地震波的层速度， 确定相邻两个 VSP道之 间的透射系数， 包括：
按照以下公式确定所述透射系数：
pk = \+ Vk ~Vk~l 其中， A表示透射系数， ^表示第 i与第 A个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
在一个实施例中， 将所述透射系数确定为指数法目标函数中的不定系数， 根据所述 指数法目标函数估计品质因子， 包括：
确定品质因子的取值范围；
在所述取值范围内通过品质因子扫描的方式， 获取使得目标函数取得最小匹配误差 的品质因子；
将取得最小匹配误差所对应的品质因子作为估计的品质因子。
在一个实施例中， 所述目标函数为：
Figure imgf000003_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道的地震波的振幅谱， 表示第 个 VSP道的地震波的振幅谱， tk表示 第 A个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， 即所述 反射， a表示第 1与第 A个 VSP道之间的品质因子， G表示匹配误差。
在一个实施例中， 基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻 两个 VSP道之间的透射系数， 包括：
对所述零偏 VSP数据进行几何扩散补偿、 和波场分离得到下行波场；
根据所述下行波场中的地震波的层速度， 确定相邻两个 VSP道之间的透射系数。 本发明实施例提供了一种基于零偏垂直地震剖面数据估计品质因子的装置， 以达到 提高品质因子的估计稳定性和估计精度的目的， 该装置包括： 确定模块， 用于基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻两 个 VSP道之间的透射系数；
估计模块， 用于将所述透射系数确定为指数法目标函数中的不定系数， 根据所述指 数法目标函数估计品质因子。
在一个实施例中， 所述确定模块具体用于按照以下公式确定所述透射系数： pk = \+ Vk ~Vk~l 其中， A表示透射系数， ^表示第 i与第 A个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
在一个实施例中， 所述估计模块包括：
确定单元， 用于确定品质因子的取值范围；
扫描单元， 用于在所述取值范围内通过品质因子扫描的方式， 获取使得目标函数取 得最小匹配误差的品质因子；
估计单元， 用于将取得最小匹配误差所对应的品质因子作为估计的品质因子。 在一个实施例中， 所述目标函数为：
Figure imgf000004_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道的地震波的振幅谱， 表示第 个 VSP道的地震波的振幅谱， tk表示 第 个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， ft表示 第 i与第 A个 VSP道之间的品质因子， G表示匹配误差。
在一个实施例中， 所述确定模块包括：
波场分离单元， 用于对所述零偏 VSP数据进行几何扩散补偿、 和波场分离得到下行 波场；
透射系数确定单元， 用于根据所述下行波场中的地震波的层速度， 确定相邻两个
VSP道之间的透射系数。
在本发明实施例中， 通过地震波的层速度来确定两个 VSP道之间的透射系数， 然后 将该透射系数作为不定系数带入到指数法的目标函数中估计品质因子， 从而解决了现有 技术中在采用指数法估计品质因子的过程中， 因为不定系数的存在而导致的估计品质因 子的稳定性和精度较低的技术问题， 达到了提高估计品质因子的稳定性和精度的技术效 果。 附图说明
此处所说明的附图用来提供对本发明的进一步理解， 构成本申请的一部分， 并不构 成对本发明的限定。 在附图中：
图 1是本发明实施例的基于零偏 VSP数据估计品质因子的方法流程图；
图 2是本发明实施例的基于零偏 VSP数据估计品质因子的装置结构框图； 图 3是本发明实施例的零偏移距 VSP观测系统示意图；
图 4是本发明实施例的衰减 VSP直达波初至记录示意图；
图 5是本发明实施例的无噪数据指数法品质因子估计结果示意图；
图 6是本发明实施例的含噪数据指数法品质因子估计结果示意图；
图 7是本发明实施例的对于无噪资料的不同的品质因子对应的指数法的匹配误差示 意图；
图 8是本发明实施例的对于无噪资料的不同品质因子对应的指数法的系数 (^示意 图；
图 9是本发明实施例的对于含噪资料的不同的品质因子对应的指数法的匹配误差示 意图；
图 10是本发明实施例的对于无噪资料的不同品质因子对应的指数法的系数 示意 图；
图 11是本发明实施例的在设定系数 = 1的情况下， 指数法品质因子估计结果示意 图；
图 12是本发明实施例的对于无噪资料的系数 已知时的指数法的稳定性分析示意 图；
图 13是本发明实施例的对于含噪资料的系数 已知时的指数法的稳定性分析示意 图；
图 14是本发明实施例的速度模型与品质因子模型示意图；
图 15是本发明实施例的对于无噪资料的谱比法、 指数法和改进指数法的品质因子估 计结果示意图；
图 16是本发明实施例的对于含噪资料的谱比法品质因子估计结果示意图； 图 17是本发明实施例的对于含噪资料的指数法品质因子估计结果示意图； 图 18是本发明实施例的对于含噪资料的改进指数法品质因子估计结果示意图。 具体实施方式
发明人对指数法进行研究后发现， 致使指数法估计的品质因子的精度出现偏差的原 因在于， 在指数法中存在未知系数， 未知系数的存在影响了指数法性能的发挥。 因此发 明人研究后发现， 可以首先求取指数法中的未知系数， 然后再估计品质因子， 从而消除 了未知系数对方法稳定性和精度的影响， 通过这种方式进行品质因子估计， 可以进一步 改善品质因子估计的稳定性。
在本发明实施例中， 提出了一种基于零偏 VSP数据估计品质因子的方法， 如图 1所 示， 包括以下步骤：
步骤 101 : 基于 VSP数据的地震波的层速度， 确定相邻两个 VSP道之间的透射系 数；
步骤 102: 将所述透射系数确定为指数法目标函数中的不定系数， 根据所述指数法 目标函数估计品质因子。
在上述实施例中， 通过地震波的层速度来确定两个 VSP道之间的透射系数， 然后将 该透射系数确定为指数法目标函数中的不定系数， 根据所述指数法目标函数估计品质因 子， 从而解决了现有技术中在采用指数法估计品质因子的过程中， 因为不定系数的存在 而导致的估计品质因子的稳定性和精度较低的技术问题， 达到了提高估计品质因子的稳 定性和精度的技术效果。
具体的， 上述步骤 101中确定透射系数， 可以按照以下公式实现： pk = \+ Vk ~Vk~l 其中， A表示透射系数， ^表示第 i与第 A个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
之所以以透射系数作为指数法中的未知系数， 主要是因为： 未知系数 一般包含两 个方面： 传播过程中的几何扩散效应 和波阻抗界面处的透射损失 。 其中， 几何扩散 效应造成的振幅衰减可以通过几何扩散校正来获得补偿， 因此求取系数 的问题就转化 为如何求解波阻抗界面造成的透射损失 。 在求取得到透射系数以后， 就可以将该透射系数代入到指数法的目标函数中求取品 质因子， 即将所述透射系数作为目标函数的不定系数代入目标函数； 通过品质因子扫描 的方式， 获取所述目标函数取得最小匹配误差时对应的品质因子； 将该最小匹配误差对 应的品质因子作为求取的品质因子。
上述目标函数为：
Figure imgf000007_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道的地震波的振幅谱， 表示第 个 VSP道的地震波的振幅谱， tk表示 第 个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， ft表示 第 i与第 A个 VSP道之间的品质因子， G表示匹配误差。
在上述各个实施例中， 上述步骤 101中， 基于零偏垂直地震剖面 VSP数据中的地震 波的层速度， 确定相邻两个 VSP道之间的透射系数， 可以包括： 对所述零偏 VSP数据 进行几何扩散补偿、 和波场分离得到下行波场； 根据所述下行波场中的地震波的层速 度， 确定相邻两个 VSP道之间的透射系数。 即， 该基于零偏 VSP数据估计品质因子的 方法的数据输入是经过几何扩散补偿的零偏 VSP数据的下行波场。
基于同一发明构思， 本发明实施例中还提供了一种基于零偏 VSP数据估计品质因子 的装置， 如下面的实施例所述。 由于基于零偏 VSP数据估计品质因子的装置解决问题的 原理与基于零偏 VSP数据估计品质因子的方法相似， 因此基于零偏 VSP数据估计品质 因子的装置的实施可以参见基于零偏 VSP数据估计品质因子的方法的实施， 重复之处不 再赘述。 以下所使用的， 术语 "单元"或者 "模块"可以实现预定功能的软件和 /或硬件 的组合。 尽管以下实施例所描述的装置较佳地以软件来实现， 但是硬件， 或者软件和硬 件的组合的实现也是可能并被构想的。 图 2是本发明实施例的基于零偏 VSP数据估计品 质因子的装置的一种结构框图， 如图 2所示， 包括： 确定模块 201和估计模块 202， 下 面对该结构进行说明。
确定模块 201， 用于基于零偏垂直地震剖面 VSP数据的地震波的层速度， 确定相邻 两个 VSP道之间的透射系数；
估计模块 202， 用于将所述透射系数确定为指数法目标函数中的不定系数， 根据所 述指数法目标函数估计品质因子。
在一个实施例中， 确定模块 201具体用于按照以下公式确定所述透射系数： P, = \ + 其中， A表示透射系数， 表示第 i与第 个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
在一个实施例中， 估计模块 202包括：
确定单元， 用于确定品质因子的取值范围；
扫描单元， 用于在所述取值范围内通过品质因子扫描的方式， 获取使得目标函数取 得最小匹配误差的品质因子；
估计单元， 用于将取得最小匹配误差所对应的品质因子作为估计的品质因子。 在一个实施例中， 所述目标函数为：
G(Qk, Ck) = f
Figure imgf000008_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道处的地震波的振幅谱， 表示第 个 VSP道处的地震波的振幅谱， tk 表示第 个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， Qk 表示第 ^ 1与第 A个 VSP道之间的品质因子， G表示匹配误差。
在一个实施例中， 确定模块 201包括： 波场分离单元， 用于对所述零偏 VSP数据进 行几何扩散补偿、 和波场分离得到下行波场； 透射系数确定单元， 用于根据所述下行波 场中的地震波的层速度， 确定相邻两个 VSP道之间的透射系数。
本发明还提供了一个具体的实施例来对上述基于零偏 VSP数据估计品质因子的方法 进行详细说明， 然而， 值得注意的是， 该具体的实施例仅是为了更好地说明本发明， 并 不构成对本发明的不当限定。
考虑地层吸收衰减效应， 传播过程中地震子波的振幅谱可以表示为：
-πί,
S,(/) = Q exp f (公式 1 ) 其中， /表示频率， 和 分别表示如图 3所示的检波点深度为 和 — , 处的地震波的振幅谱， ^表示两检波点间的地震波的传播时间， 系数 (^表示与频率无关 的量， ft表示层间品质因子， 通常假设其不依赖于频率， 值得注意的是， 此处 A在采集 数据的过程中表示检波点号， 在计算品质因子的过程中表示 VSP道的道号， 两者本质是 相同的。 基于上述公式 1， 谱比法是最常用的品质因子估计方法， 该方法是利用两个不同时 刻的地震子波的谱比对数与频率之间的线性关系来估计品质因子。 然而， 谱比对数容易 受到振幅谱陷频和噪声的影响， 导致其与频率的线性关系产生震荡， 影响了品质因子估 计的稳定性。 为此， 前人提出了一种指数法， 该方法是采用正演匹配的方式， 通过求解 方程下述方程的极小值来获得层间品质因子：
G(Qk, Ck) = f (公式 2)
Figure imgf000009_0001
其中， ;和 /2分别表示频率下限和频率上限， 为了消除上述公式 2中的未知 可以采用求导的方式， 令 dG / dCt = 0， 从而得到：
(公式 3 )
Figure imgf000009_0002
将公式 3代入公式 2， 则公式 2中仅含有一个未知量， 从而可以利用品质因子扫描 的方式获得使得公式 2取得极小值的 ft。
为了验证品质因子估计方法的性能， 可以采用正演模拟的方式得到如图 4所示的衰 减 VSP记录， 为了尽可能地消除其它因素的影响， 图 4中仅模拟了下行直达波的初至波 场， 并且相邻道的子波波形变化仅是由层间品质因子引起的。
针对如图 4所示的正演模拟记录， 引入了微弱的随机噪声 （其中， 噪声能量占有效 信号能量的 0.3%) ， 分别进行无噪资料和含噪资料的指数法品质因子估计， 估计的结果 如图 5和图 6所示， 从估计结果中可以看出， 当资料不含有噪声时， 品质因子估计结果 与理论值吻合较好 （如图 5所示） ， 含有噪声时， 品质因子估计结果与理论真实值存在 偏差 （如图 6所示） ， 其中， 在图 5和图 6中实线表示理论值， 虚线表示品质因为的估 计结果。
从图 6中可以看出， 图 4中的第 9道与第 10道之间的含噪资料的品质因子估计结果 为 206， 相对于真实值 150， 估计结果的绝对误差为 56， 相对误差为 37.3%。
为了分析指数法的误差产生的原因， 设定品质因子的扫描范围为 1~300， 分别计算 系数公式 3中的 和公式 2的匹配误差。 在无噪情况下， 最小匹配误差对应的品质因子 如图 7所示为 150， 此时的品质因子对应的系数 如图 8所示为 1， 与真实值完全吻 合； 而在含噪情况下， 最小匹配误差对应的品质因子如图 9所示为 206， 此时的品质因 子对应的系数 (^如图 10所示为 0.9948， 均偏离了真实值。
为了进一步观察其稳定性， 设定 的扫描范围为 1~300， 系数 的扫描范围为
0.8-1.35 , 分别得到了无噪资料和含噪资料的匹配误差分布图， 分析得到在真实值附近 (即， = 150， Q= l ) ， 匹配误差位于谷底， 并且趋势非常平缓， 这就造成了如果原 始资料中存在一点噪声， 最小匹配误差 （0.05595 ) 就会偏离真实的最小匹配误差 ( 0.0007557) ， 从而引起品质因子估计结果 （ft = 206) 与真实值偏离较远。
在含噪的情况下， C变成了 0.9948， 而针对图 4所示的正演记录， 实际上的理论值 应该为 1， 因此， 上述公式 3中的 的未知性是造成该方法不稳定的一个原因， 对于含 噪数据， 如果设定 为 1， 重新采用指数法计算品质因子后的结果如图 11所示， 可以发 现， 品质因子估计结果与理论值之间的吻合度获得了极大改善。
同样针对图 4中的第 9道与第 10道之间的资料， 对系数 为已知理论值时的指数 法进行了稳定性分析。 针对无噪资料和含噪资料， 设定品质因子扫描范围为 1~300， 计 算相应的匹配误差如图 12和图 13所示， 可以看到， 与图 7和图 9相比， 极小值明显位 于曲线的波谷处， 在噪声影响下， 在图 9中， 品质因子估计结果从 150变成了 155， 稳 定性获得了明显改善。
因此基于上述分析， 如果能够获得 ， 然后利用公式 2进行品质因子估计的精度就 能获得进一步改善， 下面来分析如何求取系数 。
对于零偏 VSP资料的相邻道直达波初至， 其系数 一般包含两个方面： 传播过程 中的几何扩散效应 和波阻抗界面处的透射损失 。 其中， 几何扩散效应造成的振幅衰 减可以通过几何扩散校正来获得补偿， 因此求取系数 的问题就转化为如何求解波阻抗 界面造成的透射损失 。 对于零偏 VSP资料， 根据各道直达波的到时和检波点深度， 可 以求得地表到所在检波点的平均速度， 进而将其转化为层速度， 这样就可以获得相邻两 道之间的透射系数， 即 A = l + H (公式 4) 将 带入上述公式 2中， 可以得到：
(公式 5 )
Figure imgf000010_0001
上述公式 5即为改进的指数法的目标函数， 通过 ρ扫描的方式获得的最小匹配误差 对应的 ρ值即作为所求的层间品质因子。
为了检验方法的可行性， 给定如图 14所示的速度模型与品质因子模型， 震源位于地 表， 接收点深度为 100~800m， 接收点间隔为 10m， 采用包含吸收衰减效应的传播矩阵 法进行正演数值模拟， 分别得到了零偏 VSP全波场、 下行波场和上行波场。
对于下行波场， 引入随机噪声 （其中， 噪声的能量占有效信号能量的 0.3%) ， 分别 计算无噪资料和含噪资料的品质因子。 从品质因子的估计结果发现， 对于无噪资料， 如 图 15所示， 谱比法、 指数法和本实施例中的估算方法估计结果均与真实值都吻合较好。 对于含噪资料， 当品质因子较小 （例如为： 20和 40) 时， 三种方法的估计结果均与真实 值都比较吻合， 然而当品质因子较大 （例如为： 80和 120) 时， 三种方法均存在一定程 度的估计偏差， 谱比法如图 16所示， 偏差非常大， 尤其对于品质因子为 120时， 估计结 果难以接受； 指数法如图 17所示， 估计结果相对于谱比法有了较大改善， 本例中的估计 方法如图 18所示， 又进一步改善了品质因子的估计结果， 在上述图 15至图 18中， 粗的 线代表真实值， 细的线代表估计结果。
通过上述方式估计品质因子， 可以有效提高品质因子估计的稳定性和精度。
在另外一个实施例中， 还提供了一种软件， 该软件用于执行上述实施例及优选实施 方式中描述的技术方案。
在另外一个实施例中， 还提供了一种存储介质， 该存储介质中存储有上述软件， 该 存储介质包括但不限于： 光盘、 软盘、 硬盘、 可擦写存储器等。
从以上的描述中， 可以看出， 本发明实施例实现了如下技术效果： 通过地震波的层 速度来确定两个 VSP道之间的透射系数， 然后将该透射系数确定为指数法目标函数中的 不定系数， 根据所述指数法目标函数估计品质因子， 从而解决了现有技术中在采用指数 法估计品质因子的过程中， 因为不定系数的存在而导致的估计品质因子的稳定性和精度 较低的技术问题， 达到了提高估计品质因子的稳定性和精度的技术效果。
显然， 本领域的技术人员应该明白， 上述的本发明实施例的各模块或各步骤可以用 通用的计算装置来实现， 它们可以集中在单个的计算装置上， 或者分布在多个计算装置 所组成的网络上， 可选地， 它们可以用计算装置可执行的程序代码来实现， 从而， 可以 将它们存储在存储装置中由计算装置来执行， 并且在某些情况下， 可以以不同于此处的 顺序执行所示出或描述的步骤， 或者将它们分别制作成各个集成电路模块， 或者将它们 中的多个模块或步骤制作成单个集成电路模块来实现。 这样， 本发明实施例不限制于任 何特定的硬件和软件结合。
以上所述仅为本发明的优选实施例而已， 并不用于限制本发明， 对于本领域的技术 人员来说， 本发明实施例可以有各种更改和变化。 凡在本发明的精神和原则之内， 所作 的任何修改、 等同替换、 改进等， 均应包含在本发明的保护范围之内。

Claims

权利要求书
1、 一种基于零偏垂直地震剖面数据估计品质因子的方法， 其特征在于， 包括： 基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻两个 VSP道之间 的透射系数；
将所述透射系数确定为指数法目标函数中的不定系数， 根据所述指数法目标函数估 计品质因子。
2、 如权利要求 1所述的方法， 其特征在于， 基于零偏 VSP资料中的地震波的层速 度， 确定相邻两个 VSP道之间的透射系数， 包括：
按照以下公式确定所述透射系数： pk = \+ Vk ~Vk~l 其中， A表示透射系数， ^表示第 i与第 A个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
3、 如权利要求 1所述的方法， 其特征在于， 将所述透射系数确定为指数法目标函数 中的不定系数， 根据所述指数法目标函数估计品质因子， 包括：
确定品质因子的取值范围；
在所述取值范围内通过品质因子扫描的方式， 获取使得目标函数取得最小匹配误差 的品质因子；
将取得最小匹配误差所对应的品质因子作为估计的品质因子。
4、 如权利要求 1 述的方法， 其特征在于， 所述目标函数为：
Figure imgf000013_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道的地震波的振幅谱， 表示第 个 VSP道的地震波的振幅谱， tk表示 第 个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， ft表示 第 i与第 A个 VSP道之间的品质因子， G表示匹配误差。
5、 如权利要求 1至 4中任一项所述的方法， 其特征在于， 基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻两个 VSP道之间的透射系数， 包括：
对所述零偏 VSP数据进行几何扩散补偿、 和波场分离得到下行波场； 根据所述下行波场中的地震波的层速度， 确定相邻两个 VSP道之间的透射系数。
6、 一种基于零偏垂直地震剖面数据估计品质因子的装置， 其特征在于， 包括： 确定模块， 用于基于零偏垂直地震剖面 VSP数据中的地震波的层速度， 确定相邻两 个 VSP道之间的透射系数；
估计模块， 用于将所述透射系数确定为指数法目标函数中的不定系数， 根据所述指 数法目标函数估计品质因子。
7、 如权利要求 6所述的装置， 其特征在于， 所述确定模块具体用于按照以下公式确 定所述透射系数：
pk = \+ Vk ~Vk~l 其中， A表示透射系数， ^表示第 i与第 A个 VSP道之间的层速度， 表示第 k-2与第 个 VSP道之间的层速度。
8、 如权利要求 6所述的装置， 其特征在于， 所述估计模块包括：
确定单元， 用于确定品质因子的取值范围；
扫描单元， 用于在所述取值范围内通过品质因子扫描的方式， 获取使得目标函数取 得最小匹配误差的品质因子；
估计单元， 用于将取得最小匹配误差所对应的品质因子作为估计的品质因子。
9、 6所述的装置， 其特征在于， 所述目标函数为：
Figure imgf000014_0001
其中， ;表示优势频带的频率下限， /2表示优势频带的频率上限， 表示第 A 个 VSP道的地震波的振幅谱， 表示第 个 VSP道的地震波的振幅谱， tk表示 第 A个 VSP道和第 个 VSP道之间的地震波的传播时间， 表示不定系数， ft表示 第 i与第 A个 VSP道之间的品质因子， G表示匹配误差。
10、 如权利要求 6至 9中任一项所述的装置， 其特征在于， 所述确定模块包括： 波场分离单元， 用于对所述零偏 VSP数据进行几何扩散补偿、 和波场分离得到下行 波场；
透射系数确定单元， 用于根据所述下行波场中的地震波的层速度， 确定相邻两个 VSP道之间的透射系数。
PCT/CN2013/090914 2013-12-30 2013-12-30 基于零偏垂直地震剖面数据估计品质因子的方法和装置 WO2015100544A1 (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
PCT/CN2013/090914 WO2015100544A1 (zh) 2013-12-30 2013-12-30 基于零偏垂直地震剖面数据估计品质因子的方法和装置
US14/919,053 US20160091623A1 (en) 2013-12-30 2015-10-21 Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
PCT/CN2013/090914 WO2015100544A1 (zh) 2013-12-30 2013-12-30 基于零偏垂直地震剖面数据估计品质因子的方法和装置

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US14/919,053 Continuation US20160091623A1 (en) 2013-12-30 2015-10-21 Method and Device for Estimating Quality Factor Based on Zero Offset Vertical Seismic Profile Data

Publications (1)

Publication Number Publication Date
WO2015100544A1 true WO2015100544A1 (zh) 2015-07-09

Family

ID=53492904

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2013/090914 WO2015100544A1 (zh) 2013-12-30 2013-12-30 基于零偏垂直地震剖面数据估计品质因子的方法和装置

Country Status (2)

Country Link
US (1) US20160091623A1 (zh)
WO (1) WO2015100544A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017058440A1 (en) * 2015-10-02 2017-04-06 Exxonmobil Upstream Research Company Q-compensated full wavefield inversion
CN110515127A (zh) * 2019-09-26 2019-11-29 中国石油大学(北京) 一种地震品质因子确定方法、装置、设备、介质
CN110673211A (zh) * 2019-10-13 2020-01-10 东北石油大学 一种基于测井与地震数据的品质因子建模方法
CN111596355A (zh) * 2020-06-02 2020-08-28 中油奥博(成都)科技有限公司 一种零偏移距vsp时频分析划分地层及层速度确定方法
CN111781643A (zh) * 2019-04-03 2020-10-16 中国石油天然气股份有限公司 地震层速度确定方法及装置

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107065013B (zh) * 2016-12-27 2019-01-25 中国石油大学(华东) 一种地震尺度下的层速度确定方法及装置
CN109143340B (zh) * 2018-08-20 2020-03-10 中国海洋石油集团有限公司 一种基于常q模型的粘弹介质地震波模拟方法及系统
CN112255682B (zh) * 2019-07-22 2023-10-10 中国石油化工集团有限公司 基于vsp资料的q体建模技术方法
US11391855B2 (en) 2020-03-13 2022-07-19 Saudi Arabian Oil Company Developing a three-dimensional quality factor model of a subterranean formation based on vertical seismic profiles
US11703607B2 (en) 2020-06-15 2023-07-18 Saudi Arabian Oil Company Determining a seismic quality factor for subsurface formations from a seismic source to a first VSP downhole receiver
US11859472B2 (en) 2021-03-22 2024-01-02 Saudi Arabian Oil Company Apparatus and method for milling openings in an uncemented blank pipe
US11573346B2 (en) 2021-04-15 2023-02-07 Saudi Arabian Oil Company Determining a seismic quality factor for subsurface formations for marine vertical seismic profiles
US11788377B2 (en) 2021-11-08 2023-10-17 Saudi Arabian Oil Company Downhole inflow control
US12049807B2 (en) 2021-12-02 2024-07-30 Saudi Arabian Oil Company Removing wellbore water
US12024985B2 (en) 2022-03-24 2024-07-02 Saudi Arabian Oil Company Selective inflow control device, system, and method

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120106294A1 (en) * 2010-10-28 2012-05-03 Baker Hughes Incorporated Optimization Approach to Q-Factor Estimation From VSP Data
CN103163554A (zh) * 2013-02-04 2013-06-19 西安交通大学 利用零偏vsp资料估计速度和q值的自适应波形的反演方法
CN103744114A (zh) * 2013-12-30 2014-04-23 中国石油天然气股份有限公司 基于零偏垂直地震剖面数据估计品质因子的方法和装置

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120106294A1 (en) * 2010-10-28 2012-05-03 Baker Hughes Incorporated Optimization Approach to Q-Factor Estimation From VSP Data
CN103163554A (zh) * 2013-02-04 2013-06-19 西安交通大学 利用零偏vsp资料估计速度和q值的自适应波形的反演方法
CN103744114A (zh) * 2013-12-30 2014-04-23 中国石油天然气股份有限公司 基于零偏垂直地震剖面数据估计品质因子的方法和装置

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
WANG, YANCHUN ET AL.: "computing for the viscous elastic medium zero-offset vsp data", GEOPHYSICAL PROSPECTING FOR PETROLEUM, vol. 29, no. 2, 30 June 1990 (1990-06-30) *
WU, YINTING ET AL.: "Research and modification of VSP attenuation characteristics based on Ganley's forward theory", CHINESE JOURNAL OF GEOPHYSICS, vol. 54, no. 9, 30 September 2011 (2011-09-30) *
WU, YINTING: "The Forward and Inverse of Q Value Using zero-offset VSP data", A DISSERTATION SUBMITTED FOR THE DEGREE OF MASTER OF CHANG'AN UNIVERSITY, 15 August 2009 (2009-08-15) *
YUAN, XIAOHUA ET AL.: "Forward Modeling Experiment on Attenuation Characteristics of Zero-offset VSP", INVESTIGATION SCIENCE AND TECHNOLOGY, no. 2, 20 April 2013 (2013-04-20) *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017058440A1 (en) * 2015-10-02 2017-04-06 Exxonmobil Upstream Research Company Q-compensated full wavefield inversion
CN111781643A (zh) * 2019-04-03 2020-10-16 中国石油天然气股份有限公司 地震层速度确定方法及装置
CN111781643B (zh) * 2019-04-03 2022-07-05 中国石油天然气股份有限公司 地震层速度确定方法及装置
CN110515127A (zh) * 2019-09-26 2019-11-29 中国石油大学(北京) 一种地震品质因子确定方法、装置、设备、介质
CN110673211A (zh) * 2019-10-13 2020-01-10 东北石油大学 一种基于测井与地震数据的品质因子建模方法
CN110673211B (zh) * 2019-10-13 2021-06-04 东北石油大学 一种基于测井与地震数据的品质因子建模方法
CN111596355A (zh) * 2020-06-02 2020-08-28 中油奥博(成都)科技有限公司 一种零偏移距vsp时频分析划分地层及层速度确定方法
CN111596355B (zh) * 2020-06-02 2022-04-01 中国石油集团东方地球物理勘探有限责任公司 一种零偏移距vsp时频分析划分地层及层速度确定方法

Also Published As

Publication number Publication date
US20160091623A1 (en) 2016-03-31

Similar Documents

Publication Publication Date Title
WO2015100544A1 (zh) 基于零偏垂直地震剖面数据估计品质因子的方法和装置
US10705238B2 (en) Method and apparatus for processing seismic data
Agudo et al. Acoustic full-waveform inversion in an elastic world
US8670288B2 (en) Velocity model for well time-depth conversion
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN109100786B (zh) 深度域品质因子的确定方法和装置
CN105116448B (zh) 一种转换波方位各向异性校正方法及装置
CN104820218B (zh) 一种基于频域自相关的浅海海底单参数反演方法
CN103460074A (zh) 全波场反演中小波估计和多次波预测的方法
Mordret et al. Helmholtz tomography of ambient noise surface wave data to estimate Scholte wave phase velocity at Valhall Life of the Field
CN103744114B (zh) 基于零偏垂直地震剖面数据估计品质因子的方法和装置
EA032186B1 (ru) Сейсмическая адаптивная фокусировка
CN101644782A (zh) 基于极化滤波的多波多分量地震资料的去噪方法
CN108873063B (zh) 一种微地震矩张量反演的方法及装置
CN107065013A (zh) 一种地震尺度下的层速度确定方法及装置
CN103984016B (zh) 转换波各向异性振幅随入射角变化道集抽取方法
CN109143345B (zh) 基于模拟退火的品质因子q非线性反演方法及系统
CN111337977B (zh) 品质因子的确定方法及装置
KR101693310B1 (ko) 평면파를 이용한 이방성 매질에서의 전파형 역산방법
CN112904415B (zh) 基于重构源谱一致性的q层析反演方法、系统及处理设备
CN110007349B (zh) 一种弹性参数反演方法
CN105785435A (zh) 地震物理模型实验系统及其模拟数据采集方法和装置
CN107037483B (zh) 一种基于时间慢度域的偶极反射横波三维成像方法
CN112415601B (zh) 表层品质因子q值的确定方法及装置
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 13900832

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

122 Ep: pct application non-entry in european phase

Ref document number: 13900832

Country of ref document: EP

Kind code of ref document: A1