基于量子蒙特卡罗的地震波反演方法

The primary approaches for people to understand the inner properties of the earth and the distribution of the mineral resources are mainly coming from surface geology survey and geophysical/geochemical data inversion and interpretation. The purpose of seismic inversion is to extract information of the subsurface stratum geometrical structures and the distribution of material properties from seismic wave which is used for resource prospecting, exploitation and the study for inner structure of the earth and its dynamic process. Although the study of seismic parameter inversion has achieved a lot since 1950s, some problems are still persisting when applying in real data due to their nonlinearity and ill-posedness. Most inversion methods we use to invert geophysical parameters are based on iterative inversion which depends largely on the initial model and constraint conditions. It would be difficult to obtain a believable result when taking into consideration different factors such as environmental and equipment noise that exist in seismic wave excitation, propagation and acquisition. The seismic inversion based on real data is a typical nonlinear problem, which means most of their objective functions are multi-minimum. It makes them formidable to be solved using commonly used methods such as general-linearization and quasi-linearization inversion because of local convergence. Global nonlinear search methods which do not rely heavily on the initial model seem more promising, but the amount of computation required for real data process is unacceptable. In order to solve those problems mentioned above, this paper addresses a kind of global nonlinear inversion method which brings Quantum Monte Carlo (QMC) method into geophysical inverse problems. QMC has been used as an effective numerical method to study quantum many-body system which is often governed by Schrödinger equation. This method can be categorized into zero temperature method and finite temperature method. This paper is subdivided into four parts. In the first one, we briefly review the theory of QMC method and find out the connections with geophysical nonlinear inversion, and then give the flow chart of the algorithm. In the second part, we apply four QMC inverse methods in 1D wave equation impedance inversion and generally compare their results with convergence rate and accuracy. The feasibility, stability, and anti-noise capacity of the algorithms are also discussed within this chapter. Numerical results demonstrate that it is possible to solve geophysical nonlinear inversion and other nonlinear optimization problems by means of QMC method. They are also showing that Green’s function Monte Carlo (GFMC) and diffusion Monte Carlo (DMC) are more applicable than Path Integral Monte Carlo (PIMC) and Variational Monte Carlo (VMC) in real data. The third part provides the parallel version of serial QMC algorithms which are applied in a 2D acoustic velocity inversion and real seismic data processing and further discusses these algorithms’ globality and anti-noise capacity. The inverted results show the robustness of these algorithms which make them feasible to be used in 2D inversion and real data processing. The parallel inversion algorithms in this chapter are also applicable in other optimization. Finally, some useful conclusions are obtained in the last section. The analysis and comparison of the results indicate that it is successful to bring QMC into geophysical inversion. QMC is a kind of nonlinear inversion method which guarantees stability, efficiency and anti-noise. The most appealing property is that it does not rely heavily on the initial model and can be suited to nonlinear and multi-minimum geophysical inverse problems. This method can also be used in other filed regarding nonlinear optimization.

An New Area Function for Sharp indenter Tips in Nanoindentation

A new area function is introduced and applied to a Berkovich tip in order to characterize the contact projected area between an indenter and indented material. The function can be related directly to tip-rounding, thereby having obviously physical meaning. Nanoindentation experiments are performed on a commercial Nano Indenter XPsystem. The other two area functions introduced by Oliver and Pharr and by Thurn and Cook respectively are involved in this paper for comparison. By comparison from experimental results among different area functions, the indenter tip described by the proposed area function here is very close to the experimental indenter.

Dynamic function of damage and its implications

Dynamic function of damage is the key to the problem of damage evolution of solids. In order to understand it, one must understand its mesoscopic mechanisms and macroscopic formulation. In terms of evolution equation of microdamage and damage moment, a dynamic function of damage is strictly defined. The mesoscopic mechanism underlying self-closed damage evolution law is investigated by means of double damage moments. Numerical results of damage evolution demonstrate some common features for various microdamage dynamics. Then, the dynamic function of damage is applied to inhomogeneous damage field. In this case, damage evolution rate is no longer equal to the dynamic function of damage. It is found that the criterion for damage localization is closely related to compound damage. Finally, an inversion of damage evolution to the dynamic function of damage is proposed.

Function Innovation in Combined Surface Technologies for Metals

表面技术能够显著提高材料功能而成为工程和产品设计的重要组成部分, 但前提是表面技术必须具有可设计性。为此, 需要开拓和发展表面组合加工技术, 进行创新。该技术的内涵体现了新材料与新技术、基础研究与产业化有机结合的特点, 具有重要的价值和广阔的前景。

Barium-Functionalized Multiwalled Carbon Nanotube Yarns As Low-Work-Function Thermionic Cathodes

Barium-functionalized multiwalled carbon nanotube yarns were fabricated by drawing and twisting multiwalled carbon nanotube forests through a solution containing barium nitrate. After heat activation under vacuum, the functionalized yarns were enriched in barium oxide due to the high surface-to-volume ratio of the nanotubes. The cathodes exhibited good thermionic properties, with a work function as low as 1.73-2.06 eV and thermionic current density that exceeded 185 mA/cm(2) in a field of 850 V/5 mm at 1317 K. The barium-functionalized yarns had high tensile strength of up to 420 MPa and retained strength of similar to 250 MPa after a 2 h activation process. (C) 2008 American Institute of Physics.

Semi-weight function method on computation of stress intensity factors in dissimilar materials

