含激电信息的人工源频率域电磁法正演研究——大功率人工固定源极低频电磁法激电效应预研究

ABSTRACT Base on Dias-model, this paper study the forward modeling to IP response with MT souece and CSAMT source with the author’s program，that is to provide the theoretical basis for extracting the IP information from the far field, near field and transition field in extremely low frequency electromagnetic sounding with artificial strong fixed source signal. The outline of the paper is as follows. Firstly, review the history of the complex-resistivity-model, and analyses the bases for choosing Dias model; meanwhile, effects and responses of each parameter in Dias model are analyses. Afterward, study the forward modeling to 1D layered model with MT source, numerically simulate the IP effects of classic geoelectric objects; Find the clear anomaly and relationship between the peak value of the amplitude anomaly ratio and phase anomaly difference with and without IP parameters to be considered within a frequency range. On the basis of the modeling of MT, Study the CSAMT modeling with a dipole-source, obtain the anomaly responses and the relationship between anomaly and characters of object. Base on the infinite line source, study the 2D IP effect of geoelectric objects with 2D modeling. Calculate the response with different source distance, object depth and different wall rock resistivity systemically. Finally, conclude the viewpoints and give the discussion of the result, and point out the lack of this research also.

地震P-S波时差耦合作用的斜坡崩滑效应研究

China locates between the circum-Pacific and the Mediterranean-Himalayan seismic belt. The seismic activities in our country are very frequent and so are the collapses and slides of slope triggered by earthquakes. Many collapses and slides of slope take place mainly in the west of China with many earthquakes and mountains, especially in Sichuan and Yunnan Provinces. When a strong earthquake happening, the damage especially in mountains area caused by geological hazards it triggered such as rock collapses, landslides and debris flows is heavier than that it caused directly. A conclusion which the number of lives lost caused by geological hazards triggered by a strong earthquake in mountains area often accounts for a half even more of the total one induced by the strong earthquake can be made by consulting the statistical loss of several representative earthquakes. As a result, geological hazards such as collapses and slides of slope triggered by strong earthquakes attract wide attention for their great costs. Based on field geological investigation, engineering geological exploration and material data analysis, chief conclusions have been drawn after systematic research on formation mechanism, key inducing factors, dynamic characteristics of geological hazards such as collapses and slides of slope triggered by strong earthquakes by means of engineering geomechanics comprehensive analysis, finite difference numerical simulation test, in-lab dynamic triaxial shear test of rock, discrete element numerical simulation. Based on research on a great number of collapses and landslides triggered by Wenchuan and Xiaonanhai Earthquake, two-set methods, i.e. the method for original topography recovering based on factors such as lithology and elevation comparing and the method for reconstructing collapsing and sliding process of slope based on characteristics of seism tectonic zone, structural fissure, diameter spatial distribution of slope debris mass, propagation direction and mechanical property of seismic wave, have been gotten. What is more, types, formation mechanism and dynamic characteristics of collapses and slides of slope induced by strong earthquakes are discussed comprehensively. Firstly, collapsed and slided accumulative mass is in a state of heavily even more broken. Secondly, dynamic process of slope collapsing and sliding consists of almost four stages, i.e. broken, thrown, crushed and river blocked. Thirdly, classified according to failure forms, there are usually four types which are made up of collapsing, land sliding, land sliding-debris flowing and vibrating liquefaction. Finally, as for key inducing factors in slope collapsing and sliding, they often include characteristics of seism tectonic belts, structure and construction of rock mass, terrain and physiognomy, weathering degree of rock mass and mechanical functions of seismic waves. Based on microscopic study on initial fracturing of slope caused by seismic effect, combined with two change trends which include ratio of vertical vs. horizontal peak ground acceleration corresponding to epicentral distance and enlarging effect of peak ground acceleration along slope, key inducing factor of initial slope fracturing in various area with different epicentral distance is obtained. In near-field area, i.e. epicentral distance being less than 30 km, tensile strength of rock mass is a key intrinsic factor inducing initial fracturing of slope undergoing seismic effect whereas shear strength of rock mass is the one when epicentral distance is more than 30 km. In the latter circumstance, research by means of finite difference numerical simulation test and in-lab dynamic triaxial shear test of rock shows that initial fracture begins always in the place of slope shoulder. The fact that fracture strain and shear strength which are proportional to buried depth of rock mass in the place of slope shoulder are less than other place and peak ground acceleration is enlarged in the place causes prior failure at slope shoulder. Key extrinsic factors inducing dynamic fracture of slope at different distances to epicenter have been obtained through discrete element numerical simulation on the total process of collapsing and sliding of slope triggered by Wenchuan Earthquake. Research shows that combined action of P and S seismic waves is the key factor inducing collapsing and sliding of slope at a distance less than 64 km to initial epicenter along earthquake-triggering structure. What is more, vertical tensile action of P seismic wave plays a leading role near epicenter, whereas vertical shear action of S seismic wave plays a leading role gradually with epicentral distance increasing in this range. On the other hand, single action of P seismic wave becomes the key factor inducing collapsing and sliding of slope at a distance between 64 km and 216 km to initial epicenter. Horizontal tensile action of P seismic wave becomes the key factor gradually from combined action between vertical and horizontal tensile action of P seismic wave with epicentral distance increasing in this distance range. In addition, initial failure triggered by strong earthquakes begins almost in the place of slope shoulder. However, initial failure beginning from toe of slope relates probably with gradient and rock occurrence. Finally, starting time of initial failure in slope increases usually with epicentral distance. It is perhaps that the starting time increasing is a result of attenuating of seismic wave from epicenter along earthquake-triggering structure. It is of great theoretical and practical significance for us to construct towns and infrastructure in fragile geological environment along seism tectonic belts and conduct risk management on earthquake-triggered geological hazards by referring to above conclusions.

