A finite-element algorithm for modeling variably saturated flows

A general numerical algorithm in the context of finite element scheme is developed to solve Richards’ equation, in which a mass-conservative, modified head based scheme (MHB) is proposed to approximate the governing equation, and mass-lumping techniques are used to keep the numerical simulation stable. The MHB scheme is compared with the modified Picard iteration scheme (MPI) in a ponding infiltration example. Although the MHB scheme is a little inferior to the MPI scheme in respect of mass balance, it is superior in convergence character and simplicity. Fully implicit, explicit and geometric average conductivity methods are performed and compared, the first one is superior in simulation accuracy and can use large time-step size, but the others are superior in iteration efficiency. The algorithm works well over a wide variety of problems, such as infiltration fronts, steady-state and transient water tables, and transient seepage faces, as demonstrated by its performance against published experimental data. The algorithm is presented in sufficient detail to facilitate its implementation.

Mesoscopic simulation of the impregnating process of unidirectional fibrous preform in resin transfer molding

The resin transfer molding has gained popularity in the preparation of fiber-reinforced polymer-matrix composites because of its high efficiency and low pollution. The non-uniform inter-tow and intra-tow flows are regarded as the reason of void formation in RTM. According to the process characteristics, the axisymmetric model was developed to study the interaction between the flow in the inter-tow space and that in the intra-tow space. The flow behavior inside the fiber tows was formulated using Brinkman's equation, while that in the open space around the fiber tows was formulated by Stokes' equation. The volume of fluid (VOF) method was applied to track the flow front, and the effects of filling velocity, resin viscosity, inter-tow dimension and intra-tow permeability on fluid pressure and flow front were analyzed. The results show that the flow front difference between the inter-tow and intra-tow becomes larger with the decrease of intra-tow permeability, as well as the increase of filling velocity and inter-tow dimension.

基于介电泳的电极阵列电场仿真研究

介电泳方法被广泛地应用于微纳颗粒的分离和操纵中,实现介电泳操作的关键是设计满足所需电场分布的电极阵列。针对目前在微电极阵列设计中尚缺乏简单有效的电场解析方法的现状,提出一种基于格林公式的电极阵列电场的解析方法。首先介绍了传统介电泳和行波介电泳的概念和计算模型,分析了介电泳过程与电极上所施加的交变电压的频率和幅度的关系,然后在确立电极电势的边界条件的基础上,采用基于格林公式的电场解析方法,建立了非均匀电场的解析模型,得出不同条件下的电极阵列电场分布的仿真结果,最后利用FEMLAB有限元仿真软件对解析模型进行了对比仿真,验证了该解析模型的可行性。基于格林公式的电场解析求解方法能够有效地提高电极阵列设计中的针对性以及缩短电极设计的时间。

消除采集脚印的处理方法研究

Seismic technique is in the leading position for discovering oil and gas trap and searching for reserves throughout the course of oil and gas exploration. It needs high quality of seismic processed data, not only required exact spatial position, but also the true information of amplitude and AVO attribute and velocity. Acquisition footprint has an impact on highly precision and best quality of imaging and analysis of AVO attribute and velocity. Acquisition footprint is a new conception of describing seismic noise in 3-D exploration. It is not easy to understand the acquisition footprint. This paper begins with forward modeling seismic data from the simple sound wave model, then processes it and discusses the cause for producing the acquisition footprint. It agreed that the recording geometry is the main cause which leads to the distribution asymmetry of coverage and offset and azimuth in different grid cells. It summarizes the characters and description methods and analysis acquisition footprint’s influence on data geology interpretation and the analysis of seismic attribute and velocity. The data reconstruct based on Fourier transform is the main method at present for non uniform data interpolation and extrapolate, but this method always is an inverse problem with bad condition. Tikhonov regularization strategy which includes a priori information on class of solution in search can reduce the computation difficulty duo to discrete kernel condition disadvantage and scarcity of the number of observations. The method is quiet statistical, which does not require the selection of regularization parameter; and hence it has appropriate inversion coefficient. The result of programming and tentat-ive calculation verifies the acquisition footprint can be removed through prestack data reconstruct. This paper applies migration to the processing method of removing the acquisition footprint. The fundamental principle and algorithms are surveyed, seismic traces are weighted according to the area which occupied by seismic trace in different source-receiver distances. Adopting grid method in stead of accounting the area of Voroni map can reduce difficulty of calculation the weight. The result of processing the model data and actual seismic demonstrate, incorporating a weighting scheme based on the relative area that is associated with each input trace with respect to its neighbors acts to minimize the artifacts caused by irregular acquisition geometry.

