Silin Wang, Cai Liu, Peng Li, Changle Chen, Chao Song
Seismic impedance inversion makes a significant contribution to locating hydrocarbons and interpreting seismic data. However, it suffers from non-unique solutions, and a direct linear inversion produces large errors. Global optimization methods, like simulated annealing, have been applied in seismic impedance inversion and achieved promising inversion results. Over the past decades, there has been an increasing interest in quantum computing. Due to its natural parallelism, quantum computing has a powerful computational capability and certain advantages in solving complex inverse problems. Within this article, we present a stable and efficient impedance inversion using quantum annealing with L1 norm regularization, which significantly improves the inversion accuracy compared to the traditional simulated annealing method. Tests on a one-dimensional ten-layer model with noisy data demonstrate that the new method exhibits significantly improved accuracy and stability. Additionally, we perform seismic impedance inversion for synthetic data from the Overthrust model and field data using two methods. These results demonstrate that the quantum annealing impedance inversion with L1 norm regularization dramatically enhances the accuracy and anti-noise ability.
{"title":"Stable and efficient seismic impedance inversion using quantum annealing with L1 norm regularization","authors":"Silin Wang, Cai Liu, Peng Li, Changle Chen, Chao Song","doi":"10.1093/jge/gxae003","DOIUrl":"https://doi.org/10.1093/jge/gxae003","url":null,"abstract":"\u0000 Seismic impedance inversion makes a significant contribution to locating hydrocarbons and interpreting seismic data. However, it suffers from non-unique solutions, and a direct linear inversion produces large errors. Global optimization methods, like simulated annealing, have been applied in seismic impedance inversion and achieved promising inversion results. Over the past decades, there has been an increasing interest in quantum computing. Due to its natural parallelism, quantum computing has a powerful computational capability and certain advantages in solving complex inverse problems. Within this article, we present a stable and efficient impedance inversion using quantum annealing with L1 norm regularization, which significantly improves the inversion accuracy compared to the traditional simulated annealing method. Tests on a one-dimensional ten-layer model with noisy data demonstrate that the new method exhibits significantly improved accuracy and stability. Additionally, we perform seismic impedance inversion for synthetic data from the Overthrust model and field data using two methods. These results demonstrate that the quantum annealing impedance inversion with L1 norm regularization dramatically enhances the accuracy and anti-noise ability.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2024-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139620562","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Yuyu Wu, Qiang Lai, Haitao Wang, Rong Yin, Zeyu Wang, Li Chen, Fuqiang Lai
In the Sichuan Basin, bitumen widely distributes in the carbonate reservoir, and the bitumen occurrence morphology affects the pore structure and petrophysics. Two types of bitumen, boundary-filled and centre-filled, were reconstructed using a pore morphology-based method based on dolomite carbonate digital cores from X-ray computerised tomography. The pore structures of the reconstructed bitumen-filled digital core models were determined using an extracted pore network model. The finite element method simulates the electrical resistivity of the digital core model. The investigation shows that (1) boundary-developed bitumen increases the pore size whereas centre-developed bitumen decreases the pore size, and bitumen development increases the shape factor while reducing the coordination number; (2) for electrical resistivity, boundary-filled bitumen increases the porosity exponent with increasing pore size, whereas centre-filled bitumen decreases the porosity exponent. Bitumen development caused a negative correlation between the porosity exponent and pore shape factor, a logarithmic increase between the porosity exponent and coordination number, and a power increase between the porosity exponent and porosity. Subsequently, a new porosity exponent model related to porosity is proposed based on the curve-fitting method. This new porosity exponent combines the Archie equations to calculate the saturation of the carbonate reservoir in the Sinian Dengying Formation in the Penglai gas area on the northern slope of the central Sichuan paleo-uplift, which fits the sealed core plug saturation well and provides the basis for further reservoir evaluation and fluid identification in the Penglai gas area.
在四川盆地，沥青广泛分布于碳酸盐岩储层中，沥青的赋存形态影响着孔隙结构和岩石物理。基于 X 射线计算机断层扫描的白云质碳酸盐岩数字岩心，采用基于孔隙形态的方法重建了边界充填和中心充填两种类型的沥青。利用提取的孔隙网络模型确定了重建的沥青填充数字岩心模型的孔隙结构。有限元法模拟了数字岩心模型的电阻率。研究表明：(1) 边界发育的沥青增加了孔隙尺寸，而中心发育的沥青减小了孔隙尺寸，沥青发育增加了形状系数，同时降低了配位数；(2) 对于电阻率，边界填充的沥青随着孔隙尺寸的增加而增加孔隙度指数，而中心填充的沥青则降低孔隙度指数。沥青的发展导致孔隙度指数与孔隙形状系数之间呈负相关，孔隙度指数与配位数之间呈对数增长，孔隙度指数与孔隙度之间呈幂增长。随后，基于曲线拟合方法，提出了与孔隙度相关的新孔隙度指数模型。这一新的孔隙度指数结合阿奇方程，计算出了川中古隆起北坡蓬莱气区中年登瀛层碳酸盐岩储层的饱和度，很好地拟合了密封岩心塞饱和度，为蓬莱气区进一步的储层评价和流体识别提供了依据。
{"title":"The bitumen development effect on the pore structure and electrical resistivity of carbonate reservoirs in Sinian Dengying Formation","authors":"Yuyu Wu, Qiang Lai, Haitao Wang, Rong Yin, Zeyu Wang, Li Chen, Fuqiang Lai","doi":"10.1093/jge/gxae008","DOIUrl":"https://doi.org/10.1093/jge/gxae008","url":null,"abstract":"\u0000 In the Sichuan Basin, bitumen widely distributes in the carbonate reservoir, and the bitumen occurrence morphology affects the pore structure and petrophysics. Two types of bitumen, boundary-filled and centre-filled, were reconstructed using a pore morphology-based method based on dolomite carbonate digital cores from X-ray computerised tomography. The pore structures of the reconstructed bitumen-filled digital core models were determined using an extracted pore network model. The finite element method simulates the electrical resistivity of the digital core model. The investigation shows that (1) boundary-developed bitumen increases the pore size whereas centre-developed bitumen decreases the pore size, and bitumen development increases the shape factor while reducing the coordination number; (2) for electrical resistivity, boundary-filled bitumen increases the porosity exponent with increasing pore size, whereas centre-filled bitumen decreases the porosity exponent. Bitumen development caused a negative correlation between the porosity exponent and pore shape factor, a logarithmic increase between the porosity exponent and coordination number, and a power increase between the porosity exponent and porosity. Subsequently, a new porosity exponent model related to porosity is proposed based on the curve-fitting method. This new porosity exponent combines the Archie equations to calculate the saturation of the carbonate reservoir in the Sinian Dengying Formation in the Penglai gas area on the northern slope of the central Sichuan paleo-uplift, which fits the sealed core plug saturation well and provides the basis for further reservoir evaluation and fluid identification in the Penglai gas area.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2024-01-16","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139618669","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Jian Li, Guoning Du, Dewen Qin, Wensun Yin, Jun Tan, Zhaolun Liu, Peng Song
Although the excitation-time imaging condition offers a lower memory consumption and higher computational efficiency compared to cross-correlation imaging condition, it has not been widely used in industrial applications because of the accuracy problem of traveltime calculation and the influence of low-wave-number noise. In this paper, we introduce the physics-informed neural network (PINN) algorithm to achieve a high-precision traveltime calculation of the source forward wavefield. Subsequently, we introduce a technique for high-precision wavefield decomposition of the reverse-time wavefield via the optical flow vector, enabling us to realize a correlation-weighted stacking imaging of each wavefield. Model experiments and real data processing show that the proposed traveltime calculation algorithm based on PINN offers high accuracy and good applicability in the excitation time reverse-time migration imaging of complex models, and correlation-weighted stacking imaging based on optical flow vector-based wavefield separation can significantly suppress the noise with low wavenumber and achieve high-precision imaging of complex models.
{"title":"Excitation time imaging condition reverse time migration based on physics-informed neural network traveltime calculation with wavefield decomposition using optical flow vector","authors":"Jian Li, Guoning Du, Dewen Qin, Wensun Yin, Jun Tan, Zhaolun Liu, Peng Song","doi":"10.1093/jge/gxad106","DOIUrl":"https://doi.org/10.1093/jge/gxad106","url":null,"abstract":"\u0000 Although the excitation-time imaging condition offers a lower memory consumption and higher computational efficiency compared to cross-correlation imaging condition, it has not been widely used in industrial applications because of the accuracy problem of traveltime calculation and the influence of low-wave-number noise. In this paper, we introduce the physics-informed neural network (PINN) algorithm to achieve a high-precision traveltime calculation of the source forward wavefield. Subsequently, we introduce a technique for high-precision wavefield decomposition of the reverse-time wavefield via the optical flow vector, enabling us to realize a correlation-weighted stacking imaging of each wavefield. Model experiments and real data processing show that the proposed traveltime calculation algorithm based on PINN offers high accuracy and good applicability in the excitation time reverse-time migration imaging of complex models, and correlation-weighted stacking imaging based on optical flow vector-based wavefield separation can significantly suppress the noise with low wavenumber and achieve high-precision imaging of complex models.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2024-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139534609","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Bin Wang, Jie Zhang, Haifei Lin, Jiangfei Duan, Guangliang Li, Qingsong Zhuo, Jianjun Wu, Yifeng He, Shoushi Gao
In order to study the failure mechanism of overlying strata (OS) in shallow insufficient mining areas, with a combination of such research approaches as field investigation, theoretical analysis, similarity simulation, and numerical simulation, this paper has studied the temporal and spatial evolution of the failure mechanism of overlying stratum structure in a shallow-buried interval goaf. The results show that the creep failure of temporary coal pillar (TCP) in the interval goaf is a primary reason for the failure of the basic roof. With the failure of the basic roof, stress arches in the OS of the mining section become unstable, which expands the damage range of the overlying strata. Consequently, adjacent stress arches overlap with each other, forming a ‘trapezoidal—semi-elliptical arch’ collapse shape. Thick soil layers gradually collapse to the ground surface, and the overlying strata collapse as a whole into a trapezoid-like shape. Rotary failure appears in the basic roof of mining section, forming a ‘W-shaped voussoir beam’ hinge structure. In this study, a structural model of ‘W-shaped voussoir beam’ in the OS has been established, and the mechanical characteristics of rock blocks in the basic roof of overlying strata during different mining stages of the interval goaf have been analyzed. Also, with a discrete element UDEC program, this study has performed a simulation to verify the rationality of the ‘W-shaped voussoir beam’ structural model of overlying stratum movement in the shallowly-buried interval mining section.
{"title":"Instability mechanism of overlying strata in shallow interval goaf and structural model of “W-shaped voussoir beam”","authors":"Bin Wang, Jie Zhang, Haifei Lin, Jiangfei Duan, Guangliang Li, Qingsong Zhuo, Jianjun Wu, Yifeng He, Shoushi Gao","doi":"10.1093/jge/gxae004","DOIUrl":"https://doi.org/10.1093/jge/gxae004","url":null,"abstract":"\u0000 In order to study the failure mechanism of overlying strata (OS) in shallow insufficient mining areas, with a combination of such research approaches as field investigation, theoretical analysis, similarity simulation, and numerical simulation, this paper has studied the temporal and spatial evolution of the failure mechanism of overlying stratum structure in a shallow-buried interval goaf. The results show that the creep failure of temporary coal pillar (TCP) in the interval goaf is a primary reason for the failure of the basic roof. With the failure of the basic roof, stress arches in the OS of the mining section become unstable, which expands the damage range of the overlying strata. Consequently, adjacent stress arches overlap with each other, forming a ‘trapezoidal—semi-elliptical arch’ collapse shape. Thick soil layers gradually collapse to the ground surface, and the overlying strata collapse as a whole into a trapezoid-like shape. Rotary failure appears in the basic roof of mining section, forming a ‘W-shaped voussoir beam’ hinge structure. In this study, a structural model of ‘W-shaped voussoir beam’ in the OS has been established, and the mechanical characteristics of rock blocks in the basic roof of overlying strata during different mining stages of the interval goaf have been analyzed. Also, with a discrete element UDEC program, this study has performed a simulation to verify the rationality of the ‘W-shaped voussoir beam’ structural model of overlying stratum movement in the shallowly-buried interval mining section.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2024-01-10","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139627782","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
The cross-correlation-based methods, widely used for microseismic monitoring, utilize cross-correlation to extract time differences of signals within station pairs, and subsequently use these time differences for back-projection and localization, without the need for triggering moment scanning. The selection of imaging conditions, applied to the spatial projection of all cross-correlation records determines the noise resistance and resolution of such methods. To fully exploit the azimuthal properties of station pairs for constraining the source imaging, we propose a high-order azimuth coherent imaging condition, which involves the following steps: (1) Choosing station pairs that meet specific inter-station distance criteria; (2) Combining station pairs into dual-station pairs that satisfy a certain inter-station pair angle criterion; (3) Further combining station combinations pairwise to form the final set of station pairs; (4) Multiplying the projection results of station pairs within each combination from the third step; (5) Summing the results of all combinations. This method effectively suppresses the hyperboloid in the spatial projection of a single station, enhances the coherence of seismic source imaging, and maintains noise resistance. Compared to the conventional imaging conditions, the method proposed has demonstrated superior resolution and robustness in both theoretical analysis and practical testing.
{"title":"High-order Azimuth Coherent Imaging for Microseismic Location","authors":"Ying Shi, Xuebao Guo, Youqiang Yu","doi":"10.1093/jge/gxae001","DOIUrl":"https://doi.org/10.1093/jge/gxae001","url":null,"abstract":"\u0000 The cross-correlation-based methods, widely used for microseismic monitoring, utilize cross-correlation to extract time differences of signals within station pairs, and subsequently use these time differences for back-projection and localization, without the need for triggering moment scanning. The selection of imaging conditions, applied to the spatial projection of all cross-correlation records determines the noise resistance and resolution of such methods. To fully exploit the azimuthal properties of station pairs for constraining the source imaging, we propose a high-order azimuth coherent imaging condition, which involves the following steps: (1) Choosing station pairs that meet specific inter-station distance criteria; (2) Combining station pairs into dual-station pairs that satisfy a certain inter-station pair angle criterion; (3) Further combining station combinations pairwise to form the final set of station pairs; (4) Multiplying the projection results of station pairs within each combination from the third step; (5) Summing the results of all combinations. This method effectively suppresses the hyperboloid in the spatial projection of a single station, enhances the coherence of seismic source imaging, and maintains noise resistance. Compared to the conventional imaging conditions, the method proposed has demonstrated superior resolution and robustness in both theoretical analysis and practical testing.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2024-01-05","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139383539","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Low-frequency seismic data plays a crucial role in seismic data processing and seismic wave inversion. At present, there are two methods to realize the low-frequency excitation of vibrators: one is that the low-frequency vibrators are excited by linear sweep signals, and the other is that the conventional vibrators are excited by nonlinear low-frequency sweep signals. The cost of exploration using low-frequency vibroseis is high, and it is challenging to obtain sufficient low-frequency information using traditional vibrators. To this end, this paper comparatively studies the low-frequency sweep signal characteristics and data effects of low-frequency vibrator and traditional vibrator. Therefore, three kinds of linear and nonlinear low-frequency sweep signals are designed. Theoretical analysis shows that there are certain differences between linear and nonlinear in design methods, signal shapes, etc., but after correlation calculation, the signal spectrums reflecting the seismic response and the related wavelet shapes are basically consistent. Besides, the actual force signal data shows that the linear and nonlinear harmonic distortion are basically equivalent. Finally, based on the forward simulation of three sweep signals and the comparative analysis of field test data, it can be considered that the linear and nonlinear low-frequency sweep signals of vibrator have almost the same denoising ability under the basic conditions of spectrum and wavelet. Both can achieve low-frequency excitation and obtain rich low-frequency information, and the quality of seismic data is basically the same, so they can be applied in practical production.
{"title":"Characteristic analysis and data comparative of linear and nonlinear low frequency sweep in vibroseis","authors":"Jian Zhang, Guoyong Zhao, Leiliang Xu, Xu Zhang, Xiaohan Liu, Jianggui Zhu","doi":"10.1093/jge/gxad108","DOIUrl":"https://doi.org/10.1093/jge/gxad108","url":null,"abstract":"Low-frequency seismic data plays a crucial role in seismic data processing and seismic wave inversion. At present, there are two methods to realize the low-frequency excitation of vibrators: one is that the low-frequency vibrators are excited by linear sweep signals, and the other is that the conventional vibrators are excited by nonlinear low-frequency sweep signals. The cost of exploration using low-frequency vibroseis is high, and it is challenging to obtain sufficient low-frequency information using traditional vibrators. To this end, this paper comparatively studies the low-frequency sweep signal characteristics and data effects of low-frequency vibrator and traditional vibrator. Therefore, three kinds of linear and nonlinear low-frequency sweep signals are designed. Theoretical analysis shows that there are certain differences between linear and nonlinear in design methods, signal shapes, etc., but after correlation calculation, the signal spectrums reflecting the seismic response and the related wavelet shapes are basically consistent. Besides, the actual force signal data shows that the linear and nonlinear harmonic distortion are basically equivalent. Finally, based on the forward simulation of three sweep signals and the comparative analysis of field test data, it can be considered that the linear and nonlinear low-frequency sweep signals of vibrator have almost the same denoising ability under the basic conditions of spectrum and wavelet. Both can achieve low-frequency excitation and obtain rich low-frequency information, and the quality of seismic data is basically the same, so they can be applied in practical production.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2023-12-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139145534","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Mingkun Zhang, Lingqian Wang, Hanming Chen, Hui Zhou, Peng Liu
Least-squares reverse-time migration (LSRTM) has become an advanced technique for complex structures imaging of the subsurface, as it can provide a higher resolution and more balanced amplitude migrated image than conventional reverse-time migration (RTM). However, the intrinsic attenuation of subsurface introduces amplitude attenuation and phase dispersion of seismic wavefield, which leads to the inverted image kinematically and dynamically inexactitude. Moreover, the imperfect geometry, limited bandwidth of seismic data, and inappropriate modeling kernel etc., would inevitably introduce two side-effects in migrated image, resulting in degradation of LSRTM imaging potential. To alleviate above issues, we present a data-domain sparsity constraint viscoacoustic least-squares reverse-time migration algorithm in this paper. In particular, we utilize the decoupled constant Q fractional Laplacians (DFLs) viscoacoustic wave equation as the modeling kernel to describe the attenuation effects of the subsurface, while a model constraint constructed in the misfit function via L1-2 norm is carried out to clear the migrated artefacts and boost the imaging resolution. Thanks to the excellent performance in sparsity, the drawbacks of unconstraint LSRTM can be effectively mitigated by the L1-2 norm-based regularization. In this paper, we adopt the alternating direction of multipliers method (ADMM) to iteratively address the constrained L1-2 minimization problem by implementing a proximal operator, and three synthetic examples are hired to evaluate the effectiveness and practicability of the proposed strategy. Migration results prove that the proposed scheme can effectively compensate the attenuation effects, improve the resolution, and suppress the migration artifacts of inverted images even in the complex imaging situations.
{"title":"Viscoacoustic least squares reverse-time migration using L1-2 norm sparsity constraint","authors":"Mingkun Zhang, Lingqian Wang, Hanming Chen, Hui Zhou, Peng Liu","doi":"10.1093/jge/gxad107","DOIUrl":"https://doi.org/10.1093/jge/gxad107","url":null,"abstract":"Least-squares reverse-time migration (LSRTM) has become an advanced technique for complex structures imaging of the subsurface, as it can provide a higher resolution and more balanced amplitude migrated image than conventional reverse-time migration (RTM). However, the intrinsic attenuation of subsurface introduces amplitude attenuation and phase dispersion of seismic wavefield, which leads to the inverted image kinematically and dynamically inexactitude. Moreover, the imperfect geometry, limited bandwidth of seismic data, and inappropriate modeling kernel etc., would inevitably introduce two side-effects in migrated image, resulting in degradation of LSRTM imaging potential. To alleviate above issues, we present a data-domain sparsity constraint viscoacoustic least-squares reverse-time migration algorithm in this paper. In particular, we utilize the decoupled constant Q fractional Laplacians (DFLs) viscoacoustic wave equation as the modeling kernel to describe the attenuation effects of the subsurface, while a model constraint constructed in the misfit function via L1-2 norm is carried out to clear the migrated artefacts and boost the imaging resolution. Thanks to the excellent performance in sparsity, the drawbacks of unconstraint LSRTM can be effectively mitigated by the L1-2 norm-based regularization. In this paper, we adopt the alternating direction of multipliers method (ADMM) to iteratively address the constrained L1-2 minimization problem by implementing a proximal operator, and three synthetic examples are hired to evaluate the effectiveness and practicability of the proposed strategy. Migration results prove that the proposed scheme can effectively compensate the attenuation effects, improve the resolution, and suppress the migration artifacts of inverted images even in the complex imaging situations.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2023-12-26","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139157326","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Haihong Xu, Jianshe Wei, Xiaofeng Han, Fei Zhao, Junlin Zhou, Ting Jiang, Jizhong Shi, Wei Xu, Bo Song, Baowen Wang
This research aims to explore petroleum resource distribution characteristics in the Yingen-Ejinaqi Basin and provide a basis for assessing the prospect of petroleum potential in the region. Based on recent and previous gravity and aeromagnetic data obtained in the basin, anomalies including regional and residual gravity and magnetic anomalies, various derivative anomalies, and anomalies in structure-enhancing filtering were determined using multiple techniques. By analysing gravity, magnetic, and seismic anomaly characteristics of existing typical oil and gas wells, the petroleum resource distribution was studied by combining such data with exploration results arising from wells where oil and gas have been discovered. The research shows that the distribution ranges of the sags revealed by seismic profiles correspond to those of negative residual gravity anomalies; the distribution ranges of rocks identified in the sags along the seismic profiles show some correspondence to the positive residual magnetic anomalies. Analysis indicates that various sags corresponding to areas where negative residual gravity anomalies appear in the north and south of the basin are the primary and secondary favourable areas for oil and gas exploration in Mesozoic-Upper Paleozoic strata. Therefore, breakthroughs in oil and gas exploration should be made by focusing on these sags in the region. The positive residual magnetic anomalies in the basin can reflect the distribution of concealed rocks in the region, while concealed rocks in sags are an important area for oil and gas exploration in igneous rocks in the basin.
{"title":"Characteristics of gravity and magnetic anomalies and their petroleum geological significance in the Yingen-Ejinaqi Basin, Inner Mongolia, China","authors":"Haihong Xu, Jianshe Wei, Xiaofeng Han, Fei Zhao, Junlin Zhou, Ting Jiang, Jizhong Shi, Wei Xu, Bo Song, Baowen Wang","doi":"10.1093/jge/gxad102","DOIUrl":"https://doi.org/10.1093/jge/gxad102","url":null,"abstract":"This research aims to explore petroleum resource distribution characteristics in the Yingen-Ejinaqi Basin and provide a basis for assessing the prospect of petroleum potential in the region. Based on recent and previous gravity and aeromagnetic data obtained in the basin, anomalies including regional and residual gravity and magnetic anomalies, various derivative anomalies, and anomalies in structure-enhancing filtering were determined using multiple techniques. By analysing gravity, magnetic, and seismic anomaly characteristics of existing typical oil and gas wells, the petroleum resource distribution was studied by combining such data with exploration results arising from wells where oil and gas have been discovered. The research shows that the distribution ranges of the sags revealed by seismic profiles correspond to those of negative residual gravity anomalies; the distribution ranges of rocks identified in the sags along the seismic profiles show some correspondence to the positive residual magnetic anomalies. Analysis indicates that various sags corresponding to areas where negative residual gravity anomalies appear in the north and south of the basin are the primary and secondary favourable areas for oil and gas exploration in Mesozoic-Upper Paleozoic strata. Therefore, breakthroughs in oil and gas exploration should be made by focusing on these sags in the region. The positive residual magnetic anomalies in the basin can reflect the distribution of concealed rocks in the region, while concealed rocks in sags are an important area for oil and gas exploration in igneous rocks in the basin.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2023-12-23","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139161966","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
C. Maschio, Gilson M. Silva Neto, A. Davolio, Vinicius de Souza Rios, D. Schiozer
The use of 4D seismic (4DS) (or time-lapse seismic, TLS) in data assimilation (DA) makes the process more complex due to the higher amount of data to be assimilated, requiring more robust methods and better computational resources (processing capacity and memory). The development and application of permanent seismic monitoring technologies have increased in the last years, improving the overall 4D seismic quality, in terms of signal resolution and repeatability. However, a massive amount of data is generated from the multiple monitors, making the incorporation of 4DS data in the DA process much more complex. Therefore, robust DA methods capable of dealing with huge amount of data effectively and efficiently are essential. This paper aims to assess the performance of an iterative ensemble smoother method, named Subspace Ensemble Randomized Maximum Likelihood with a local analysis (SEnRML-LA), to assimilate a big data set. The method was applied in a challenging pre-salt-like benchmark case with eight seismic surveys, one base, and seven monitors. The 4DS data are the impedance ratios (between two consecutive monitors) in 15 seismic horizons, totaling 105 maps to be assimilated. To our best knowledge, this is state of the art in terms of practical applications in data assimilation. It was possible to assimilate all the data simultaneously: the 105 horizons for the 4DS data and the wells’ production and pressure data. The data assimilation was successful in terms of results quality and method performance. We also ran a case assimilating only well data for comparison purposes.
{"title":"Performance assessment of an iterative ensemble smoother with local analysis to assimilate big 4D seismic dataset applied to a complex pre-salt-like benchmark case","authors":"C. Maschio, Gilson M. Silva Neto, A. Davolio, Vinicius de Souza Rios, D. Schiozer","doi":"10.1093/jge/gxad099","DOIUrl":"https://doi.org/10.1093/jge/gxad099","url":null,"abstract":"The use of 4D seismic (4DS) (or time-lapse seismic, TLS) in data assimilation (DA) makes the process more complex due to the higher amount of data to be assimilated, requiring more robust methods and better computational resources (processing capacity and memory). The development and application of permanent seismic monitoring technologies have increased in the last years, improving the overall 4D seismic quality, in terms of signal resolution and repeatability. However, a massive amount of data is generated from the multiple monitors, making the incorporation of 4DS data in the DA process much more complex. Therefore, robust DA methods capable of dealing with huge amount of data effectively and efficiently are essential. This paper aims to assess the performance of an iterative ensemble smoother method, named Subspace Ensemble Randomized Maximum Likelihood with a local analysis (SEnRML-LA), to assimilate a big data set. The method was applied in a challenging pre-salt-like benchmark case with eight seismic surveys, one base, and seven monitors. The 4DS data are the impedance ratios (between two consecutive monitors) in 15 seismic horizons, totaling 105 maps to be assimilated. To our best knowledge, this is state of the art in terms of practical applications in data assimilation. It was possible to assimilate all the data simultaneously: the 105 horizons for the 4DS data and the wells’ production and pressure data. The data assimilation was successful in terms of results quality and method performance. We also ran a case assimilating only well data for comparison purposes.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2023-11-30","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139197913","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
Surface offset gathers (SOGs) are crucial in velocity updating in seismic data migration and muting the stretched waveform in shallow parts of migrated seismic data. Kirchhoff-based migration methods can output SOGs because it takes every seismic trace as the input. As a result, the imaging results can be rearranged according to the offset value of the input data to obtain SOGs. However, regarding more accurate wave equation methods (such as reverse time migration (RTM)), the SOGs cannot be obtained directly with a single migration calculation. Attribute migration is an easy-to-implement method to output SOGs of wave equation migration methods. It calculates the offset value of each imaging point by dividing the migration result modulated by the offset attribute with the conventional migration result. However, the division error may cause the calculation result outside the given offset value. This paper proposes a constrained least square inversion method for stable calculation of RTM offset range. The method ensures that the calculated results are within the given offset range. We tested the method against the direct division and least square methods without constraints using the Marmousi model and real 3D dataset. We show that the SOGs obtained by the constrained least squares inversion attribute migration method had the same characteristics as those derived using the division-based attribute migration method, but also has higher vertical resolution, better energy distribution, and improved lateral continuity. In further study, we expect to develop a stable and direct method to efficiently calculate SOGs for RTM and an iterative closed loop for RTM velocity updating.
{"title":"Reverse time migration surface offset gathers by attribute migration and constrained least squares inversion","authors":"Yingdong Zheng, Guofeng Liu","doi":"10.1093/jge/gxad097","DOIUrl":"https://doi.org/10.1093/jge/gxad097","url":null,"abstract":"Surface offset gathers (SOGs) are crucial in velocity updating in seismic data migration and muting the stretched waveform in shallow parts of migrated seismic data. Kirchhoff-based migration methods can output SOGs because it takes every seismic trace as the input. As a result, the imaging results can be rearranged according to the offset value of the input data to obtain SOGs. However, regarding more accurate wave equation methods (such as reverse time migration (RTM)), the SOGs cannot be obtained directly with a single migration calculation. Attribute migration is an easy-to-implement method to output SOGs of wave equation migration methods. It calculates the offset value of each imaging point by dividing the migration result modulated by the offset attribute with the conventional migration result. However, the division error may cause the calculation result outside the given offset value. This paper proposes a constrained least square inversion method for stable calculation of RTM offset range. The method ensures that the calculated results are within the given offset range. We tested the method against the direct division and least square methods without constraints using the Marmousi model and real 3D dataset. We show that the SOGs obtained by the constrained least squares inversion attribute migration method had the same characteristics as those derived using the division-based attribute migration method, but also has higher vertical resolution, better energy distribution, and improved lateral continuity. In further study, we expect to develop a stable and direct method to efficiently calculate SOGs for RTM and an iterative closed loop for RTM velocity updating.","PeriodicalId":54820,"journal":{"name":"Journal of Geophysics and Engineering","volume":null,"pages":null},"PeriodicalIF":1.4,"publicationDate":"2023-11-29","publicationTypes":"Journal Article","fieldsOfStudy":null,"isOpenAccess":false,"openAccessPdf":"","citationCount":null,"resultStr":null,"platform":"Semanticscholar","paperid":"139211502","PeriodicalName":null,"FirstCategoryId":null,"ListUrlMain":null,"RegionNum":3,"RegionCategory":"地球科学","ArticlePicture":[],"TitleCN":null,"AbstractTextCN":null,"PMCID":"","EPubDate":null,"PubModel":null,"JCR":null,"JCRName":null,"Score":null,"Total":0}
京公网安备 11010802042870号
京ICP备2023020795号-1
扫码关注我们
求助内容：
应助结果提醒方式：