INTRODUCTION

The ground motion attenuation relationship is influenced by the source effect，path effect and site effect，which can be used to calculate the ground motion of an engineering site.It has been widely used in seismic zoning map compilation and seismic safety evaluation.It is also one of the important research directions of engineering seismology.The methods of ground motion estimation are divided into the empirical method and theoretical method （Jiang Hui， 2005）.The empirical method has a direct statistical method and conversion method.The direct statistical method is statistical regression according to an attenuation relationship formula with the existing target strong motion records.We've considered both the source effect and path effect in this statistical relationship.The method is fit for the countries and regions which have a lot of strong ground motion records，such as Japan and the Western United States.The theoretical method is also called the mapping method， that is to say， we choose a region which has a lot of strong motion records and also has similar regional tectonic conditions as a reference area，then we assumed that any point in the reference area has a matching point in the research area.If the two same epicentral distance points have the same earthquake intensity，then they have the same ground motion parameters.Although the physical conception of this method is not clear，it is a way of dealing with a lack of strong motion records.The theoretical method uses the dislocation and rise time to describe the source，and uses the seismic wave theory to simulate shear wave propagation from the source to site in homogeneous semi-infinite space.The theoretical method includes the theoretical Greens function method，the experiential Greens function method and the numerical Greens function method.The theoretical Greens function method can obtain a theoretical vibration map based on mathematics and physics.This method is suitable for a horizontal stratified medium，and also needs accurate seismic wave propagation path information and local site threedimensional shear wave velocity structure information.In fact，it is difficult to calculate the theoretical Greens function （Li Qicheng et al.， 2010）.The experiential Greens function method is a half-theoretical method，which regards large earthquakes as a series of small earthquakes with the same focal mechanism.The large earthquake ground motion time history can be obtained through superposition of the Greens function of small earthquakes.The result contains source effect， path effect and site effect as this method is based on small earthquake records， however the section of high frequency seismic waves is unstable，in which an artificial cycle exists.

McGuire et al.（1980）， Hanks et al.（1981） used Parsevel theorem of random vibration theory and statistical property of maximum random process value based on the ω2point source model to estimate the attenuation relationship between RMS of ground motion acceleration and peak ground motion acceleration.Boore （1983） compared the results between the random vibration method and point source method to synthesize the ground motion time history，and both methods have good consistency.Based on the limited duration，limited bandwidth white noise random vibration theory， Atkinson （1984） used Hanks'（Hanks T.C.et al.， 1981） method to calculate the ground motion acceleration attenuation relationship in Eastern Canada，where there was a lack of strong motion records.Atkinson et al.（1992） used the regression and inverse method to research the ground motion attenuation relationship of Canada，and determined the geometric attenuation and inelastic attenuation.Boore （2002） combined the regional parameters，the geometric attenuation parameters and the crust structure study in North America to develop the procedure of the stochastic method.Wang Guoxin et al. （2001） used the two-step regression method to build the ground motion acceleration attenuation relationship of Northern China based on the improved source spectrum model.Tao Zhengru （2010） used the last 10 records from the Northeastern Japan F-net stations and selected the Fourier spectrum as the objective function to invert the stress drop Δσ，quality factors Q0and η， using the results to build the peak ground motion acceleration attenuation relationship by the random vibration method.Cui Anping （2013）selected 147 records of 82 small earthquakes in Sichuan Province and 863 records of 154 small earthquakes in Yunnan Province and took the acceleration Fourier spectra as the objective function to invert the seismic source parameters and the crustal media parameters， such as stress drop Δσ，quality factors Q0and η， geometric attenuation parameters R1and R2， building the peak ground motion acceleration attenuation relationship of the Yunan-Sichuan region with the results.The author（Li Wenqian，2014） selected 592 records of 33 small earthquakes in Lanzhou region to invert the peak ground motion acceleration attenuation relationship of the Lanzhou region，and applied this attenuation relationship to the Lanzhou region zoning map and found that most of the regions are in good agreement， while some areas are lower.At the same time， the author used 1995 records of 28 small earthquakes in Northern China，and adopted Kappa filtering to invert the peak ground motion acceleration attenuation relationship.