南海北部陆缘深水区构造－热演化与烃源岩成熟度史研究

Abstract In order to provide basic data for evaluation of the petroleum potential in the deep water area of the northern margin of the South China Sea (SCS), present-day thermal regime and basin tectonothermal evolution are reconstructed and the maturation history of the Cenozoic major source rocks in the study area is derived. The present-day geothermal regime in the deep water area of the northern margin of SCS is defined according to the geothermal gradient, thermal properties and heat flow data. Tectonic subsidence history is reconstructed based on borehole and seismic data, and accordingly the stretching episodes are determined from the subsidence pattern. Heat flow history in the deep water area of the northern margin of SCS is estimated on a finite time, laterally non-uniform and multi-episode stretching model. Maturation history of the main source rocks in the study area is estimated through EASYRo% kinetic model and thermal history, and the potential of petroleum in the deep water area of the northern margin of SCS is evaluated based on the data above. The results show that the present-day geothermal regime in the deep water area of the northern margin of SCS is characterized by “hot basin” with high geothermal gradient (39.1±7.4℃/km) and high heat flow (77.5±14.8 mW/m2), and that the Qiongdongnan Basin (QDNB) underwent three stretching episodes and consequently suffered three heating episodes (Eocene, Oligocene and Pliocene time) with highest paleo-heat flow of 65～90 mW/m2 at the end of the Pliocene, that the Pearl River Mouth Basin (PRMB) two stretching and two heating episodes (Eocene, Oligocene time) with highest paleo-heat flow of 60～70 mW/m2 at the end of the Oligocene, and that the source rocks matured drastically responding to the heating episodes. There are four hydrocarbon generation kitchens in the deep water area of the northern margin of SCS which are favor of its bright petroleum perspective. Tectonothermal analysis indicates that the present-day geothermal regime which is characterized with “hot basin” in the deep water area of the PRMB resulted mainly from the Cenozoic stretching as well as faulting and magmatic activities during the Neotectonic period, and that the Pliocene heating episode of the QDNB is coupled with the transition from sinistral to dextral gliding of the Red Rive fault, and that the deep water basins in the northern margin of SCS are typical of multiple rifting which caused multi-episode heating process.

复杂介质中地震波场数值模拟研究——基于广义正交多项式迭积微分算子法

In this paper, we propose a new numerical modeling method – Convolutional Forsyte Polynomial Differentiator (CFPD), aimed at simulating seismic wave propagation in complex media with high efficiency and accuracy individually owned by short-scheme finite differentiator and general convolutional polynomial method. By adjusting the operator length and optimizing the operator coefficient, both global and local informations can be easily incorporated into the wavefield which is important to invert the undersurface geological structure. The key issue in this paper is to introduce the convolutional differentiator based on Forsyte generalized orthogonal polynomial in mathematics into the spatial differentiation of the first velocity-stress equation. To match the high accuracy of the spatial differentiator, this method in the time coordinate adopts staggered grid finite difference instead of conventional finite difference to model seismic wave propagation in heterogeneous media. To attenuate the reflection artifacts caused by artificial boundary, Perfectly Matched Layer (PML) absorbing boundary is also being considered in the method to deal with boundary problem due to its advantage of automatically handling large-angle emission. The PML formula for acoustic equation and first-order velocity-stress equation are also derived in this paper. There is little difference to implement the PML boundary condition in all kind of wave equations, but in Biot media, special attenuation factors should be taken. Numerical results demonstrate that the PML boundary condition is better than Cerjan absorbing boundary condition which makes it more suitable to hand the artificial boundary reflection. Based on the theories of anisotropy, Biot two-phase media and viscous-elasticity, this paper constructs the constitutive relationship for viscous-elastic and two-phase media, and further derives the first-order velocity-stress equation for 3D viscous-elastic and two-phase media. Numerical modeling using CFPD method is carried out in the above-mentioned media. The results modeled in the viscous-elastic media and the anisotropic pore elastic media can better explain wave phenomena of the true earth media, and can also prove that CFPD is a useful numerical tool to study the wave propagation in complex media.

