Full 3-D numerical modeling of borehole electric image logging and the evaluation model of fracture

A full 3-D finite element method numerical modeling program is written based on the principle and technical specification of borehole electric image well logging tool. The response of well logging is computed in the formation media model with a single fracture. The effect of changing fracture aperture and resistivity ratio to the logging response is discussed. The identification ability for two parallel fractures is also present. A quantitative evaluation formula of fracture aperture from borehole electric image logging data is set up. A case study of the model well is done to verify the accuracy of the for-mula. The result indicates that the formula is more accurate than the foreign one.     

Regional gravity field modeling based on rectangular harmonic analysis

Regional gravity field modeling with high-precision and high-resolution is one of the most important scientific objectives in geodesy,and can provide fundamental information for geophysics,geodynamics,seismology,and mineral exploration.Rectangular harmonic analysis(RHA)is proposed for regional gravity field modeling in this paper.By solving the Laplace’s equation of gravitational potential in local Cartesian coordinate system,the rectangular harmonic expansions of disturbing potential,gravity anomaly,gravity disturbance,geoid undulation and deflection of the vertical are derived,and so are the formula for signal degree variance and error degree variance of the rectangular harmonic coefficients(RHC).We also present the mathematical model and detailed algorithm for the solution of RHC using RHA from gravity observations.In order to reduce the edge effects caused by periodic continuation in RHA,we propose the strategy of extending the size of computation domain.The RHA-based modeling method is validated by conducting numerical experiments based on simulated ground and airborne gravity data that are generated from geopotential model EGM2008 and contaminated by Gauss white noise with standard deviation of 2 mGal.The accuracy of the 2.5′×2.5′geoid undulations computed from ground and airborne gravity data is 1 and 1.4cm,respectively.The standard error of the gravity disturbances that downward continued from the flight height of 4 km to the geoid is only 3.1 mGal.Numerical results confirm that RHA is able to provide a reliable and accurate regional gravity field model,which may be a new option for the representation of the fine structure of regional gravity field.     

Prediction of sedimentation in reservoirs by combining catchment based model and stream based model with limited data

Estimation of sedimentation in reservoirs helps in the management and design of the reservoir's useful capacity. This research was done on the Awash River basin at the Koka Dam Reservoir in Ethiopia. The method applied was the loose integration of the Soil and Water Assessment Tool(SWAT) model and Hydrologic Engineering Center-River Analysis System(HEC-RAS) model for the estimation of the sediment load reaching the reservoir. The SWAT model was used for the estimation of erosion at the catchment level,and the HEC-RAS model was applied to estimate the sediment transport in the river channel. The implemented method allows sedimentation in the floodplains and bed shear stress to be considered in the sediment modeling, which cannot be considered in the SWAT model. In addition, the river cross sectional properties and the hydrodynamic processes in the rivers were considered in the modeling process. The data used in this study are a combination of i) observed data collected by government agencies, ii) data available online in data repositories, and iii) data extracted from remote sensing in the Shuttle Radar Topographic Mission(SRTM) Digital Elevation Model(DEM). The calibration and validation of the SWAT model was done by using Sequential Uncertainty Fitting(SUIF-2) calibration and validation tools. The HEC-RAS model was calibrated by adjusting the roughness factor. The output from the integrated approaches gives better estimates of flow and sediment near the inlet to the reservoir, with coefficients of determination of 0.85 and 0.67, respectively, and Nash Sutcliffe coefficients of model fit efficiency of 0.90 and 0.62, respectively, for daily simulations.     

太湖风生流数值模式及其对总磷分布的影响

A new up-winding finite element numerical model, which is two-dimensional and suitable for modeling lake current and the distribution of total phosphorus(TP) in shallow lakes, is derived. Moreover, it is used in the study of wind-driven current in Taihu Lake, and the impact of lake current field on the distribution of TP is also discussed.     

A laboratory study of rheological properties of mudflows in Hangzhou Bay,China

