The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field ...The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are deriv...This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.展开更多
On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharp...On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocenters were precisely located. Using 8 aftershocks as empirical Green(s function (EGF) we calculated the source time function (STF) of the ML=4.1 event respectively, and by superposing these results, we obtained the average time function which has better signal to noise ratio. The final result shows that the ML=4.1 earthquake consisted of two events, one of which was stronger, the other was weaker, STFs from each station demonstrate an obvious (seismic Doppler effect(. Through a trial-and-error method, the following results were drawn: for the 1st event, the rupture length is 0.44 km, rupture velocity is 4.0 km/s, while for the 2nd event those parameters are 0.25 km, 3.0 km/s, the two events ruptured in the direction of 140° and 90° respectively (measured from the strike of the fault, anticlockwise positive), the two events were separated by 0.57 km in space, and 0.09 s in time. Using single EGF, we also retrieved the source time function of other 13 events with ML(2.1 in the Huailai earthquake sequence. These results were reached: each of the 13 events consisted of single event, with the duration of the source time functions ranging between 0.05～0.16 s. Using far-field ground displacement measurement method, the scalar seismic moment, stress drop and rupture radii of 25 events with 0.95(ML(4.1 were obtained. The trend that scalar seismic moment and stress drop increase with magnitude was obvious. In the magnitude range 0.9(ML(2.4, the obtained radii show no correlativity with magnitude, we tend to believe, in such small magnitude range, precise measurement for radii of such small events is beyond the data used in this study.展开更多
This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional a...This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.展开更多
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...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 aft展开更多
Analytical and numerical investigation is made of the source potential for-floating structure with forward speed in waves. A particular form is selected for numerical applications, where the double integral of the Gre...Analytical and numerical investigation is made of the source potential for-floating structure with forward speed in waves. A particular form is selected for numerical applications, where the double integral of the Green function is transformed into the single one and the oscillation characteristics for integrands in the specific computation domain are treated numerically. A comparison of calculated examples with published data is given and it shows that the numerical simulation is satisfactory and the accuracy is adequate to engineering application.展开更多
A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surfa...A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare.展开更多
In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the rou...In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the routine methods for arrival picking are insufficiently efficient due to these strong influences. Here, we modified the empirical green function method and applied it to extract the green functions and the radiation coefficients of microseismic events with the largest energy in clusters.Multiple-channel records were used to estimate an "average" source spectrum and then the "average" source wavelet was removed from the records by deconvolution. We applied this method to the real data,and the result indicated clear improvement in extracting the dominant event of the clusters.展开更多
基金supported by the National Natural Science Foundation of China (Grant No. 50879090)
文摘The singularities, oscillatory performances and the contributing factors to the 3-＇D translating-pulsating source Green function of deep-water Havelock form which consists of a local disturbance part and a far-field wave-like part, are analyzed systematically. Relative numerical integral methods about the two parts are presented in this paper. An improved method based on LOBATTO rule is used to eliminate singularities caused respectively by infinite discontinuity and jump discontinuous node from the local disturbance part function, which makes the improvement of calculation efficiency and accuracy possible. And variable substitution is applied to remove the singularity existing at the end of the integral interval of the far-field wave-like part function. Two auxiliary techniques such as valid interval calculation and local refinement of integral steps technique in narrow zones near false singularities are applied so as to avoid unnecessary integration of invalid interval and improve integral accordance. Numerical test results have proved the efficiency and accuracy in these integral methods that thus can be applied to calculate hydrodynamic performance of floating structures moving in waves.
基金Program for New Century Excellent Talents in University Under Grant No. NCET-05-0248the Key Program for Applied Basic Research of Tianjin Municipality Under Grant No. 07JCZDJC10100
文摘This paper presents an indirect boundary integration equation method for diffraction of plane SV waves by a 2-D cavity in a poroelastic half-space.The Green＇s functions of compressive and shear wave sources are derived based on Biot＇s theory. The scattered waves are constructed using fictitious wave sources close to the boundary of the cavity, and their magnitudes are determined by the boundary conditions. Verification of the accuracy is performed by： （1） checking the satisfaction extent of the boundary conditions, （2） comparing the degenerated solutions of a single-phased case with well- known solutions, and （3） examining the numerical stability of the solutions. The nature of diffraction of plane SV waves around a cavity in a poroelastic half-space is investigated by numerical examples.
基金National Key Basic Research Project "Mechanism and Prediction of Continental Strong Earthquakes".
文摘On July 20, 1995 an M_L= 4.1 earthquake occurred in Huailai basin, northwest of Beijing, with epicenter coordinates 40.326°N, 115.448°E and focal depth 5.5 km. Following the main shock, microseismicity sharply increased in the basin. This earthquake sequence was recorded by Sino-European Cooperative Huailai Digital Seismograph Network (HDSN) and the hypocenters were precisely located. Using 8 aftershocks as empirical Green(s function (EGF) we calculated the source time function (STF) of the ML=4.1 event respectively, and by superposing these results, we obtained the average time function which has better signal to noise ratio. The final result shows that the ML=4.1 earthquake consisted of two events, one of which was stronger, the other was weaker, STFs from each station demonstrate an obvious (seismic Doppler effect(. Through a trial-and-error method, the following results were drawn: for the 1st event, the rupture length is 0.44 km, rupture velocity is 4.0 km/s, while for the 2nd event those parameters are 0.25 km, 3.0 km/s, the two events ruptured in the direction of 140° and 90° respectively (measured from the strike of the fault, anticlockwise positive), the two events were separated by 0.57 km in space, and 0.09 s in time. Using single EGF, we also retrieved the source time function of other 13 events with ML(2.1 in the Huailai earthquake sequence. These results were reached: each of the 13 events consisted of single event, with the duration of the source time functions ranging between 0.05～0.16 s. Using far-field ground displacement measurement method, the scalar seismic moment, stress drop and rupture radii of 25 events with 0.95(ML(4.1 were obtained. The trend that scalar seismic moment and stress drop increase with magnitude was obvious. In the magnitude range 0.9(ML(2.4, the obtained radii show no correlativity with magnitude, we tend to believe, in such small magnitude range, precise measurement for radii of such small events is beyond the data used in this study.
基金support from the Program for New Century Excellent Talents in University (NCET-05-0248)the Key Program for Applied Basic Research of Tianjin Municipality (07JCZDJC10100)
文摘This paper presents an indirect boundary integration equation method for diffraction of plane P waves by a two-dimensional canyon of arbitrary shape in poroelastic half-space. The Green＇s functions of compressional and shear wave sources in poroelastic half-space are derived based on Biot＇s theory. The scattered waves are constructed using the fictitious wave sources close to the boundary of the canyon, and magnitude of the fictitious wave sources are determined by the boundary conditions. The precision of the method is verified by the satisfaction extent of boundary conditions, the comparison between the degenerated solutions of single-phased half-space and the well-known solutions, and the numerical stability of the method.
文摘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 aft
文摘Analytical and numerical investigation is made of the source potential for-floating structure with forward speed in waves. A particular form is selected for numerical applications, where the double integral of the Green function is transformed into the single one and the oscillation characteristics for integrands in the specific computation domain are treated numerically. A comparison of calculated examples with published data is given and it shows that the numerical simulation is satisfactory and the accuracy is adequate to engineering application.
基金financially supported by the National Basic Research Program of China(973 Program,Grant No.2014CB046203)
文摘A domain decomposition and matching method in the time-domain is outlined for simulating the motions of ships advancing in waves. The flow field is decomposed into inner and outer domains by an imaginary control surface, and the Rankine source method is applied to the inner domain while the transient Green function method is used in the outer domain. Two initial boundary value problems are matched on the control surface. The corresponding numerical codes are developed, and the added masses, wave exciting forces and ship motions advancing in head sea for Series 60 ship and S175 containership, are presented and verified. A good agreement has been obtained when the numerical results are compared with the experimental data and other references. It shows that the present method is more efficient because of the panel discretization only in the inner domain during the numerical calculation, and good numerical stability is proved to avoid divergence problem regarding ships with flare.
基金supported by the National Natural Science Foundation of China(Grant:41404038)the Natural Science Foundation of Jiangsu Province(Grant:BK20130570)the National Key R&D Program of China(2016YFC0600302)
文摘In mining, reservoir impoundment and hydraulic fracturing exploitation, the induced microseismic clusters may present complex waveforms at receivers since the individual events generally arrive very close. And the routine methods for arrival picking are insufficiently efficient due to these strong influences. Here, we modified the empirical green function method and applied it to extract the green functions and the radiation coefficients of microseismic events with the largest energy in clusters.Multiple-channel records were used to estimate an "average" source spectrum and then the "average" source wavelet was removed from the records by deconvolution. We applied this method to the real data,and the result indicated clear improvement in extracting the dominant event of the clusters.
