弹性波的三维有限元模拟

本文讨论了弹性波瞬态传播问题的三维有限元计算方法及当前存在的实际困难.针对要求计算机内存大和计算时间长的问题,采取了改进措施.由于采用了集中质量矩阵和修正的中心差分时间积分显格式相结合的方法,可以使计算机内存和计算时间大为减少;由于采用结点定位法,最适合用于目前发展的并行计算机系统,可使计算速度大大增快;还采用了有效激发震源法,有效激发区是随时间步进的增加而逐步增大,这不仅能节省计算时间,而且使波场的传播过程一目了然,本文计算了由两种介质组成的三维楔形问题,得到若干典型剖面的瞬时波场图及随时间变化的合成地震图.     

Effect of some additives on synthesis of zeolite from coal fly ash

Hydrothermal conversion of fly ash into zeolites was conducted and the effects of the addition of sodium halide and waste solutions produced after zeolitization of fly ash, as well as the adjustment of the Si/Al ratio prior to synthesis process on the formation and cation exchange capacity (CEC) of zeolite product were evaluated. Both the addition of NaCl and NaF ameliorated the crystallinity and CEC of synthesized zeolite, but NaF had a better improvement effect. Na⁺ was considered to enhance the crystallization of zeolite, while F⁻ favored the dissolution of fly ash. The type of zeolite formed depended on the Si/Al ratio of the starting material prior to the nucleation and crystallization of zeolite. The adjustment of the Si/Al ratio of fly ash by addition of Na₂SiO₄ and Al(OH)₃ changed the type and CEC of zeolite. Waste solutions contained large amount of Si and little Al due to the formation of a zeolite named NaP1 in zeolite terminology with the Joint Committee of Powder Diffraction Standard (JCPDS) code of 39-0219. The alkalinity decreased largely. As a result, the CEC value of zeolite products synthesized with waste solution as alkali source decreased. The supplementation of new alkali to adjust the alkalinity of waste solution could enhance the CEC of synthesized product. It was concluded that: (1) addition of sodium halide and adjustment of the Si/Al ratio prior to synthesis can improve the quality of zeolite; (2) waste solutions produced following the zeolitization of fly ash can be reused as an alkali source in the activation of fly ash; zero-emission of waste solution in the synthesis of zeolite from fly ash is possible.     

A novel strategy using biodegradable EDDS for the chemically enhanced phytoextraction of soils contaminated with heavy metals

For the sake of cost and potential environmental risk, it is necessary to minimize the amount of chelates used in chemically-enhanced phytoextraction. In the present study, a biodegradable chelating agent, EDDS was added in a hot solution at 90℃ to the soil in which garland chrysanthemum （Chrysanthemum coronarium L.） and beans （Phaseolus vulgaris L., white bean） were growing. The application of hot chelate solutions was much more efficient than the application of normal chelate solutions （25℃） in improving the uptake of heavy metals by plants. When 1 mmol kg1 of EDDS as a hot solution was applied to soil, the concentrations ofCu, Zn and Cd and the total phytoextraction by the shoots of the two plant species exceeded or approximated those in the shoots of plants treated with 5 mmol kg^-1 of normal EDTA solution. The concentrations of metals in the shoots of beans were significantly correlated with the relative electrolyte leakage rate of root cells, indicating that the root damage resulting from the hot solution might play an important role in the process of chelate-enhanced metal uptake. The soil leaching study demonstrated that decreasing the dosage of chelate resulted in decreased concentrations of soluble metals in soils. On the 28th day following the application of chelate, the concentrations of soluble metals in the EDDS treated soil were not significantly different from the concentrations in the control soil to which chelates had not been applied.     

Deep structure of the Nojima Fault, southwest Japan, estimated from borehole observations of fault-zone trapped waves

