“拒地”释义补说

韩愈《忆昨行和张十一》诗中的“拒地”是中古汉语习见词语,前人释义不一,清代学者陈景云的“同‘據地’,以脚據持地”说略近其实。细考“拒地”(又作“距地”)与“據地”,其含义有细微差别:“據地”自先秦以来承袭已久,侧重“用手按(據持)着地(跪坐着)”之义,主语一般都是手,唐代以前未见例外;而“拒(距)地”唐代始见,强调“用力抵(撑蹬)地”之义,施事主语可以是手,也可以是足或其他事物。宋代以后以“據地”混同于“拒地”,当属不规范表达。

Optimizing data for electrical resistivity tomography in hardened sites through the ratio method

To reduce the shielding effect of hardened layers on electrical resistivity tomography,a ratio method based on the distortion correction principle and the isolation coefficient is proposed.The effects of the resistivity and thickness of hardened concrete layers on the detection of target objects are explored.Both numerical simulations and indoor tank tests indicate that when the ratio method is employed to correct the original collected data,the maximum allowable error for the isolation coefficient should not exceed 1%.Notably,when the ratio of hardened layer thickness to electrode spacing does not exceed 1,correction through this method significantly enhances the recognition capability of target objects.However,when the hardened layer thickness is greater than the electrode spacing by a factor of 2 or more,the ratio method cannot achieve satisfactory results.The case study of flood control engineering detection in the Zhangxi section of the Huangpen River in Dongzhi County demonstrates that the detection effect after correction by the ratio method is comparable to that for the adjacent unhardened pavement,and the influence of the hardened layer is obviously weakened,resulting in more reliable results.

Irregular seismic data reconstruction based on exponential threshold model of POCS method

Irregular seismic data causes problems with multi-trace processing algorithms and degrades processing quality. We introduce the Projection onto Convex Sets （POCS） based image restoration method into the seismic data reconstruction field to interpolate irregularly missing traces. For entire dead traces, we transfer the POCS iteration reconstruction process from the time to frequency domain to save computational cost because forward and reverse Fourier time transforms are not needed. In each iteration, the selection threshold parameter is important for reconstruction efficiency. In this paper, we designed two types of threshold models to reconstruct irregularly missing seismic data. The experimental results show that an exponential threshold can greatly reduce iterations and improve reconstruction efficiency compared to a linear threshold for the same reconstruction result. We also analyze the anti- noise and anti-alias ability of the POCS reconstruction method. Finally, theoretical model tests and real data examples indicate that the proposed method is efficient and applicable.

The algorithm of 3D multi-scale volumetric curvature and its application

To fully extract and mine the multi-scale features of reservoirs and geologic structures in time/depth and space dimensions, a new 3D multi-scale volumetric curvature （MSVC） methodology is presented in this paper. We also propose a fast algorithm for computing 3D volumetric curvature. In comparison to conventional volumetric curvature attributes, its main improvements and key algorithms introduce multi-frequency components expansion in time-frequency domain and the corresponding multi-scale adaptive differential operator in the wavenumber domain, into the volumetric curvature calculation. This methodology can simultaneously depict seismic multi-scale features in both time and space. Additionally, we use data fusion of volumetric curvatures at various scales to take full advantage of the geologic features and anomalies extracted by curvature measurements at different scales. The 3D MSVC can highlight geologic anomalies and reduce noise at the same time. Thus, it improves the interpretation efficiency of curvature attributes analysis. The 3D MSVC is applied to both land and marine 3D seismic data. The results demonstrate that it can indicate the spatial distribution of reservoirs, detect faults and fracture zones, and identify their multi-scale properties.

Research on velocity and acceleration geophones and their acquired information

Based on the research of two geophone types （10 Hz moving-coil velocity and piezoelectric acceleration） and their velocity and acceleration parameters, frequency response functions have been determined, as well as the differences between them. Also, shock- vibration tests have been accomplished, not only to explain the two shock response signal differences, but also to analyze the response signal characteristics and its ability to carry information. In addition, seismic data acquisition experiments have been carried out under comparable conditions in the field. A contrast analysis of shot gathers and stack profiles acquired with the two geophone types is given in this paper. The results show that the acceleration signal from the acceleration geophone has a better advantage in terms of high signal-to-noise ratio, high accuracy, high resolution, and quantity of information to better meet current and future requirements for seismic exploration.

An Overview of Coilia ectenes in Jiangsu Section of the Yangtze River