粘弹性介质中的地震波传播与波形反演研究

The real earth is far away from an ideal elastic ball. The movement of structures or fluid and scattering of thin-layer would inevitably affect seismic wave propagation, which is demonstrated mainly as energy nongeometrical attenuation. Today, most of theoretical researches and applications take the assumption that all media studied are fully elastic. Ignoring the viscoelastic property would, in some circumstances, lead to amplitude and phase distortion, which will indirectly affect extraction of traveltime and waveform we use in imaging and inversion. In order to investigate the response of seismic wave propagation and improve the imaging and inversion quality in complex media, we need not only consider into attenuation of the real media but also implement it by means of efficient numerical methods and imaging techniques. As for numerical modeling, most widely used methods, such as finite difference, finite element and pseudospectral algorithms, have difficulty in dealing with problem of simultaneously improving accuracy and efficiency in computation. To partially overcome this difficulty, this paper devises a matrix differentiator method and an optimal convolutional differentiator method based on staggered-grid Fourier pseudospectral differentiation, and a staggered-grid optimal Shannon singular kernel convolutional differentiator by function distribution theory, which then are used to study seismic wave propagation in viscoelastic media. Results through comparisons and accuracy analysis demonstrate that optimal convolutional differentiator methods can solve well the incompatibility between accuracy and efficiency, and are almost twice more accurate than the same-length finite difference. They can efficiently reduce dispersion and provide high-precision waveform data. On the basis of frequency-domain wavefield modeling, we discuss how to directly solve linear equations and point out that when compared to the time-domain methods, frequency-domain methods would be more convenient to handle the multi-source problem and be much easier to incorporate medium attenuation. We also prove the equivalence of the time- and frequency-domain methods by using numerical tests when assumptions with non-relaxation modulus and quality factor are made, and analyze the reason that causes waveform difference. In frequency-domain waveform inversion, experiments have been conducted with transmission, crosshole and reflection data. By using the relation between media scales and characteristic frequencies, we analyze the capacity of the frequency-domain sequential inversion method in anti-noising and dealing with non-uniqueness of nonlinear optimization. In crosshole experiments, we find the main sources of inversion error and figure out how incorrect quality factor would affect inverted results. When dealing with surface reflection data, several frequencies have been chosen with optimal frequency selection strategy, with which we use to carry out sequential and simultaneous inversions to verify how important low frequency data are to the inverted results and the functionality of simultaneous inversion in anti-noising. Finally, I come with some conclusions about the whole work I have done in this dissertation and discuss detailly the existing and would-be problems in it. I also point out the possible directions and theories we should go and deepen, which, to some extent, would provide a helpful reference to researchers who are interested in seismic wave propagation and imaging in complex media.

复杂裂隙岩体统计断裂本构关系与变形参数研究

