Quantitative characterization of tight gas sandstone reservoirs using seismic data via an integrated rock-physics-based framework

Seismic characterizing of tight gas sandstone (TGS) reservoirs is essential for identifying promising gas-bearing regions. However, exploring the petrophysical significance of seismic-inverted elastic properties is challenging due to the complex microstructures in TGSs. Meanwhile, interbedded structures of sandstone and mudstone intensify the difficulty in accurately extracting the crucial tight sandstone properties. An integrated rock-physics-based framework is proposed to estimate the reservoir quality of TGSs from seismic data. TGSs with complex pore structures are modeled using the double-porosity model, providing a practical tool to compute rock physics templates for reservoir parameter estimation. The VP/VS ratio is utilized to predict the cumulative thickness of the TGS reservoirs within the target range via the threshold value evaluated from wireline logs for lithology discrimination. This approach also facilitates better capturing the elastic properties of the TGSs for quantitative seismic interpretation. Total porosity is estimated from P-wave impedance using the correlation obtained based on wireline log analysis. After that, the three-dimensional rock-physics templates integrated with the estimated total porosity are constructed to interpret microfracture porosity and gas saturation from velocity ratio and bulk modulus. The integrated framework can optimally estimate the parameters dominating the reservoir quality. The results of the indicator proposed based on the obtained parameters are in good agreement with the gas productions and can be utilized to predict promising TGS reservoirs. Moreover, the results suggest that considering microfracture porosity allows a more accurate prediction of high-quality reservoirs, further validating the applicability of the proposed method in the studied region.

Correlations between trace elements in pyrite and gold mineralization of gold deposits on the North China platform

By studying both the microscopic physical and chemical typomorphic characteristics of typical mineral pyrite samples associated with representative gold deposits on the north-central margin of the North China Platform,this paper seeks to identify macroscopic metallogenic mechanisms of gold deposits and to reveal the formation mechanism of lattice gold in pyrite.Typomorphic characteristics of pyrite reveal that pyrite grain size has a negative correlation with gold content.Cubic pyrite,as the dominant crystal form,contains more gold than pentagonal dodecahedral pyrite.Both pyrite crystal forms and chemical compositions indicate that the replacement style of gold deposit formed in a low saturability,low sulfur fugacity,and at temperatures either much higher or much lower than its best forming temperature;comparatively,that of the quartz vein style of gold deposit occurred under conditions with the best temperature,rich in sulfur,and with high sulfur fugacity.The Au/Ag ratios of the pyrites show that both the replacement and quartz vein styles of deposits are mesothermal and hypothermal,while the Co/Ni ratios of the pyrites indicate that the quartz vein style is of magmatic-hydrothermal origin.The X-ray diffraction intensity of pyrite rich in gold is lower than that of pyrite poor in gold at the quartz vein style.In general,with an increase in gold content in pyrite,the total sum intensityΣI decreases.The pyroelectricity coefficient has a negative correlation trend with the values of(Co+Ni+Se+Te)-As and(Co+Ni+Se+Te)/As.The pyrite pyroelectricity of the replacement style is N-type,indicating that it formed under low sulfur fugacity,while that of the quartz vein style is a mixture of P-N types,indicating that it formed under high sulfur fugacity.On the pyroelectricity-temperature diagram,pyrite of the replacement style is mainly distributed between 200 and 270°C,while that of the quartz vein style varies between 90–118 and274–386°C,demonstrating a multistage forming process.In contrast to previous researchers'conclusions,the authors confirm the existence of lattice gold in pyrites through the use of an electron paramagnetic resonance(EPR)test.Au in the form of Au~+,entering pyrite as an isomorph and producing electron–hole centers,makes the centers produce spin resonance absorption and results in EPR absorption peak II.The intensity of auriferous pyrite absorption peak II has certain direct positive correlations with pyrite gold content.The#I and#III absorption peaks of pyrites possibly result from the existence of Ni^(2+)and/or Cu^(2+).γ1,γ2,andγ3 are the strongest and most typical absorption peaks of the infrared spectra of the pyrites.Generally,with the increase in gold content in the pyrite samples,γ1,γ2,andγ3 tend to shift to higher wavenumbers,and the gold content in the pyrite samples has a positive correlation with their relative absorbance.

Research on the generation method of seawater sound velocity model based on Perlin noise