人工源频率域大地电磁法正反演研究

The CSAMT method is playing an important role in the exploration of geothermal and the pre-exploration in tunnel construction project recently. In order to instruct the interpretation technique for the field data, the forward method from ID to 3D and inversion method in ID and 2D are developed in this paper for the artificial source magnetotelluric in frequency domain. In general, the artificial source data are inverted only after the near field is corrected on the basis of the assumption of half-homogeneous space; however, this method is not suitable for the complex structure because the assumption is not valid any more. Recently the new idea about inversion scheme without near field correction is published in order to avoid the near field correction error. We try to discuss different inversion scheme in ID and 2D using the data without near field correction.The numerical integration method is used to do the forward modeling in ID CSAMT method o The infinite line source is used in the 2D finite-element forward modeling, where the near-field effect is occurred as in the CSAMT method because of using artificial source. The pseudo-delta function is used to modeling the source distribution, which reduces the singularity when solving the finite-element equations. The effect on the exploration area is discussed when anomalous body exists under the source or between the source and exploration area; A series of digital test show the 2D finite element method are correct, the results of modeling has important significant for CSAMT data interpretation. For 3D finite-element forward modeling, the finite-element equation is derived by Galerkin method and the divergence condition is add forcedly to the forward equation, the forward modeling result of the half homogeneous space model is correct.The new inversion idea without near field correction is followed to develop new inversion methods in ID and 2D in the paper. All of the inversion schemes use the data without near field correction, which avoid introducing errors caused by near field correction. The modified grid parameter method and the layer-by-layer inversion method are joined in the ID inversion scheme. The RRI method with artificial source are developed and finite-element inversion method are used in 2D inversion scheme. The inversion results using digital data and the field data are accordant to the model and the known geology data separately, which means the inversion without near field correction is accessible. The feasibility to invert the data only in exploration area is discussed when the anomalous body exists between the source and the exploration area.

Experimental investigation of an unstable ring resonator with 90-deg beam rotation for a chemical oxygen iodine laser

We report the experimental results of an unstable ring resonator with 90-deg beam rotation for a kilowatt class chemical oxygen iodine laser (COIL). The distributions of near-field phase and far-field intensity were measured. A beam quality of 1.6 was achieved when the COIL average output power was approximately 5 kW. (C) 1999 Optical Society of America.

Large strain field near a crack tip in a rubber sheet

The distribution of stress-strain near a crack tip in a rubber sheet is investigated by employing the constitutive relation given by Gao (1997). It is shown that the crack tip field is composed of two shrinking sectors and one expanding sector. The stress state near the crack tip is in uniaxial tension. The analytical solutions are obtained for both expanding and shrinking sectors.

DPIV and Interferometry Technique for Velocity Field and Concentration Fields Measurement in Crystal Growth

The property of crystal depends seriously on the solution concentration distribution near the growth surface of a crystal. However, the concentration distributions are affected by the diffusion and convection of the solution. In the present experiment, the two methods of optical measurement are used to obtained velocity field and concentration field of NaClO₃ solution. The convection patterns in sodium chlorate (NaClO₃) crystal growth are measured by Digital Particle image Velocimetry (DPIV) technology. The 2-dimentional velocity distributions in the solution of NaClO3 are obtained from experiments. And concentration field are obtained by a Mach-Zehnder interferometer with a phase shift servo system. Interference patterns were recorded directly by a computer via a CCD camera. The evolution of velocity field and concentration field from dissolution to crystallization are visualized clearly. The structures of velocity fields were compared with that of concentration field.

Effect of non-uniform magnetic field on crystal growth by floating-zone method in microgravity

The magnetic damping effect of the non-uniform magnetic field on the floating-zone crystal growth process in microgravity is studied by numerical simulation. The results show that the non-uniform magnetic field with designed configuration can effectively reduce the flow near the free surface and then in the melt zone. At the same time, the designed magnetic field can improve the impurity concentration non-uniformity along the solidification interface. The primary principles of the magnetic field configuration design are also discussed.

Numerical simulation of magnetic field design for damping thermocapillary convection in a floating half-zone

The magnetic fields produced by electrical coils are designed for damping the the thermocapillary convection in a floating half-zone in microgravity. The fields are designed specially to reduce the flow near the free surface and then in the melt zone by adjusting the longitudinal coil positions close to the melt zone. The effects of the designed magnetic fields on reducing the flow velocity and temperature distribution non-uniformity in the melt zone are stronger than those of the case of an uniform longitudinal magnetic field obtained by numerical simulation, particularly at the melt-rod interface. It brings fundamental insights into the heat and mass transfer control at the solidification interface by the magnetic field design for crystal growth by the floating full-zone method.