Pishan County （37.6°N， 78.2°E）， located in the Hotan Region of Xinjiang Uygur Autonomous Region experienced the MS6.5 earthquake on July 3，2015.It is a kind of thrust type seismic event（Li Jin et al.，2016） occurred in the west Kunlun fault zone.By the end of July 27，2015，the Xinjiang Strong Motion Network had recorded more than 120 records including 39 main shock records， and the Xinjiang Regional Digital Seismic Network had recorded 2，173 aftershocks including 27 earthquakes above MS4.0 and 1 earthquake greater than MS5.0.In this paper， we used Pishan earthquake sequences waveform data， adopted micro-genetic algorithm to invert 5 medium parameters of the focal area.We use these parameters to calculate the Fourier spectrum and combine them with the random phase spectrum to get the ground motion time history，and build the strong ground motion acceleration attenuation relationship.It is a kind of half-experience and half-theoretical method， which considers the effects of source， path and site.At the same time，this method based on seismological small earthquake data can solve the current absence of strong motion records to some degree.It is also suitable for regions which have abundant small earthquake records but lack strong motion records.The purpose of this paper is to examine the applicability of this method in the Xinjiang region.

1 BUILDING THE GROUND MOTION ATTENUATION RELATIONSHIP

We assuming that the far field acceleration is an elastic half-space with finite duration and finite bandwidth white noise，we use seismological methods to build a point source Fourier spectrum equation.We use the source spectrum to express the source effect，use geometric attenuation and inelastic attenuation to express the path effect，and use surface amplitude amplification factor and high cutoff term to express the site effect.We can get the velocity Fourier spectrum by fast Fourier transform with small earthquake records，then choose the velocity Fourier spectrum as the target curve with the micro-genetic algorithm method to get the five medium parameters of the source region.We use these parameters to calculate the Fourier spectrum and combine them with the random phase spectrum to get the ground motion time history.In order to avoid the error caused by randomness， we took 50 times'peak ground motion accelerations to get the average value，and built the ground motion acceleration attenuation relationship.

程凤萍介绍，80年代并不是每个人都能够坐飞机的，能够坐飞机的一般都是领导干部或者改革开放初期最先富裕起来的人，那时候乘机还需要开介绍信。那时的飞机也很小，她执飞的机型是安—24，容纳乘客48人。“当时的飞机上没有水箱、没有烤箱，也不配餐，就是登机的时候每个人发一个纪念品和一盒饮料，纪念品是一个有飞机标志的钥匙扣或者是一把折扇，饮料是纸盒装的冬瓜茶，那时算是非常高档的饮品了。”程凤萍说：“那时在飞机上根本就不需要做什么服务工作，旅客上机后都在睡觉，一直到降落。”

In this equation，R is epicentral distance，expressed as R ＝D is the shortest distance on the surface of the vertical projection from the strong motion stations to the rupture surface which is the empirical results；R1and R2are distance subsection points reflecting the change of seismic wave in the propagation process.The seismic wave is taken as a shear wave when R is less than

狗可以不懂事，人必须明事理。依法养犬，文明养犬，是每个养犬人的义务。从现实来看，大多数养狗之人也是爱狗之人，日常能够做到文明养犬。但也有一些养狗人只图自己开心方便，无视公共利益和他人权益。他们或认为自家爱犬听话、不会咬人，或认为拴狗绳对狗健康不利，或认为犬只随地便溺反正有环卫工人清扫，何必污了自己的手。如此一来，势必将自己的爱犬置于矛盾的焦点，不仅是对爱犬的不负责任，也会影响所有爱狗人士的形象。

The ground motion Fourier spectrum caused by point source can be expressed as（Boore D.M.， 1983）：

辛弃疾的早年，是充满传奇色彩的壮烈生涯，从二十余岁离开山东老家投奔南宋朝廷后，一直任职地方，奔走于江南各地，时刻未忘故国故乡。但是，“南共北，正分裂”的局势，使他有家不得归，备受亲人失散的哀伤。因此，他创作了抒写离情别绪而陷于苦闷忧伤的“代人赋”闺怨词，写得是婉曲妩媚，格调高洁。既反映南宋离乱动荡的社会生活，又是他自身感情的写照。

In this equation， S（M0， f） is source spectrum， G（R） and D（R， f） are propagation path effects，A（f） and P（f） are site condition effects， and I（f） is the spectrum type parameter.From equation （1）， the ground motion Fourier spectrum caused by the point source is related to the seismic moment M0，frequency f and distance R.

1.1 The Source

where fmaxis high cut-off frequency，taken as 5-10Hz in general.

Rθφ reflects the source radiation pattern and station orientation effect， taken as 0.6 generally； F is the free surface amplification effect， taken as 2.0； V is the horizontal component coefficient of seismic energy，taken as 1； R0is the reference distance when select the ρSand βS， taken as 1.0km. ρSand βSis the medium density and shear wave velocity near the source， taken as 2.9g／cm3and 3.4km／s respectively （Zhao Cuiping et al.， 2005）.