In the processing of conventional marine seismic data,seawater is often assumed to have a constant velocity model.However,due to static pressure,temperature difference and other factors,random disturbances may often frequently in seawater bodies.The impact of such disturbances on data processing results is a topic of theoretical research.Since seawater sound velocity is a difficult physical quantity to measure,there is a need for a method that can generate models conforming to seawater characteristics.This article will combine the Munk model and Perlin noise to propose a two-dimensional dynamic seawater sound velocity model generation method,a method that can generate a dynamic,continuous,random seawater sound velocity model with some regularity at large scales.Moreover,the paper discusses the influence of the inhomogeneity characteristics of seawater on wave field propagation and imaging.The results show that the seawater sound velocity model with random disturbance will have a significant influence on the wave field simulation and imaging results.

Influence of Temperature and Frequency on Microwave Dielectric Properties of Lunar Regolith Simulant

The dielectric constant of the lunar regolith can directly influence the reflection coefficient and the trans-mission coefficient of the Moon′s surface, and plays an important role in the Moon research. In order to study the di-electric properties of the lunar regolith, the lunar regolith simulant was made according to the making procedure of the CAS-1 simulant made by Chinese Academy of Sciences. Then the dielectric constants of the lunar regolith simulant were measured with 85070E Aiglent Microwave Network Analyzer in the frequency ranging from 0.2 GHz to 20.0 GHz and at temperature of 25.1℃, 17.7℃, 13.1℃, 11.5℃, 9.6℃, 8.0℃, 4.1℃, -0.3℃, -4.7℃, -9.5℃, -18.7℃, -27.7℃, and -32.6℃, respectively. The Odelevsky model was employed to remove the influence of water in the air on the final effective dielectric constants. The results indicate that frequency and temperature have apparent influences on the dielectric constants of the lunar regolith simulant. The real parts of the dielectric constants increase fast over the range of 0.2 GHz to 3.0 GHz, but decrease slowly over the range of 4.0 GHz to 20.0 GHz. The opposite phenomenon occurs in the imaginary parts. The influences of the frequency and temperature on the brightness temperature were also estimated based on the radiative transfer equation. The result shows that the variation of the frequency and temperature results in great changes of the microwave brightness temperature emitting from the lunar regolith.

Acousto-electric well logging by eccentric source and extraction of shear wave

The nonaxisymmetric acousto-electric field excited by an eccentric acoustic source in the borehole based on Pride seismoelectric theory is considered. It is shown that the acoustic field inside the borehole, converted electric and magnetic fields and coupled fields outside the borehole are composed of an infinitude of multipole fields with different orders. The numerical results show that both the electromagnetic waves and the seismoelectric field in the borehole, and the three components of both electric field and magnetic field can be detected. Measurements on the borehole axis will be of advantage to determining shear velocity information. The components of the symmetric and nonsymmetric acoustic and electromagnetic fields can be strengthened or weakened by adding or subtracting the two full waveforms logged in some azimuths. It may be a new method of directly measuring the shear wave velocity by using the borehole seismoelectric effect.

Retrieving volume FeO and TiO_2 abundances of lunar regolith with CE-2 CELMS data using BPNN method

The volume FeO and TiO_2 abundances(FTAs) of lunar regolith can be more important for understanding the geological evolution of the Moon compared to the optical and gamma-ray results. In this paper, the volume FTAs are retrieved with microwave sounder(CELMS) data from the Chang'E-2 satellite using the back propagation neural network(BPNN) method. Firstly, a three-layered BPNN network with five-dimensional input is constructed by taking nonlinearity into account. Then, the brightness temperature(TB) and surface slope are set as the inputs and the volume FTAs are set as the outputs of the BPNN network.Thereafter, the BPNN network is trained with the corresponding parameters collected from Apollo, Luna,and Surveyor missions. Finally, the volume FTAs are retrieved with the trained BPNN network using the four-channel TBderived from the CELMS data and the surface slope estimated from Lunar Orbiter Laser Altimeter(LOLA) data. The rationality of the retrieved FTAs is verified by comparing with the Clementine UV-VIS results and Lunar Prospector(LP) GRS results. The retrieved volume FTAs enable us to re-evaluate the geological features of the lunar surface. Several important results are as follows. Firstly, very-low-Ti(<1.5 wt.%) basalts are the most spatially abundant, and the surfaces with TiO_2> 5 wt.% constitute less than 10% of the maria. Also, two linear relationships occur between the FeO abundance(FA) and the TiO_2 abundance before and after the threshold, 16 wt.% for FA. Secondly, a new perspective on mare volcanism is derived with the volume FTAs in several important mare basins, although this conclusion should be verified with more sources of data. Thirdly, FTAs in the lunar regolith change with depth to the uppermost surface,and the change is complex over the lunar surface. Finally, the distribution of volume FTAs hints that the highlands crust is probably homogeneous, at least in terms of the microwave thermophysical parameters.

