Dissipation source function and an improvement to LAGFD-WAM model

In this paper, a new theoretical expression of dissipation term is presented on the basis of statistical model of breaking wave, which is an improvement to LAGFD-WAM wave model. The computational results in three typical wind fields show a good improvement to LAGFD-WAM model and a better accuracy in comparison with the observed data in the South China Sea.

The optimal semi-analytical modeling for the infinite-conductivity horizontal well performance under rectangular bounded reservoir based on a new instantaneous source function

The main objective of this study is to develop the optimal semi-analytical modeling for the infiniteconductivity horizontal well performance under rectangular bounded reservoir based on a new instantaneous source function.The available semi-analytical infinite-conductivity models(ICMs)for horizontal well under rectangular bounded reservoir in literature were developed by applying superposition of pressures in space(SPS).A new instantaneous source function(i.e.,instantaneous uniform-flux segmentary source function under bounded reservoir)is derived to be used instead of SPS to develop the optimal semi-analytical ICM.The new semi-analytical ICM is verified with ICM of Schlumberger[1]and with previous semi-analytical ICMs in terms of bottom hole pressure(BHP)profile and inflow rate distribution along the wellbore.The model is also validated with real horizontal wells in terms of inflow rate distribution along the wellbore.The results show that the developed model gives the optimal semi-analytical modeling for the infinite-conductivity horizontal well performance under rectangular bounded reservoir.Besides that,high computationalefficiency and high-resolution of wellbore discretization have been achieved(i.e.,wellbore segment number could be tens of hundreds depending on solution requirement).The results also show that at pseudosteady state(PSS)flow regime,inflow rate distribution along the wellbore by previous semi-analytical ICMs is stabilized U-shaped as performance of inflow rate distribution at late radial flow regime.Therefore,the previous semi-analytical ICMs are incorrectly modeling inflow rate distribution at PSS flow regime due to the negative influence of applying SPS.The optimal semi-analytical ICM is in a general form and real time domain,and can be applicable for 3D horizontal well and 2D vertical fracture well under infinite and rectangular bounded reservoirs,of uniform-flux and infinite-conductivity wellbore conditions at any time of well life.

Source function extracted from single-event HBT correlation functions of hydrodynamical sources

We investigate the single-event two-pion correlation functions for the hydrodynamic particle-emitting sources with the fluctuating initial conditions generated by the Heavy Ion Jet Interaction Generator （HIJING）. Using a three-dimension fast Fourier transform （FFT）, we further extract the source functions from the single-event correlation functions. It is found that the inhomogeneity of the hydrodynamic sources with the fluctuating initial conditions lead to event-by-event fluctuations of the correlation functions and source functions.

Source time functions of the Gonghe，China earthquake retrieved from long－period digital waveform data using empirical Green＇s function technique