To estimate the deep structure of the southern part of the Nojima Fault, southwest Japan without the influence of near-surface structures, we analyzed the Love-wave-type fault-zone trapped waves (LTWs) recorded by a borehole seismometer at 1800 m depth. We examined the polarization, dispersion, and dominant frequency of the wavetrain following the direct S-wave in each seismogram to identify the LTW. We selected eight candidates for typical LTWs from 462 records. Because the duration of the LTW increases with hypocentral distance, we infer that the low velocity fault-zone of the Nojima Fault continues towards the seismogenic depth. In addition, since the duration of the LTW increases nonlinearly with hypocentral distance, we infer that the S-wave velocity of the fault-zone increases with depth. The location of events showing the LTW indicates that the fault-zone dips to the southeast at 75° and continues to a depth of approximately 10 km. We assumed a uniform low-velocity waveguide to estimate the average structure of the fault-zone. We estimated the average width, S-wave velocity, and Q_s of the fault-zone by comparing an analytical solution of the LTW with measured data. The average width, S-wave velocity, and Q_s of the fault-zone are 150 to 290 m, 2.5 to 3.2 km/s, and 40 to 90, respectively. Hence the fault-zone structure with a larger width and smaller velocity reduction than the fault-zone model estimated by previous surface observation is more suitable to represent the average fault-zone structure of the Nojima fault. The present study also indicated that the shallow layers and/or a shallow fault-zone structure drastically changes the characteristics of the LTW recorded at the surface, and therefore cause a discrepancy in the fault-zone model between the borehole observation and surface observation.     

冲击凿岩钎杆的瞬态动力学分析

以ANSYS为平台，按照真实工作状态对凿岩钎杆进行了瞬态动力学分析，并对不同活塞形状下分析结果数据进行了比较，为冲击凿岩系统设计提供了理论依据。     

Modeling the seismic response of basin-edge structure by a sequential waveform method

A sequential waveform method is developed to simulate the seismic response of basin-edge structure excited by a plane incident P-wave. The full procedure involves: (a) a previous parameterization of the investigated model using the seismic wave velocities and depths of the sedimentary stratifications; (b) an input motion determined from the records at stations installed on hard rock; (c) forward computation of the P-SV elastic wave field by means of a two-dimensional finite difference (FD) method; (d) the optimization of the model vector using simulated annealing technique and comparing the simulated seismic response of the tested structure with the observed wave field; (e) the correction of the initial model by trial-and-error by testing the differences between synthetics and observed data, and (f) the final solution obtained by iteration using the conjugate gradient algorithm. The search of an optimal basin-edge model has been parallel processed by varying the shapes and velocities of strata on the basis of the fitting of relative timing, amplitude and phase between the output and the observed data. The input motion and sensitivity have been checked and the validity of the method has been demonstrated by numeric analysis. Using the teleseismic records generated by 7 earthquakes recorded at 26 broadband seismic stations, we have studied the seismic velocity structure of the southern edge of the Jiyang depression located in the Bohai Bay basin, northern China. Two cross sections show an agreement between the velocity results and the geological sections available in the region. In addition, we obtain evidence of three hidden faults under the sections and features that suggest major extensions at the Paleogene.     

Estimating the minimum in-stream flow requirements via wetted perimeter method based on curvature and slope techniques

Under the assumptions of triangular cross section channel and uniform stable flow, an analytical solution of the minimum ecological in-stream flow requirement (MEIFR) is deduced. Based on the analytical solution, the uncertainty of the wetted perimeter method is analyzed by comparing the two techniques for the determination of the critical point on the relationship curve between wetted perimeter, P and discharge, Q. It is clearly shown that the results of MEIFR based on curvature technique (corresponding to the maximum curvature) and slope technique (slope being 1) are significantly different. On the P-Q curve, the slope of the critical point with the maximum curvature is 0.39 and the MEIFR varied prominently with the change of the slope threshold. This indicates that if a certain value of the slope threshold is not available for slope technique, curvature technique may be a better choice. By applying the analytical solution of MEIFR in the losing rivers of the Western Route South-to-North Water Transfer Project in China, the MEIFR value via curvature technique is 2.5%-23.7% of the multi-year average annual discharge, while that for slope technique is 11%-105.7%. General conclusions would rely on the more detailed research for all kinds of cross-sections.