Review on retrieval of lunar regolith thickness by active and passive microwave measurements

It is one of the important methods to retrieve lunar regolith thickness using active and passive microwave techniques.The retrieval of lunar regolith thickness is based on microwave radiation transfer process simulation in the regolith media.The lunar regolith model is first introduced,and the features of the involved physical parameters are indicated thereafter,such as dielectric constants,surface roughness,particle size and thermal grads of the lunar regolith.The time delay and the migration of the radar echoes from the different interfaces is the key problem for active microwave measurement.And the simulation of the microwave radiative transfer in the regolith media is the important technique for the passive microwave measurement.The important parameters and the physical mechanism for the two measurements are also presented.

Predicting Surface Roughness and Moisture of Bare Soils Using Multi- band Spectral Reflectance Under Field Conditions

Soil surface roughness, denoted by the root mean square height(RMSH), and soil moisture(SM) are critical factors that affect the accuracy of quantitative remote sensing research due to their combined influence on spectral reflectance(SR). In regards to this issue, three SM levels and four RMSH levels were artificially designed in this study; a total of 12 plots was used, each plot had a size of 3 m × 3 m. Eight spectral observations were conducted from 14 to 30 October 2017 to investigate the correlation between RMSH, SM, and SR. On this basis, 6 commonly used bands of optical satellite sensors were selected in this study, which are red(675 nm), green(555 nm), blue(485 nm), near infrared(845 nm), shortwave infrared 1(1600 nm), and shortwave infrared 2(2200 nm). A negative correlation was found between SR and RMSH, and between SR and SM. The bands with higher coefficient of determination R^2 values were selected for stepwise multiple nonlinear regression analysis. Four characterized bands(i.e., blue, green, near infrared, and shortwave infrared 2) were chosen as the independent variables to estimate SM with R^2 and root mean square error(RMSE) values equal to 0.62 and 2.6%, respectively. Similarly, the four bands(green, red, near infrared, and shortwave infrared 1) were used to estimate RMSH with R^2 and RMSE values equal to 0.48 and 0.69 cm, respectively. These results indicate that the method used is not only suitable for estimating SM but can also be extended to the prediction of RMSH. Finally, the evaluation approach presented in this paper highly restores the real situation of the natural farmland surface on the one hand, and obtains high precision values of SM and RMSH on the other. The method can be further applied to the prediction of farmland SM and RMSH based on satellite and unmanned aerial vehicle(UAV) optical imagery.

In situ electrical resistance and activation energy of solid C_(60) under high pressure

The in situ electrical resistance and transport activation energies of solid C60 fullerene have been measured under high pressure up to 25 GPa in the temperature range of 300-423 K by using a designed diamond anvil cell. In the experiment, four parts of boron-doped diamond films fabricated on one anvil were used as electrical measurement probes and a W-Ta thin film thermocouple which was integrated on the other diamond anvil was used to measure the temperature. The current results indicate that the measured high-pressure resistances are bigger than those reported before at the same pressure and there is no pressure-independent resistance increase before 8 GPa. From the temperature dependence of the resistivity, the C60 behaviors as a semiconductor and the activation energies of the cubic C60 fullerene are 0,49, 0.43, and 0.36 eV at 13, 15, and 19 GPa, respectively.

Necessary Maximum Principle of Stochastic Optimal Control with Delay and Jump Diffusion

In this paper, we have studied the necessary maximum principle of stochastic optimal control problem with delay and jump diffusion.

Methods of identifying gas-bearing reservoirs by logs and applications of volcanic formations

The authors generalized the methods how to identify and evaluate gas zones using logs,and put forward the methods of crossplots and overlays of porosity logs on the identification of volcanic gas-bearing reservoirs in the northern Songliao Basin with good results.This study provides technical clues in deep formations and offers references for other areas to identify oil and gas layers.

A Cloud Framework for High Spatial Resolution Soil Moisture Mapping from Radar and Optical Satellite Imageries