We choose Wang Guoxin's et al. （Wang Guoxin et al.， 2001）point source spectrum model as our source spectrum S（M0， f） which is based on the improved result of Brune's ω2source spectrum model.It can express the change of corner frequency with the rupture area.It not only approaches the ω2source spectrum in the high frequency but the characteristics of source spectrum amplitude will also not quickly increase with the magnitude increase in the low frequency.It is expressed by：

We selected 63 small earthquakes in the earthquake sequence.The magnitudes of those small events were between MW3.0 and MW3.5 from July 3，2015 to July 17，2015 and the focal depths distributed in 0-11km.The 23 digital seismic stations recorded 63 earthquakes （1，179 records in total）within the scope of an epicentral distance of about 300km.We used these small earthquake records to calculate the Pishan area strong ground motion attenuation.We also collected 39 peak ground motion accelerations which are the larger between the EW component and NS components of the main shock to comparatively analyze and calculate the result.These strong ground motion stations are soil sites mostly distributed in the northwest of the epicenter.Fig.1 shows the distribution of the epicenter and the stations.All the records with the epicentral distance from 188km to 432km are distributed in the Kashi-Wuqia intersection area as the Hotan region has no strong ground motion station.Selibuya station was the closest station which recorded 30.4gal of the peak acceleration with an epicenter 191.8km，and the Qiongkuerqiake station recorded 62.5gal of the peak acceleration with the epicenter 192km，which was the largest peak acceleration of this earthquake.The strong ground motion stations which obtained the records are both soil sites except for Wuheshalu station，which is a rock site.The information of the strong motion stations and the seismological stations is listed in Tables 1 and 2.

We selected 30 seconds of waveform data which includes the S-wave to take the Fourier transform，thus we could get the Fourier spectral envelope scattering through envelope discretization.We adopted the micro-genetic algorithm to invert the regional parameters.Combined with original genetic method， this method canceled the mutation， and at the same time took the best individual preservation strategy.We used the known magnitude，the known epicentral distance and scattering of the Fourier spectral envelope to invert the stress drop Δσ，quality factors （Q0，η） and geometric attenuation parameters （R1， R2）.The range of the inversion parameters were summarized from the Xinjiang regional seismological research（Meng Lingyuan et al.，2014； Zhou Yunhao et al.，2004； Zhao Cuiping et al.，2011； Li Zhihai et al.，2010； Ayixiangu·Maimaiti et al.，2015； Pan Zhensheng et al.，2010； Wang Ji et al.，2010），and the inversion results are listed in Table 3.

The magnitude-epicenter distance distribution and the magnitude-focal depth distribution of the records were shown in the Fig.2 and Fig.3.In the magnitude-epicenter distance distribution，the magnitudes of the seismic sequence we adopted were between MW3.0 and MW3.3 except for several earthquakes for which the magnitudes were larger than MW3.3.Most of the epicenter distances were between 20km to 120km or 150km to 300km.Epicenter distances from 130km to 150km only have several records.The focal depths were from 6km to 10km.

1.2 Path

G（R） and D（R，f） are geometrical attenuation and inelastic attenuation.G（R） is related to the regional crustal velocity structure，and the seismic wave has different components at different distances.The main ingredients of the seismic wave are the shear wave in the near field，and the surface wave in the middle field.The geometrical attenuation is expressed by the three-stage geometric attenuation model（Boore D.M.， 1983）：

显然，这种聚合轴上的“选择”与“被选”揭示的是底本与述本之间潜在的且无法穷尽的替代之可能性。无论是底本中材料的取舍还是再现方式的“选择”，所最终呈现出来的述本都只露出冰山之一角。如果将底本与述本之间的“双轴关系”结合可能世界理论，则可以发现底本存在于无边无界的“可能世界”。而述本则是“三界”(实在世界、虚构世界、可能世界)混合所形成的文本形态。因此，从“跨世界通达”这一理论视角来看，“底本”指涉的是可能显现(但不一定显现)的可能世界之素材与再现方式集合。而“述本”指涉的是显现在“三界”之中的叙述文本。

