Test of frequency characteristic for well and its response to seismic waves

The permeation parameters have been calculated by forefathers on the basis of permeation theory by means of the Slug test (Yin, Zheng, 1992) and the restoration curves of well level. We are interested in oscillation of the well level when we make Slug test. Both the permeation parameters and frequency parameters, i.e., natural period and damping coefficients of well aquifer, have been calculated on the basis of vibration theory by means of the oscillation curves. Not only this has given a new method, but also the different response of well level to seismic waves has been explained by it in theory.     

Using Grey System Theory for Earthquake Forecast

By combining conventional grey correlation analysis,grey clustering method and grey forecasting methods with our multi-goal forecast thoughts and the techniques of grey time series processing,we develop six different grey earthquake forecast models in this paper,Using the record of major earthquakes in Japan from 1872 to 1995,we forecast future earthquakes in Japan.We develop an earthquake forecast model.By using the major earthquakes in Japan from 1872 to 1984,we forecast earthquakes from 1985 to 1995 and check the precision of the grey earthquake models.We find that the grey system theory can be applied to earthquake forecast.We introduce the above analysis methods and give a real example to evaluate and forecast.We also further discuss the problems of how to improve the precision of earthquake forecast and how to strengthen the forecast models in future research.     

Seismicity acceleration model and its application to several earthquake regions in China

With the theory of subcritical crack growth, we can deduce the fundamental equation of regional seismicity acceleration model. Applying this model to intraplate earthquake regions, we select three earthquake subplates: North China Subplate, Chuan-Dian Block and Xinjiang Subplate, and divide the three subplates into seven researched regions by the difference of seismicity and tectonic conditions. With the modified equation given by Sornette and Sammis (1995), we analysis the seismicity of each region. To those strong earthquakes already occurred in these region, the model can give close fitting of magnitude and occurrence time, and the result in this article indicates that the seismicity acceleration model can also be used for describing the seismicity of intraplate. In the article, we give the magnitude and occurrence time of possible strong earthquakes in Shanxi, Ordos, Bole-Tuokexun, Ayinke-Wuqia earthquake regions. In the same subplate or block, the earthquake periods for each earthquake region are similar in time interval. The constant αin model can be used to describe the intensity of regional seismicity, and for the Chinese Mainland, α is 0.4 generally. To the seismicity in Taiwan and other regions with complex tectonic conditions, the model does not fit well at present.     

Preface for the Special Issue of 2000 International Conference for CT & 3D Imaging, Beijing

In the last 20-years, CT theory and applications in medical diagnosis have developed fast in Asia &China, and extended into the fields of seismology, geophysics, bio-medicine, industry and engineering,and got a lot of achievements and beneficial results. A lot of academic worldwide exchanges of CT theoryand application are important for human. Hence, we have gotten the corresponding official permissiondocument from China Ministry of Science and Technology and China Association for Science…     

Theory and Method for Identifying Well Water Level Anomalies in a Groundwater Overdraft Area

The overexploitation of underground water leads to the continuous drawdown of groundwater levels, change of water quality and dry-up in dynamic water level observation wells. Due to land subsidence, the well pipes uplift and the observation piping systems are damaged. These environmental geology problems can present serious difficulties for the identification of earthquake anomalies by groundwater level observation. Basied on hydrogeological theories and methods, the paper analyzes the relations of the water balance state of aquifers with stress-strain conditions and the water level regime, and then discusses preliminarily the theory and method for identifying well water level anomalies in a groundwater overdraft area. The result shows that we can accurately judge the nature of the anomaly according to the diffusion character of the drawdown funnel in the well area in combination with the aforementioned theory and method and multi-year variation patterns obtained from existing data. The results of the research are helpful for distinguishing the influence of single centralized water pumping from the long-term overdraft of water on the water level, correctly recognizing water level anomalies in the groundwater overdraft area and increasing the level of earthquake analysis and prediction.     

Sensitivity studies of a high accuracy surface modeling method

