Control mechanisms for gas hydrate production by depressurization in different scale hydrate reservoirs

The methane hydrate was formed in a pressure vessel 38 mm in id and 500 mm in length. Experimental works on gas production from the hydrate-bearing core by depressurization to 0.1, 0.93, and 1.93 MPa have been carried out. The hydrate reservoir simulator TOUGH-Fx/Hydrate was used to simulate the experimental gas production behavior, and the intrinsic hydration dissociation constant (K-0) fitted for the experimental data was on the order of 104 mol m(-2) Pa-1 s(-1), which was one order lower than that of the bulk hydrate dissociation. The sensitivity analyses based on the simulator have been carried out, and the results suggested that the hydrate dissociation kinetics had a great effect on the gas production behavior for the laboratory-scale hydrate-bearing core. However for a field-scale hydrate reservoir, the flow ability dominated the gas production behavior and the effect of hydrate dissociation kinetics on the gas production behavior could be neglected.

Tunable parametric fluorescence using a single crystal

In this paper, we investigate the mechanism of tunable parametric superfluorescence (PS) based on the second harmonic generation and parametric processes taking place in the same nonlinear crystal (BBO). The tunable spectra of PS has been generated between 480 nm and 530 nm, which is pumped by the second-harmonic from the high-power Ti: sapphire laser system at 1 kHz repetition rate. We present the generation mechanism of PS theoretically and simulate the process of PS ring using the amplification transfer function. The experiment and the theory show that PS will appear when the phase matching angle for second-harmonic generation is close to the optimal pump angle for optical parametric generation, and then the tunable spectra of PS are generated by slightly adjusting the crystal angle. The result provides a theoretical basis for controlling the generation of PS and quantum entanglement states, which is of great significance for the development of quantum imaging, quantum communications and other applieations.

Effects of shock waves on spatial distribution of proton beams in ultrashort laser-foil interactions

The characteristics of proton beam generated in the interaction of an ultrashort laser pulse with a large prepulse with solid foils are experimentally investigated. It is found that the proton beam emitted from the rear surface is not well collimated, and a "ring-like" structure with some "burst-like" angular modulation is presented in the spatial distribution. The divergence of the proton beam reduces significantly when the laser intensity is decreased. The "burst-like" modulation gradually fades out for the thicker target. It is believed that the large divergence angle and the modulated ring structure are caused by the shock wave induced by the large laser prepulse. A one-dimensional hydrodynamic code, MED103, is used to simulate the behavior of the shock wave produced by the prepulse. The simulation indicates that the rear surface of the foil target is significantly modified by the shock wave, consequently resulting in the experimental observations. (c) 2006 American Institute of Physics.

Optimization of multiple filamentation of femtosecond laser pulses in air using a pinhole

The robustness and prolongation of multiple filamentation (MF) for femtosecond laser propagation in air are investigated experimentally and numerically. It is shown that the number, pattern, propagation distance, and spatial stability of MF can be controlled by a variable-aperture on-axis pinhole. The random MF pattern can be optimized to a deterministic pattern. In our numerical simulations, we configured double filaments to principlly simulate the experimental MF interactions. It is experimentally and numerically demonstrated that the pinhole can reduce the modulational instability of MF and is favorable for a more stable MF evolution. (c) 2007 Optical Society of America.

紫花苜蓿不同根系分布模式的土壤水分模拟和验证

根系分布影响着土壤水分养分吸收,实测根系分布费时费力,经验根系分布函数参数简单,应用方便。该研究在田间采用苜蓿栽培土柱试验,测定根系分布,并将其和不同经验根系分布函数分别应用于Hydrus-1D对土壤水分进行动态模拟,通过土壤水分实测值和模拟值比较,验证分析了经验根系分布函数的适用性以及对土壤水分动态变化的影响。结果表明:拟合的根系分布、Prasad分布、Hoffman和van Genuchten分布3种根系分布函数的根长密度模拟值与36cm以下的根长密度实测值较为吻合,Raats根系分布模拟值与实测值及其他分布函数则差别较大。不同根系分布下土壤水分模拟差别不大,平均相对均方根误差在3.5%以下。非胁迫生长条件下,Prasad根系分布、Hoffman和vanGenuchten根系分布都可描述紫花苜蓿实际根系分布状况。

黄土高原不同降水类型区旱作玉米田土壤干燥化效应与土壤水分承载力模拟研究