Anchovies are migratory fish and adult individuals migrate in group from the sea into the river every year because of reproductive needs. TheYangtze River is the largest river in China, with rich food and fish species, hence, the Yangtze River would be an ideal habitat for growth and reproduction of anchovy. In this pa- per, based on fishing log data of anchovy in Jiangsu section of the Yangtze River, realted analysis was conducted in digital map using statistical software and GIS. From 2008 to 2011, the number of anchovy fishing vessels and fishing days in Jiangsu section remained stable, but the fishing volume and economic profit of single-vessel fluctuated constantly. This study provides references and advices for plan and management of anchovy fishery resource to administrative department.

Controlled illumination and seismic acquisition geometry for target-oriented imaging

The conventional method of seismic data acquisition geometry design is based on the assumption of horizontal subsurface reflectors, which often is not suitable for complex structure. We start from a controlled illumination analysis and put forward a method of seismic survey geometry design for target-oriented imaging. The method needs a velocity model obtained by a preliminary seismic interpretation. The one-way Fourier finite-difference wave propagator is used to extrapolate plane wave sources on the target layer to the surface. By analyzing the wave energy distribution at the surface extrapolated from the target layer, the shot or receiver locations needed for target layer imaging can be determined. Numerical tests using the SEG-EAGE salt model suggest that this method is useful for confirming the special seismic acquisition geometry layout for target-oriented imaging.

Integrated application of 3D seismic and microseismic data in the development of tight gas reservoirs

The development of unconventional resources, such as shale gas and tight sane gas, requires the integration of multi-disciplinary knowledge to resolve many engineering problems in order to achieve economic production levels. The reservoir heterogeneit3 revealed by different data sets, such as 3D seismic and microseismic data, can more full3 reflect the reservoir properties and is helpful to optimize the drilling and completioT programs. First, we predict the local stress direction and open or close status of the natura fractures in tight sand reservoirs based on seismic curvature, an attribute that reveals reservoi heterogeneity and geomechanical properties. Meanwhile, the reservoir fracture network is predicted using an ant-tracking cube and the potential fracture barriers which can affec hydraulic fracture propagation are predicted by integrating the seismic curvature attribute anc ant-tracking cube. Second, we use this information, derived from 3D seismic data, to assis in designing the fracture program and adjusting stimulation parameters. Finally, we interpre the reason why sand plugs will occur during the stimulation process by the integration of 3E seismic interpretation and microseismic imaging results, which further explain the hydraulic fracure propagation controlling factors and open or closed state of natural fractures in tigh sand reservoirs.

Key parameter optimization and analysis of stochastic seismic inversion

Stochastic seismic inversion is the combination of geostatistics and seismic inversion technology which integrates information from seismic records, well logs, and geostatistics into a posterior probability density function （PDF） of subsurface models. The Markov chain Monte Carlo （MCMC） method is used to sample the posterior PDF and the subsurface model characteristics can be inferred by analyzing a set of the posterior PDF samples. In this paper, we first introduce the stochastic seismic inversion theory, discuss and analyze the four key parameters： seismic data signal-to-noise ratio （S/N）, variogram, the posterior PDF sample number, and well density, and propose the optimum selection of these parameters. The analysis results show that seismic data S/N adjusts the compromise between the influence of the seismic data and geostatistics on the inversion results, the variogram controls the smoothness of the inversion results, the posterior PDF sample number determines the reliability of the statistical characteristics derived from the samples, and well density influences the inversion uncertainty. Finally, the comparison between the stochastic seismic inversion and the deterministic model based seismic inversion indicates that the stochastic seismic inversion can provide more reliable information of the subsurface character.

A study on the temporal and spatial characteristics of droughts following earthquakes

According to the ＂jacking-up＂ theory, which relates the cause of earthquakes to outer core convection ascension bodies, the crust will gradually recover after an earthquake. In such cases, the crust is stretched, the underground temperature is reduced, precipitation decreases, and drought occurs. In this paper, precipitation is compared with ground temperature and seismic data to determine the spatial and temporal relationship between earthquakes and subsequent droughts. Our objective is to develop a new method of drought prediction. With a few exceptions in location, the analysis of the first drought to occur after the Ms 〉 7 earthquakes in China's Mainland and the adjacent areas since 1950 shows that droughts tended to occur in regions near earthquake epicenters and in the eastern regions of the epicenters at the same latitude within six months after the earthquakes. In addition, and the differences between the starting time of the earthquakes and the droughts nearly share the same probability of 0 to 6 months. After careful analysis of 34 Ms 〉 6.5 earthquakes occurring in western China from 1980 to 2011, we determined that a second drought tends to occur approximately six months following the first drought, indicating a quasi-half-year period. Moreover, the duration of the quasi-half-year fluctuation increases with the magnitude of earthquake, at approximately 2.5 years for Ms 6.5 earthquake and approximately 5 years for Ms 8 earthquake.