A laboratory study of the rheology of mudflows in Hangzhou Bay, China, is reported in this paper. Both the steady and oscillatory （dynamic） rheological properties are studied using RMS-605 rheometer. A Dual-Bingham model is proposed for analyzing flow curves and compared with Worrall-Tuliani model. It is found that Dual-Bingham plastic rheological model is easier to implement than Worrall-Tuliani model and can provide satisfactory representations of the steady mudflows in Hangzhou Bay and other published data. The dependence of the yield stress and viscosity on sediment concentration is discussed based on the data from Hangzhou Bay mud and other published data. For the dynamic rheologieal properties of Hangzhou Bay mud, empirical expressions for elastic modulus and dynamic viscosity are provided in the form of exponential functions of sediment volume concentration, and comparisons with other published data also discussed.     

Suppression of strong random noise in seismic data by using time-frequency peak filtering

Time-frequency peak filtering (TFPF) is highly efficient in suppressing random noise in seismic data. Although the hypothesis of stationary Gaussian white noise cannot be fulfilled in practical seismic data, TFPF can effectively suppress white and colored random noise with different intensities, as can be theoretically demonstrated by detecting such noise in synthetic seismic data. However, a "zero-drift" effect is observed in the filtered signal and is independent of the average power and variance of the random noise, but related to its mean value. Furthermore, we consider the situation where the local linearization of the seismic data cannot be satisfied absolutely and study the "distortion" characteristics of the filtered signal using TFPF on a triangular wave. We found that over-compensation is possible in the frequency band for the triangular wave. In addition, it is nonsymmetrical and has a relationship to the time-varying curvature of the seismic wavelet. The results also present an improved approach for TFPF.     

KINETIC MODELING OF CONSTITUTIVE RELATIONS FOR PARTICLE MOTION IN LOW TO MODERATELY CONCENTRATED FLOWS

Formulating underlying mechanisms of concentrated solid-liquid flows is essential for simulation of various industrial processes and natural phenomena. A generalized constitutive model for particle motion in flows with low to moderate solids concentrations is developed. This generalized model facilitates characterization of inelastic collisions, particle-fluid interactions, and shearing effects. Moderately concentrated simple shear flows of a sand-water mixture are analyzed, and comparisons of model predictions and experimental data are in good agreement. This model exhibits sound performance in characterizing particle motion for wide ranges of concentration and shear rate, and may supply a reasonable and competent alternative to previous models developed for dilute and rapid-granular flows when applied to moderately concentrated situations. The concentration approaches zero （C --~ O） asymptote is observed at a relatively high shear rate in model predictions. Assumption of low collisional dissipation of the particle phase as C → 0 is more reasonable for this observation, compared to that without the interstitial fluid effect. Accurately modeling energy dissipation is important for characterizing the stability of dilute simple shear flows of solid-liquid mixtures. Incorporating friction forces will also facilitate improvement of the applicability of this generalized model to flows at extremely high concentrations.     

New Understanding of Variation Characteristics of Well Water Level in Majiagou Mine, Tangshan City, Hebei, China

Based on the well water level data in Majiagou Mine in Tangshan City since 2000, it is suggested that the water level rise from 2008-2014 is possibly associated with reduced infiltration caused by less precipitation, together with the effect of less groundwater exploitation. However, the water level rise from 1992-1998 cannot be interpreted properly if we extend the data source back to 1991. After comparing the data with long-term water level data of five wells with different hydro-geological units in Tangshan and Qinhuangdao, the study indicates that the long-term water level variation cannot be only attributed to the effect of less precipitation and the control of groundwater exploitation, but also with the influence of regional stress field change. With the support of regional geological and tectonic information, combined with the comprehensive analysis of CBM data obtained from the coal field, new understandings of the dynamic characteristics of annual variation are obtained.     

Seismic inversion for underground fractures detection based on effective anisotropy and fluid substitution