Q is the quality factor，the high-frequency of the ground motion is attenuated lower than the part of the low-frequency in the same situation.Q value becomes larger，the attenuation will be slower when distance R and frequency f are determined.The quality factor shows the regional characteristic from a great deal of studies（Boore D.M.et al.，1984； Chen P.S.et al.，1984），different areas have different Q values.Q ＝ Q0fη， Q0and η are regional parameters， Q0is the quality factor when f＝ 1Hz.We take Q0and η as source regional inversion parameters.

1.3 Site

A（f） and P（f） are the amplification factor and attenuation factor.A（f） shows different frequency amplitude variation near the surface caused by difference of crustal velocity gradient.It can be approximately determined by the quarter-wave method according to the regional crustal velocity structure transfer function （Boore D.M.et al.， 1997）.P（f） is high cut-off frequency item.Papageorgiou et al. （1983） suggested that high frequency attenuation is caused by the source effect， and Hanks T.C. （1982） considered that it is the result of the site effect， or the combined action of the source effect and the site effect.We adopted the fmaxfilter，expressed as（Hanks T.C.， 1982）：R1， or taken as a surface wave when R is greater than R2， otherwise the seismic wave is taken as a mixture of shear wave and Moho reflection wave.R1and R2reflect the characteristics of wave propagation in the crust medium.It is related to the regional rather than the size of the magnitude.We took R1and R2as source regional inversion parameters.Seismic wave amplitude attenuate present index form with distance increase caused by dissipation of the energy（Atkinson et al.， 1995）， expressed by：

Proportionality coefficient C can be expressed as（Boore D.M.， 1983）：

2 THE GROUND MOTION ACCELERATION ATTENUATION RELATIONSHIP OF THE PISHAN AREA

现在，这种做法估计很难玩下去了。从有关报道来看，这种张冠李戴的公文类不正文风，已成为检查的关注点之一。6月6日，《中国纪检监察报》便有这样的报道：江西省德安县人社局在出台《德安县人力资源和社会保障系统精准扶贫行动实施方案》时，因为照搬照抄上级文件，影响扶贫工作成效，引起了九江市委巡察组的注意，相关责任人因此受到处分。而在今年新修订的《中国共产党纪律处分条例》中，已明确工作中的形式主义、官僚主义行为属于违纪行为。这是一个强有力的信号：做虚假文章也是有风险的！

Fig.1 The distribution of stations and the epicenter

Table 1 Information of strong motion stations from which the data were recorded

Station name Station code Instrument type Site type Selibuya Station 65SBY ETNA（Inner set） Soil Qiongku'erqiake Station 65QQK ETNA（Inner set） Soil Ahonglukumu Station 65ALM ETNA（Inner set） Soil Maola Station 65MLA ETNA（Inner set） Soil 42 Regiment Station 65SRT ETNA（Inner set） Soil Yuepuhu Station 65YPH ETNA（Inner set） Soil Yekexianbazha Station 65YBZ ETNA（Inner set） Soil Jiashizongchang Station 65JZC ETNA（Inner set） Soil Akekashi Station 65AKS ETNA（Inner set） Soil Guleluke Station 65GLK ETNA（Inner set） Soil Hannanlike Station 65HLK ETNA（Inner set） Soil Heiziwei Station 65HZW ETNA（Inner set） Soil Mushi Station 65MUS ETNA（Inner set） Soil Shang′atushi Station 65STS ETNA（Inner set） Soil Shufu Station 65SUF ETNA（Inner set） Soil Tashimilike Station 65TLK ETNA（Inner set） Soil Wolituohulage Station 65WLG ETNA（Inner set） Soil Wulukeqiati Station 65WQT ETNA（Inner set） Soil Xike′er Station 65XKR ETNA（Inner set） Soil Xiaputule Station 65XTL ETNA（Inner set） Soil Aoyitake Station 65YTK ETNA（Inner set） Soil

Table 1