Using 4C OBS to reveal the distribution and velocity attributes of gas hydrates at the northern continental slope of South China Sea

To investigate the distribution and velocity attributes of gas hydrates in the northern continental slope of South China Sea, Guangzhou Marine Geological Survey conducted four-component （4C） ocean-bottom seismometer （OBS） surveys. A case study is presented to show the results of acquiring and processing OBS data for detecting gas hydrates. Key processing steps such as repositioning, reorientation, PZ summation, and mirror imaging are discussed. Repositioning and reorientation find the correct location and direction of nodes. PZ summation matches P- and Z-components and sums them to separate upgoing and downgoing waves. Upgoing waves are used in conventional imaging, whereas downgoing waves are used in mirror imaging. Mirror imaging uses the energy of the receiver ghost reflection to improve the illumination of shallow structures, where gas hydrates and the associated bottom-simulating reflections （BSRs） are located. We developed a new method of velocity analysis using mirror imaging. The proposed method is based on velocity scanning and iterative prestack time migration. The final imaging results are promising. When combined with the derived velocity field, we can characterize the BSR and shallow structures; hence, we conclude that using 4C OBS can reveal the distribution and velocity attributes of gas hydrates.

Seismic data reconstruction based on CS and Fourier theory

Traditional seismic data sampling follows the Nyquist sampling theorem. In this paper, we introduce the theory of compressive sensing （CS）, breaking through the limitations of the traditional Nyquist sampling theorem, rendering the coherent aliases of regular undersampling into harmless incoherent random noise using random undersampling, and effectively turning the reconstruction problem into a much simpler denoising problem. We introduce the projections onto convex sets （POCS） algorithm in the data reconstruction process, apply the exponential decay threshold parameter in the iterations, and modify the traditional reconstruction process that performs forward and reverse transforms in the time and space domain. We propose a new method that uses forward and reverse transforms in the space domain. The proposed method uses less computer memory and improves computational speed. We also analyze the antinoise and anti-aliasing ability of the proposed method, and compare the 2D and 3D data reconstruction. Theoretical models and real data show that the proposed method is effective and of practical importance, as it can reconstruct missing traces and reduce the exploration cost of complex data acquisition.

An improved predictive deconvolution based on maximization of non-Gaussianity

The predictive deconvolution algorithm （PD）, which is based on second-order statistics, assumes that the primaries and the multiples are implicitly orthogonal. However, the seismic data usually do not satisfy this assumption in practice. Since the seismic data （primaries and multiples） have a non-Gaussian distribution, in this paper we present an improved predictive deconvolution algorithm （IPD） by maximizing the non-Gaussianity of the recovered primaries. Applications of the IPD method on synthetic and real seismic datasets show that the proposed method obtains promising results.

Experimental analysis and application of sparsity constrained deconvolution

Sparsity constrained deconvolution can improve the resolution of band-limited seismic data compared to conventional deconvolution. However, such deconvolution methods result in nonunique solutions and suppress weak reflections. The Cauchy function, modified Cauchy function, and Huber function are commonly used constraint criteria in sparse deconvolution. We used numerical experiments to analyze the ability of sparsity constrained deconvolution to restore reflectivity sequences and protect weak reflections under different constraint criteria. The experimental results demonstrate that the performance of sparsity constrained deconvolution depends on the agreement between the constraint criteria and the probability distribution of the reflectivity sequences; furthermore, the modified Cauchy- constrained criterion protects the weak reflections better than the other criteria. Based on the model experiments, the probability distribution of the reflectivity sequences of carbonate and clastic formations is statistically analyzed by using well-logging data and then the modified Cauchy-constrained deconvolution is applied to real seismic data much improving the resolution.

Seismic data analysis based on spatial subsets

There are some limitations when we apply conventional methods to analyze the massive amounts of seismic data acquired with high-density spatial sampling since processors usually obtain the properties of raw data from common shot gathers or other datasets located at certain points or along lines. We propose a novel method in this paper to observe seismic data on time slices from spatial subsets. The composition of a spatial subset and the unique character of orthogonal or oblique subsets are described and pre-stack subsets are shown by 3D visualization. In seismic data processing, spatial subsets can be used for the following aspects： （1） to check the trace distribution uniformity and regularity; （2） to observe the main features of ground-roll and linear noise; （3） to find abnormal traces from slices of datasets; and （4） to QC the results of pre-stack noise attenuation. The field data application shows that seismic data analysis in spatial subsets is an effective method that may lead to a better discrimination among various wavefields and help us obtain more information.