An earthquake of Ms= 6, 9 occurred at the Gonghe, Qinghai Province, China on April 26, 1990. Three larger aftershocks took place at the same region, Ms= 5. 0 on May 7, 1990, Ms= 6. 0 on Jan. 3, 1994 and Ms= 5. 7on Feb. 16, 1994. The long-period recordings of the main shock from China Digital Seismograph Network (CDSN) are deconvolved for the source time functions by the correspondent0 recordings of the three aftershocks asempirical Green's functions (EGFs). No matter which aftershock is taken as EGF, the relative source time functions (RSTFs) Obtained are nearly identical. The RSTFs suggest the Ms= 6. 9 event consists of at least two subevents with approximately equal size whose occurrence times are about 30 s apart, the first one has a duration of 12 s and a rise time of about 5 s, and the second one has a duration of 17 s and a rise time of about & s. COmParing the RSTFs obtained from P- and SH-phases respectively, we notice that those from SH-phases are a slightly more complex than those from p-phases, implying other finer subevents exist during the process of the main shock. It is interesting that the results from the EGF deconvolution of long-Period way form data are in good agreement with the results from the moment tensor inversion and from the EGF deconvolution of broadband waveform data. Additionally, the two larger aftershocks are deconvolved for their RSTFs. The deconvolution results show that the processes of the Ms= 6. 0 event on Jan. 3, 1994 and the Ms= 5. 7 event on Feb. 16,1994 are quite simple, both RSTFs are single impulses.The RSTFs of the Ms= 6. 9 main shock obtained from different stations are noticed to be azimuthally dependent, whose shapes are a slightly different with different stations. However, the RSTFs of the two smaller aftershocks are not azimuthally dependent. The integrations of RSTFs over the processes are quite close to each other, i. e., the scalar seismic moments estimated from different stations are in good agreement. Finally the scalar seismic moments of the three aftershocks are compared. The relative scalar seismic moment Of the three aftershocks deduced from the relative scalar seismic moments of the Ms=6. 9 main shock are very close to those inverted directly from the EGF deconvolution. The relative scalar seismic moment of the Ms =6. 9 main shock calculated using the three aftershocks as EGF are 22 (the Ms= 6. 0 aftershock being EGF), 26 (the Ms= 5. 7 aftershock being EGF) and 66 (the Ms= 5. 5 aftershock being EGF), respectively. Deducingfrom those results, the relative scalar sesimic moments of the Ms= 6. 0 to the Ms= 5. 7 events, the Ms= 6. 0 tothe Ms= 5. 5 events and the Ms= 5. 7 to the Ms= 5. 5 events are 1. 18, 3. 00 and 2. 54, respectively. The correspondent relative scalar seismic moments calculated directly from the waveform recordings are 1. 15, 3. 43, and 3. 05.

Source time functions of the 1999, Jiji (Chi-Chi) earthquake from GDSN long period waveform data using aftershocks as empirical Green's functions

A large earthquake (Mw=7.6) occurred in Jiji (Chi-Chi), Taiwan, China on September 20, 1999, and was followed by many moderate-size shocks in the following days. Two of the largest aftershocks with the magnitudes of Mw=6.1 and Mw=6.2, respectively, were used as empirical Green's functions (EGFs) to obtain the source time functions (STFs) of the main shock from long-period waveform data of the Global Digital Seismograph Network (GDSN) including IRIS, GEOSCOPE and CDSN. For the Mw=6.1 aftershock of September 22, there were 97 pairs of phases clear enough from 78 recordings of 26 stations; for the Mw=6.2 aftershock of September 25, there were 81 pairs of phases clear enough from 72 recordings of 24 stations. For each station, 2 types of STFs were retrieved, which are called P-STF and S-STF due to being from P and S phases, respectively. Totally, 178 STF individuals were obtained for source-process analysis of the main shock. It was noticed that, in general, STFs from most of the stations had similarities except that those in special azimuths looked different or odd due to the mechanism difference between the main shock and the aftershocks; and in detail, the shapes of the STFs varied with azimuth. Both of them reflected the stability and reliability of the retrieved STFs. The comprehensive analysis of those STFs suggested that this event consisted of two sub-events, the total duration time was about 26 s, and on the average, the second event was about 7 s later than the first one, and the moment-rate amplitude of the first event was about 15% larger than that of the second one.

Source parameters determination for earthquakes in Kushiro,Japan considering source time function

This paper applies a new formulation to do moment tensor inversion for earthquakes in the Kushiro area of Japan. Comparing with conventional moment tensor inversion method, the new one takes the effect of source time function into consideration. For the inversion, best solution is obtained by minimizing the difference between the observed seismograms and the synthetic ones. And the best-fitting focal depth is determined from the variance reduction. The results indicate that half duration of source time function is proportional to the magnitude of earthquakes. Large earthquakes have long half duration, whereas that of moderate-small earthquakes is comparatively shorter. The focal mechanisms of all three earthquakes are of thrust fault type, which is mainly ascribed to the collision of the North American plate with the Eurasia plate in the late Cretaceous or Paleogene.