S t a t i o n n a m e S t a t i o n c o d e I n s t r u m e n t t y p e S i t e t y p e A l a l i S t a t i o n 6 5 A L L E T N A（I n n e r s e t） S o i l J i a s h i S t a t i o n 6 5 J A S E T N A（I n n e r s e t） S o i l K a s h i F i a n c i a l S c h o o l S t a t i o n 6 5 K C X E T N A（I n n e r s e t） S o i l S h u l e S t a t i o n 6 5 S U L E T N A（I n n e r s e t） S o i l B u l a k e s u S t a t i o n 6 5 B K S E T N A（I n n e r s e t） S o i l A t u s h i S t a t i o n 6 5 A T S E T N A（I n n e r s e t） S o i l H o n g q i F a r m S t a t i o n 6 5 H Q C E T N A（I n n e r s e t） S o i l H a l a j u n S t a t i o n 6 5 H L J E T N A（I n n e r s e t） S o i l Z h o n g y a n g c h a n g S t a t i o n 6 5 Z Y C E T N A（I n n e r s e t） S o i l W u q i a S t a t i o n 6 5 W U Q E T N A（I n n e r s e t） S o i l W u h e s h a l u S t a t i o n 6 5 W S L E T N A（I n n e r s e t） B e d r o c k J i g e n S t a t i o n 6 5 J I G E T N A（I n n e r s e t） S o i l B a y i n g k u l u t i S t a t i o n 6 5 B L T E T N A（I n n e r s e t） S o i l T u o y u n S t a t i o n 6 5 T O Y E T N A（I n n e r s e t） S o i l A k t o S t a t i o n 6 5 A K T E T N A（I n n e r s e t） S o i l B a i s h i k e r a n m u S t a t i o n 6 5 B R M E T N A（I n n e r s e t） S o i l B o s i t a n t i e l i e k e S t a t i o n 6 5 B L K E T N A（I n n e r s e t） S o i l Y a m a n y a S t a t i o n 6 5 Y M Y E T N A（I n n e r s e t） S o i l

Table 2 The information of seismic stations from which the data was recorded

Station name Station code Instrument type Site type Bachu station BCH CTS-1E Bedrock Bapanshuimo station BPM CMG-3ESPC Bedrock Hotan station HTA JCZ-1 Bedrock Hotan seismic array 0 station HTTZ0 CTS-1E Bedrock Hotan seismic array 1 station HTTZ1 CMG-40T-1 Bedrock Hotan seismic array 2 station HTTZ2 CMG-40T-1 Bedrock Hotan seismic array 3 station HTTZ3 CMG-40T-1 Bedrock Hotan seismic array 4 station HTTZ4 CMG-40T-1 Bedrock Hotan seismic array 5 station HTTZ5 CMG-40T-1 Bedrock Hotan seismic array 6 station HTTZ6 CMG-40T-1 Bedrock Hotan seismic array 7 station HTTZ7 CMG-40T-1 Bedrock Hotan seismic array 8 station HTTZ8 CMG-40T-1 Bedrock Hotan seismic array 9 station HTTZ9 CMG-40T-1 Bedrock Kashi station KSH CTS-1E Bedrock Kashizhongji station KSZ CMG-3ESPC Bedrock Tashiku'ergan station TAG CMG-3ESPC Bedrock Xike'er station XKR BBVS-60 Bedrock Yecheng station YCH CMG-3ESPC Bedrock Yingjisha station YJS CMG-3ESPC Bedrock Yuepuhu station YPH FSS-3DBH Bedrock Floating seismic station 1 DL1 CMG-40T Bedrock Floating seismic station 2 DL2 CMG-40T Bedrock

Fig.2 The magnitude-epicentral distance distribution of siesmological records

Fig.3 The magnitude-focal depth distribution of siesmological records

We adopt the viewpoint of some seismologists（Abercrombie R.E.， 1995； Kanamori H.et al.，1975； Shearer P.M.， 2009）， they considered that the stress drop of small earthquakes is the same as large earthquakes which is basically constant，and therefore we take the stress drop as the source regional inversion parameter.

控制变量中营业毛利率与企业盈利能力有显著的正相关关系。即企业营业毛利率越高，其盈利能力就越强。具体而言，企业营业毛利润每上升一个单位，企业的盈利能力会增加0.70%。而资产负债率、企业规模与企业盈利能力为负相关，即生产性服务业上市公司的资产负债率越高、企业规模越大，其企业盈利能力可能就会越低。但是，资产负债率与企业盈利能力的这种负相关关系显著，而企业规模与企业盈利能力的负相关关系并不显著。

M0is seismic moment， f is frequency， factors are a ＝3.05-0.33MWand b＝2.0／a， f0is corner frequency， the relationship between f0and Δσ expressed as：

Table 3 The result of parameters inversion

Parameter Parameter scope Inversion result Δσ 40-200 bars 75.97 bars Q0 100-400 186.33 η 0.1-1.0 0.26 R1 50km-100km 72.18km R2 100km-150km 139.70km