在模拟精度验证基础上,应用WinEPIC模型模拟研究了黄土高原不同降水类型区1960～2000年期间不同肥力水平下连作春玉米的产量变化和土壤水分效应。模拟结果表明:(1)洛川、延安、寿阳、榆林连作春玉米产量地区间差异显著,其年平均值分别为4.40、3.75、2.50、1.75t/hm2。施肥增产率和水分利用效率地区间差异均显著,均表现为洛川>延安>寿阳>榆林;同一地区水分利用效率随施肥量的增加而提高,不同地区施肥增产率最大值范围不同。(2)0～7m土层逐月土壤有效含水量地区间差异显著,洛川>延安>寿阳>榆林;同一地区内肥力水平越高,春玉米耗水量越大,逐月土壤有效含水量平均值越低。在春玉米生育期内0～7m土层年度土壤有效水分增减量地区间差异显著,洛川<延安<寿阳<榆林;不同肥力处理间的差异在模拟前期为显著,在模拟后期为不显著。(3)所有地区的春玉米在连作期间0～7m土层土壤湿度剖面分布,在经历了土壤湿度逐年降低、土壤干层逐年加厚的干燥化过程后,均出现了稳定的土壤干层。土壤干燥化速度地区间差异显著,洛川<延安<寿阳<榆林;不同肥力处理间随着施肥水平的提高土壤干燥化速度加快。(4)洛川、延安、寿阳、榆林连作春玉米的...

黄土旱塬Hybrid-Maize模型适应性及春玉米生产潜力估算

通过对Hybrid-Maize玉米高产模型进行田间验证,应用该模型对黄土旱塬春玉米生产潜力进行初步估算。结果表明,Hybrid-Maize模型在黄土旱塬表现出较好模拟效果,总生物量、秸秆生物量和籽粒产量模拟值与实测值间具有极显著线性相关性,其决定系数分别为0.9469、0.8164和0.9650,回归系数分别为1.0198、0.9787和1.1844,接近于1。黄土旱塬区多年光温生产潜力和气候生产潜力因品种不同有所差别,对多年平均光温籽粒和总生物量生产潜力,紧凑型玉米品种分别为13.25和22.45t/hm2,平展型玉米品种分别为12.32和20.62t/hm2,年际变化小;对多年平均气候籽粒和总生物量生产潜力,紧凑型玉米品种分别为11.97和19.94t/hm2,平展型玉米品种分别为11.37和18.63t/hm2,年际波动大。在黄土旱塬区,玉米产量潜力挖掘的主要途径应集中在提高密度和水分限制条件下,Hybrid-Maize玉米模型在指导玉米高产栽培上具有较好应用。

施肥水平对冬小麦田土壤水分影响的模拟

在模拟精度验证基础上,应用WinEPIC模型模拟研究了黄土高原长武旱塬地1957～2001年期间不同肥力水平下连作冬小麦田的土壤干燥化效应。结果表明:4种肥力处理下连作冬小麦逐月土壤有效含水量均呈现波动性降低趋势,平均每年分别减少8.6、8.6、11.1和13.7mm,无肥和低肥与中肥和高肥处理间土壤有效含水量差异十分显著,但无肥和低肥处理间无明显差异;在模拟初期出现了土壤湿度逐年降低、土壤干层逐年加厚的干燥化过程,在连作第5～8年以后均出现了稳定的土壤干层,无肥和低肥处理土壤干层均分布于2～3m土层,中肥和高肥肥处理分别分布于2～4m和2～5m土层,随肥力和作物产量水平的提高,土壤干层厚度增厚;4种肥力处理麦田生产年度耗水量接近且呈现波动性缓慢降低趋势,但前1～4年中肥和高肥处理麦田共计较无肥处理麦田多耗水91和203mm,长武旱塬地麦田土壤水分承载力为1.422～2.405t/hm2,相应的肥力水平为N45～90kg/hm2和P15～30kg/hm2。

分布式水文模型的误差分析

分布式水文模型以其具有明确物理意义的参数结构和对空间分异性的全面反映,能够准确详尽地描述和模拟流域内真实的降水径流过程而被广泛需求和关注。在模拟土地利用、土地覆盖、水土流失等各种变化过程的水文响应,面源污染、陆面过程、气候变化影响评价等诸多领域都有广泛的应用。模型的预报精度和误差至关重要,决定了模型的应用和推广。在分析分布式水文模型建立和验证过程的基础上,提出了模型的4类误差来源:被排除在外的因素引起的误差,实测历史记录资料的随机或系统误差,参数误差和模型结构误差,讨论了各类误差的分析与计算方法,为模型的发展和成长提供了依据。