Geophysical data sparse reconstruction based on L0-norm minimization

Missing data are a problem in geophysical surveys, and interpolation and reconstruction of missing data is part of the data processing and interpretation. Based on the sparseness of the geophysical data or the transform domain, we can improve the accuracy and stability of the reconstruction by transforming it to a sparse optimization problem. In this paper, we propose a mathematical model for the sparse reconstruction of data based on the LO-norm minimization. Furthermore, we discuss two types of the approximation algorithm for the LO- norm minimization according to the size and characteristics of the geophysical data： namely, the iteratively reweighted least-squares algorithm and the fast iterative hard thresholding algorithm. Theoretical and numerical analysis showed that applying the iteratively reweighted least-squares algorithm to the reconstruction of potential field data exploits its fast convergence rate, short calculation time, and high precision, whereas the fast iterative hard thresholding algorithm is more suitable for processing seismic data, moreover, its computational efficiency is better than that of the traditional iterative hard thresholding algorithm.

Research on Evaluation and Optimization of Regional Land Use Structure Efficiency

[Objective] The purpose of this study is to judge the rationality of land use structure from new perspectives and method. [Method] The theory of data en- velopment analysis （DEA） has been employed in this article. A case study of evalu- ation of land use structure for Wujin district in 2012 is provided to illustrate the ap- plication of this research model. At the micro perspective, land use efficiency has been analyzed by use of data envelopment analysis model. The optimization and adjustment direction in land use structure has been pointed out from the view point of land use efficiency. [Result] The result has shown that the research method based on information entropy and data envelopment analysis model can effectively overcome the drawback of traditional evaluation techniques. Compared with the tra- ditional methods, the new optimization model based on a structure of the multi-crite- ria factors and objective weighting method can evaluate the rationality of land use structure more comprehensively.

3D electrical resistivity inversion using prior spatial shape constraints

To minimize the number of solutions in 3D resistivity inversion, an inherent problem in inversion, the amount of data considered have to be large and prior constraints need to be applied. Geological and geophysical data regarding the extent of a geological anomaly are important prior information. We propose the use of shape constraints in 3D electrical resistivity inversion, Three weighted orthogonal vectors （a normal and two tangent vectors） were used to control the resistivity differences at the boundaries of the anomaly. The spatial shape of the anomaly and the constraints on the boundaries of the anomaly are thus established. We incorporated the spatial shape constraints in the objective function of the 3D resistivity inversion and constructed the 3D resistivity inversion equation with spatial shape constraints. Subsequently, we used numerical modeling based on prior spatial shape data to constrain the direction vectors and weights of the 3D resistivity inversion. We established a reasonable range between the direction vectors and weights, and verified the feasibility and effectiveness of using spatial shape prior constraints in reducing excessive structures and the number of solutions. We applied the prior spatially shape-constrained inversion method to locate the aquifer at the Guangzhou subway. The spatial shape constraints were taken from ground penetrating radar data. The inversion results for the location and shape of the aquifer agree well with drilling data, and the number of inversion solutions is significantly reduced.

The anti-correlation method for removing harmonic distortion in vibroseis slip-sweep data

The slip-sweep technique is one of the high-efficiency, high-fidelity, and environmental vibroseis seismic prospecting techniques which consists of a vibrator group sweeping without waiting for the previous group＇s sweep to terminate. The cycle time can be reduced drastically and hence the production efficiency can be increased significantly but harmonic distortion of one sweep will leak into the record of the other sweep. In this paper, we propose an anti-correlation method for removing harmonic distortion in vibroseis data. This method is based on decomposition of the ground force signal into fundamental and harmonic components. Then the corresponding anti-correlation operator can be computed to estimate the energy of each harmonic after correlating the vibroseis data with the corresponding harmonic component. Finally, the vibroseis harmonic noise to be removed can be obtained by subtracting the extracted harmonic noise from the traces of the previous group＇s sweep. The advantage of the proposed method is that it can process both uncorrelated and correlated vibroseis seismic data. Moreover, the algorithm is simple, stable, and computationally fast. Especially, the significant contribution of this method is a considerable reduction in the harmonic without any alteration of the desired signals. The method was tested on both synthetic and field data sets to validate the good harmonic noise suppression results.

Separation of P. and SV-wavefields from multi-componen seismic data

In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.