Tomographic inversion for microseismic source parameters in mining

We propose a new method for inverting source function of microseismic event induced in mining. The observed data from microseismic monitoring during mining are represented by a wave equation in a spherical coordinate system and then the data are transformed from the time-space domain to the time-slowness domain based on tomographic principle, from whichwe can obtain the signals related to the source in the time-slowness domain. Through analyzing the relationship between the signal located at the maximum energy and the source function, we derive the tomographic equations to compute the source function from the signals and to calculate the effective radiated energy based on the source function. Moreover, we fit the real amplitude spectrum of the source function computed from the observed data into the co-2 model based on the least squares principle and determine the zero-frequency level spectrum and the corner frequency, finally, the source rupture radius of the event is calculated and The synthetic and field examples demonstrate that the proposed tomographic inversion methods are reliable and efficient

Source Characteristics of the 2012 Ahar-Varzaghan Earthquake

The August 11, 2012 Mw 6.4 earthquakes in northwestern Iran occurred as a result of oblique strike-slip faulting in the shallow crust of the Eurasia plate, approximately 300 kilometer east of the plate boundary between the Eurasia and Arabia plates. The seismotectonics of this region is controlled by the collision of the Arabia and Eurasia plates;at the latitude of the earthquakes, the Arabia plate moves almost due north with respect to the Eurasia plate at a rate of approximately 26 millimeter per year. Over the past forty years, seven earthquakes of Mw 6 or greater have occurred within 300 kilometer of today’s events. The nearest was a Mw 6.1 earthquake in February of 1997, approximately 100 kilometer to the east, which caused 1100 fatalities. The studied area encourage the authors to determine the focal mechanism, source time function and sub events’ features which can guide us to reliable judges.

Energy budget of surface waves in the global ocean

Mechanical energy input from atmosphere and losses from wave-breaking dissipation of sea surface waves are estimated by a direct scheme. This scheme is based on the integration in the wavenumber space of the wind input and breaking dissipation source functions of the MASNUM wave model. The global amount of wind energy input, averaged in 2005, is about 57 TW, and the wave-breaking dissipation summed in deep-water is about 33 TW, over a half of the wind energy input. The residual may be dissipated by beach processes. Global distributions of the energy input and breaking dissipation concentrate in the westerlies of the Southern Hemisphere.

THE MESHLESS VIRTUAL BOUNDARY METHOD AND ITS APPLICATIONS TO 2D ELASTICITY PROBLEMS

A novel numerical method for eliminating the singular integral and boundary effect is processed. In the proposed method, the virtual boundaries corresponding to the numbers of the true boundary arguments are chosen to be as simple as possible. An indirect radial basis function network （IRBFN） constructed by functions resulting from the indeterminate integral is used to construct the approaching virtual source functions distributed along the virtual boundaries. By using the linear superposition method, the governing equations presented in the boundaries integral equations （BIE） can be established while the fundamental solutions to the problems are introduced. The singular value decomposition （SVD） method is used to solve the governing equations since an optimal solution in the least squares sense to the system equations is available. In addition, no elements are required, and the boundary conditions can be imposed easily because of the Kronecker delta function properties of the approaching functions. Three classical 2D elasticity problems have been examined to verify the performance of the method proposed. The results show that this method has faster convergence and higher accuracy than the conventional boundary type numerical methods.

Tempo-spatial rupture process of the 1997Mani, Xizang(Tibet), China earthquake of Ms=7.9