The sensitivities of the initial value and the sampling information to the accuracy of a high accuracy surface modeling(HASM) are investigated and the implementations of this new modeling method are modified and enhanced. Based on the fundamental theorem of surface theory, HASM is developed to correct the error produced in geographical information system and ecological modeling process. However, the earlier version of HASM is theoretically incomplete and its initial value must be produced by other surface modeling methods, such as spline, which limit its promotion. In other words, we must use other interpolators to drive HASM. According to the fundamental theorem of surface theory, we modify HASM, namely HASM.MOD, by adding another important nonlinear equation to make it independent of other methods and, at the same time, have a complete and solid theory foundation. Two mathematic surfaces and monthly mean temperature of 1951–2010 are used to validate the effectiveness of the new method. Experiments show that the modified version of HASM is insensitive to the selection of initial value which is particular important for HASM. We analyze the sensitivities of sampling error and sampling ratio to the simulation accuracy of HASM.MOD. It is found that sampling information plays an important role in the simulation accuracy of HASM.MOD. Another feature of the modified version of HASM is that it is theoretically perfect as it considers the third equation of the surface theory which reflects the local warping of the surface. The modified HASM may be useful with a wide range of spatial interpolation as it would no longer rely on other interpolation methods.     

Study on the method for comprehensive discrimination of small earthquakes

IntroductionThe present-day digital seismic observation technique provides us a large amount of accurate seismic records of high-precision, wide-band, and large dynamic range and conditions for higher-solution and higher-speed discrimination of earthquakes. But the problem we face is how to use the modern theory and method for digital signal processing in fast and precise extraction of the characteristic parameters from the large amount of digital seismic records. A basis or the important prob…     

A wave propagation model with the Biot and the fractional viscoelastic mechanisms

Energy loss in porous media containing fluids is typically caused by a variety of dynamic mechanisms.In the Biot theory,energy loss only includes the frictional dissipation between the solid phase and the fluid phase,resulting in underestimation of the dispersion and attenuation of the waves in the low frequency range.To develop a dynamic model that can predict the high dispersion and strong attenuation of waves at the seismic band,we introduce viscoelasticity into the Biot model and use fractional derivatives to describe the viscoelastic mechanism,and finally propose a new wave propagation model.Unlike the Biot model,the proposed model includes the intrinsic dissipation of the solid frame.We investigate the effects of the fractional order parameters on the dispersion and attenuation of the P-and S-waves using several numerical experiments.Furthermore,we use several groups of experimental data from different fluid-saturated rocks to testify the validity of the new model.The results demonstrate that the new model provides more accurate predictions of high dispersion and strong attenuation of different waves in the low frequency range.     

Characteristics and Implication of the Earthquake Swarm Occurred in Fuzhou in September 1999

on Septmeber 23,1999,an earthquake swarm occured in Fuzhou,Because the swarm occurred in the region where earthquaks occurred scarcely before and very close to the center of the city as well as shortly after the Jiji earthquake with Ms7.6 in Taiwan,September 21,1999,has aroused interest broadly.In this paper,we analyzed the characteristics of spatial and temporal distribution of the earthquake swarm and validated magnitude-number constituent of the swarm is special.In present theory,the earthquake swarm means that a small scale macro original rupture has formed in the layer of the crust in Fuzhou region where moderately strong earthquake risk exists.     

Analytical estimation of mixing coefficient induced by surface wave-generated turbulence based on the equilibrium solution of the second-order turbulence closure model

Based on the equations of motion and the assumption that ocean turbulence is of isotropy or quasi-isotropy, we derived the closure equations of the second-order moments and the variation equations for characteristic quantities, which describe the mechanisms of advection transport and shear instability by the sum of wave-like and eddy-like motions and circulation. Given that ocean turbulence generated by wave breaking is dominant at the ocean surface, we presented the boundary conditions of the turbulence kinetic energy and its dissipation rate, which are determined by energy loss from wave breaking and entrainment depth respectively. According to the equilibrium solution of the variation equations and available data of the dissipation rate, we obtained an analytical estimation of the characteristic quantities of surface-wave-generated turbulence in the upper ocean and its related mixing coefficient. The derived kinetic dissipation rate was validated by field measurements qualitatively and quantitatively, and the mixing coefficient had fairly good consistency with previous results based on the Prandtl mixing length theory.     