Rock mass is widely recognized as a kind of geologic body which consists of rock blocks and discontinuities. The deformation and failure of rock mass is not only determined by rock block，but also by discontinuity which is virtually more important. Mutual cutting and combination of discontinuities controlled mechanical property of rock mass. The complex cutting of discontinuities determine the intense anisotropy on mechanical property of rock mass，especially under the effect of ground stress. Engineering practice has show that the brittle failure of hard rock always occurs when its working stress is far lower than the yield strength and compressive strength，the failure always directly related to the fracture propagation of discontinuities. Fracture propagation of discontinuities is the virtue of hard rock’s failure. We can research the rock mass discontinuous mechanical properties precisely by the methods of statistical analysis of discontinuities and Fracture Mechanics. According to Superposition Principle in Fracture Mechanics，A Problem or C Problem could be chosen to research. Problem A mainly calculates the crack-tip stress field and displacement field on internal discontinuities by numerical method. Problem C calculate the crack-tip stress field and displacement field under the assumption of that the mainly rock mass stress field has been known. So the Problem C avoid the complex mutual interference of stress fields of discontinuities，which is called crack system problem in Fracture Mechanics. To solve Problem C, field test on stress field in the rock mass is needed. The linear Superposition of discontinuities strain energies are Scientific and Rational. The difference of Fracture Mechanics between rock mass and other materials can mostly expression as：other materials Fracture Mechanics mostly face the problem A，and can’t avoid multi-crack puzzle, while the Rock mass Fracture Mechanics answer to the Problem C. Problem C can avoid multi-discontinuities mutual interference puzzle via the ground stress test. On the basis of Problem C, Fracture Mechanics could be used conveniently in rock mass. The rock mass statistics fracture constitutive relations, which introduced in this article, are based on the Problem C and the Discontinuity Strain Energy linear superposition. This constitutive relation has several merits: first, it is physical constitutive relation rather than empirical; second, it is very fit to describe the rock mass anisotropy properties; third, it elaborates the exogenous factors such as ground stress. The rock mass statistics fracture constitutive relation is the available approach to answer to the physical, anisotropic and ground stress impacted rock mass problems. This article stand on the foundation of predecessor’s statistics fractures constitutive relation, and improved the discontinuity distributive function. This article had derived the limitation of negative exponential distribution in the course of regression analysis, and advocated to using the two parameter negative exponential distribution for instead. In order to solve the problems of two-dimension stability on engineering key cross-sectional view in rock mass, this article derived the rock mass planar flexibility tensor, and established rock mass two-dimension penetrate statistics fracture constitutive relation on the basis of penetrate fracture mechanics. Based on the crack tip plasticity research production of penetrate fracture, for example the Irwin plasticity equifinality crack, this article established the way to deal with the discontinuity stress singularity and plastic yielding problem at discontinuity tip. The research on deformation parameters is always the high light region of rock mass mechanics field. After the dam foundation excavation of XiaoWan hydroelectric power station, dam foundation rock mass upgrowthed a great deal of unload cracks, rock mass mechanical property gotten intricacy and strong anisotropy. The dam foundation rock mass mostly upgrowthed three group discontinuities： the decantation discontinuity, the steep pitch discontinuity, and the schistosity plane. Most of the discontinuities have got partial unload looseness. In accordance with ground stress field data, the dam foundation stress field greatly non-uniform, which felled under the great impaction of tectonic stress field, self-weight stress field, excavation geometric boundary condition, and excavation, unload. The discontinuity complexity and stress field heterogeneity, created the rock mass mechanical property of dam foundation intricacy and levity. The research on the rock mass mechanics, if not take every respected influencing factor into consideration as best as we can, major errors likely to be created. This article calculated the rock mass elastic modulus that after Xiao Wan hydroelectric power station dam foundation gutter excavation finished. The calculation region covered possession monolith of Xiao Wan concrete double-curvature arch dam. Different monolith were adopted the penetrate fracture statistics constitutive relation or bury fracture statistics constitutive relation selectively. Statistics fracture constitutive relation is fit for the intensity anisotropy and heterogeneity rock mass of Xiao Wan hydroelectric power station dam foundation. This article had contrastive analysis the statistics fracture constitutive relation result with the inclined plane load test actual measurement elastic modulus and RMR method estimated elastic modulus, and find that the three methods elastic modulus have got greatly comparability. So, the statistics fracture constitutive relations are qualified for trust. Generally speaking，this article had finished following works based on predecessors job: “Argumentation the C Problems of superposition principle in Fracture Mechanics, establish two-dimension penetrate statistics fracture constitutive relation of rock mass, argue the negative exponential distribution limitation and improve it, improve of the three-dimension berry statistics fracture constitutive relation of rock mass, discontinuity-tip plastic zone isoeffect calculation, calculate the rock mass elastic modulus on two-dimension cross-sectional view”. The whole research clue of this article inherited from the “statistics rock mass mechanics” of Wu Faquan（1992）.

各向异性与非线性介质中弹性波场传播特征研究