We can get the acceleration Fourier amplitude spectrum with equation 1 by giving a pair of magnitude and distance according to the source parameters which we inverted.The calculated Fourier spectrum is combined with random phase spectrum to get the ground motion time history.We picked the peak ground motion acceleration of each time history and built the strong ground motion acceleration attenuation relationship.We added the earthquakes within epicenter 300km，the range of magnitude of（6.5 ±0.5）， including the Wuqia MS6.9 earthquake on October 5，2008 （26 records）， the Akto MS6.2 earthquake on October 6， 2008 （3 records） and the Lop MS6.0 earthquake on March 9， 2012 （23 records）， as most epicentral distance of the Pishan earthquake records are far less than 200km.Because the Kunlun Mountains and Altun Mountains have few strong ground motion stations，the location of the recorded stations of those three earthquakes were in accord with the Pishan earthquake recorded stations.We compared the calculated result and empirical attenuation relationships （Zhang Zhenbin et al.， 2010； Yu Yanxiang et al.，2013） with strong motion records in Fig.4.

The calculation result is expressed by a series of scatter points which were applied to the mid-far field as the adoption point source spectrum model，and epicentral distance is from 10km to 300km in order to consider engineering significance.The result is similar to a straight line because we use three-stage geometry attenuation in 70km-130km.From Fig.4 we can see the epicentral distance of Wuqia earthquake is distributed in 40km-300km，and the epicentral distance of the Lop earthquake is distributed in 200km-300km.The short axis results of Yu Yanxiang et al. （2013） indicates that it goes through the Wuqia earthquake records preferable but lower than the Lop and Pishan earthquake records； the result of Zhang Zhenbin et al. （2010）is consistent with the Wuqia and Lop earthquake records within the epicentral distance of 200km-300km，and is higher than the Wuqia earthquake records and lower than the Pishan earthquake records within 40km-200km.The result was higher than the Wuqia earthquake records and lower than the Pishan earthquake records，but can go through the Lop earthquake records.Pishan earthquake records are higher than the result of calculation and empirical attenuation relationships due to the Pishan earthquake occurring in the Hotan region of the Kunlun fault zone，and most of the strong ground motion stations which recorded records are distributed in the Wuqia area，where the intersecting area of western part of southern Tianshan Mountains and Kunlun Mountains is.The Pishan earthquake records have the geological characteristics of the Tianshan Mountains region.The seismic wave energy has different dissipation due to passing through two different crustal media in the propagation process，and finally expressed by the Pishan earthquake records higher than the calculate result and empirical attenuation relationships.We can indicate that the attenuation of Hotan region is partly faster than the Wuqia region.At the same time，most of the observation peak ground accelerations are soil site conditions，while empirical attenuation relationships are established by the rock site.The Lop earthquake records have similar characteristics to the Pishan earthquake records as Lop County is located 166km east of Pishan，so the Lop earthquake records are higher than the empirical attenuation relationships as well.Wuqia earthquake records are lower than the calculated result which is mainly caused by inversion parameter selection.There are differences in the quality factor parameters，because we used the seismological small earthquakes records to calculate the ground motion attenuation relationship，but adopted the parameters'scope of the South Tianshan Mountain fault（Kashi-Wuqia intersection area） to inversion.Meanwhile， even though we added three earthquakes to conduct comparison，but have only 91 records distributed 100km away.It reflects the lack of strong motion records and also indicates that verifying the empirical attenuation relationships by using the limited strong motion data is unreliable.

元宵节是我国的传统节日之一。元宵之夜，大街小巷张灯结彩，人们点起万盏花灯，携亲伴友出门赏灯，早已成为我们的传统习俗。时代虽然在快速地发展，但这些传统文化依旧是人们心中割舍不断的情愫。

Fig.4 The comparison of empirical attenuation relationship，our study result with strong motion records

3 CONCLUSION AND DISCUSSION

We took the Pishan MS6.5 earthquake as the research object，used the seismic wave data with magnitudes between 3.0 and 3.5 recorded by the broadband digital seismic network.We could get 5 medium parameters using the micro genetic algorithm.After many Fourier transforms with the time function，we can get the peak ground motion acceleration attenuation relationship.We compared the calculated result and empirical results with the strong motion records including the Pishan MS6.5 earthquake on July 3， 2015， Wuqia MS6.9 earthquake on October 5， 2008，the Akto MS6.2 earthquake on October 6， 2008 and the Lop MS6.0 earthquake on March 9，2012.The Pishan and Lop earthquakes occurred in the Hotan region of Kunlun Fault Zone，and most of the strong ground motion stations which recorded records are distributed in the Wuqia area which is the intersecting area of the western part of the Southern Tianshan Mountains and Kunlun Moutains.The seismic wave energy has different dissipation as it passes through two different crustal media in the propagation process.The records have the geological characteristics of the Tianshan Mountains region，and finally expressed by the Pishan and Lop earthquake records being higher than the calculated result and empirical attenuation relationships.To a certain extent，the attenuation of the Hotan region is faster than the Wuqia region.The Wuqia earthquake and Akto earthquake records did not show this characteristic，while it is lower than the calculated result.

