Bayesian prestack seismic inversion with a self-adaptive Huber-Markov random-field edge protection scheme

Seismic inversion is a highly ill-posed problem, due to many factors such as the limited seismic frequency bandwidth and inappropriate forward modeling. To obtain a unique solution, some smoothing constraints, e.g., the Tikhonov regularization are usually applied. The Tikhonov method can maintain a global smooth solution, but cause a fuzzy structure edge. In this paper we use Huber-Markov random-field edge protection method in the procedure of inverting three parameters, P-velocity, S-velocity and density. The method can avoid blurring the structure edge and resist noise. For the parameter to be inverted, the Huber- Markov random-field constructs a neighborhood system, which further acts as the vertical and lateral constraints. We use a quadratic Huber edge penalty function within the layer to suppress noise and a linear one on the edges to avoid a fuzzy result. The effectiveness of our method is proved by inverting the synthetic data without and with noises. The relationship between the adopted constraints and the inversion results is analyzed as well.     

Evolution of isotopic composition and deuterium excess of brines in the Sichuan Basin, China

The isotopic composition and parameters for deuterium excess of brines, which were sampled in the Si-chuan Basin, show obvious regularities of distribution. The brine isotopic composition shows distinct two systems of marine and terrestrial deposits, with the Middle Triassic strata as the boundary. Brine hydrogen isotopic composition of marine deposits is lower while oxygen isotopic composition is higher than that of the SMOW, respectively, indicating that the brines were derived from seawater with different evaporating degrees at different times. From the Sinian strata, up to the Cambrian, Permian Maokou Formation and the Triassic Jialingjiang Formation, the δD values of brines tend to become relatively positive with the strata becoming younger. Brines of terrestrial deposits are considered to have been derived from precipitation and their isotopic composition is close to the globe meteoric water line (GMWL). Brines of transitional deposits between marine and terrestrial ones (the Upper Triassic Xujiahe Formation) have δD and δ18O values falling between the two end members of marine deposit brines and precipitation, indicating that the brines are a mixture of precipitation and vaporing seawater. Water samples from the brine-bearing strata of different ages show various deuterium excesses (d) with an evident decreasing trend as the age of strata gets older and older. Brine-bearing strata of the Triassic Leikoupo-Jialingjiang Formation, the Permian Maokou Formation, the Cambrian and Sinian strata are all carbonate rocks which have experienced intensive water/rock reaction and the deuterium excess essentially changes with time. All brine-bearing-strata surrounding the basin or faults, as well as those brine wells exploited for resources, have been obviously influenced by the precipitation supply. Therefore, the deuterium excesses of their brines have increased to different extents, depending on the amount of involvement of meteoric water. The variation and distribution of d values of the brines from different Triassic strata are related to the embedded depth of the strata. The deuterium excesses of brines become lower with increasing burial depth of the strata.     

Seasonal variability of zonal heat advection in the mixed layer of the tropical Pacific

Zonal heat advection （ZHA） plays an important role in the variability of the thermal structure in the tropical Pacific Ocean, especially in the western Pacific warm pool （WPWP）. Using the Simple Ocean Data Assimilation （SODA） Version 2.02/4 for the period 1958-2007, this paper presents a detailed analysis of the climatological and seasonal ZHA in the tropical Pacific Ocean. Climatologically, ZHA shows a zonal- band spatial pattern associated with equatorial currents and contributes to forming the irregular eastern boundary of the WPWP （EBWP）. Seasonal variation of ZHA with a positive peak from February to July is most prominent in the Nifio3.4 region, where the EBWP is located. The physical mechanism of the seasonal cycle in this region is examined. The mean advection of anomalous temperature, anomalous advection of mean temperature and eddy advection account for 31%, 51%, and 18% of the total seasonal variations, respectively. This suggests that seasonal changes of the South Equatorial Current induced by variability of the trade winds are the dominant contributor to the anomalous advection of mean temperature and hence, the seasonality of ZHA. Heat budget analysis shows that ZHA and surface heat flux make comparable contributions to the seasonal heat variation in the Nifio3.4 region, and that ZHA cools the upper ocean throughout the calendar year except in late boreal spring. The connection between ZHA and EBWP is further explored and a statistical relationship between EBWP, ZHA and surface heat flux is established based on least squares fitting.     

华南近岸海域水质评价及其影响因素

