瞬变

无论在生活中，还是在工作中，我们都会遇到各种改变。股市会发生恐慌，企业会有兼并，即便聚焦到个人，我们也无时无刻不面对着各种改变：大多数人都会告别单身最终步入婚姻的围城之中，拥有他们的家庭和孩子。所有这些，都是改变。然而，对生活或者工作中的一些改变，我们却常抱以抗拒的态度。老师们常感叹顽皮的学生“冥顽不化”，

Synthetic analysis of meteorological elements for the sea fogs occurred in northwest the Yellow Sea under low pressure control and comparison with sea fogs occurred under high pressure control in summer

Based on the principle of transient perturbation analysis,in this paper,a method to objectively determine the weather pattern formed by sea fog is provided.On the basis of the classification results,the circulation situation,divergence and vertical velocity field,and the vertical profile of temperature and humidity are synthesized and analyzed.The basic characteristics of the circulation and physical field of sea fog under low pressure control(L type sea fog)are obtained,and the results are compared with the sea fog under the control of high pressure(H type sea fog):a)L type sea fogs potential height anomaly disturbance is mainly manifested in the low layer,and its average value is-65.66 gpm,gradually weakening upward;b)L type sea fogs inversion structure is weaker than H type sea fogs when it occurs,the fog layer is thicker and the high relative humidity level is high over the fog layer,while the H type sea fogs fog layer has a relatively obvious dry layer;c)L sea fog has three layers of structure at the vertical direction.The first layer 1000-950 hPa is convergence accompanied by weak rise and subsidence,the second layer 950-850 hPa is divergence accompanied by weak subsidence,and the third layer 850 to 500hPa is gradually strengthened.While there are two layer structures of the H type sea fog.1000 hPa is divergence accompanied by weak rising and sinking movement,950-500 hPa is a uniform subsidence movement.d)Probability density statistical analysis further quantified the vertical movement of L and H type sea fog and the distribution of relative humidity in each layer.These conclusions provide an important reference for forecasting the sea fog in the northwest of the Yellow Sea under the condition of low pressure circulation in summer.

Cone-shaped source characteristics and inductance effect of transient electromagnetic method

Small multi-turn coil devices are used with the transient electromagnetic method （TEM） in areas with limited space, particularly in underground environments such as coal mines roadways and engineering tunnels, and for detecting shallow geological targets in environmental and engineering fields. However, the equipment involved has strong mutual inductance coupling, which causes a lengthy turn-off time and a deep “blind zone”. This study proposes a new transmitter device with a conical-shape source and derives the radius formula of each coil and the mutual inductance coefficient of the cone. According to primary field characteristics, results of the two fields created, calculation of the conical-shaped source in a uniform medium using theoretical analysis, and a comparison of the inductance of the new device with that of the multi-turn coil, show that inductance of the multi-turn coil is nine times greater than that of the conical source with the same equivalent magnetic moment of 926.1 A·m2. This indicates that the new source leads to a much shallower “blind zone.” Furthermore, increasing the bottom radius and turn of the cone creates a larger mutual inductance but increasing the cone height results in a lower mutual inductance. Using the superposition principle, the primary and secondary magnetic fields for a conical source in a homogeneous medium are calculated; results indicate that the magnetic behavior of the cone is the same as that of the multi-turn coils, but the transient responses of the secondary field and the total field are more stronger than those of the multi-turn coils. To study the transient response characteristics using a cone-shaped source in a layered earth, a numerical filtering algorithm is then developed using the fast Hankel transform and the improved cosine transform, again using the superposition principle. During development, an average apparent resistivity inverted from the induced electromotive force using each coil is defined to represent the comprehensive resistivity of the conical source. To verify the forward calculation method, the transient responses of H type models and KH type models are calculated, and data are inverted using a “smoke ring” inversion. The results of inversion have good agreement with original models and show that the forward calculation method is effective. The results of this study provide an option for solving the problem of a deep “blind zone” and also provide a theoretical indicator for further research.

Topographic effects on long offset transient electromagnetic response