ACKNOWLEDGEMENT

The seismographic records were acquired from the Monitoring Center of the Earthquake Agency of Xinjiang Uygur Autonomous Region.Thanks to assistant researcher Li Jin，and specially thanks for the guidance of research professors Tao Xiaxin and Tao Zhengru.

REFERENCES

Abercrombie R.E.Earthquake source scaling relationships from-1 to 5 MLusing seismograms recorded at 2.5-km depth ［J］.Journal of Geophysical Research， 1995， 100（B12）： 24015-24036.

Atkinson G.M.Attenuation of strong ground motion in Canada from a random vibrations approach ［J］.Bulletin of the Seismological Society of America， 1984， 74（6）： 2629-2653.

Atkinson G.M.， Mereu R.F.The shape of ground motion attenuation curves in Southeastern Canada ［J］.Bulletin of the Seismological Society of America， 1992， 82（5）： 2014-2031.

Atkinson G.M.， Boore D.M.Ground-motion relations for Eastern North America ［J］.Bulletin of the Seismological Society of America， 1995， 85（1）： 17-30.

Ayixiangu·Maimaiti， Nilupa'er·Maimaitusun， Zulipiya·Bake， Xia Aiguo， Jiang Yuan.Research on quality factor Q of small and moderate earthquakes in the western part of the southern Tianshan ［J］.Plateau Earthquake Research， 2015， 27（1）： 23-30 （in Chinese with English abstract）.

Boore D.M.Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra ［J］.Bulletin of the Seismological Society of America， 1983， 73（6A）： 1865-1894.

Boore D.M.， Boatwright J.Average body-wave radiation coefficients［J］.Bulletin of the Seismological Society of America， 1984， 74（5）： 1615-1621.

Boore D.M.， Joyner W.B.Site amplifications for generic rock sites［J］.Bulletin of the Seismological Society of America， 1997， 87（2）： 327-341.

Boore D.M.SMSIM-Fortran programs for simulating ground motions from earthquakes：version 2.0-a Revision of OFR 96-80-A［R］.U.S.Geological Survey Open-File Report， 2002.

Boore D.M.Simulation of ground motion using the stochastic method ［J］.Pure and Applied Geophysics， 2003，160（3）： 635-676.

Brune J.N.Tectonic stress and the spectra of seismic shear waves from earthquakes ［J］.Journal of Geophysical Research， 1970， 75（26）： 4997-5009.

Brune J.N.Correction to“tectonic stress and the spectra， of seismic shear waves from earthquakes”［J］.Journal of Geophysical Research， 1971， 76（20）： 5002.

Chen P.S.，Nuttli O.W.Estimates of magnitudes and short-period wave attenuation of Chinese earthquakes from modified Mercalli intensity data ［J］.Bulletin of the Seismological Society of America， 1984， 74（3）： 957-968.

Cui Anping.Strong Ground Motion Attenuation Relationship for Sichuan-Yunnan Region from Small Earthquake Records by China Earthquake Networks［D］.Master's thesis.Harbin： Institute of Engineering Mechanics，China Earthquake Administration， 2013 （in Chinese with English abstract）.

Hanks T.C.， McGuire R.K.The character of high-frequency strong ground motion ［J］.Bulletin of the Seismological Society of America， 1981， 71（6）： 2071-2095.

Hanks T.C.fmax［J］.Bulletin of the Seismological Society of America， 1982， 72（6A）： 1867-1879.

Jiang Hui.Site Effect Research in Ground Motion Random Simulation ［D］.Doctoral thesis.Beijing： Institute of Geophysics， China Earthquake Administration， 2005 （in Chinese）.

Kanamori H.， Anderson D.L.Theoretical basis of some empirical relations in seismology ［J］.Bulletin of the Seismological Society of America， 1975， 65（5）： 1073-1095.

Li Jin， Wang Qiong， Wu Chuanyong， Xiang Yuan.Preliminary study for seismogenic structure of the Pishan MS6.5 earthquake of July 3， 2015 ［J］.Chinese Journal of Geophysics， 2016， 59（8）： 2859-2870 （in