The topic of this study is about the propagation features of elastic waves in the anisotropic and nonlinear media by numerical methods with high accuracy and stability. The main achievements of this paper are as followings: Firstly, basing on the third order elastic energy formula, principle of energy conservation and circumvolved matrix method, we firstly reported the equations of non-linear elastic waves with two dimensions and three components in VTI media. Secondly, several conclusions about some numerical methods have been obtained in this paper. Namely, the minimum suitable sample stepth in space is about 1/8-1/12 of the main wavelength in order to distinctly reduce the numerical dispersion resulted from the numerical mehtod, at the same time, the higher order conventional finite difference (CFD) schemes will give little contribution to avoid the numerical solutions error accumulating with time. To get the similar accuracy with the fourth order center finite difference method, the half truncation length of SFFT should be no less than 7. The FDFCT method can present with the numerical solutions without obvious dispersion when the paprameters of FCT is suitable (we think they should be in the scope from 0.0001 to 0.07). Fortunately, the NADM method not only can reported us with the higher order accuracy solutions (higher than that of the fourth order finite difference method and lower than that of the sixth order finite difference method), but also can distinctly reduce the numerical dispersion. Thirdly, basing on the numerial and theoretical analysis, we reported such nonlinear response accumulating with time as waveform aberration, harmonic generation and resonant peak shift shown by the propagation of one- and two-dimensional non-linear elasticwaves in this paper. And then, we drew the conclusion that these nonlinear responses are controlled by the product between nonlinear strength (SN) and the amplitude of the source. At last, the modified FDFCT numerical method presented by this paper is used to model the two-dimensional non-linear elastic waves propagating in VTI media. Subsequently, the wavelet analysis and polarization are adopted to investigate and understand the numerical results. And then, we found the following principles (attention: the nonlinear strength presented by this paper is weak, the thickness of the -nonlinear media is thin (200m), the initial energy of the source is weak and the anisotropy of the media is weak too): The non-linear response shown by the elastic waves in VTI media is anisotropic too; The instantaneous main frequency sections of seismic records resulted from the media with a non-linear layer have about 1/4 to 1/2 changes of the initial main frequency of source with that resulted from the media without non-linear layer; The responses shown by the elasic waves about the anisotropy and nonlinearity have obvious mutual reformation, namely, the non-linear response will be stronger in some directions because of the anisotropy and the anisotropic strength shown by the elastic waves will be stronger when the media is nonlinear.

网络成瘾者心率变异性频谱特征及人格相关因素的研究

The aim of this research is to explore heart rate variability frequency field characteristics and personality influential factors of internet addicts by experimentation and questionnaire. Two studies were carried out: study 1 was to explore the pathological and mental mechanism of internet addicts by heart rate variability physiological index and EPQ and internet addiction scale. Study 2 was to compare the personality and mental characteristics between internet addicts and Non-addicts. The testees were 30 internet addict schoolboys who were in-patient and 43 schoolboys who accorded with qualification from grade 2 in senior high school and sophomore. It is found that: 1、 Internet addicts have obviously lower HFNU than Non-addicts， but have obviously higher LFNU and LF/HF ratio than Non-addicts. Internet addicts have dysfunction in their sympathetic and parasympathetic system; 2、 Internet addicts and Non-addicts have no significant difference in their EPQ except their lying point, but the nervous characteristics of EPQ of internet addicts have influence on their equipoise of parasympathetic system, that is to say when the score of nervous characteristics of EPQ become higher, equipoise of parasympathetic system become worse and worse. However the EPQ personality characteristics of non-addicts have no influence on their sympathetic and parasympathetic system. 3、 The extent of internet addiction of the internet addicts is independent of their equipoise of parasympathetic system, but non-addicts use internet more time, their equipoise of parasympathetic system become lower, there is the significant difference in the neurophysiology between internet addicts and non-addicts. That is to say internet environment is safe to most adolescent, they can make use of internet environment accurately, but internet addicts, their physiological and psychological level has changed to a certain, need to be treated.； 4、 Serious internet game addicts have the metal characteristics of low social support, low purpose in life and low adventure; 5、 The objective support、support utilization of social support questionnaire and sensation seeking characteristics have prognosticative function for internet addiction degree; 6、 Serious internet game addicts have the metal characteristics of low social support, low lying particularity, they are inclined to self abandonment. Serious internet game addicts of low N characteristic have high sensation seeking characteristic, and at the same E personality foundation, they also seek new stimulus with higher intension. There is a prompt that we should pay attention to internet addicts’ personality so that obtain better curative effect for internet addiction therapy.