以华南沿海32个地市为研究对象,首先基于单一要素评价和综合水质指数法对其近岸海域的水环境质量进行分析,再采用典范对应分析和逐步线性回归探究了各类影响因素与水质指标间的相关性.结果表明:2017年华南近岸海域的主要污染要素为无机氮(N)、活性磷酸盐(P)和悬浮物(SS),严重超标现象主要出现在福建厦门、宁德、福州,广西钦...     

Numerical well model for non-Darcy flow through isotropic porous media[*]This work was supported in part by the EPA under grant # R 825207-01-1, by the State of Texas under ARP/ATP grant # 010366-168 and by a gift grant from Mobil Oil Corp.

Numerical simulation of fluid flow in a hydrocarbon reservoir has to account for the presence of wells. The pressure of a grid cell containing a well is different from the average pressure in that cell and different from the bottom-hole pressure for the well [17]. This paper presents a study of grid pressures obtained from the simulation of single phase flow through an isotropic porous medium using different numerical methods. Well equations are proposed for Darcy flow with Galerkin finite elements and mixed finite elements. Furthermore, high velocity (non-Darcy) flow well equations are developed for cell-centered finite difference, Galerkin finite element and mixed finite element techniques.     

Pulse-length-tolerant features and detectors for sector-scan sonar imagery

This paper presents a neural-network-based system to detect small man-made objects in sequences of sector-scan sonar images created using signals of various pulse lengths. The detection of such objects is considered out to ranges of 150 m by using an experimental sector-scan sonar system mounted on a vessel. The sonar system considered in this investigation has three modes of operation to create images over ranges of 200, 400, and 800 m from the vessel using acoustic pulses of a different duration for each mode. After an initial cleaning operation performed by compensating for the motion of the vessel, the imagery is segmented to extract objects for analysis. A set of 31 features extracted from each object is examined. These features consist of basic object size and contrast features, shape moment-based features, moment invariants, and features extracted from the second-order histogram of each object. Optimal sets of 15 features are then selected for each mode and over all modes using sequential forward selection (SFS) and sequential backward selection (SBS). These features are then used to train neural networks to detect man-made objects in each sonar mode. By the addition of a feature describing the sonar's mode of operation, a neural network is trained to detect man-made objects in any of the three sonar modes. The multimode detector is shown to perform very well when compared with detectors trained specifically for each sonar mode setting. The proposed detector is also shown to perform well when compared to a number of statistical detectors based on the same set of features. The proposed detector achieves a 92.4% probability of detection at a mean false-alarm rate of 10 per image, averaged over all sonar mode settings.     

雷击大地密度Ng数值研究

雷击大地密度Ng作为一个重要的气象参数,直接影响雷电灾害风险评估的科学计算。Ng数值的获取有两种方式:一种是由地闪定位网络系统提供,一种是按Ng≈0.1T_d(年均雷暴日)估算。在利用后一种方法获取Ng值时,简单地取0.1T_d并非科学。因为T_d是较大区域范围的平均值,未必是雷电灾害风险评估项目所处区域的真实值;况且Ng值不仅受T_d的影响,还与雷暴时间、雷暴路径、闪电特征、雷灾频次与强度有关。因此,利用后一种方法获取Ng值时,需计入其他气象因素的影响。基于对影响Ng的各种气象因素的分析,将焦作市按行政区划定为65个区域,根据各个区域的气象环境,分析出各个区域的雷电灾害风险等级;根据风险等级,设定风险系数,并结合周边区域气象因素对该区域的影响,计算出Ng的修正系数,修正Ng的取值,使修正后的Ng值更接近真实值。     

利用粉晶衍射数据精测点阵参数计算程序

本文着重讨论了利用粉晶数据进行晶体点阵参数的计算问题(PDOLC)。计算程序采用外推最小二乘拟合技术,外推函数任选,作了详细的系统消光处理,分别估计了六个点阵参数的精度,适用于用各种粉晶方法获得的衍射数据处理。     

古田水口库区ML4.8地震成因探讨

介绍水口库区地震地质构造特征与地震活动特点,重点分析2008年3月6日古田4.8级地震序列以及震时福建地区应力场特征。分析认为：1) 2008年3月水口库区处于明显的应力积累区域;2）古田4.8级地震为构造型水库地震, NE向湾口断裂（F2）是可能的发震构造之一;3）水口水库前两次4级地震发震时间相隔近12年,表明水库水位的升降与发震之间的相关性并不明显,区域构造应力的变化是导致水口库区地震的决定因素。