信息的不可修改性与当代科学某些问题

本文指出了"信息确定性"所涉及的信息的不可 修改性、信息"保真性"是认识观念的重大变革。 以取消数量稳定性限制条件,得出了非线性不稳 定具有特殊性信息及其周期性,并可解释翁文波 的"可公度"法是运用了特殊性信息的特殊规律, 在物质演化问题上不能将数量的一般规律拓广至 特殊性信息中;由此涉及了当代科学没有解决物 质演化的"时间"、特殊性、非惯性系及其数量分 析伪造信息等问题。信息的客观存在性和不可修 改性,已实质上表明了将近代的"信息科学"的 信息分析等同于当代科学的数量处理是认识上的 误解。即"信息"不能等同于"数量"。     

Relationships between Atmospheric Precipitation,Evaporation, Temperature,and Production of Aquatic Ecosystems

The average production capacity of aquatic ecosystems is shown to be proportional to the coefficient of variations of substance concentrations in atmospheric precipitation and the sum of effective temperatures. The values of ichthyomass and primary production of phytoplankton per 1°C of the sum of effective temperatures multiplied by the coefficient of variations of substance concentration in atmospheric precipitation are proposed as ecological thermometers to be used in water bodies in different geographic zones of Europe and North America.     

Temperature and heat flux spectra in the turbulent buoyancy subrange

The physical nature of motions with scales intermediate between approximately isotropic turbulence and quasi-linear internal gravity waves is not understood at the present time. Such motions play an important role in the energetics of small scales processes, both in the ocean and in the atmosphere, and in vertical transport of heat and constituents. This scale range is currently interpreted either as a saturated gravity waves field or as a buoyancy range of turbulence.We first discuss some distinctive predictions of the classical (Lumley, Phillips) buoyancy range theory, recently improved (Weinstock, Dalaudier and Sidi) to describe potential energy associated with temperature fluctuations. This theory predicts the existence of a spectral gap in the temperature spectra and of an upward mass flux (downward buoyancy and heat fluxes), strongly increasing towards large scales. These predictions are contrasted with an alternate theory, assuming energetically insignificant buoyancy flux, proposed by Holloway.Then we present experimental evidences of such characteristic features obtained in the lower stratosphere with an instrumented balloon. Spectra of temperature, vertical velocity, and cospectra of both, obtained in homogeneous, weakly turbulent regions, are compared with theoretical predictions. These results are strongly consistent with the improved classical buoyancy range theory and support the existence of a significant downward heat flux in the buoyancy range.The theoretical implications of the understanding of this scale range are discussed. Many experimental evidences consistently show the need for an anisotropic theory of the buoyancy range of turbulence.     

电磁辐射"EMOLS"仪观测结果原理及震例

“九五”期间研制的“EMAOS”电磁辐射仪,其接收原理是按大地电磁测深的原理改进而研制的。选择接收频段为0．01～10Hz,该频段探测深度在0～20km,是浅源地震发生的层位。记录的电磁扰动是在应力场作用下,动电效应产生的电磁场。电场和磁场同时记录增强了信息可信度的判断能力,记录震例表明,异常规律与已有模拟结果基本相同。因此,已总结的模拟记录预报指标等可在数字记录结果中使用。     

地震观测标准子体系研究

根据《地震行业标准体系表》(DB／T1-2000)和2002年及2004年国家科技部项目《数字化地震前兆观测技术标准》实施过程中得到的认识,运用标准体系的原理和建立标准体系的原则与方法,研究和分析了地震观测标准子体系的结构和层次划分以及地震观测专业通用标准和地震观测门类通用标准,为进一步修订地震行业标准体系表提出了建议性的结论。     

