日本的天然气水合物地质调查工作

介绍了日本天然气水合物研究概况、研究计划和对日本南海海槽天然气水合物的地质调查工作。

3D seismic forward modeling from the multiphysical inversion at the Ketzin CO_(2) storage site

From June 2008 to August 2013,approximately 67 kt of CO_(2) was injected into a deep saline formation at the Ketzin pilot CO_(2) storage site.During injection,3D seismic surveys have been performed to monitor the migration of sequestered CO_(2).Seismic monitoring results are limited by the acquisition and signal-to-noise ratio of the acquired data.The multiphysical reservoir simulation provides information regarding the CO_(2) fluid behavior,and the approximated model should be calibrated with the monitoring results.In this work,property models are delivered from the multiphysical model during 3D repeated seismic surveys.The simulated seismic data based on the models are compared with the real data,and the results validate the effectiveness of the multiphysical inversion method.Time-lapse analysis shows the trend of CO_(2) migration during and after injection.

Tectonic framework of Qilian orogen: reveal from an aeromagnetic anomaly feature

The Qilian Orogenic belt is one of the typical orogenic belts globally and a natural laboratory for studying plate tectonics.Many researchers have studied the ophiolite and high pressure and ultra-high pressure metamorphic rocks in the Qilian orogen and obtained valuable achievements.However,a hot debate exists on the basement property,the distribution of ophiolite,and the boundaries of tectonic units.Large-scale high-precision aeromagnetic surveys have recently been conducted in the Qilian Orogenic belt and adjacent areas.In this study,we are trying to analysis the tectonic framework of the Qilian Orogen using 1:500,000 aeromagnetic data.The results provide geophysical perspectives for studying the structural framework and deformation of this area.According to the aeromagnetic∆T anomaly map,the central and Southern Qilian have the same magnetic anomaly feature that noticeably differs from the North Qilian Orogenic belt and the Qaidam Block.This result indicates that the central and Southern Qilian have a unified magnetic basement and differ from the North Qilian orogenic belt and Qaidam Block.The map shows the distribution of ophiolite in the North Qilian orogenic belt.Linear magnetic anomalies represent the ophiolites because the mafic–ultramafic rocks usually have high magnetic susceptibility.The ophiolite belts are continuously distributed in the western part of North Qilian orogenic belt and have a large scale.However,the scale of the ophiolite belt and the outcropping of mafic–ultramafic rocks reduces when they pass through Qilian County to the east.The results indicate differences in the evolution process between the eastern and western parts of North Qilian,with Qilian County as the transition zone.This study also systematically defines the geophysical boundaries of the Qaidam Block,Qilian Block,North Qilian Orogenic belt,and Alxa block.It is proposed that the sinistral displacement of the Altun Fault is adjusted and absorbed by the series of NE-trending faults in the Qilian orogen and merge into the Longshoushan–Gushi Fault.The extension of the North Qilian Orogenic belt is strengthened by the neotectonics movement along the shearing direction,which separated the North Qilian Orogenic belt into several segments and formed a series of northeast-trending faults.

A Hybrid Dung Beetle Optimization Algorithm with Simulated Annealing for the Numerical Modeling of Asymmetric Wave Equations

In the generalized continuum mechanics(GCM)theory framework,asymmetric wave equations encompass the characteristic scale parameters of the medium,accounting for microstructure interactions.This study integrates two theoretical branches of the GCM,the modified couple stress theory(M-CST)and the one-parameter second-strain-gradient theory,to form a novel asymmetric wave equation in a unified framework.Numerical modeling of the asymmetric wave equation in a unified framework accurately describes subsurface structures with vital implications for subsequent seismic wave inversion and imaging endeavors.However,employing finite-difference(FD)methods for numerical modeling may introduce numerical dispersion,adversely affecting the accuracy of numerical modeling.The design of an optimal FD operator is crucial for enhancing the accuracy of numerical modeling and emphasizing the scale effects.Therefore,this study devises a hybrid scheme called the dung beetle optimization(DBO)algorithm with a simulated annealing(SA)algorithm,denoted as the SA-based hybrid DBO(SDBO)algorithm.An FD operator optimization method under the SDBO algorithm was developed and applied to the numerical modeling of asymmetric wave equations in a unified framework.Integrating the DBO and SA algorithms mitigates the risk of convergence to a local extreme.The numerical dispersion outcomes underscore that the proposed SDBO algorithm yields FD operators with precision errors constrained to 0.5‱while encompassing a broader spectrum coverage.This result confirms the efficacy of the SDBO algorithm.Ultimately,the numerical modeling results demonstrate that the new FD method based on the SDBO algorithm effectively suppresses numerical dispersion and enhances the accuracy of elastic wave numerical modeling,thereby accentuating scale effects.This result is significant for extracting wavefield perturbations induced by complex microstructures in the medium and the analysis of scale effects.

生物岩石学的定义