施肥水平对长武旱塬地冬小麦产量影响的模拟

为了在实时气象条件下确定旱塬区适宜的肥力水平和产量水平,在模拟精度验证基础上,应用WinEPIC模型模拟研究了黄土高原长武旱塬地1957～2001年期间不同肥力水平下连作冬小麦田产量效应。模型验证结果表明CK、N、NP处理产量模拟值与观测值之间的相关系数分别为0.740、0.764和0.740,均达到显著水平,WinEPIC模型对不同肥力处理下冬小麦籽粒产量模拟较为准确。模拟结果表明,无肥、低肥、中肥和高肥处理下连作冬小麦的产量均呈现波动性降低趋势,其平均值分别为1.055、1.422、2.405和3.170t/hm2,不同肥力处理间差异显著,以中肥和高肥处理增产效果最好。为黄土高原南部沟壑区旱地小麦生产的可持续发展提供了科学依据。

MATRIX DESCRIPTION OF DYNAMIC HOLZ DIFFRACTION TESTED ON THE STRAINED LAYER SUPERLATTICE SI/GESI

A matrix formulation has been developed and applied to simulate large-angle convergent-beam electron diffraction (LACBED) patterns from the Si/GexSi1-x strained layer superlattice (SLS). Good quantitative agreement has been achieved between experimental and simulated patterns. By utilizing dynamical HOLZ line patterns, we demonstrate that an accuracy of better than 0.1% can be achieved in the determination of the averaged lattice constant of a SLS, and the averaged number of layers of atoms within one period of SLS can be determined up to a single monolayer.

2-DELTA-FUNCTION GENERALIZED PLASMA POLE MODEL FOR 1ST PRINCIPLE CALCULATIONS OF QUASI-PARTICLE ENERGIES IN MANY-ELECTRON SYSTEMS

To evaluate the dynamical effects of the screened interaction in the calculations of quasiparticle energies in many-electron systems a two-delta-function generalized plasma pole model (GPP) is introduced to simulate the dynamical dielectric function. The usual single delta-function GPP model has the drawback of over simplifications and for the crystals without the center of symmetry is inappropriate to describe the finite frequency behavior for dielectric function matrices. The discrete frequency summation method requires too much computation to achieve converged results since ab initio calculations of dielectric function matrices are to be carried out for many different frequencies. The two-delta GPP model is an optimization of the two approaches. We analyze the two-delta GPP model and propose a method to determine from the first principle calculations the amplitudes and effective frequencies of these delta-functions. Analytical solutions are found for the second order equations for the parameter matrices entering the model. This enables realistic applications of the method to the first principle quasiparticle calculations and makes the calculations truly adjustable parameter free.

Effect of Si-Ge interdiffusion on the waveguide properties of SiGe-Si MQW photodetector

Because of Si-Ge interdiffusion in the Si-SiGe interface during the growth process, the square-wave refractive index distribution of a SiGe-Si multiple-quantum-web (MQW) will become smooth. In order to simulate the actual refractive index profile, a staircase approximation is applied. Based on this approach, the dispersion equation of the MQW waveguide is obtained by using a transfer matrix method, The effects of index changes caused by the interdiffusion on the optical field and the characteristics of the photodetector are evaluated by solving the dispersion equation, It is shown that the Si-Ge interdiffusion can result in a reduction of the effective absorption coefficient and the quantum efficiency.

Gain Switch of an AlGaInP Red Light Semiconductor Laser Diode

We solve the single mode coupled rate equations by computer, simulate the behavior of a gain switch of an AlGalnP red light semiconductor laser diode, and find the characteristic of FWHM of pulses changing with the amplitude of modulation signal, the bias current, and the modulated frequency. On this basis, we conduct experiments. The experiment results accord with the simulations well.

Applicatioo of Pade Approximation in Simulating Photonic Crystals

To save finite-difference time-domain(FDTD) computing time, several methods are proposed to convert the time domain FDTD output into frequency domain. The Padé approximation with Baker's algorithm and the program are introduced to simulate photonic crystal structures. For a simple pole system with frequency 160THz and quality factor of 5000,the intensity spectrum obtained by the Padé approximation from a 28-item sequence output is more exact than that obtained by fast Fourier transformation from a 220-item sequence output. The mode frequencies and quality factors are calculated at different wave vectors for the photonic crystal slab from a much shorter FDTD output than that required by the FFT method,and then the band diagrams are obatined. In addition,mode frequencies and Q-factors are calculated for photonic crystal microcavity.