Soil moisture plays an important role in crop yield estimation,irrigation management,etc.Remote sensing technology has potential for large-scale and high spatial soil moisture mapping.However,offline remote sensing data processing is time-consuming and resource-intensive,and significantly hampers the efficiency and timeliness of soil moisture mapping.Due to the high-speed computing capabilities of remote sensing cloud platforms,a High Spatial Resolution Soil Moisture Estimation Framework(HSRSMEF)based on the Google Earth Engine(GEE)platform was developed in this study.The functions of the HSRSMEF include research area and input datasets customization,radar speckle noise filtering,optical-radar image spatio-temporal matching,soil moisture retrieving,soil moisture visualization and exporting.This paper tested the performance of HSRSMEF by combining Sentinel-1,Sentinel-2 images and insitu soil moisture data in the central farmland area of Jilin Province,China.Reconstructed Normalized Difference Vegetation Index(NDVI)based on the Savitzky-Golay algorithm conforms to the crop growth cycle,and its correlation with the original NDVI is about 0.99(P<0.001).The soil moisture accuracy of the random forest model(R 2=0.942,RMSE=0.013 m3/m3)is better than that of the water cloud model(R 2=0.334,RMSE=0.091 m3/m3).HSRSMEF transfers time-consuming offline operations to cloud computing platforms,achieving rapid and simplified high spatial resolution soil moisture mapping.

Gravity compression forward modeling and multiscale inversion based on wavelet transform

在三维的严肃倒置的主要问题是解决方案的非唯一和大数据集合的高计算的费用。最小化高计算的费用，我们建议一个新排序方法在使用小浪变换以前减少变化和敏感矩阵的高频率。因而，敏感矩阵的稀少和压缩比率象向前建模的精确性一样被改进。而且，记忆存储要求被减少，向前建模与未压缩的前面的建模相比被加速。向前当模特儿的结果建议敏感矩阵的压缩比率能是超过 300。而且，多尺度的倒置基于小浪变换被用于严肃倒置。由把严肃倒置分解成不同规模的 subproblems，当多尺度的数据被考虑，严肃倒置的非唯一和稳定性被改进。最后，我们在模仿并且测量的数据上使用了常规集中的倒置和多尺度的倒置表明建议严肃倒置方法的有效性。

Soil and Vegetation Spectral Coupling Difference (SVSCD) for Minerals Extraction from Hyperion Data in Vegetation Covered Area

Remote sensing data have been widely applied to extract minerals in geologic exploration, however, in areas covered by vegetation, extracted mineral information has mostly been small targets bearing little information. In this paper, we present a new method for mineral extraction aimed at solving the difficulty of mineral identification in vegetation covered areas. The method selected six sets of spectral difference coupling between soil and plant(SVSCD). These sets have the same vegetation spectra reflectance and a maximum different reflectance of soil and mineral spectra from Hyperion image based on spectral reflectance characteristics of measured spectra. The central wavelengths of the six, selected band pairs were 2314 and 701 nm, 1699 and 721 nm, 1336 and 742 nm, 2203 and 681 nm, 2183 and 671 nm, and 2072 and 548 nm. Each data set's reflectance was used to calculate the difference value. After band difference calculation, vegetation information was suppressed and mineral abnormal information was enhanced compared to the scatter plot of original band. Six spectral difference couplings, after vegetation inhibition, were arranged in a new data set that requires two components that have the largest eigenvalue difference from principal component analysis(PCA). The spatial geometric structure features of PC1 and PC2 was used to identify altered minerals by spectral feature fitting(SFF). The collecting rocks from the 10 points that were selected in the concentration of mineral extraction were analyzed under a high-resolution microscope to identify metal minerals and nonmetallic minerals. Results indicated that the extracted minerals were well matched with the verified samples, especially with the sample 2, 4, 5 and 8. It demonstrated that the method can effectively detect altered minerals in vegetation covered area in Hyperion image.

Simulation of Electromagnetic Wave Logging Response in Deviated Wells Based on Vector Finite Element Method

The vector finite element method of tetrahedral elements is used to model 3D electromagnetic wave logging response. The tangential component of the vector field at the mesh edges is used as a degree of freedom to overcome the shortcomings of node-based finite element methods. The algorithm can simulate inhomogeneous media with arbitrary distribution of conductivity and magnetic permeability. The electromagnetic response of well logging tools are studied in dipping bed layers with the borehole and invasion included. In order to simulate realistic logging tools, we take the transmitter antennas consisting of circular wire loops instead of magnetic dipoles. We also investigate the apparent resistivity of inhomogeneous formation for different dip angles.

Fhe application of electromagnetic-induction on the measurement of sea ice thickness in the Antarctic