An earthquake of Ms=7.4 occurred in Mani, Xizang (Tibet), China on November 8, 1997. The moment tensor ofthis earthquake was inverted using the long period body wave form data from China Digital Seismograph Network(CDSN). The apparent source time functions (AS TFs) were retrieved from P and S waves, respectively, using thedeconvolution technique in frequency domain, and the tempo-spatial rupture process on the fault plane was imagedby inverting the azimuth dependent AS TFs from different stations. The result of the moment tensor inversionindicates that the P and T axes of earthquake-generating stress field were nearly horizontal, with the P axis in theNNE direction (29), the T axis in the SEE direction (122) and that the NEE-SWW striking nodal plane andNNW-SSE striking nodal plane are mainly left-lateral and right-lateral strike-slip, respectively; that this earthquakehad a scalar seismic moment of 3.4xl02o N. .m, and a moment magnitude of Mw=7.6. Taking the aftershock distribution into account, we proposed that the earthquake rupture occurred in the fault plane with the strike of 250,the dip of 88 and the rake of 19. On the basis of the result of the moment tensor inversion, the theoretical seismograms were synthesized, and then the AS T Fs were retrieved by deconvoving the synthetic seismograms fromthe observed seismograms. The A S T Fs retrieved from the P and S waves of different stations identically suggestedthat this earthquake was of a simple time history, whose ASTF can be approximated with a sine function with thehalf period of about 10 s. Inverting the azimuth dependent A S T Fs from P and S waveforms led to the imageshowing the tempo-spatial distribution of the rupture on the fault plane. From the 'remembering' snap-shots, therupture initiated at the western end of the fault, and then propagated eastward and downward, indicating an overallunilateral rupture. However, the slip distribution is non-uniform, being made up of three sub-areas, one in thewestern end, about 10 km deep ('western area'), another about 55 kin away from the western end and about 35 Iondeep ('eastern area'), the third about 30 km away from the western end and around 40 km deep ('central area').The total rupture area was around 70 km long and 60 km wide. From the 'forgetting' snap-shots, the rupturingappeared quite complex, with the slip occurring in different position at different time, and the earthquake being ofthe characteristics of 'healing pulse'. Another point we have to stress is that the locations in which the ruptureinitiated and terminated were not where the main rupture took place. Eventually, the static slip distribution wascalculated, and the largest slip values of the three sub-areas were 956 cm, 743 cm and 1 060 cm, for the western.eastern and central areas, respectively. From the slip distribution, the rupture mainly distributed in the fault about70 km eastern to the epicenter; from the aftershock distribution. however, the aftershocks were very sparse in thewest to the epicenter while densely clustered in the east to the epicenter It indicated that the Maul Ms=7.9 earthquake was resulted from the nearly eastward extension of the NEE-SWW to nearly E-W striking fault in thenorthwestern Tibetan plateau.

Chemical composition and source apportionment of PM10 and PM2.5 in different functional areas of Lanzhou,China

To elucidate the air pollution characteristics of northern China, airborne PM_10（atmospheric dynamic equivalent diameter ≤ 10 μm） and PM_（2.5）（atmospheric dynamic equivalent diameter ≤ 2.5 μm） were sampled in three different functional areas（Yuzhong County,Xigu District and Chengguan District） of Lanzhou, and their chemical composition（elements, ions, carbonaceous species） was analyzed. The results demonstrated that the highest seasonal mean concentrations of PM_10（369.48 μg/m^3） and PM_（2.5）（295.42 μg/m^3） were detected in Xigu District in the winter, the lowest concentration of PM_（2.5）（53.15 μg/m^3） was observed in Yuzhong District in the fall and PM_10（89.60 μg/m^3） in Xigu District in the fall.The overall average OC/EC（organic carbon/elemental carbon） value was close to the representative OC/EC ratio for coal consumption, implying that the pollution of Lanzhou could be attributed to the burning of coal. The content of SNA（the sum of sulfate, nitrate,ammonium, SNA） in PM_（2.5）in Yuzhong County was generally lower than that at other sites in all seasons. The content of SNA in PM_（2.5）and PM_10 in Yuzhong County was generally lower than that at other sites in all seasons（0.24–0.38）, indicating that the conversion ratios from precursors to secondary aerosols in the low concentration area was slower than in the area with high and intense pollutants. Six primary particulate matter sources were chosen based on positive matrix factorization（PMF） analysis, and emissions from dust, secondary aerosols, and coal burning were identified to be the primary sources responsible for the particle pollution in Lanzhou.