The forward modeling procedure used in this article is formulated with the volume integral equation based on the tensor Green＇s function. The electromagnetic components responses are first calculated in the frequency domain and then transformed to the time domain by digital filtering. The valley and hill topography with a layered earth is stimulated by a horizontal electric dipole （HED） transmitter, which is common in field surveys, and the TEM responses are calculated at the transmitter and receivers. The topography effects on the long offset electromagnetic transient （LOTEM） responses are discussed in detail. The results show that both valley and hill topography has significant effect on the LOTEM measurement. If the HED is located in the bottom of a valley, the distortion of the observed anomalous field at distance is severe. A valley at the receiver locations show a strong effect but are localized in space and time. In general, hill-shaped topography shows smaller effects no matter where its located. When the topography is located between source and receivers, the influence is negligible. We conclude that the location of the source is much more important than the receivers and it is critical to put the transmitter in an open flat area in the field survey.

The analysis on IP signals in TEM response based on SVD

During transient electromagnetic method （TEM） exploration of a copper mine, we detected the late-channel TEM signal reversal phenomenon （a voltage change from positive to negative） caused by the influence of the induced polarization （IP） effect, which affects the depth and precision of the TEM detection. The conventional inversion method is inefficient because it is difficult to process the data. In this paper, the Cole-Cole model is adopted to analyze the effect of Dc resistivity, chargeability, time constant, and frequency exponent on the TEM response in an homogeneous half space model. Singular Value Decomposition （SVD） is used to invert the measured TEM data, and the Dc resistivity, chargeability, time constant and frequency exponent were extracted from the measured TEM data in the mine area. The extracted parameters are used for interpreting the detection result as a supplement. This reveals why the TEM data acquired in the area has a low resolution. It was found that the DC resistivity and time constant do not significantly change the results, however, the chargeability and frequency exponent have a significant effect. Because of these influences, the SVD method is more accurate than the conventional method in the apparent resistivity profile. The area of the copper mine is confined accurately based on the SVD inverted data. The conclusion has been verified by drill and is identical to the practical geological situation.

Numerical analysis of multicomponent responses of surface-hole transient electromagnetic method

We calculate the multicomponent responses of surface-hole transient electromagnetic method. The methods and models are unsuitable as geoelectric models of conductive surrounding rocks because they are based on regular local targets. We also propose a calculation and analysis scheme based on numerical simulations of the subsurface transient electromagnetic fields. In the modeling of the electromagnetic fields, the forward modeling simulations are performed by using the finite-difference time-domain method and the discrete image method, which combines the Gaver–Stehfest inverse Laplace transform with the Prony method to solve the initial electromagnetic fields. The precision in the iterative computations is ensured by using the transmission boundary conditions. For the response analysis, we customize geoelectric models consisting of near-borehole targets and conductive wall rocks and implement forward modeling simulations. The observed electric fields are converted into induced electromotive force responses using multicomponent observation devices. By comparing the transient electric fields and multicomponent responses under different conditions, we suggest that the multicomponent-induced electromotive force responses are related to the horizontal and vertical gradient variations of the transient electric field at different times. The characteristics of the response are determined by the varying the subsurface transient electromagnetic fields, i.e., diffusion, attenuation and distortion, under different conditions as well as the electromagnetic fields at the observation positions. The calculation and analysis scheme of the response consider the surrounding rocks and the anomalous field of the local targets. It therefore can account for the geological data better than conventional transient field response analysis of local targets.

A study on the discrete image method for calculation of transient electromagnetic fields in geological media

We conducted a study on the numerical calculation and response analysis of a transient electromagnetic field generated by a ground source in geological media. One solution method, the traditional discrete image method, involves complex operation, and its digital filtering algorithm requires a large number of calculations. To solve these problems, we proposed an improved discrete image method, where the following are realized： the real number of the electromagnetic field solution based on the Gaver-Stehfest algorithm for approximate inversion, the exponential approximation of the objective kernel function using the Prony method, the transient electromagnetic field according to discrete image theory, and closed-form solution of the approximate coefficients. To verify the method, we tentatively calculated the transient electromagnetic field in a homogeneous model and compared it with the results obtained from the Hankel transform digital filtering method. The results show that the method has considerable accuracy and good applicability. We then used this method to calculate the transient electromagnetic field generated by a ground magnetic dipole source in a typical geoelectric model and analyzed the horizontal component response of the induced magnetic field obtained from the ＂ground excitation-stratum measurement method. We reached the conclusion that the horizontal component response of a transient field is related to the geoelectric structure, observation time, spatial location, and others. The horizontal component response of the induced magnetic field reflects the eddy current field distribution and its vertical gradient variation. During the detection of abnormal objects, positions with a zero or comparatively large offset were selected for the drill- hole measurements or a comparatively long observation delay was adopted to reduce the influence of the ambient field on the survey results. The discrete image method and forward calculation results in this paper can be used as references for relevant research.