偶极子假设引起的大回线源瞬变电磁响应偏差分析

偶极子假设下的场点近似解与场点的瞬变电磁精确解存在一定的误差,文中对偶极子假设引起的误差进行了探索性研究.首先分析了由恒定电流偶极子公式“比拟”出的谐变偶极子位函数解.分别计算了静态场中磁偶极子和电偶极子的近似解与各自对应的未做偶极子假设的电流环和载流导线的精确解之间的误差值,然后分析了谐变偶极子微元与点电荷微元的误差...     

Displacements of ridges and subduction zones in the models of mantle convection with lithospheric plates

In the existing kinematic theory of the tectonics of lithospheric plates, the position and the parameters of plates are assigned a priori in the initial and boundary conditions. However, in the self-consistent dynamic theory, the properties of oceanic plates (as well as the structure of the mantle’s convection) should appear automatically as the solution of the differential equations of energy, mass, and momentum transfer for a viscous fluid. In this case, the viscosity of the mantle’s substance as a function of temperature, pressure, shear stress, and chemical composition must be taken from the data of laboratory measurements. In the present work, the results of the numerical solution of the equations of convection are presented in the problem formulation mentioned above on a simple model of heated viscous fluid with properties that correspond to the mantle’s substance. In this case, to reveal the main reason for the generation of plates and their influence on the convection, a number of simplifications are introduced; in particular, temperature variations in the viscosity in the mantle are disregarded. In spite of the undertaken simplifications, the models show how the rigid lithosphere can be split into separate plates immersed in the mantle, how in the course of evolution the sizes of plates and their number can change, and how in this case the ridges and subduction zones can be displaced.     

Quasilinear heating of electrons in the earth's plasmasphere

In the framework of a quasilinear theory we examine the interaction between thermal electrons and ion-cyclotron waves (ICW) in the outer plasmasphere of the earth. For this type of wave-particle interaction, a simplification of the quasilinear diffusion integral in a magneto-active plasma under plasmaspheric conditions is given. Under the assumption of a Maxwellian distribution of electrons we have calculated the collision frequency and the heating source as the electrons are scattered by ICW. The obtained values of intrinsic parameters of the outer plasmasphere may exceed greatly, accordingly, the Coulomb frequency of collisions and the heating source due to suprathermal electrons. ICW-heating causes the ionosphereward thermal flux to increase, and this must lead to an increase in electron temperature in lower-lying plasmaspheric regions and in the subauroral ionosphere. A quantitative estimation of the electron temperature for the hot zone, made in this paper, is consistent with available experimental data.     

俯冲初始时板块分界面形状对俯冲过程的影响

自板块理论建立以来,俯冲一直是学者们关心的热点问题.前人结合地质、地球物理、实验室物理实验和数值模拟等多种手段对这一问题进行了大量的研究.以往的研究更为关注俯冲过程中板块的作用、地幔流动的规律和物质的迁移与相变等问题,却常忽视了俯冲是如何开始的这一基本问题.同时,由于相关数据资料较为有限,更限制了俯冲启动的相关研究.因此,本文选取俯冲启动问题中板块分界面形状对俯冲过程的影响这一问题,使用有限元的方法进行了数值模拟.我们选择针对倾斜型、垂直型和弯曲型三种不同形状的板块分界面建立对比模型,比较它们演化至10Ma的过程我们发现:分界面几何形状的不同的确会对俯冲板块演化和海沟的深度产生影响.倾斜型模型的俯冲角度最大,海沟深度最深,俯冲深度最深;垂直型模型的俯冲角度和海沟深度仅次于倾斜型模型,俯冲深度最浅;弯曲型模型的俯冲角度最小,海沟深度最小,俯冲的深度介于倾斜型和垂直型之间.结合以上结论不难看出,俯冲角和海沟深度变化具有一定的相关性,俯冲角度越大,相应的海沟深度越大.     

A triple porosity scheme for fluid/solid substitution: theory and experiment

Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.