A fast integration method for translating-pulsating source Green＇s function in Bessho form

The singularities and oscillatory performance of translating-pulsating source Green＇s function in Bessho form were analyzed. Relative numerical integration methods such as Gaussian quadrature rule, variable substitution method （VSM）, and steepest descent integration method （SDIM） were used to evaluate this type of Green＇s function. For SDIM, the complex domain was restricted only on the 0-plane. Meanwhile, the integral along the real axis was computed by use of the VSM to avoid the complication of a numerical search of the steepest descent line. Furthermore, the steepest descent line was represented by the B-spline function. Based on this representation, a new self-compatible integration method corresponding to parametric t was established. The numerical method was validated through comparison with other existing results, and was shown to be efficient and reliable in the calculation of the velocity potentials for the 3D seakeeping and hydrodynamic performance of floating struc- tures moving in waves.

Spatial and temporal rupture process of the January 26, 2001, Gujarat, India, M_S=7.8 earthquake

The source parameters, such as moment tensor, focal mechanism, source time function (STF) and temporal-spatial rupture process, were obtained for the January 26, 2001, India, MS=7.8 earthquake by inverting waveform data of 27 GDSN stations with epicentral distances less than 90? Firstly, combining the moment tensor inversion, the spatial distribution of intensity, disaster and aftershocks and the orientation of the fault where the earthquake lies, the strike, dip and rake of the seismogenic fault were determined to be 92? 58?and 62? respectively. That is, this earthquake was a mainly thrust faulting with the strike of near west-east and the dipping direction to south. The seismic moment released was 3.51020 Nm, accordingly, the moment magnitude MW was calculated to be 7.6. And then, 27 P-STFs, 22 S-STFs and the averaged STFs of them were determined respectively using the technique of spectra division in frequency domain and the synthetic seismogram as Greens functions. The analysis of the STFs suggested that the earthquake was a continuous event with the duration time of 19 s, starting rapidly and ending slowly. Finally, the temporal-spatial distribution of the slip on the fault plane was imaged from the obtained P-STFs and S-STFs using an time domain inversion technique. The maximum slip amplitude on the fault plane was about 7 m. The maximum stress drop was 30 MPa, and the average one over the whole rupture area was 7 MPa. The rupture area was about 85 km long in the strike direction and about 60 km wide in the down-dip direction, which, equally, was 51 km deep in the depth direction. The rupture propagated 50 km eastwards and 35 km westwards. The main portion of the rupture area, which has the slip amplitude greater than 0.5 m, was of the shape of an ellipse, its major axis oriented in the slip direction of the fault, which indicated that the rupture propagation direction was in accordance with the fault slip direction. This phenomenon is popular for strike-slip faulting, but rather rare for thrust faulting. The eastern portion of the rupture area above the initiation point was larger than the western portion below the initiation point, which was indicative of the asymmetrical rupture. In other words, the rupturing was kind of unilateral from west to east and from down to up. From the snapshots of the slip-rate variation with time and space, the slip rate reached the largest at the 4th second, that was 0.2 m/s, and the rupture in this period occurred only around the initiation point. At the 6th second, the rupture around the initiation point nearly stopped, and started moving outwards. The velocity of the westward rupture was smaller than that of the eastward rupture. Such rupture behavior like a circle mostly stopped near the 15th second. After the 16th second, only some patches of rupture distributed in the outer region. From the snapshots of the slip variation with time and space, the rupture started at the initiation point and propagated outwards. The main rupture on the area with the slip amplitude greater than 5 m extended unilaterally from west to east and from down to up between the 6th and the 10th seconds, and the western segment extended a bit westwards and downwards between the 11th and the 13th seconds. The whole process lasted about 19 s. The rupture velocity over the whole rupture process was estimated to be 3.3 km/s.

