Nonstationary Maximum Likelihood Method for Sonic Logging Dispersion

Resolve Bottlenecks,
Find Innovative Solutions
Generate Solutions

Solution Overview

Problem

Current methods for estimating dispersion spectra in full waveform sonic (FWS) logging, such as the stationary maximum likelihood method, face challenges with over-estimation and inaccurate representations due to the need for partial linear moveout correction and pseudo-ensemble averaging, leading to high computational costs and reduced resolution.

Innovation Solution

A nonstationary maximum likelihood method is introduced, which uses nonstationary predictive error filtering within the Burg algorithm to estimate local matrices, eliminating the need for partial linear moveout correction and pseudo-ensemble averaging, thereby reducing computational costs and improving resolution.

Engineering Contradictions & Design Principles

VSEngineering Contradiction Analysis

1Reliability

If partial linear moveout correction and pseudo-ensemble averaging are used in stationary maximum likelihood method, then computational stability is improved, but computational cost increases and resolution decreases

Engineering Contradiction:
Improvecomputational stabilityVSAvoidcomputational cost
Core Design Contradiction:
ReliabilityVSProductivity

Solution Approach 1:

The patent transitions from stationary to nonstationary maximum likelihood method, changing the fundamental assumption about signal stationarity. This parameter change allows the method to adapt to local variations in the data without requiring pseudo-ensemble averaging, thereby reducing computational cost while maintaining stability through the nonstationary framework's ability to handle local covariance structures.

Inventive Principle:
Principle #35Parameter changes

Solution Approach 2:

The nonstationary method segments the signal analysis into local regions with their own covariance matrices, rather than applying a global stationary assumption. This segmentation allows computational stability to be achieved locally without the need for pseudo-ensemble averaging across the entire dataset, reducing overall computational cost.

Inventive Principle:
Principle #1Segmentation

2Reliability

If partial linear moveout correction and pseudo-ensemble averaging are used in stationary maximum likelihood method, then computational stability is improved, but measurement precision deteriorates

Engineering Contradiction:
Improvecomputational stabilityVSAvoiddispersion spectra estimation accuracy
Core Design Contradiction:
ReliabilityVSMeasurement precision

Solution Approach 1:

By changing from stationary to nonstationary assumptions, the method accurately captures local variations in dispersion characteristics without the blurring effect of pseudo-ensemble averaging. This improves measurement precision while the nonstationary framework's rigorous local covariance estimation maintains computational stability.

Inventive Principle:
Principle #35Parameter changes

Solution Approach 2:

The nonstationary method applies local quality by estimating separate covariance matrices for different segments of the data, allowing each local region to be analyzed with its own characteristics. This eliminates the need for pseudo-ensemble averaging and improves dispersion spectra estimation accuracy while maintaining stability through local adaptation.

Inventive Principle:
Principle #3Local quality

3Ease of operation

If conventional stationary method is used, then processing simplicity is maintained, but resolution in characterizing velocity dispersion decreases

Engineering Contradiction:
Improveprocessing simplicityVSAvoidvelocity dispersion resolution
Core Design Contradiction:
Ease of operationVSMeasurement precision

Solution Approach 1:

The transition to nonstationary method changes the processing framework to accommodate local variations, which improves velocity dispersion resolution. While the processing becomes more complex, the patent implements efficient algorithms that compute local covariance matrices and apply nonstationary maximum likelihood estimation, achieving high resolution with manageable computational effort.

Inventive Principle:
Principle #35Parameter changes

4Measurement precision

If nonstationary method is used, then resolution in characterizing velocity dispersion is improved, but device complexity increases

Engineering Contradiction:
Improvevelocity dispersion resolutionVSAvoidprocessing algorithm complexity
Core Design Contradiction:
Measurement precisionVSDevice complexity

Solution Approach 1:

The nonstationary method improves resolution by changing the stationarity parameter, but the patent manages complexity through efficient implementation strategies including localized covariance matrix computation and systematic application of the nonstationary maximum likelihood framework, making the increased complexity tractable.

Inventive Principle:
Principle #35Parameter changes

Solution Approach 2:

By segmenting the data into local regions for separate covariance estimation, the patent manages algorithmic complexity through modular processing. Each segment is handled independently with its own covariance matrix, making the overall nonstationary analysis more manageable while maintaining high resolution.

Inventive Principle:
Principle #1Segmentation

Data Source

PatentEP3938811B1Nonstationary maximum likelihood method to estimate dispersion spectra for full waveform sonic logging
Publication Date: 2023.12.27 SAUDI ARABIAN OIL CO
  • EP3938811B1 patent drawingFigure 1
  • EP3938811B1 patent drawingFigure 2
  • EP3938811B1 patent drawingFigure 3

AI summary

The present disclosure describes methods and systems for estimating dispersion spectra for full waveform sonic (FWS) logging. One computer-implemented method includes receiving FWS data, performing frequency-spatial (FX) transform on the FWS data, using a nonstationary predictive error filtering (PEF) inversion on the transformed FWS data to estimate local matrix L and matrix P, calculating an inverse covariance matrix based on the estimated local matrix L and matrix P, and obtaining a nonstationary maximum likelihood method (MLM) spectra based on the inverse covariance matrix.