Numerical modeling of the 2D time-domain transient electromagnetic secondary field of the line source of the current excitation

To effectively minimize the electromagnetic field response in the total field solution, we propose a numerical modeling method for the two-dimensional （2D） time- domain transient electromagnetic secondary field of the line source based on the DuFort- Frankel finite-difference method. In the proposed method, we included the treatment of the earth-air boundary conductivity, calculated the normalized partial derivative of the induced electromotive force （Emf）, and determined the forward time step. By extending upward the earth-air interface to the air grid nodes and the zero-value boundary conditions, not only we have a method that is more efficient but also simpler than the total field solution. We computed and analyzed the homogeneous half-space model and the fiat layered model with high precision--the maximum relative error is less than 0.01% between our method and the analytical method--and the solution speed is roughly three times faster than the total-field solution. Lastly, we used the model of a thin body embedded in a homogeneous half-space at different delay times to depict the downward and upward spreading characteristics of the induced eddy current, and the physical interaction processes between the electromagnetic field and the underground low-resistivity body.

“Abnormal transient creep” in fine-grained Al-5356 alloy observed at low strain rates by high-resolution strain measurement

Transient creep at very low strain rates （less than 10-10 s-1） is still unclear. The traditional uniaxial creep testing is less useful due to unsatisfied resolution strain （~10-6）. To study transient creep behavior at such low strain rates, a high-resolution strain measurement using the helicoid spring specimen technique was employed in a fine-grained Al-5356 alloy at temperatures ranging from 0.47Tm to 0.74Tm （Tm： melting point）. To clarify transient creep mechanism at such low strain rates, transmission electron microscopy （TEM） was used in microstructure observation of crept specimens. The abnormal transient creep, high temperature strengthening at T〉Tp （Tp： the phase transformation temperature, 0.58Tm） or intermediate temperature softening at 0.4Tm〈T￡Tp and double-normal type （creep curves including double work-hardening stages） at T=Tp, were firstly observed. The substructure observation in a crept specimen at T=0.58Tm and e=1×10-4 shows pile-up dislocations including many small jogs with equal interval, and dislocations emitted from grain boundaries. The b-Al3Mg2 phase dissolves under the condition of testing temperatures higher than 523 K, which causes solid-solution quantity of Mg atoms to increase. Therefore, the “abnormal transient creep” may be related to the difference of solid solution strengthening caused by phase change during the creep tests.

Multi-field dynamic modeling and numerical simulation of aluminum alloy resistance spot welding

In order to explore the influence of welding parameters and to investigate the Al alloy (AA) nugget formation process, a comprehensive model involving electrical-thermal-mechanical and metallurgical analysis was established to numerically display the resistance spot welding (RSW) process within multiple fields and understand the AA-RSW physics. A multi-disciplinary finite element method (FEM) framework and a empirical sub-model were built to analyze the affecting factors on weld nugget and the underlying nature of welding physics with dynamic simulation procedure. Specifically, a counter-intuitive phenomenon of the resistance time-variation caused by the transient inverse virtual variation (TIVV) effect was highlighted and analyzed on the basis of welding current and temperature distribution simulation. The empirical model describing the TIVV phenomenon was used for modifying the dynamic resistance simulation during the AA spot welding process. The numerical and experimental results show that the proposed multi-field FEM model agrees with the measured AA welding feature, and the modified dynamic resistance model captures the physics of nugget growth and the electrical-thermal behavior under varying welding current and fluctuating heat input.

Integrated identification method of rheological model of sandstone in Sanmenxia bauxite

Based on the uniaxial compression creep experiments conducted on bauxite sandstone obtained from Sanmenxia,typical creep experiment curves were obtained.From the characteristics of strain component of creep curves,the creep strain is composed of instantaneous elastic strain,ε（me）,instantaneous plastic strain,ε（mp）,viscoelastic strain,ε（ce）,and viscoplastic strain,ε（cp）.Based on the characteristics of instantaneous plastic strain,a new element of instantaneous plastic rheology was introduced,instantaneous plastic modulus was defined,and the modified Burgers model was established.Then identification of direct screening method in this model was completed.According to the mechanical properties of rheological elements,one- and three-dimensional creep equations in different stress levels were obtained.One-dimensional model parameters were identified by the method of least squares,and in the process of computation,Gauss-Newton iteration method was applied.Finally,by fitting the experimental curves,the correctness of direct method model was verified,then the examination of posterior exclusive method of the model was accomplished.The results showed that in the improved Burgers models,the rheological characteristics of sandstone are embodied properly,microscopic analysis of creep curves is also achieved,and the correctness of comprehensive identification method of rheological model is verified.