生物岩石学是研究生物岩(即生物成因岩石)的特征、形成机制、形成环境及其与矿产资源关系的一门新兴交叉学科,其研究内容至少包括生物矿化、现代生物礁、古代生物礁、现代微生物岩石和古代微生物岩石5个方面。生物矿化作用和现代生物岩研究为古代生物岩研究提供了认识基础。由于现代生物圈和环境不同于古代,现代生物矿化作用和生物岩的研究成果并不能全部直接应用于古代生物岩研究。古代生物岩和生物矿化作用的类型比现代丰富得多,不可能全部从现代生物岩和生物矿化作用中找到参照,但可以为现代生物矿化实验研究提供设计思路。生物岩石学的相关学科包括生物学、微生物学、古生物学、古微生物学、沉积学、沉积岩石学、矿物学、地球化学、地质微生物学等。生物岩石学研究需要应用这些学科的知识,也会反哺这些学科。

浅谈大气科学与地质学的学科交叉

随着新一轮的科技革命蓬勃兴起,大气科学学科正步入地球系统科学的新时代,学科交叉必然产生新的增长点.大气科学的发展经历了观测-理论-模型的三个阶段,很好地践行了“数据-模式驱使科学”的研究范式.然而,地质学具有更深远的时空复杂性,需要更长时间的数据积累.目前,地质学正面临着研究范式由观测向理论和模型的转变.大气科学与地质学的交叉将为这一转变提供经验和启示.同时,大气科学与地质学的交叉,需要研究气候系统上边界(大气顶的太阳辐射)和下边界(固体地球形态)变化导致的大气和海洋环流的响应.研究这些问题,将成为发展同时统辖“分钟、小时”直至“地史”的时间尺度的大气科学理论的关键,也是未来地球系统模式发展的重要方向之一.中国地质大学(武汉)的大气科学专业,作为推动大气科学与地质学交叉的排头兵,任重而道远.

Continuous TDEM for monitoring shale hydraulic fracturing

Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic methods, are expensive, time- consuming, and do not show the changes in the formation with time. The resistivities of hydraulic fracturing fluid and reservoir rocks were measured. The results suggest that the injection fluid and consequently the injected reservoir are characterized by very low resistivity and high chargeability. This allows using of the controlled-source electromagnetic method （CSEM） to monitor shale gas hydraulic fracturing. Based on the geoelectrical model which was proposed according to the well-log and seismic data in the test area the change rule of the reacted electrical field was studied to account for the change of shale resistivity, and then the normalized residual resistivity method for time lapse processing was given. The time-domain electromagnetic method （TDEM） was used to continuously monitor the shale gas fracturing at the Fulin shale gas field in southern China. A high-power transmitter and multi-channel transient electromagnetic receiver array were adopted. 9 h time series of Ex component of 224 sites which were laid out on the surface and over three fracturing stages of a horizontal well at 2800 m depth was recorded. After data processing and calculation of the normalized resistivity residuals, the changes in the Ex signal were determined and a dynamic 3D image of the change in resistivity was constructed. This allows modeling the spatial distribution of the fracturing fluid. The model results suggest that TDEM is promising for monitoring hydraulic fracturing of shale.

Feasibility of central loop TEM method for prospecting multilayer water-fi lled goaf

With deep mining of coal mines, prospecting multilayer water-filled goaf has become a new content that results from geophysical exploration in coalfields. The central loop transient electromagnetic （TEM） method is favorable for prospecting conductive layers because of the coupling relationship between its field structure and formation. However, the shielding effect of conductive overburden would not only require a longer observation time when prospecting the same depth but also weaken the anomalous response of underlying layers. Through direct time domain numerical simulation and horizontal layered earth forward modeling, this paper estimates the length of observation time required to prospect the target, and the distinguishable criterion of multilayer water-filled goal is presented with observation error according to the effect of noise on observation data. The observed emf curves from Dazigou Coal Mine, Shanxi Province can distinguish multilayer water-filled goaf. In quantitative inversion interpretation of observed curves, using electric logging data as initial parameters restrains the equivalence caused by coal formation thin layers. The deduced three-layer and two-layer water-filled goals are confirmed by the drilling hole. The result suggests that when observation time is long enough and with the anomalous situation of underlying layers being greater than the observation error, the use of the central loop TEM method to orosoect a multilaver water-filled goaf is feasible.

济阳孔店转型期识别标志及盆地特征

济阳地区孔店组与上覆地层的分界线,一直是困扰着孔店组构造层序石油地质研究与勘探的关键问题,笔者通过区域地震大剖面的研究,首次从地震剖面找到了孔店组顶、底不整合界面,为研究孔店组原型盆地奠定了基础.研究发现东营凹陷孔店组盆地与沙河街组盆地同受北部控盆断层控制,控盆断层东部为北西向,向西转为近东西向,东部可见小角度的角度不整合;而惠民凹陷孔店组盆地受北掉的南部北东向断层系控制,形成南断北超的半地堑,与上覆沙河街组北断南超半地堑叠置,孔店组与沙四段之间有明显跷跷板现象,沙四段向南翘倾作用使得孔店组根部地层大幅度剥蚀,产生孔店组与上覆地层之间明显角度不整合.东营与惠民凹陷的孔店组盆地通过林樊家～青城构造调节带衔接,济阳地区北部沾化凹陷、车镇凹陷孔店期为零星小盆地,孔店组的油气勘探重点应放在东营与惠民,由于东营孔店组盆地所受的改造相对较少,应为首选目标.