Underground fractures play an important role in the storage and movement of hydrocarbon fluid. Fracture rock physics has been the useful bridge between fracture parameters and seismic response. In this paper, we aim to use seismic data to predict subsurface fractures based on rock physics. We begin with the construction of fracture rock physics model. Using the model, we may estimate P-wave velocity, S-wave velocity and fracture rock physics parameters. Then we derive a new approximate formula for the analysis of the relationship between fracture rock physics parameters and seismic response, and we also propose the method which uses seismic data to invert the elastic and rock physics parameters of fractured rock. We end with the method verification, which includes using well-logging data to confirm the reliability of fracture rock physics effective model and utilizing real seismic data to validate the applicability of the inversion method. Tests show that the fracture rock physics effective model may be used to estimate velocities and fracture rock physics parameters reliably, and the inversion method is resultful even when the seismic data is added with random noise. Real data test also indicates the inversion method can be applied into the estimation of the elastic and fracture weaknesses parameters in the target area.     

3D Modeling and Born approximation inversion for the borehole surface electromagnetic method

We present a 3D approach to numerical modeling of the borehole-surface electromagnetic （BSEM） method. The 3D electromagnetic response created by a vertical line current source in a layered medium is modeled using the 3D integral equation method. The modeling results are consistent with analytical solutions. 3D Born approximation inversion of BSEM data is also conducted for reservoir delineation. The inversion method is verified by a synthetic reservoir model.     

Three-dimensional inversion of borehole-surface electrical data based on quasi-analytical approximation

3D inversion of borehole-surface electrical data for complex geo-electrical models is still a challenging problem in geophysical exploration. We have developed a program for 3D inversion to borehole-surface electrical data based on the quasi-analytical approximation （QA） and re-weighted regularized conjugate gradient method （RRCG） algorithms using Visual Fortran 6.5. Application of the QA approximation to forward modeling and Frechet derivative computations speeds up the calculation dramatically. The trial calculation for synthetic data of theoretical model showed that the program is fast and highly precise.     

A holistic computational model for prediction of clay suspension structure

The formation of clay suspensions involves complex interactions among clay particles subjected to the geochemical environment during the sedimentation process. The structural characteristics have a major influence on the physical and mechanical behavior of the suspension. A modeling framework involving a Discrete Element Method (DEM) model with customized particle mechanical interactions is proposed in this paper for the holistic prediction of the physical structure of clay suspensions. The particle interaction force model is based on the Derjaguin, Landau, Verwey, and Overbeek (DLVO) theory that accounts for electrostatic repulsion, van der Waals attraction, contact repulsion, etc. Kaolinite is used as the model clay to demonstrate the model performance. The surface charge density of kaolinite is obtained through Atomic Force Microscope (AFM) force measurement and is implemented in the particle interaction force model in the subsequent simulations. Influencing factors, such as centrifugal acceleration, ionic concentration, platy structure of particles, and particle size, on the formation of kaolinite suspensions are studied with the numerical model and compare favorably well with the experimental data. This work lays down a unique framework consisting of computational modeling and microscale characterization of clay particles to holistically predict the characteristics of clay suspensions, which paves the basis to model and predict their bulk physical and mechanical behavior.     

Co-seismic fault geometry and slip distribution of the 26 December 2004, giant Sumatra–Andaman earthquake constrained by GPS,coral reef,and remote sensing data

We analyze co-seismic displacement field of the 26 December 2004, giant Sumatra–Andaman earthquake derived from Global Position System observations,geological vertical measurement of coral head, and pivot line observed through remote sensing. Using the co-seismic displacement field and AK135 spherical layered Earth model, we invert co-seismic slip distribution along the seismic fault. We also search the best fault geometry model to fit the observed data. Assuming that the dip angle linearly increases in downward direction, the postfit residual variation of the inversed geometry model with dip angles linearly changing along fault strike are plotted. The geometry model with local minimum misfits is the one with dip angle linearly increasing along strike from 4.3oin top southernmost patch to 4.5oin top northernmost path and dip angle linearly increased. By using the fault shape and geodetic co-seismic data, we estimate the slip distribution on the curved fault. Our result shows that the earthquake ruptured *200-km width down to a depth of about 60 km.0.5–12.5 m of thrust slip is resolved with the largest slip centered around the central section of the rupture zone78N–108N in latitude. The estimated seismic moment is8.2 9 1022 N m, which is larger than estimation from the centroid moment magnitude(4.0 9 1022 N m), and smaller than estimation from normal-mode oscillation data modeling(1.0 9 1023 N m).     