Magnetic field design for floating zone crystal growth

The magnetic fields produced by electrical coils are designed for P-doped Si crystal growth in a floating full zone in microgravity environment. The fields are designed specially to reduce the how near the free surface and then in the melt zone by adjusting the coil positions near the melt zone. The effects of the designed magnetic fields on reducing the Row velocity and the non-uniformity of the concentration distribution in the melt zone are better than those of the case of a uniform longitudinal magnetic field, obtained by numerical simulation. It is expected to improve the radial macro-segregation and reduce the convection in the crystal growth at the same time by using the designed magnetic field.

Crack Tip Field and J-Integral with Strain Gradient Effect

The mode I plane strain crack tip field with strain gradient effects is presented in this paper based on a simplified strain gradient theory within the framework proposed by Acharya and Bassani. The theory retains the essential structure of the incremental version of the conventional J_2 deformation theory No higher-order stress is introduced and no extra boundary value conditions beyond the conventional ones are required. The strain gradient effects are considered in the constitutive relation only through the instantaneous tangent modulus. The strain gradient measures are included into the tangent modulus as internal parameters. Therefore the boundary value problem is the same as that in the conventional theory Two typical crack Problems are studied: (a) the crack tip field under the small scale yielding condition induced by a linear elastic mode-I K-field and (b) the complete field for a compact tension specimen. The calculated results clearly show that the stress level near the crack tip with strain gradient effects is considerable higher than that in the classical theory The singularity of the strain field near the crack tip is nearly equal to the square-root singularity and the singularity of the stress field is slightly greater than it. Consequently, the J-integral is no longer path independent and increases monotonically as the radius of the calculated circular contour decreases.

High-order asymptotic field of tensile plane-strain nonlinear crack problems

In this paper, the governing equations and the analytical method of the secondorder asymptotic field for the plane-straln crack problems of mode I have been presented. The numerical calculation has been carried out. The amplitude coefficients k2 of the second term of the asymptotic field have been determined after comparing the present results with some fine results of the finite element calculation. The variation of coefficients k2 with changes of specimen geometry and developments of plastic zone have been analysed. It is shown that the second term of the asymptotic field has significant influence on the near-crack-tip field. Therefore, we may reasonably argue that both the J-integral and the coefficient k2 can beeome two characterizing parameters for crack initiation.

Mixed-mode fracture mechanisms near the fatigue threshold of aisi-316 stainless-steel

Near threshold, mixed mode (I and II), fatigue crack growth occurs mainly by two mechanisms, coplanar (or shear) mode and branch (or tensile) mode. For a constant ratio of ΔK_I/ΔK_II the shear mode growth shows a self-arrest character and it would only start again when ΔK_I and ΔK_II are increased. Both shear crack growth and the early stages of tensile crack growth, are of a crystallographic nature; the fatigue crack proceeds along slip planes or grain boundaries. The appearance of the fracture surfaces suggest that the mechanism of crack extension is by developing slip band microcracks which join up to form a macrocrack. This process is thought to be assisted by the nature of the plastic deformation within the reversed plastic zone where high back stresses are set up by dislocation pile-ups against grain boundaries. The interaction of the crack tip stress field with that of the dislocation pile-ups leads to the formation of slip band microcracks and subsequent crack extension. The change from shear mode to tensile mode growth probably occurs when the maximum tensile stress and the microcrack density in the maximum tensile plane direction attain critical values.

The asymptotic near-tip solution for mode-iii crack in steady growth in power hardening media

Two local solutions, one perpendicular and one parallel to the direction of solar gravitational field, are discussed. The influence of gravity on the gas-dynamical process driven by the piston is discussed in terms of characteristic theory, and the flow field is given quantitatively. For a typical piston trajectory similar to the one for an eruptive prominence, the velocity of the shock front which locates ahead the transient front is nearly constant or slightly accelerated, and the width of the compressed flow region may be kept nearly constant or increased linearly, depending on the velocity distribution of the piston. Based on these results, the major features of the transient may be explained. Some of the fine structure of the transient is also shown, which may be compared in detail with observations.

Near dipole-dipole effects on the propagation of few-cycle pulse in a dense two-level medium

The propagation behaviors, which include the carrier-envelope phase, the area evolution and the solitary pulse number of few-cycle pulses in a dense two-level medium, are investigated based on full-wave Maxwell-Bloch equations by taking Lorentz local field correction (LFC) into account. Several novel features are found: the difference of the carrier-envelope phase between the cases with and without LFC can go up to pi at some location; although the area of ultrashort solitary pulses is lager than 2 pi, the area of the effective Rabi frequency, which equals to that the Rabi frequency pluses the product of the strength of the near dipole-dipole (NDD) interaction and the polarization, is consistent with the standard area theorem and keeps 2 pi; the large area pulse penetrating into the medium produces several solitary pulses as usual, but the number of solitary pulses changes at certain condition. (C) 2005 Optical Society of America.