Semi-weight function method is developed to solve the plane problem of two bonded dissimilar materials containing a crack along the bond. From equilibrium equation, stress and strain relationship, conditions of continuity across interface and free crack surface, the stress and displacement fields were obtained. The eigenvalue of these fields is lambda. Semi-weight functions were obtained as virtual displacement and stress fields with eigenvalue-lambda. Integral expression of fracture parameters, K-I and K-II, were obtained from reciprocal work theorem with semi-weight functions and approximate displacement and stress values on any integral path around crack tip. The calculation results of applications show that the semi-weight function method is a simple, convenient and high precision calculation method.

Synchronization of modified Chua's circuit with x vertical bar x vertical bar function

This paper considers the chaos synchronization of the modified Chua's circuit with x vertical bar x vertical bar function. We firstly show that a couple of the modified Chua systems with different parameters and initial conditions can be synchronized using active control when the values of parameters both in drive system and response system are known aforehand. Furthermore, based on Lyapunov stability theory we propose an adaptive active control approach to make the states of two identical Chua systems with unknown constant parameters asymptotically synchronized. Moreover the designed controller is independent of those unknown parameters. Numerical simulations are given to validate the proposed synchronization approach.

Self-consistent kinetic modeling of electron distribution function in a low-pressure inductively coupled plasma

Based upon the spatially inhomogeneous Boltzmann equation in two-term approximation coupled with electromagnetic and fluid model analysis for the recently developed inductively coupled plasma sources, a self-consistent electron kinetic model is developed. The electron distribution function, spatial distributions of the electron density and ionization rate are calculated and discussed.

Optimization of tried wave-function in quantum monte-carlo method

A method for optimizing tried wave functions in quantum Monte Carlo method has been found and used to calculate the energies of molecules, such as H-2, Li-2, H-3+, H-3 and H-4. Good results were obtained.

Statistical-theory of multiphoton excitation of polyatomic-molecules based on the wigner function

This paper is aimed at establishing a statistical theory of rotational and vibrational excitation of polyatomic molecules by an intense IR laser. Starting from the Wigner function of quantum statistical mechanics, we treat the rotational motion in the classical approximation; the vibrational modes are classified into active ones which are coupled directly with the laser and the background modes which are not coupled with the laser. The reduced Wigner function, i.e., the Wigner function integrated over all background coordinates should satisfy an integro-differential equation. We introduce the idea of ``viscous damping'' to handle the interaction between the active modes and the background. The damping coefficient can be calculated with the aid of the well-known Schwartz–Slawsky–Herzfeld theory. The resulting equation is solved by the method of moment equations. There is only one adjustable parameter in our scheme; it is introduced due to the lack of precise knowledge about the molecular potential. The theory developed in this paper explains satisfactorily the recent absorption experiments of SF6 irradiated by a short pulse CO2 laser, which are in sharp contradiction with the prevailing quasi-continuum theory. We also refined the density of energy levels which is responsible for the muliphoton excitation of polyatomic molecules.

The gated integration technique for the accurate measurement of the autocorrelation function of speckle intensities scattered from random phase screens

A new approach based on the gated integration technique is proposed for the accurate measurement of the autocorrelation function of speckle intensities scattered from a random phase screen. The Boxcar used for this technique in the acquisition of the speckle intensity data integrates the photoelectric signal during its sampling gate open, and it repeats the sampling by a preset number, in. The average analog of the in samplings output by the Boxcar enhances the signal-to-noise ratio by root m, because the repeated sampling and the average make the useful speckle signals stable, while the randomly varied photoelectric noise is suppressed by 1/ root m. In the experiment, we use an analog-to-digital converter module to synchronize all the actions such as the stepped movement of the phase screen, the repeated sampling, the readout of the averaged output of the Boxcar, etc. The experimental results show that speckle signals are better recovered from contaminated signals, and the autocorrelation function with the secondary maximum is obtained, indicating that the accuracy of the measurement of the autocorrelation function is greatly improved by the gated integration technique. (C) 2006 Elsevier Ltd. All rights reserved.

Optical interconnection for neural networks by use of a self-imaging function

A grating-lens combination unit is developed to form a scaling self-transform function that can self-image on scale. Then an array of many such grating-lens units is used for the optical interconnection of a two-dimensional neural network, and experiments are carried out. We find that our idea is feasible, the optical interconnection system is simple, and optical adjustment is easy. (C) 1998 Optical Society of America.

Computation of the ambiguity function for circularly symmetric pupils

A method of computing the ambiguity function (AF) for a circularly symmetric pupil function is presented. The AFs of a clear aperture and two shaded apertures are considered in detail and an explicit expression for the first of these AFs is given. We explain these results in the context of the well-known optical transfer function theory and show a primary application of these computations. A good analytic approximation is also introduced, providing an alternative method for calculating the AF, in a simpler way.

Analysis of optical systems with extended depth of field using the Wigner distribution function

We describe the use of a Wigner distribution function approach for exploring the problem of extending the depth of field in a hybrid imaging system. The Wigner distribution function, in connection with the phase-space curve that formulates a joint phase-space description of an optical field, is employed as a tool to display and characterize the evolving behavior of the amplitude point spread function as a wave propagating along the optical axis. It provides a comprehensive exhibition of the characteristics for the hybrid imaging system in extending the depth of field from both wave optics and geometrical optics. We use it to analyze several well-known optical designs in extending the depth of field from a new viewpoint. The relationships between this approach and the earlier ambiguity function approach are also briefly investigated. (c) 2006 Optical Society of America.