3-D physical model in strong ground motion numerical simulation： A case study of Kunming basin

Seismic hazard assessment based on urban active faults can provide scientific bases for city planning and projectconstruction,while numerical simulation of strong ground motion is an important method for seismic hazard pre-diction and assessment.A 3-D physical model in conformity with real strata configuration of(mainly)the Quater-nary is a prerequisite to ensure the reliability of the simulation results.In this paper,we give a detailed account ofthe technical scheme and process for creating a 3-D physical model in Kunming basin.The data used are synthe-sized from seismogeological data,borehole data,topographic data,digital elevation mode(DEM)data,seismicexploration results and wave velocity measurements.Stratigraphic division is based mainly on shear wave velocity,with strata sequence taken into consideration.The model construction is finally accomplished with ArcGIS andmany relevant programming techniques via layer-by-layer stacking(in depth direction)of the adjacent mediuminterfaces(meshes).Meanwhile,a database of 3-D physical models is set up,which provides model data and pa-rameters for strong ground motion simulation.Some processing methods and significant issues are also addressedin the paper in accordance with different types of exploration and experimental data.     

Vulnerability of ordinary moment resistant concrete frames

The effects of structural modeling (bar slip in lap splice), ground motion selection process (epsilon effect) and size of a structure (number of bays and stories) on the fragility of reinforced concrete ordinary moment resisting frames are investigated. An analytical model is developed to account for bar slip in lap splice, which exhibits good correlation with existing experimental data. Then, incremental dynamic analysis is used to derive the fragility curves for four model structures. The model structures simulate frames with a different number of bays and stories. Finally, the fragility curves are corrected for the epsilon effect. The results show that slip in the lap splice can significantly increase the failure probability, especially for the collapse prevention limit state. At the same time, it is found that spectral shape has a significant impact on the fragility curves. It is also found that accounting for or ignoring bar slip or epsilon effects increases the probability of failure for larger structures. This indicates an unconservative bias in the safety of larger structures.     

Hierarchy modeling of subsurface palaeochannel reservoir architecture

The studies on fluvial reservoir architecture are mainly aimed at outcrop and modern deposition,but rarely at the subsurface reservoir,so there are few effective methods to predict the distribution of subsurface reservoir architectures. In this paper,taking the meandering river reservoir of Guantao formation Gudao Oilfield,Jiyang depression,Baohai Gulf Basin,East China as an example,the archi-tectural modeling method of complex meandering belt reservoir is proposed,that is hierarchy con-straint,pattern fitting and multi-dimensional interaction. Architectures of meandering river reservoir can be divided into three hierarchies: meandering channel sandbody,point bar and lateral accretion body. Different hierarchies of the quantitative architecture pattern are fitted to subsurface well data (including dynamic monitoring data) in different hierarchies through one-dimensional hole,2D profiles and plane and 3D space,which are verified by each other. And then 3D model in different hierarchies is established. At the same time,the quantificational relationship between width of active river and the scale of point bar is set up,and the scale of lateral accretion sand body and shale beddings is con-firmed with horizontal well data. The study not only has significant meaning on the development of geology,but also can improve the oilfield exploitation greatly.     

Preliminary Discussion on the Implication of Determination of Source Depth of Earthquakes

In an earthquake catalog,the source depth of earthquakes is one of the most basic parameters; however it is also one of the most difficult to measure.The source depths in the PDE bulletin,the centroid depths from CMT inversion of USGS and Harvard University,and those obtained by depth phases of broadband waveform data are not equal to each other,and they are also have different meanings.Also,the depth obtained from broadband waveform modeling is another kind of depth of different a rupture "core." This paper preliminarily discusses the physical implication of the source depth of earthquakes with different measuring methods.     

大地电磁法三维共轭梯度反演研究

Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.