同位网格摄动有限体积格式求解浮力驱动方腔流

利用对流扩散方程的摄动有限体积格式,在Rayleigh数从10~3到10~8的范围内对浮力驱动方腔流动问题作了数值模拟.对流扩散方程的摄动有限体积格式具有一阶迎风格式的简洁形式,使用相同的基点,重构近似精度高,特别是两相邻控制体中心到公共界面的距离相等或不相等,PFV格式公式相同等优点.在数值模拟中,无论均匀网格还是非均匀网格均获得与DSC方法、自适应有限元法、多重网格法等Benchmark解相符较好的数值结果,证明UPFV格式对高Rayleigh数对流传热问题的适用性和有效性.

Influence of Local Magnetic Fields on P-Doped Si Floating Zone Melting Crystal Growth in Microgravity

A two-dimensional axisymmetric numerical model is presented to study the influence of local magnetic fields on P-doped Si floating zone melting crystal growth in microgravity. The model is developed based on the finite difference method in a boundary-fitted curvilinear coordinate system. Extensive numerical simulations are carried out, and parameters studied include the curved growth interface shape and the magnetic field configurations. Computed results show that the local magnetic field is more effective in reducing the impurity concentration nonuniformity at the growth interface in comparison with the longitudinal magnetic field. Moreover, the curved growth interface causes more serious impurity concentration nonuniformity at the growth interface than the case with a planar growth interface.

The Liquefaction and Displacement of Highly Saturated Sand under Water Pressure Oscillation

In order to investigate the characteristics of water wave induced liquefaction in highly saturated sand in vertical direction, a one-dimensional model of highly saturated sand to water pressure oscillation is presented based oil the two-phase continuous media theory. The development of the effective stresses and the liquefaction thickness are analyzed. It is shown that water pressure oscillating loading affects liquefaction severely and the developing rate of liquefaction increases with the decreasing of the sand strength or the increasing of the loading strength. It is shown also that there is obvious phase lag in the sand Column. If the sand permeability is non-uniform, the pore pressure and the strain rise sharply at which the smallest permeability occurs. This solution may explain Why the fracture occurs in the sand column in some conditions.

Experimental Studies of the Material Properties of the Forewing of Cicada (Homoptera, Cicadidae)

Detailed investigations on the structural and mechanical properties of the forewing of the cicada were carried out. Measurement of the structures of the wings showed that the thickness of the membrane of each cell and the diameter of each vein were non-uniform in both the longitudinal and transverse directions, and their means were approximately 12.2 and 133.3 mum, respectively. However, the aspect ratios of the wings and the bodies were quite uniform and were approximately equal to 2.98 and 2.13, respectively. Based on the measured thickness, mass and area of the membranes of the cells, the mean density and the mean area density of the wing were approximately 2.3 g cm(-3) and 2.8 x 10(-3) g cm(-2), respectively. In addition, the diameters of the veins of the wings, including the diameters of the holes in the vein of the leading edge, were examined. The mechanical properties of the wing were investigated separately by nanoindentation and tensile testing. The results indicated that the mean Young's modulus, hardness and yield stress of the membranes of the wings were approximately 3.7 Gpa, 0.2 Gpa and 29 Mpa, respectively, and the mean Young's modulus and strength of the veins along the direction of the venation of wings were approximately 1.9 Gpa and 52 Mpa, respectively. Finally, the relevant results were briefly analyzed and discussed, providing a guideline to the biomimetic design of the aerofoil materials of micro air vehicles.

高超声速喷管非平衡尺度效应的数值分析

采用轴对称热化学非平衡全N-S方程,数值分析了高超声速喷管流场中非平衡引起的尺度效应.并给出计算条件下喷管轴线上的物理量和化学组元分布.计算结果表明,在喉部下游不远处流动趋近于化学冻结和热力冻结状态.流场(特别是核心无粘流区)的特性不仅取决于流动尺度的相对值,也取决于其绝对值.