Chinese with English abstract）.

Li Qicheng， Jing Liping.Discussion about theory of empirical Green's function ［J］.Earthquake Engineering and Engineering Vibration， 2010， 30（3）： 39-44 （in Chinese with English abstract）.

Li Wenqian.Strong Ground Motion Attenuation Relationships for Lanzhou and North China from Small Earthquake Records by China Earthquake Networks［D］.Master's thesis.Harbin： Institute of Engineering Mechanics，China Earthquake Administration， 2014 （in Chinese with English abstract）.

Li Zhihai， Zhu Chengying， Xia Aiguo， Feng Lei， Tang Lanlan.Study on non-elasticity attenuation， site response and source parameters in the Jiashi region， Xinjiang ［J］.Earthquake Research in China， 2010，26（3）： 340-348. （in Chinese with English abstract）.

McGuire R.K.，Hanks T.C.RMS accelerations and spectral amplitudes of strong ground motion during the San Fernando， California earthquake ［J］.Bulletin of the Seismological Society of America， 1980， 70（5）： 1907-1919.

Meng Lingyuan， Zhou Longquan， Zhang Xiaotao， Yang Wen， Yao Qi， Song Jin.Study the characteristics of the Yutian， Xinjiang MS7.3 earthquake， February 12， 2014［J］.Earthquake Research in China， 2014， 30（2）： 159-167 （in Chinese with English abstract）.

Pan Zhensheng， Liu Hui， Gu Meiju， Li Xiaodong， Shen Qi.Preliminary study on source parameters of earthquakes occurred in Keping block， Xinjiang ［J］.Northwestern Seismological Journal， 2010， 32（4）：357-362 （in Chinese with English abstract）.

Papageorgiou A.S.，Aki K.A specific barrier model for the quantitative description of inhomogeneous faulting and the prediction of strong ground motion.I.Description of the model［J］.Bulletin of the Seismological Society of America， 1983， 73（3）： 693-722.

Shearer P.M.Introduction to Seismology［M］.2nd ed.Cambridge： Cambridge University Press，2009.

Tao Zhengru.Ground Motion Attenuation Relationship Based on Stochastic Method ［D］.Doctoral thesis.Harbin：National Graduate Institute for Policy Studies and Building Research Institute，2010 （in Chinese with English abstract）.

Wang Guoxin.Study on Ground Motion Attenuation ［D］.Doctoral thesis.Harbin： Institute of Engineering Mechanics， China Earthquake Administration， 2001 （in Chinese with English abstract）.

Wang Guoxin， Tao Xiaxin.A new two-steps method for fitting ground motion attenuation relationship ［J］.Earthquake Engineering and Engineering Vibration， 2001， 20（1）： 24-28 （in Chinese with English abstract）.

Wang Ji， Yu Yanxiang， Gong Fei， Wang Suyun， Ding Li， Li Liyuan.Suitability analysis of ground motion attenuation in West China—Some of Xinjiang earthquakes as an example ［J］.Earthquake Research in China， 2008， 24（2）： 126-133 （in Chinese with English abstract）.

Yu Yanxiang， Li Shanyou， Xiao Liang.Development of ground motion attenuation relations for the new seismic hazard map of China ［J］.Technology for Earthquake Disaster Prevention， 2013， 8（1）： 24-33 （in Chinese with English abstract）.

Zhang Zhenbin，Tang Lihua.The attenuation relationship of peak ground motion acceleration on soil site at the south Tianshan region in Xinjiang ［J］.Earthquake Research in China， 2010， 26（4）： 401-406 （in Chinese with English abstract）.

Zhao Cuiping， Xia Aiguo， Zheng Sihua， Chen Zhangli.Study on source parameters of earthquakes occurred on the Middle East area of Northern Tianshan in Xinjiang ［J］.Earthquake Research in China，2005， 21（1）： 61-69 （in Chinese with English abstract）.

Zhao Cuiping， Chen Zhangli， Hua Wei， Wang Qincai， Li Zhixiong， Zheng Sihua.Study on source parameters of small to moderate earthquakes in the main seismic active regions， China mainland ［J］.Chinese Journal of Geophysics， 2011， 54（6）： 1478-1489 （in Chinese with English abstract）.

Zhou Yunhao， Chen Zhangli， Miao Fajun.Source process of the 14 November 2001 western Kunlun mountain MS8.1 earthquake ［J］.Acta Seismologica Sinica， 2004， 26（S1）： 9-20 （in Chinese with English abstract）.