Seismic attribute extraction based on HHT and its application in a marine carbonate area

The Hilbert-Huang transform（HHT） is a new analysis method suitable for nonlinear and non-stationary signals.It is very appropriate to seismic signals because they show both non-stationary and nonlinear characteristics.We first introduce the realization of HHT empirical mode decomposition（EMD） and then comparatively analyze three instantaneous frequency algorithms based on intrinsic mode functions（IMF） resulting from EMD,of which one uses the average instantaneous frequency of two sample intervals having higher resolution which can determine that the signal frequency components change with time.The method is used with 3-D poststack migrated seismic data of marine carbonate strata in southern China to effectively extract the three instantaneous attributes.The instantaneous phase attributes of the second intrinsic mode functions（IMF2） better describe the reef facies of the platform margin and the IMF2 instantaneous frequency attribute has better zoning.Combining analysis of the three IMF2 instantaneous seismic attributes and drilling data can identify the distribution of sedimentary facies well.

Instantaneous rotation speed measurement and error analysis for variable speed hydraulic system

In order to monitor the working state of piston motor and measure its instantaneous rotation speed accurately, the measuring principle and method of instantaneous rotation speed based on industrial personal computer and data acquisition card are introduced, and the major error source, influence mechanism and processing method of data quantization error are dis- cussed. By means of hybrid programming approach of LabVIEW and MATLAB, the instantaneous rotation speed measurement system for the piston motor in variable speed hydraulic system is designed. The simulation and experimental results show that the designed instantaneous speed measurement system is feasible. Furthermore, the sampling frequency has an important influ- ence on the instantaneous rotation speed measurement of piston motor and higher sampling frequency can lower quantization er- ror and improve measurement accuracy.

Light-Induced Changes in a-Si∶H Films Studied by Transient Photoconductivity

Light induced changes in a-Si∶H films are investigated by transient photoconductivity.The transient photoconductivity decay data can neither be fit well by common power-law for transient photocurrent in amorphous semiconductors,nor by stretched exponential rule for transient decay from the steady state in photoconductivity.Instead,the data are fit fairly well with a sum of two exponential functions.The results show that the long time decay is governed by deep traps rather than band tail states,and two different traps locating separately at 0.52 and 0.59eV below E _c are responsible for the two exponential functions.They are designated as negatively charged dangling bond D - centers.The light-induced changes in photoconductivity are attributed mainly to the decrease in electron lifetime caused by the increase of recombination centers after light soaking.

THE BOUNDARY ELEMENT ANALYSIS PROGRAM FOR 3-D TRANSIENT DYNAMIC FIELD

The matrix expression for the 3 D transient dynamic boundary integral equation in Laplace transform space is obtained and the degenerative element method has been implemented to treat the kernel function over the singular element. In the computer program BEMTDY the Koizumi′s numerical inversion method is used and three examples of the 3 D vibrated foundation under harmonic forces and the influence with both adjacent foundations are studied.

Three-dimensional numerical modeling of fullspace transient electromagnetic responses of water in goaf

The full-space transient electromagnetic response of water-filled goaves in coal mines were numerically modeled. Traditional numerical modeling methods cannot be used to simulate the underground full-space transient electromagnetic field. We used multiple transmitting loops instead of the traditional single transmitting loop to load the transmitting loop into Cartesian grids. We improved the method for calculating the z-component of the magnetic field based on the characteristics of full space. Then, we established the full- space 3D geoelectrical model using geological data for coalmines. In addition, the transient electromagnetic responses of water-filled goaves of variable shape at different locations were simulated by using the finite-difference time-domain （FDTD） method. Moreover, we evaluated the apparent resistivity results. The numerical modeling results suggested that the resistivity differences between the coal seam and its roof and floor greatly affect the distribution of apparent resistivity, resulting in nearly circular contours with the roadway head at the center. The actual distribution of apparent resistivity for different geoelectrical models of water in goaves was consistent with the models. However, when the goal water was located in one side, a false low-resistivity anomaly would appear on the other side owing to the full-space effect but the response was much weaker. Finally, the modeling results were subsequently confirmed by drilling, suggesting that the proposed method was effective.

