Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuit...Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.展开更多
The low-wavenumber components in the gradient of full waveform inversion(FWI)play a vital role in the stability of the inversion.However,when FWI is implemented in some high frequencies and current models are not far ...The low-wavenumber components in the gradient of full waveform inversion(FWI)play a vital role in the stability of the inversion.However,when FWI is implemented in some high frequencies and current models are not far away from the real velocity model,an excessive number of low-wavenumber components in the gradient will also reduce the convergence rate and inversion accuracy.To solve this problem,this paper firstly derives a formula of scattering angle weighted gradient in FWI,then proposes a hybrid gradient.The hybrid gradient combines the conventional gradient of FWI with the scattering angle weighted gradient in each inversion frequency band based on an empirical formula derived herein.Using weighted hybrid mode,we can retain some low-wavenumber components in the initial lowfrequency inversion to ensure the stability of the inversion,and use more high-wavenumber components in the high-frequency inversion to improve the convergence rate.The results of synthetic data experiment demonstrate that compared to the conventional FWI,the FWI based on the proposed hybrid gradient can effectively reduce the low-wavenumber components in the gradient under the premise of ensuring inversion stability.It also greatly enhances the convergence rate and inversion accuracy,especially in the deep part of the model.And the field marine seismic data experiment also illustrates that the FWI based on hybrid gradient(HGFWI)has good stability and adaptability.展开更多
Full-waveform inversion(FWI) uses the full information of seismic data to obtain a quantitative estimation of subsurface physical parameters. Anisotropic FWI has the potential to recover high-resolution velocity and a...Full-waveform inversion(FWI) uses the full information of seismic data to obtain a quantitative estimation of subsurface physical parameters. Anisotropic FWI has the potential to recover high-resolution velocity and anisotropy parameter models, which are critical for imaging the long-offset and wideazimuth data. We develop an acoustic anisotropic FWI method based on a simplified pure quasi P-wave(qP-wave) equation, which can be solved efficiently and is beneficial for the subsequent inversion.Using the inverse Hessian operator to precondition the functional gradients helps to reduce the parameter tradeoff in the multi-parameter inversion. To balance the accuracy and efficiency, we extend the truncated Gauss-Newton(TGN) method into FWI of pure qP-waves in vertical transverse isotropic(VTI) media. The inversion is performed in a nested way: a linear inner loop and a nonlinear outer loop.We derive the formulation of Hessian-vector products for pure qP-waves in VTI media based on the Lagrange multiplier method and compute the model update by solving a Gauss-Newton linear system via a matrix-free conjugate gradient method. A suitable preconditioner and the Eisenstat and Walker stopping criterion for the inner iterations are used to accelerate the convergence and avoid prohibitive computational cost. We test the proposed FWI method on several synthetic data sets. Inversion results reveal that the pure acoustic VTI FWI exhibits greater accuracy than the conventional pseudoacoustic VTI FWI. Additionally, the TGN method proves effective in mitigating the parameter crosstalk and increasing the accuracy of anisotropy parameters.展开更多
The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods fo...The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.展开更多
Millimeter-wave(mmWave)radar communication has emerged as an important technique for future wireless systems.However,the interference between the radar signal and communication data is the main issue that should be co...Millimeter-wave(mmWave)radar communication has emerged as an important technique for future wireless systems.However,the interference between the radar signal and communication data is the main issue that should be considered for the joint radar communication system.In this paper,a co-sharing waveform(CSW)is proposed to achieve communication and radar sensing simultaneously.To eliminate the co-interference between the communication and sensing signal,signal splitting and processing methods for communication data demodulation and radar signal processing are given respectively.Simulation results show that the bit error rate(BER)of CSW is close to that of the pure communication waveform.Moreover,the proposed CSW can achieve better performance than the existing waveforms in terms of range and velocity estimation.展开更多
Underwater pulse waveform recognition is an important method for underwater object detection.Most existing works focus on the application of traditional pattern recognition methods,which ignore the time-and space-vary...Underwater pulse waveform recognition is an important method for underwater object detection.Most existing works focus on the application of traditional pattern recognition methods,which ignore the time-and space-varying characteristics in sound propagation channels and cannot easily extract valuable waveform features.Sound propagation channels in seawater are time-and space-varying convolutional channels.In the extraction of the waveform features of underwater acoustic signals,the effect of high-accuracy underwater acoustic signal recognition is identified by eliminating the influence of time-and space-varying convolutional channels to the greatest extent possible.We propose a hash aggregate discriminative network(HADN),which combines hash learning and deep learning to minimize the time-and space-varying effects on convolutional channels and adaptively learns effective underwater waveform features to achieve high-accuracy underwater pulse waveform recognition.In the extraction of the hash features of acoustic signals,a discrete constraint between clusters within a hash feature class is introduced.This constraint can ensure that the influence of convolutional channels on hash features is minimized.In addition,we design a new loss function called aggregate discriminative loss(AD-loss).The use of AD-loss and softmax-loss can increase the discriminativeness of the learned hash features.Experimental results show that on pool and ocean datasets,which were collected in pools and oceans,respectively,by using acoustic collectors,the proposed HADN performs better than other comparative models in terms of accuracy and mAP.展开更多
In this paper,we formulate the precoding problem of integrated sensing and communication(ISAC)waveform as a non-convex quadratically constrained quadratic programming(QCQP),in which the weighted sum of communication m...In this paper,we formulate the precoding problem of integrated sensing and communication(ISAC)waveform as a non-convex quadratically constrained quadratic programming(QCQP),in which the weighted sum of communication multi-user interference(MUI)and the gap between dual-use waveform and ideal radar waveform is minimized with peak-toaverage power ratio(PAPR)constraints.We propose an efficient algorithm based on alternating direction method of multipliers(ADMM),which is able to decouple multiple variables and provide a closed-form solution for each subproblem.In addition,to improve the sensing performance in both spatial and temporal domains,we propose a new criteria to design the ideal radar waveform,in which the beam pattern is made similar to the ideal one and the integrated sidelobe level of the ambiguity function in each target direction is minimized in the region of interest.The limited memory Broyden-Fletcher-Goldfarb-Shanno(LBFGS)algorithm is applied to the design of the ideal radar waveform which works as a reference in the design of the dual-function waveform.Numerical results indicate that the designed dual-function waveform is capable of offering good communication quality of service(QoS)and sensing performance.展开更多
Objective:To observe and compare the clinical effects of different electroacupuncture waveforms on primary dysmenorrhea.Methods: This was a prospective,randomized,three-group,parallel-controlled trial.Participants wit...Objective:To observe and compare the clinical effects of different electroacupuncture waveforms on primary dysmenorrhea.Methods: This was a prospective,randomized,three-group,parallel-controlled trial.Participants with primary dysmenorrhea were randomly divided into dense-sparse wave,continuous wave,and discontinuous wave groups in a 1:1:1 ratio.Two lateral Ciliao(BL 32)points were used.All three groups started treatment 3–5 days before menstruation,once a day for six sessions per course of treatment,one course of treatment per menstrual cycle,and three menstrual cycles.The primary outcome measure was the proportion with an average visual analog scale(VAS)score reduction of≥50%from baseline for dysmenorrhea in the third menstrual cycle during treatment.The secondary outcome measures included changes in dysmenorrhea VAS scores,Cox Menstrual Symptom Scale scores and the proportion of patients taking analgesic drugs.Results: The proportion of cases where the average VAS score for dysmenorrhea decreased by≥50%from baseline in the third menstrual cycle was not statistically significant(P>.05).Precisely 30 min after acupuncture and regarding immediate analgesia on the most severe day of dysmenorrhea,there was a statistically significant difference in the dense-sparse wave group compared with the other two groups during the third menstrual cycle(P<.05).Additionally,there was a statistically significant difference between the dense-sparse wave and discontinuous wave groups 24 h after acupuncture(P<.05).Conclusions: Waveform electroacupuncture can alleviate primary dysmenorrhea and its related symptoms in patients.The three groups showed similar results in terms of short-and long-term analgesic efficacy and a reduction in the number of patients taking analgesic drugs.Regarding achieving immediate analgesia,the dense-sparse wave group was slightly better than the other two groups.展开更多
Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with un...Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.展开更多
文摘Seismic migration and inversion are closely related techniques to portray subsurface images and identify hydrocarbon reservoirs.Seismic migration aims at obtaining structural images of subsurface geologic discontinuities.More specifically,seismic migration estimates the reflectivity function(stacked average reflectivity or pre-stack angle-dependent reflectivity)from seismic reflection data.On the other hand,seismic inversion quantitatively estimates the intrinsic rock properties of subsurface formulations.Such seismic inversion methods are applicable to detect hydrocarbon reservoirs that may exhibit lateral variations in the inverted parameters.Although there exist many differences,pre-stack seismic migration is similar with the first iteration of the general linearized seismic inversion.Usually,seismic migration and inversion techniques assume an acoustic or isotropic elastic medium.Unconventional reservoirs such as shale and tight sand formation have notable anisotropic property.We present a linearized waveform inversion(LWI)scheme for weakly anisotropic elastic media with vertical transversely isotropic(VTI)symmetry.It is based on two-way anisotropic elastic wave equation and simultaneously inverts for the localized perturbations(ΔVp_(0)/Vp_(0)/Vs_(0)/Vs_(0)/,Δ∈,Δδ)from the long-wavelength reference model.Our proposed VTI-elastic LWI is an iterative method that requires a forward and an adjoint operator acting on vectors in each iteration.We derive the forward Born approximation operator by perturbation theory and adjoint operator via adjoint-state method.The inversion has improved the quality of the images and reduces the multi-parameter crosstalk comparing with the adjoint-based images.We have observed that the multi-parameter crosstalk problem is more prominent in the inversion images for Thomsen anisotropy parameters.Especially,the Thomsen parameter is the most difficult to resolve.We also analyze the multi-parameter crosstalk using scattering radiation patterns.The linearized waveform inversion for VTI-elastic media presented in this article provides quantitative information of the rock properties that has the potential to help identify hydrocarbon reservoirs.
基金jointly supported by Young Scientists Cultivation Fund Project of Harbin Engineering University(79000013/003)the Mount Taishan Industrial Leading Talent Project+1 种基金the Great and Special Project under Grant KJGG-2022-0104 of CNOOC Limitedthe National Natural Science Foundation of China(42006064,42106070,42074138)。
文摘The low-wavenumber components in the gradient of full waveform inversion(FWI)play a vital role in the stability of the inversion.However,when FWI is implemented in some high frequencies and current models are not far away from the real velocity model,an excessive number of low-wavenumber components in the gradient will also reduce the convergence rate and inversion accuracy.To solve this problem,this paper firstly derives a formula of scattering angle weighted gradient in FWI,then proposes a hybrid gradient.The hybrid gradient combines the conventional gradient of FWI with the scattering angle weighted gradient in each inversion frequency band based on an empirical formula derived herein.Using weighted hybrid mode,we can retain some low-wavenumber components in the initial lowfrequency inversion to ensure the stability of the inversion,and use more high-wavenumber components in the high-frequency inversion to improve the convergence rate.The results of synthetic data experiment demonstrate that compared to the conventional FWI,the FWI based on the proposed hybrid gradient can effectively reduce the low-wavenumber components in the gradient under the premise of ensuring inversion stability.It also greatly enhances the convergence rate and inversion accuracy,especially in the deep part of the model.And the field marine seismic data experiment also illustrates that the FWI based on hybrid gradient(HGFWI)has good stability and adaptability.
基金supported by the National Natural Science Foundation of China (grant No. 42174156)the Young Science and Technology Star Project of Shaanxi Province (grant No. 2023KJXX-021)the Fundamental Research Funds for the Central Universities, CHD (grant Nos. 300102263401 and 300102264204)。
文摘Full-waveform inversion(FWI) uses the full information of seismic data to obtain a quantitative estimation of subsurface physical parameters. Anisotropic FWI has the potential to recover high-resolution velocity and anisotropy parameter models, which are critical for imaging the long-offset and wideazimuth data. We develop an acoustic anisotropic FWI method based on a simplified pure quasi P-wave(qP-wave) equation, which can be solved efficiently and is beneficial for the subsequent inversion.Using the inverse Hessian operator to precondition the functional gradients helps to reduce the parameter tradeoff in the multi-parameter inversion. To balance the accuracy and efficiency, we extend the truncated Gauss-Newton(TGN) method into FWI of pure qP-waves in vertical transverse isotropic(VTI) media. The inversion is performed in a nested way: a linear inner loop and a nonlinear outer loop.We derive the formulation of Hessian-vector products for pure qP-waves in VTI media based on the Lagrange multiplier method and compute the model update by solving a Gauss-Newton linear system via a matrix-free conjugate gradient method. A suitable preconditioner and the Eisenstat and Walker stopping criterion for the inner iterations are used to accelerate the convergence and avoid prohibitive computational cost. We test the proposed FWI method on several synthetic data sets. Inversion results reveal that the pure acoustic VTI FWI exhibits greater accuracy than the conventional pseudoacoustic VTI FWI. Additionally, the TGN method proves effective in mitigating the parameter crosstalk and increasing the accuracy of anisotropy parameters.
基金supported by the Major Research Plan on West-Pacific Earth System Multispheric Interactions (Nos.91858215,91958206)the National Natural Science Foundation of China (NSFC)Shiptime Sharing Project (No.41949581)the Key Research and Development Program of Shandong Province (No.2019GHY112019)。
文摘The precision and reliability of first-arrival picking are crucial for determining the accuracy of geological structure inversion using active source ocean bottom seismometer(OBS)refraction data.Traditional methods for first-arrival picking based on sample points are characterized by theoretical errors,especially in low-sampling-frequency OBS data because the travel time of seismic waves is not an integer multiple of the sampling interval.In this paper,a first-arrival picking method that utilizes the spatial waveform variation characteristics of active source OBS data is presented.First,the distribution law of theoretical error is examined;adjacent traces exhibit variation characteristics in their waveforms.Second,a label cross-correlation superposition method for extracting highfrequency signals is presented to enhance the first-arrival picking precision.Results from synthetic and field data verify that the proposed approach is robust,successfully overcomes the limitations of low sampling frequency,and achieves precise outcomes that are comparable with those of high-sampling-frequency data.
基金supported by the National Natural Science Foundation of China(No.62171052 and No.61971054)the Fundamental Research Funds for the Central Universities(No.24820232023YQTD01).
文摘Millimeter-wave(mmWave)radar communication has emerged as an important technique for future wireless systems.However,the interference between the radar signal and communication data is the main issue that should be considered for the joint radar communication system.In this paper,a co-sharing waveform(CSW)is proposed to achieve communication and radar sensing simultaneously.To eliminate the co-interference between the communication and sensing signal,signal splitting and processing methods for communication data demodulation and radar signal processing are given respectively.Simulation results show that the bit error rate(BER)of CSW is close to that of the pure communication waveform.Moreover,the proposed CSW can achieve better performance than the existing waveforms in terms of range and velocity estimation.
基金partially supported by the National Key Research and Development Program of China(No.2018 AAA0100400)the Natural Science Foundation of Shandong Province(Nos.ZR2020MF131 and ZR2021ZD19)the Science and Technology Program of Qingdao(No.21-1-4-ny-19-nsh).
文摘Underwater pulse waveform recognition is an important method for underwater object detection.Most existing works focus on the application of traditional pattern recognition methods,which ignore the time-and space-varying characteristics in sound propagation channels and cannot easily extract valuable waveform features.Sound propagation channels in seawater are time-and space-varying convolutional channels.In the extraction of the waveform features of underwater acoustic signals,the effect of high-accuracy underwater acoustic signal recognition is identified by eliminating the influence of time-and space-varying convolutional channels to the greatest extent possible.We propose a hash aggregate discriminative network(HADN),which combines hash learning and deep learning to minimize the time-and space-varying effects on convolutional channels and adaptively learns effective underwater waveform features to achieve high-accuracy underwater pulse waveform recognition.In the extraction of the hash features of acoustic signals,a discrete constraint between clusters within a hash feature class is introduced.This constraint can ensure that the influence of convolutional channels on hash features is minimized.In addition,we design a new loss function called aggregate discriminative loss(AD-loss).The use of AD-loss and softmax-loss can increase the discriminativeness of the learned hash features.Experimental results show that on pool and ocean datasets,which were collected in pools and oceans,respectively,by using acoustic collectors,the proposed HADN performs better than other comparative models in terms of accuracy and mAP.
基金supported in part by the National Natural Science Foundation of China under Grant 62271142in part by the Key Research and Development Program of Jiangsu Province BE2023021+2 种基金in part by the Jiangsu Key Research and Development Program Project under Grant BE2023011-2in part by the Young Scholar Funding of Southeast Universityin part by the Fundamental Research Funds for the Central Universities 2242022k60001。
文摘In this paper,we formulate the precoding problem of integrated sensing and communication(ISAC)waveform as a non-convex quadratically constrained quadratic programming(QCQP),in which the weighted sum of communication multi-user interference(MUI)and the gap between dual-use waveform and ideal radar waveform is minimized with peak-toaverage power ratio(PAPR)constraints.We propose an efficient algorithm based on alternating direction method of multipliers(ADMM),which is able to decouple multiple variables and provide a closed-form solution for each subproblem.In addition,to improve the sensing performance in both spatial and temporal domains,we propose a new criteria to design the ideal radar waveform,in which the beam pattern is made similar to the ideal one and the integrated sidelobe level of the ambiguity function in each target direction is minimized in the region of interest.The limited memory Broyden-Fletcher-Goldfarb-Shanno(LBFGS)algorithm is applied to the design of the ideal radar waveform which works as a reference in the design of the dual-function waveform.Numerical results indicate that the designed dual-function waveform is capable of offering good communication quality of service(QoS)and sensing performance.
基金supported by Technology Innovation Special Project of Dongzhimen Hospital affiliated to Beijing University of Chinese Medicine.
文摘Objective:To observe and compare the clinical effects of different electroacupuncture waveforms on primary dysmenorrhea.Methods: This was a prospective,randomized,three-group,parallel-controlled trial.Participants with primary dysmenorrhea were randomly divided into dense-sparse wave,continuous wave,and discontinuous wave groups in a 1:1:1 ratio.Two lateral Ciliao(BL 32)points were used.All three groups started treatment 3–5 days before menstruation,once a day for six sessions per course of treatment,one course of treatment per menstrual cycle,and three menstrual cycles.The primary outcome measure was the proportion with an average visual analog scale(VAS)score reduction of≥50%from baseline for dysmenorrhea in the third menstrual cycle during treatment.The secondary outcome measures included changes in dysmenorrhea VAS scores,Cox Menstrual Symptom Scale scores and the proportion of patients taking analgesic drugs.Results: The proportion of cases where the average VAS score for dysmenorrhea decreased by≥50%from baseline in the third menstrual cycle was not statistically significant(P>.05).Precisely 30 min after acupuncture and regarding immediate analgesia on the most severe day of dysmenorrhea,there was a statistically significant difference in the dense-sparse wave group compared with the other two groups during the third menstrual cycle(P<.05).Additionally,there was a statistically significant difference between the dense-sparse wave and discontinuous wave groups 24 h after acupuncture(P<.05).Conclusions: Waveform electroacupuncture can alleviate primary dysmenorrhea and its related symptoms in patients.The three groups showed similar results in terms of short-and long-term analgesic efficacy and a reduction in the number of patients taking analgesic drugs.Regarding achieving immediate analgesia,the dense-sparse wave group was slightly better than the other two groups.
文摘Simultaneous waveform inversion was used to predict lithofacies and fluid type across the field. Very often, characterizing reservoirs in terms of lithology and fluid type using conventional methods is replete with uncertainties, especially in marginal fields. An approach is employed in this study that integrated rock physics and waveform inverse modelling for lithology and fluid-type characterization to appropriately identify potential hydrocarbon saturated zones and their corresponding lithology. Seismic and well-log data were analyzed using Hampson Russel software. The method adopted includes lithofacies and fluid content analysis using rock physics parameters and seismic simultaneous inverse modelling. Rock physics analysis identified 2 broad reservoirs namely: HDZ1 and HDZ2 reservoirs. Results from the inverse modelling showed that low values of acoustic impedance from 19,743 to 20,487 (ft/s)(g/cc) reflect hydrocarbon-bearing reservoirs while medium to high values shows brine and shale respectively, with brine zone ranging from 20,487 to 22,531 (ft/s)(g/cc) and shale above 22,531 (ft/s)(g/cc). Two lithofacies were identified from inversion analysis of Vp/Vs and Mu-Rho, namely: sand and shale with VpVs 1.95 values respectively. Mu-Rho > 12.29 (GPa)(g/cc) and <12.29 (GPa) (g/cc) represent sand and shale respectively. From 3D volume, it was observed that a high accumulation of hydrocarbon was observed to be saturated at the north to the eastern part of the field forming a meandering channel. Sands were mainly distributed around the northeastern to the southwestern part of the field, that tends to be away from Well 029. This was also validated by the volume of rigidity modulus (Mu-Rho) showing high values indicating sands fall within the northeastern part of the field.