Boundary-refl ected waves and ultrasonic coda waves in rock physics experiments

Ultrasonic coda waves are widely usea to stuay hign-trequency scattering, however, ultrasonic coda waves are strongly affected by interference from by boundary-reflected waves. To understand the effect of boundary-reflected waves, we performed ultrasonic experiments using aluminum and shale samples, and the rotating staggered-mesh finite-difference method to simulate the wavefield. We analyzed the wavefield characteristics at the different receiving points and the interference characteristics of the boundary-reflected waves with the ultrasonic coda wave, and the effect of sample geometry on the ultrasonic coda waves. The increase in the aspect ratio of the samples delays the interference effect of the laterally reflected waves and reduces the effect on the ultrasonic coda waves. The main waves interfering with the ultrasonic coda waves are laterally reflected PP-, PS-, PPP-, and PPS-waves. The scattering and attenuation of the high-frequency energy in actual rocks can weaken the interference of laterally reflected waves with the ultrasonic coda waves.

A Bregman adaptive sparse-spike deconvolution method in the frequency domain

To improve the anti-noise performance of the time-domain Bregman iterative algorithm,an adaptive frequency-domain Bregman sparse-spike deconvolution algorithm is proposed.By solving the Bregman algorithm in the frequency domain,the influence of Gaussian as well as outlier noise on the convergence of the algorithm is effectively avoided.In other words,the proposed algorithm avoids data noise effects by implementing the calculations in the frequency domain.Moreover,the computational efficiency is greatly improved compared with the conventional method.Generalized cross validation is introduced in the solving process to optimize the regularization parameter and thus the algorithm is equipped with strong self-adaptation.Different theoretical models are built and solved using the algorithms in both time and frequency domains.Finally,the proposed and the conventional methods are both used to process actual seismic data.The comparison of the results confirms the superiority of the proposed algorithm due to its noise resistance and self-adaptation capability.

Normal moveout for long offset in isotropic media using the Padé approximation

The nomaal moveout correction is important to long-offset observations, especially deep layers. For isotropic media, the conventional two-term approximation of the normal moveout function assumes a small offset-to-depth ratio and thus fails at large offset-to-depth ratios. We approximate the long-offset moveout using the Pade approximation. This method is superior to typical methods and flattens the seismic gathers over a wide range of offsets in multilayered media. For a four-layer model, traditional methods show traveltime errors of about 5 ms for offset-to-depth ratio of 2 and greater than 10 ms for offset-to-depth ratio of 3; in contrast, the maximum traveltime error for the [3, 3]-order Pade approximation is no more than 5 ms at offset-to-depth ratio of 3. For the Cooper Basin model, the maximum oft＇set-to-depth ratio for the [3, 3]-order Pade approximation is typically double of those in typical methods. The [7, 7]-order Pade approximation performs better than the [3.3]-order Pade armroximation.

Fast and accurate calculation of seismic wave travel time in 3D TTI media

The Tilted tilted transversely isotropic(TTI)media,a kind of anisotropic medium,widely exists within the earth.For faster calculation of travel times in the TTI anisotropic media,we modifi ed a minimum traveltime tree algorithm with high effi ciency by dynamical modifi cation of the secondary wave propagation region during the spread of seismic waves.To manage the wavefront points in the modified version,we used a novel minimum heap sorting technique to reduce the time spent on selecting secondary waves points.In this study,seismic group velocities were obtained from analytical solutions in terms of phase angle,and the corresponding phase angles were determined by binary search rather than approximate equations for weakly anisotropic media.For the most time-consuming part of the secondary wave traveltime calculation,the parallel computation was initially performed using multiple cores and threads.Numerical examples showed that the improved method can calculate seismic travel times and ray paths faster and accurately in a 3D TTI medium.For four cores and eight threads,the computing speed increased by six times when compared to the conventional method.

Removal of the airwave effect by main-part decomposition of the anomalous field of MCSEM data

The airwave effect greatly influences the observational data from controlledsource electromagnetic exploration in shallow seas, which obscures the abnormal effects generated by exploration targets and, hence, affects the accuracy of the late exploration data interpretation. In this study, we propose a method to separate the main part from the anomalous field of marine controlled-source electromagnetic method （MCSEM） data based on Stratton-Chu integral transforms to eliminate the airwave effect, which dominates observed electromagnetic （EM） response in shallow seawater. This method of separating the main part from the anomalous field is a type of finite impulse response filter based on a discrete data set. Theoretical analysis proved that the method is stable and able to effectively depress noise. A numerical test indicated that the method could successfully eliminate the airwave effect from the observed EM signals generated by an air water interface and a seawater layer. This technique is applicable for seawater models with either flat or rough seabeds.

Super-resolution least-squares prestack Kirchhoff depth migration using the L_0-norm

Least-squares migration （LSM） is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L2-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L0-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.