Analysis on rupture feature of a great complicated earthquake in Kamchatka

The M S=7.3 earthquake of June, 8 1993, off the eastern coast of Kamchatka was very complicated in the rupture history. The rupture feature of this event was discussed by the broadband waveform modelling method as well as the combining analysis on the subevent stack and the quasi time difference. The results suggest that the rupture propagation of the event was in a strong unidirection and its main rupture processes can be expressed as: rupture nucleation→NEE→near east by north→near east by south→stop, from deep to shallow.

A head-related transfer function model for fast synthesizing multiple virtual sound sources

A head-related transfer function （HRTF） model for fast and real-time synthesizing multiple virtual sound sources is proposed. A head-related impulse response （HRIR, time- domain version of HRTF） is first decomposed by a two-level wavelet packet and then represented by a model composed of subband filters and reconstruction filters. The coefficients of the subband filters are the zero interpolation of the wavelet coefficients of the HRIR. The coefficients of the reconstruction filters can be calculated from the wavelet function. The model is simplified by applying a threshold method to reduce the wavelet coefficients. The calculated results indicate that for a model with 30 wavelet coefficients, the error of reconstructed HRIR is about 1%. And the result of a psychoacoustic test shows that a model with 35 wavelet coefficients is perceptually indistinguishable from the original HRIR. When multiple virtual sound sources are synthesized simultaneously, the computational cost of the proposed model is much less than the traditional HRTF filters.

Design and validation on a multiple sound source fast-measurement system of near-field head-related transfer functions

Near-field head-related transfer functions （HRTFs） are essential to scientific re- searches of binaural hearing and practical applications of virtual auditory display. High ef- ficiency, accuracy and repeatability are required in a near-field HRTF measurement. Hence, there is no reference which intents on solving the measuring difficulties of near-field HRTF for human subjects. In present work, an efficient near-field HRTF measurement system based on computer control is designed and implemented, and a fast calibration method for the system is proposed to first solve the measurement of near-field HRTF for human subjects. The efficiency of measurement is enhanced by a comprehensive design on the acoustic, electronic and mechanical parts of the system. And the accuracy and repeatability of the measurement are greatly im- proved by carefully calibrating the positions of sound source, subject and binaural microphones. This system is suitable for near-field HRTF measurement at various source distances within 1.0 m, for both human subject and artificial head. The time costs of HRTF measurement at a single sound source distance and full directions has been reduced to less than 20 minutes. The measurement results indicate that the accuracy of the system satisfies the actual requirements. The system is applicable to scientific research and can be used to establish an individualized near-field HRTF database for human subjects.

THE COMPUTATION OF PARTIAL COHERENCE FUNCTION APPLIED FOR SOUND SOURCE IDENTIFICATION

In this paper, a new method is applied to get the computation formula of partial coherence function. The main attention is paid to the computation formula of the partial coherence function with three and four signals. The advantages of the method discussed in the paper are clear in physical meaning and easy to compute at the end of the paper,the application of the method to the identification of an air compressor noise source is presented and the results are satisfactory.

Adhesion Improvement of Zirconium Coating on Polyurethane Modified by Plasmas

In order to improve the adhesion of the middle frequency magnetic sputtered zirconium coating on a polyurethane film,an anode layer source was used to pretreat the polyurethane film with nitrogen and oxygen ions.SEMs and AFM roughness profiles of treated samples and the contrast groups were obtained.Besides,XPS survey spectrums and high resolution spectrums were also investigated.The adhesion test revealed that ion bombardment could improve the adhesion to the polyurethane coating substrate.A better etching result of oxygen ions versus nitrogen predicts a higher bonding strength of zirconium coating on polyurethane and,indeed,the highest bonding strengths are for oxygen ion bombardment upto 13.3 MPa.As demonstrated in X-ray photoelectron spectroscopy,the oxygen ion also helps to introduce more active groups,and,therefore,it achieves a high value of adhesion strength.