Continuous TDEM for monitoring shale hydraulic fracturing

Monitoring and delineating the spatial distribution of shale fracturing is fundamentally important to shale gas production. Standard monitoring methods, such as time-lapse seismic, cross-well seismic and micro-seismic methods, are expensive, time- consuming, and do not show the changes in the formation with time. The resistivities of hydraulic fracturing fluid and reservoir rocks were measured. The results suggest that the injection fluid and consequently the injected reservoir are characterized by very low resistivity and high chargeability. This allows using of the controlled-source electromagnetic method （CSEM） to monitor shale gas hydraulic fracturing. Based on the geoelectrical model which was proposed according to the well-log and seismic data in the test area the change rule of the reacted electrical field was studied to account for the change of shale resistivity, and then the normalized residual resistivity method for time lapse processing was given. The time-domain electromagnetic method （TDEM） was used to continuously monitor the shale gas fracturing at the Fulin shale gas field in southern China. A high-power transmitter and multi-channel transient electromagnetic receiver array were adopted. 9 h time series of Ex component of 224 sites which were laid out on the surface and over three fracturing stages of a horizontal well at 2800 m depth was recorded. After data processing and calculation of the normalized resistivity residuals, the changes in the Ex signal were determined and a dynamic 3D image of the change in resistivity was constructed. This allows modeling the spatial distribution of the fracturing fluid. The model results suggest that TDEM is promising for monitoring hydraulic fracturing of shale.

Experiment in metal disturbance during advanced detection using a transient electromagnetic method in coal mines

Tunneling machines, or excavators, are large and good conductors and affect the reliability of data gathering and interpretation in advanced detection using transient electromagnetic methods. In our experiment, we used a coincident-loop and central loop type of configuration, where the coil plane l） vertical to and 2） parallel to the working face. A SIROTEM instrument at different locations was used to observe the transient electromagnetic responses of the excavator and to analyze the response amplitudes. The result shows that the tunneling machine affects the advanced detection data and is related to the way the coil is coupled. When the excavator is 6 m from the observatory, the interference of tunneling machine can be ignored.

Aerodynamic effects on a railway tunnel with partially changed cross-sectional area

The present study numerically explored the aerodynamic performance of a novel railway tunnel with a partially reduced cross-section.The impact of the reduction rate of the tunnel cross-section on wave transmissions was analyzed based on the three-dimensional,unsteady,compressible,and RNG k-εturbulence model.The results highlight that the reduction rate(S)most affects pressure configurations at the middle tunnel segment,followed by the enlarged segments near access,and finally the exit.The strength of the newly generated compression wave at the tunnel junction where the cross-section abruptly changes increases exponentially with the decrease of the cross-sectional area.The maximum peak-to-peak pressureΔP on the tunnel and train surface for non-uniform tunnels is reduced by 10.7%and 13.8%,respectively,compared with those of equivalent uniform tunnels.Overall,the economic analysis suggests that the aerodynamic performance of the developed tunnel prototype surpasses those conventional tunnels based on the same excavated volume.

Design of ground-airborne TEM three-component air-core coil sensor

Transient electromagnetic method (TEM),as a non-seismic geophysical exploration mainstream electromagnetic method,is widely used in oil,gas,mineral and other underground resources exploration areas. The coil sensor is generally used to collect data. In view of the problems of incomplete information of the abnormal body and the data loss in the existing TEM single-component coil sensor,a three-component TEM coil sensor is designed. By analyzing the relationship between sensor sensitivity and coil structure parameters,the coil structure and turns are designed. By analyzing the frequency response characteristics of the TEM magnetic field sensor,the signal distortion is reduced by using the under-damped matching mode. By analyzing the distribution of various noise sources of the magnetic sensor,the appropriate amplifier is selected to reduce the background noise. Finally,a three-component TEM induction magnetic field sensor is designed. The weight of the sensor is controlled at 3.2 kg and the working frequency is 10 mHz-10 kHz. The background noises of X and Y components probably keep in 1.5×10^-8 V/ Hz and sensitivities are 8.4 and 9.8 nT/s,respectively,the background noise of vertical component is 2.1× 10^-7 V/ Hz and sensitivity is 18.5 nT/s. Compared with the existing single-component TEM receiving magnetic field sensor,the designed sensor realizes the signal acquisition of three components. Without too much increase in volume and total weight,it improves the sensitivity of the sensor and reduces the background noise,thus the signal-to-noise ratio (SNR) of the signal is improved.