Mô phỏng băng gia tốc bằng phương pháp ngẫu nhiên

ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - HÀ VĨNH LONG MÔ PHỎNG BĂNG GIA TỐC BẰNG PHƢƠNG PHÁP NGẪU NHIÊN LUẬN VĂN THẠC SĨ KHOA HỌC Hà Nội – 2015 ĐẠI HỌC QUỐC GIA HÀ NỘI TRƢỜNG ĐẠI HỌC KHOA HỌC TỰ NHIÊN - HÀ VĨNH LONG MÔ PHỎNG BĂNG GIA TỐC BẰNG PHƢƠNG PHÁP NGẪU NHIÊN Chuyên ngành: Vật lý Địa cầu Mã số: 60440111 LUẬN VĂN THẠC SỸ KHOA HỌC NGƢỜI HƢỚNG DẪN KHOA HỌC: TS NGUYỄN LÊ MINH HÀ NỘI – 2015 LỜI CÁM ƠN Trong q trình nghiên cứu hồn thành luận văn , ngồi nỡ lực thân, nhận đƣợc quan tâm, giúp đỡ nhiệt tình nhiều quan, tổ chức, cá nhân, bạn bè đồng nghiệp Nhân dịp xin gửi lời cảm ơn chân thành quan tâm q báu Trƣớc hết tơi xin trân trọng cảm ơn thầy, cô giáo Khoa Vật lý, Trƣờng Đại học Khoa học Tự nhiên - Đại học Quốc gia Hà Nội dạy dỗ, giúp đỡ tạo điều kiện thuận lợi cho suốt thời gian học tập trƣờng Tơi xin bày tỏ lòng biết ơn sâu sắc đế n TS Nguyễn Lê Minh hƣớng dẫn, bảo phƣơng pháp làm việc, nghiên cứu giúp đỡ tơi hồn thành luận văn Xin trân trọng cảm ban lãnh đạo Viện Vật lý Địa cầu, Phòng Quan sát Động đất, bạn bè, đồng nghiệp ủng hộ suốt thời gian học tập Tôi xin bày tổ giúp đỡ khoa học tài của: Nhiệm vụ hợp tác quốc tế khoa học công nghệ cấp Viện Hàn lâm Khoa học Công nghệ Việt Nam “Cấu trúc địa chấn phía biển Đơng Việt Nam ý nghĩa với hoạt động kiến tạo khu vực Đông Nam Á” thời gian thực luận văn Cuối cùng, xin chân thành cảm ơn gia đình, bạn bè ghiệp ln tạo điều kiện giúp đỡ, động viên suốt thời gian nghiên cứu hoàn thiện luận văn./ Xin trân trọng cảm ơn! Hà Nội, ngày tháng 12 năm 2015 Tác giả Hà Vĩnh Long LỜI CAM ĐOAN Tôi xin cam đoan cơng trình nghiên cứu riêng Các số liệu, kết nêu luận văn trung thực chƣa đƣợc công bố cơng trình khác Tơi xin cam đoan thơng tin trích dẫn luận văn có xuất xứ, nguồn gốc cụ thể Việc sử dụng thơng tin q trình nghiên cứu hoàn toàn hợp lệ Hà Nội, ngày tháng 12 năm 2015 Tác giả Hà Vĩnh Long MỤC LỤC DANH MỤC CHỮ VIẾT TẮT DANH MỤC HÌNH DANH MỤC BẢNG, BIỂU MỞ ĐẦU CHƢƠNG TỔNG QUAN 1.1 Phƣơng pháp tất định 1.2 Phƣơng pháp hàm Green thực nghiệm 1.3 Phƣơng pháp ngẫu nhiên 10 CHƢƠNG PHƢƠNG PHÁP NGẪU NHIÊN VÀ CHƢƠNG TRÌNH SMSIM SỬ DỤNG PHƢƠNG PHÁP NGẪU NHIÊN…………… ………………………… 12 2.1 Lịch sử phát triển phƣơng pháp 12 2.2.Cơ phƣơng pháp ngẫu nhiên 13 2.2.1 Thuyết dao động ngẫu nhiên 13 2.2.2 Phương pháp ngẫu nhiên tốn mơ dao động mạnh 16 2.3 Chƣơng trình SMSIM sử dụng phƣơng pháp ngẫu nhiên 27 2.3.1 Tham số xác định mô hình 30 2.3.2 Tham số đầu vào ứng với kịch 30 2.3.3 Tham số đầu chương trình 30 CHƢƠNG ÁP DỤNG PHƢƠNG PHÁP NGẪU NHIÊN CHO ĐIỀU KIỆN ĐỘNG ĐẤT THỰC TẾ 38 3.1 Động đất khu vực California 38 Tham số sử dụng để mô động đất khu vực California 42 3.2 Chuỗi động đất Điện Biên 2001 44 Tham số sử dụng để mô động đất Điện Biên 46 3.3 Kết thảo luận 49 3.3.1 Động đất khu vực California 49 3.3.2 Chuỗi động đất Điện Biên 2001 61 KẾT LUẬN 69 TÀI LIỆU THAM KHẢO 70 PHỤ LỤC 75 DANH MỤC CHỮ VIẾT TẮT PGA: Peak Ground Acceleration PGV: Peak Ground Verlocity RVT: Random Vibration Theory FAS: Fourier Acceleration Spectrum SDOF: Single Degree of freedom Oscillator transfer Function GMT: Greenwich Mean Time RMS: Root-mean-square nnk: Những ngƣời khác DANH MỤC HÌNH Hình 1.1 Phổ gia tốc đƣợc thành lập từ tập hợp băng ghi gia tốc………………….5 Hình 1.2 Hàm lọc F(t) sử dụng để hiệu chỉnh dành cho hàm phụ thuộc thời gian vận tốc (…)…………………………………………………………………………9 Hình 1.3 Diện tích đứt gãy trận động đất lớn đƣợc coi nhƣ hệ đứt gãy thứ cấp kiện động đất nhỏ (…)…………………………………………………… 10 Hình 2.1 Mơ tả cách xác định Ne, Nz tín hiệu thuộc tần số khác (Trái) Chuỗi thời gian có tần số Hz; (…)……………………………………………….16 Hình 2.2 Mơ tả dao động thành phần tổng hợp hiệu ứng của:Nguồn, đƣờng truyền,hiệu ứng nền……………………………………………………………………17 Hình 2.3 Các băng ghi mơ hình giả định lớp vỏ lớp miền Trung Hoa Kỳ đƣợc sử dụng tính tốn hiệu ứng đƣờng truyền………………………………… 20 Hình 2.4 Mơ tả hàm lan truyền hình học Z(R)cho miền Trung Tây Bắc Hoa Kỳ 21 Hình 2.5 Mơ tả tốn tử suy giảm Q dành cho vùng Đông Bắc Hoa Kỳ………………22 Hình 2.6: Tổng hợp hiệu ứng khuếch đại suy giảm dành cho chuỗi tham số suy giảm κ0……………………………………………………………………………26 Hình 2.7: Tổ hợp khuếch đại tham số suy giảm κ0 ………………………26 Hình 2.8 Sơ đồ phác thảo hoạt động chƣơng trình SMSIM………………………29 Hình 3.1 Bản đồ california với đứt gãy (trái), đồ chấn tâm trận động đất sở liệu dao động mạnh California (phải) .39 Hình 3.2 Bản đồ chấn tâm trạm ghi số liệu dao động mạnh mạnh Hector Mine, California Thời gian 1999-10-16 12:46:45 GMT……… …………………….41 Hình 3.3 Bản đồ dao động (trái) đồ đƣờng đẳng chấn (phải) trận động đất Hector Mine……………………………………………………………………………… 42 Hình 3.4 So sánh giá trị thực giá trị mơ cho chƣơng trình đƣợc xây dựng Atkinson Silva (2000)…………………………………………………………… 44 Hình 3.5 Động đất Điện Biên ngày 19/2/2001(…)…………………… 46 Hình 3.6 Đƣờng cong suy giảm cho miền Bắc Việt Nam số khu vực khác… 47 Hình 3.7 (a) Đƣờng cong thể giá trị gia tốc đỉnh tƣơng ứng với khoảng cách; (b) phân bố mạng trạm so với chấn tâm trận động đất Hector Mine M=7.1 .53 Hình 3.8 So sánh kết băng gia tốc mô băng gia tốc quan sát tƣơng ứng với động đất Hector(…)….……………………………………………………….54 Hình 3.9 So sánh kết băng gia tốc mô băng gia tốc quan sát tƣơng ứng với động đất Hector(…)….……………………………………………………….55 Hình 3.10 So sánh kết băng gia tốc mô băng gia tốc quan sát tƣơng ứng với động đất Hector(…)….……………………………………………………….56 Hình 3.11 Đƣờng cong thể giá trị gia tốc đỉnh tƣơng ứng với khoảng cách trận động đất Parkfield M=6.0 ngày 28/12/2004…………………………………… 58 Hình 3.12 (a) Bản đồ dịch chuyển đƣợc công bố; (b) phân bố mạng trạm số liệu so với chấn tâm động đất trận động đất Parkfield M=6.0 ngày 28/12/2004 58 Hình 3.13 So sánh kết băng gia tốc mô băng gia tốc quan sát tƣơng ứng với động đất Parkfield(…)……………………………………………………….59 Hình 3.14 So sánh kết băng gia tốc mô băng gia tốc quan sát tƣơng ứng với động đất Parkfield(…)……………………………………………………….60 Hình 3.15 So sánh kết băng gia tốc mô (∆σ=150 bar ;Q=200f 0.3) băng gia tốc quan sát ( ) .63 Hình 3.16 So sánh kết băng gia tốc mô (∆σ=150 bar ;Q=200f 0.3) băng gia tốc quan sát ( )…………………………………………………………… 64 Hình 3.17 Phổ biên độ Fourier phổ gia tốc băng động đất quan sát băng động đất mô tƣơng ứng với ∆σ=150, 100 50 bar ( ) 65 Hình 3.18 Phổ biên độ Fourier phổ gia tốc băng động đất quan sát băng động đất mô tƣơng ứng với ∆σ=150, 100 50 bar ( ) 66 (a) (b) (c) (d) Hình 3.16 So sánh kết băng gia tốc mô (∆σ=50 bar ;Q=200f0.3) băng gia tốc quan sát trân động đất tương ứng với động đất ngày 19/2/2001 có M=5.3, khoảng cách R= 66 km ghi trạm Tuần Giáo (a) So sánh băng sóng; (b) So sánh phổ biên độ Fourier; (c) So sánh cường độ Arias;(d) So sánh phổ gia tốc với độ suy giảm 5% 64 (a1) (b1) (a2) (b2) Hình 3.17 Phổ biên độ Fourier phổ gia tốc băng động đất quan sát băng động đất mô tương ứng với ∆σ=150, 100 50 bar trận động đất ngày 19/2/2001 có M=5.3, khoảng cách R= 19 km ghi trạm Điện Biên (a1) Phổ biên độ Fourier với Q=180f0.5; (b1) Phổ gia tốc với Q=200f0.3; (a2) Phổ biên độ Fourier với Q=180f0.5;(b2) Phổ gia tốc với Q=200f0.3 65 (a1) (b1) (a2) (b2) Hình 3.18 Phổ biên độ Fourier phổ gia tốc băng động đất quan sát băng động đất mô tương ứng với ∆σ=150, 100 50 bar trận động đất ngày 19/2/2001 có M=5.3, khoảng cách R= 166 km ghi trạm Tuần Giáo (a1) Phổ biên độ Fourier với Q=180f0.5; (b1) Phổ gia tốc với Q=200f0.3; (a2) Phổ biên độ Fourier với Q=180f0.5;(b2) Phổ gia tốc với Q=200f0.3 66 Việc thử nghiệm mô chuỗi động đất Điện Biên 2001 với tham số ƣớc lƣợng cho thấy giá trị đƣợc tính tốn với hàm Q = 180f0.5 Q = 200f0.3 có kết tƣơng đồng với băng ghi thực Giá trị mô trạm Điện Biên lại cho thấy phù hợp với độ sụt giảm ứng suất nguồn 150 bar với M=5.3 100 bar với M=4.8 4.7 Điều cho thấy mối tƣơng quan sụt giảm ứng suất độ lớn động đất Từ hình 3.14 đến 3.17 cho thấy kết mô trạm Điện Biên có phù hợp với kết quan sát Điều thể qua phù hợp giá trị PGA, hình dạng băng sóng, phổ biên độ Fourier, phổ gia tốc cƣờng độ Arias Ngƣợc lại, giá trị tính tốn trạm Tuần Giáo lại cho thấy phù hợp với giá trị 50 bar Kết so sánh trạm Tuần Giáo cho thấy phù hợp không cao 150 100 bar, khoảng cách dƣới 70 km Ta thấy rõ cƣờng độ Arias hình 3.15 phổ biên độ phổ gia tốc hình 3.17 Điều giải thích thông tin điều kiện trạm Tuần Giáo chƣa rõ ràng (nền trạm Tuần Giáo xảy động đất M=5.3 nên đá) Từ kết áp dụng cho hai trạm cho thấy việc sử dụng nguồn điểm cho thấy tiện dụng điều kiện vùng nguồn chƣa đƣợc xác định cách chi tiết xác Ngồi ra, việc tính tốn đến tham số quan trọng nguồn nhƣ độ sụt giảm ứng suất, hoạt động phá hủy vùng nguồn v.v quan trọng Ngoài ra, quan điểm khác xạ sóng địa chấn, lan truyền sóng địa chấn cần đƣợc xem xét đến cách cụ thể Do đó, để có đƣợc kết mơ tin cậy cần thiết phải có kết nghiên cứu chi tiết nguồn, mơi trƣờng lan truyền sóng, điều kiện mang tính địa phƣơng Qua thảo luận cho thấy chƣơng trình SMSIM có khả việc xây dựng sở liệu băng ghi gia tốc động đất phạm vi lớn khoảng cách không gần nguồn Tuy vậy, thực tế việc có đƣợc giá trị khu vực gần nguồn động đất lớn thƣờng khó xảy Do đó, 67 đánh giá chƣơng trình SMSIM phƣơng pháp ngẫu nhiên đáng tin cậy cho việc xây dựng sở liệu gia tốc điều kiện mạng trạm quan trắc không đầy đủ 68 KẾT LUẬN Trong luận văn này, học viên đạt đƣợc số kết sau: - Tìm hiểu sở lý thuyết ngẫu nhiên, phƣơng pháp mô băng gia tốc sử dụng phƣơng pháp ngẫu nhiên Tìm hiểu chƣơng trình SMSIM, tham số mơ hình chƣơng trình mô dao động mạnh sử dụng phƣơng pháp ngẫu nhiên - Áp dụng mô thành công băng ghi gia tốc tần số lớn hớn Hz cho khu vực California dựa thông số đƣợc cơng bố Theo đó, kết so sánh hình dạng băng ghi, giá trị gia tốc cực đại, phổ biên độ Fourier, phổ gia tốc, cƣờng độ Arias tƣơng đối phù hợp Bƣớc đầu thử nghiệm mô băng gia tốc cho động đất Điện Biên với độ lớn M=5.3, M=4.8 M=4.7 Các kết thu đƣợc ban đầu cho thấy tiềm áp dụng phƣơng pháp việc xây dựng sở liệu dao động mạnh dành cho Việt Nam Từ trình nghiên cứu áp dụng, học viên nhận thấy: Để thực tốt việc áp dụng phƣơng pháp ngẫu nhiên cho điều kiện thực tế Việt Nam, cần có nghiên cứu cụ thể tham số đƣợc sử dụng trình mô nhƣ hoạt động nguồn động đất, hàm lan truyền sóng Z(R), hàm chất lƣợng Q(f), hệ số suy giảm 𝜅, hàm khuếch đại A(f) Các hàm đƣợc xác định thông quan nghiên cứu mang tính thực nghiệm Phƣơng pháp ngẫu nhiên chƣa thể xây dựng đƣợc băng ghi gia tốc giả định dải tần số dƣới Hz Do đó, để có đƣợc sở liệu dao động mạnh cách hiệu xác nên tiến hành sử dụng phƣơng pháp lai (hay phƣơng pháp tổ hợp) phƣơng pháp ngẫu nhiên phƣơng pháp tất định, phƣơng pháp ngẫu nhiên nhiều phƣơng pháp khác để có đƣợc tham số cách hoàn thiện nhiều giải tần số khác mà nhà thiết kế kháng chấn quan tâm./ 69 TÀI LIỆU THAM KHẢO Tài liệu tiếng Việt Bộ xây dựng (2006), TCXDVN 375 : 2006, thiết kế cơng trình chịu động đất, Nhà xuất Xây dựng Lê Tử Sơn, Nguyễn Quốc Dũng (2003), Kết quan sát gia tốc Việt Nam Tạp trí khoa học Trái đất T.25, 1, 78-85 Lê Tử Sơn (2004), Động đất Điện Biên Ms=5.3, ngày 19-02-2001 Tạp chí khoa học Trái đất T.26, 1, 112-121 Nguyễn Ánh Dƣơng (2004), Mô chuyển động mạnh phương pháp hàm Green thực nghiệm, Luận văn thạc sỹ khoa học, Trƣờng đại học Khoa học tự nhiên, Đại học Quốc gia Hà Nội Phạm Văn Thục (2007), Địa chấn học động đất Việt Nam, Nhà xuất khoa học tự nhiên công nghệ Trầ n Thi My ̣ ̃ Thành (2008), Mô phỏng băng gia tố c dao động nề n cho thiế t kế kháng chân Việt Nam , Tuyể n tâ ̣p các công trin ̀ h nghiên cƣ́u Vâ ̣t lý Điạ cầ u 2008, Viê ̣n Vâ ̣t lý Điạ cầ u – Viê ̣n Khoa ho ̣c và Công nghê ̣ Viê ̣t Nam Tiếng Anh Aki, K (1966), Generation and Propagation of G Waves from the Niigata Earthquake of June 16, 1964 Part Estimation of Earthquake Moment, Released Energy, and Stress-strain Drop from the G-Wave Spectrum, Bull Earthq Res Inst 44, 73–88 Aki, K (1967), Scaling law of seismic spectrum, J Geophys Res, 72, 12171231 Aki, K (1980), Attenuation of Shear-waves in the Lithosphere for Frequencies from 0.05 to 25 Hz, Phys.Earth Planet Inter 21, 50–60 70 10 Anderson, J G and Hough, S E (1984), A Model for the Shape of the Fourier Amplitude Spectrum of Acceleration at High Frequencies, BSSA 74, 1969– 1993 11 Arias, A A (1970), Measure of Earthquake Intensity.In Seismic Design for Nuclear Power Plants (Robert J.Hansen, ed.), (The M.I.T Press, Cambridge, Mass 1970) pp 438–483 12 Atkinson, G M and D M Boore (1995) Ground-motion relations for eastern North America, BSSA 85, 17-30 13 Atkinson, G M and Hanks, T C (1995), A High-frequency Magnitude Scale, BSSA 85, 825–833 14 Atkinson, G M and Silva, W (2000), Stochastic Modeling of California Ground Motions, BSSA 90, 255–274 15 Boore, D M (1983), Stochastic simulation of high-frequency ground motions based on seismological models of the radiated spectra, BSSA.73, 1865–1894 16 Boore, D M (1984), Use of Seismoscope Records to Determine ML and Peak Velocities, BSSA 74, 315–324 17 Boore, D M.and Boatwright, J (1984), Average Body-wave Radiation Coefficients, BSSA 74, 1615–1621 18 Boore, D M and Joyner, W B (1984), A Note on the Use of Random Vibration Theory to Predict Peak Amplitudes of Transient Signals, BSSA 74, 2035–2039 19 Boore, D M and Joyner, (1997), Site amplifications for Generic Rock Sites, BSSA 87, 327–341 20 Boore, D M (2003) Simulation of ground motion using the stochasticmethod, Pure Appl Geophys 160 635–675 71 21 Boore, D M (2005) SMSIM—Fortran programs for simulating ground motions from earthquakes: Version 2.3—A revision of OFR 96-80-A, U.S Geol Surv Open-File Rept 00-509, 55 pp 22 Brune, J N (1970) Tectonic stress and the spectr a of seismic shear waves from earthquakes, J Geophys Res., 75, 4997-5009 23 Brune, J N (1971) Correction, J Geophys Res 76, 5002 24 Cong, et al, (2003), Distribution of Lg coda Q in the Chinese continent and its adjacent region, Science in China (series D) 46(6), 530-539 25 Guzman, R.A and P.C Jennings (1976) Design spectra for nuclear power plants, Journal of the Power Division, ASCE, 102, No P02, Proc Paper 12521, 165-178 26 Hanks, T C (1982), fmax, BSSA 72, 1867–1879 27 Hanks, T C and Boore, D M (1984), Moment-Magnitude Relations in Theory and Practice, J Geophys.Res 89, 6229–6235 28 Hanks, T C., and R K McGuire (1981) The character of high-frequency strong ground motion, BSSA 71, 2071-2095 29 Hanks, T C and Kanamori, H (1979), A Moment Magnitude Scale, J Geophys Res 84, 2348–2350 30 Haskell, N.A (1966), Part II A statistical source model, BSSA 56, 125-140 31 Hartzell, S H (1978), Earthquake Aftershocks as Green’s Functions, Geophis Res Lett., 5, 32 Heaton, T.H (1990), Overview of Seismological Methods for the Synthesis of Strong Ground Motion, Proceedings of the EPRI/Stanford/USGS workshop on modeling ground motion at close distances, held at Palo Alto, CA, Sept 33 Irikura, K (1983), Semi-empirical estimation of strong ground motions during large earthquake, BDPRI, Kyoto Univ, 33 Part 2, No 298, 63-104 72 34 Irikura, K (1986), Prediction of Strong Acceleration Motion using Empirical Green’s Function, Proc th Japan Earthquake Engineering Symposium, 151156 35 Irikura, K and H Miyake (2011), Lecture Note on Strong Motion Seismology, IISEE-UNESCO Lecture Notes Archive, S1-100-2010, 74pp 36 Kanamori, H (1977), The energy release in great earthquakes, J Geophys Res 82 (20): 2981–2876 37 Kanamori, H and Anderson, D.L (1975) Amplitude of the earth’s free oscillations and long period characteristics of the earthquake source J Geophys Res 82, 1075-1078 38 Le Minh Nguyen et al (2011) The first ML scale for North of Vietnam, Journal of Asian Earth Sciences 40 279–286 39 McGuire, R K et al (1984) Spectral estimates of seismic shear waves, BSSA 74, 1427-1440 40 Papageprgiou, A S and Aki, K (1983b), A Specific Barrier Model for the Quantitative Description of Inhomogeneous Faulting and the Prediction of Strong Ground Motion Part II Applications of the Model, BSSA 73, 953– 978 41 Silva, W.J and Lee, K (1987), WES RASCAL code for Synthesizing Earthquake Ground Motions, State-of-the-Art for Assessing Earthquake Hazards in the United States, Report 24, U.S Army Engineers Waterways Experiment Station, Misc Paper S-73-1 42 Silva, W J Et al (1997), Description and Validation of the Stochastic Ground Motion Model, Final Report Brookhaven National Laboratory, Associated Universities, Inc Upton, New York 73 43 Tran Viet Hung and Osamu Kiyomiya (2013) Source parameter estimation and stochastic ground motion based on recorded accelerograms in northwestern Vietnam, Journal of Earthquake Engineering, 17:3, 304-332 44 Wells, D and K Coppersmith (1994) New empirical relationships among magnitude, rupture length, rupture width, rupture area, and surface displacement, BSSA 84, 974–1002 45 Xu, Y et al, (2010) Preliminary high-frequency ground-motion scaling in Yunnan and Southern Sichuan, China BSSA 100(5B), 2508–2517 46 Youngs, R R and Silva, W J (1992), Fitting the ω-2 Brune SourceModel to California Empirical Strong Motion Data (abs.), Seism Res Lett 63, 34 Trang Web 47 Số liệu dao động mạnh California http://www.strongmotioncenter.org/ 48 Số liệu băng ghi động đất toàn cầu: http://ds.iris.edu/wilber3/ 49 Trang thơng tin động đất cục địa chấn hoa kì: http://www.usgs.gov/ 50 Cơ sở liệu danh mục động đất giới ISC: http://www.isc.ac.uk/iscbulletin/search/catalogue/ 51 Trang web cá nhân D M Boore: http://daveboore.com/ 52 Trang chủ bang California: http://www.conservation.ca.gov/ 53 Tiêu chuẩn xây dựng cơng trình Hoa Kỳ: http://peer.berkeley.edu/ngawest2/databases/ 74 PHỤ LỤC FILE THAM SỐ ĐẦU VÀO CỦA CHƢƠNG TRÌNH SMSIM A-1 AS00.PARAMS Params for AS2000 (Atkinson and Silva, BSSA 90, 255 274) source, applied to WNA !rho, beta, prtitn, radpat, fs: 2.8 3.5 0.707 0.55 2.0 !spectral shape: source number, pf_a, pd_a, pf_b, pd_b 2.0 1.0 0.0 0.0 !spectral scaling: ! stressc, dlsdm, fbdfa, amagc, c1_fa, c2_fa, amagc4fa, c1_eps, c2_eps, amagc4eps 100.0 0.0 4.0 7.0 0.0 0.0 0.0 0.0 0.0 0.0 !finite_fault factor specification: 0 0.0 0.0 0.0 0.0 0.0 ! No f_ff !Geometrical spreading option: !Parameters for the frequency dependent gsprd: 10.0 0.0 0.0 0.0 0.0 0.0 ! Use these for option !gsprd: r_ref, nsegs, (rlow(i), a_s, b_s, m_s(i)) (Usually set r_ref = 1.0 km) 1.0 1.0 -1.0 0.0 6.5 40.0 -0.5 0.0 6.5 !q: fr1, Qr1, s1, ft1, ft2, fr2, qr2, s2, c_q 1.0 180 0.45 1.0 1.0 1.0 180 0.45 3.5 !source duration: weights of 1/fa, 1/fb 0.5 0.5 !path duration: nknots, (rdur(i), dur(i), slope of last segment 0.0 0.0 0.05 !crustal amplification, from the source to the site (note that this can include ! local site amplification): namps, (famp(i), amp(i)) 11 75 0.01 1.00 0.09 1.10 0.16 1.18 0.51 1.42 0.84 1.58 1.25 1.74 2.26 2.06 3.17 2.25 6.05 2.58 16.6 3.13 61.2 4.00 !site diminution parameters: fmax, kappa, dkappadmag, amagkref 0.0 0.03 0.0 0.0 !low-cut filter parameters: fcut, nslope (=4, 8, 12, etc) 0.04 !rv params: zup, eps_int (integration accuracy), amp_cutoff (for fup), osc_crrctn(0=no correction; ! 1=BJ84;2=LP99; 3=BT Drms/Dex, file 1; 4=BT Drms/Dex, file 2; 5=average of BT file & file 2) 10.0 0.00001 0.001 !window params: idxwnd(0=box,1=exp), tapr(