海冰作为冰雪圈中的重要组成部分,对气候变化十分敏感,海冰物理过程研究需要得到精确的海冰厚度。本文介绍一种能够高效探测海冰厚度的电磁感应方法及其在南极考察中的成功应用。该方法针对海冰和海水的电特性,利用电磁场原理精确探测仪器至冰水交界面的距离,以实现海冰厚度的测定。通过电磁感应(EM)仪观测的视电导率与同点位钻孔测量数据对比分析,获得视电导率与海冰厚度的转换关系式,通过对该关系式计算出的海冰厚度进行验证,表明电磁感应技术能够获得可靠的海冰厚度数据,相对平均误差仅为 5.5%。对探测剖面的海冰厚度统计表明,2005-2006年夏季南极内拉峡湾内海冰以平整冰为主,其厚度范围在0.8m到1.4m之间。

Time correction of the ocean bottom seismometers deployed at the southwest Indian ridge using ambient noise cross-correlation

Seismic monitoring using ocean bottom seismometers（OBS） is an efficient method for investigating earthquakes in mid-ocean ridge far away from land. Clock synchronization among the OBSs is difficult without direct communication because electromagnetic signals cannot propagate efficiently in water. Time correction can be estimated through global positioning system（GPS） synchronization if clock drift is linear before and after the deployment. However, some OBSs in the experiments at the southwest Indian ridge（SWIR） on the Chinese DY125-34 cruise had not been re-synchronized from GPS after recovery. So we attempted to estimate clock drift between each station pairs using time symmetry analysis（TSA） based on ambient noise cross-correlation. We tested the feasibility of the TSA method by analyzing daily noise cross-correlation functions（NCFs） that extract from the data of another OBS experiment on the Chinese DY125-40 cruise with known clock drift and the same deployment site. The results suggest that the NCFs＇ travel time of surface wave between any two stations are symmetrical and have an opposite growing direction with the date. The influence of different band-pass filters,different components and different normalized methods was discussed. The TSA method appeared to be optimal for the hydrophone data within the period band of 2–5 s in dozens of km-scale interstation distances. A significant clock drift of ~2 s was estimated between OBSs sets through linear regression during a 108-d deployment on the Chinese cruise DY125-34. Time correction of the OBS by the ambient noise cross-correlation was demonstrated as a practical approach with the appropriate parameters in case of no GPS re-synchronization.

Fast local wavenumber （FLW） method for the inversion of magnetic source parameters

为磁性的异例的解释的当前的本地 wavenumber 方法由解决复杂矩阵计算地质的身体的地点。一会儿，如此的方法要求知道结构的索引，它是代表来源类型的一个参数。结构的索引是难的在真实数据知道；因而，当前的方法的精确是低的。我们在场快本地 wavenumber (FLW ) 方法，并且把水平、垂直的本地 wavenumbers 的摆平的和定义为累积本地 wavenumber。FLW 方法是 umulative 本地人 wavenumberand 的线性联合另外的 wavenumbers，并且没有一个 priori 信息和矩阵答案，被用来计算地点和来源的结构的索引。我们把 FLW 方法用于合成磁性的异例，并且结果建议 FLW 方法对背景和倾斜的磁化感觉迟钝。下次，我们把 FLW 方法用于真实磁性的数据获得地点和来源的结构的索引。

Numerical Simulation of Shear-Horizontal-Wave-Induced Electromagnetic Field in Borehole Surrounded by Porous Formation

Seismoelectric fieM excited by purely torsional loading applied directly to the borehole wall is considered. A brief formulation and some computed waveforms show the advantage of using shear-horizontal （SH） transverseelectric （TE） seismoelectric waves logging to measure shear velocity in a fluid-saturated porous formation. By assuming that the acoustic field is not influenced by its induced electromagnetic field due to seismoeleetric effect, the coupling governing equations for electromagnetic field are reduced to Maxwell equations with a propagation current source. It is shown that this simplification is valid and the borehole seismoelectric conversion efficient is mainly dependent on the electrokinetic coupling coefficient. The receivers to detect the conversion electromagnetic field and to obtain shear velocity can be set in the borehole fluid in the SH-TE seismoelectric wave log.

Log evaluation of fractured igneous reservoirs in Songliao Basin

The petrophysical parameters are bad in most igneous reservoirs of Songliao Basin because of the very low porosity and permeability.The evaluation of igneous reservoirs has not been fully studied so far.The current technique of formation evaluation and interpretation used in sedimentary formations face a series of problems and difficulties.In this study,The PCA was used to identify lithology,a multi-mineral model "QAPM" was proposed."Surface effect" must be considered when evaluating saturation.A software "SIMPLE" was developed and was used to deal with the logging data in over 70 wells with good results were achieved.