এই পেপারটি ETAS মডেল ব্যবহার করে ভূমিকম্পন পূর্বাভাস নিয়ে গবেষণা করে। ঐতিহ্যবাহী MCMC পদ্ধতির গণনা ব্যয়বহুলতা এবং শক্তিশালী পরামিতি সম্পর্কের সমস্যা সমাধানের জন্য, এটি INLA-ভিত্তিক "inlabru" পদ্ধতি ব্যবহার করে বেইজিয়ান অনুমান প্রস্তাব করে। পেপারটি তিনটি পরীক্ষার মাধ্যমে পদ্ধতির কার্যকারিতা যাচাই করে: (১) স্থির পরামিতিগুলি অন্যান্য পরামিতির পশ্চাৎ বিতরণে কীভাবে প্রভাব ফেলে তা অন্বেষণ করে, α এবং K পরামিতির মধ্যে স্পষ্ট পারস্পরিক প্রভাব আবিষ্কার করে; (२) বাস্তব ডেটা এবং সংশ্লেষিত ক্যাটালগের স্বাভাবিকীকৃত ঘটনা-মধ্যবর্তী সময় বিতরণ তুলনা করে; (३) প্রধান কম্পনের আগের ঘটনাগুলি ব্যবহার করে প্রধান কম্পনের পরে দশ সপ্তাহের মধ্যে ভূমিকম্পনের পূর্বাভাস দেয়, গণনা পরীক্ষা এবং ক্রমাগত র্যাঙ্কিং সম্ভাব্যতা স্কোর (CRPS) ব্যবহার করে পূর্বাভাসের নির্ভুলতা মূল্যায়ন করে।
ভূমিকম্পন মানুষের জীবন এবং সম্পত্তির নিরাপত্তার জন্য হুমকিস্বরূপ প্রাকৃতিক দুর্যোগের একটি, যা হঠাৎ ঘটে এবং সময়মত পালানো কঠিন। যদিও বর্তমানে প্রধান কম্পন পূর্বাভাস এই গবেষণার সর্বোচ্চ আদর্শ অর্জন করা সম্ভব নয়, তবে পরবর্তী কম্পন পূর্বাভাসে ভাল অগ্রগতি হয়েছে। নির্ভরযোগ্য পরবর্তী কম্পন পূর্বাভাস সরকার এবং উদ্ধার কর্মীদের কম্পন-পরবর্তী উদ্ধার কাজ আরও ভালভাবে পরিকল্পনা করতে সাহায্য করতে পারে।
পরবর্তী কম্পন পূর্বাসের নির্ভুলতা সরাসরি সম্পর্কিত: ১. ক্ষতিগ্রস্ত ব্যক্তিদের জন্য আরও বেশি বেঁচে থাকার সুযোগ প্রদান করা २. উদ্ধার কর্মীদের জীবন সুরক্ষা আরও ভালভাবে রক্ষা করা ३. কম্পন-পরবর্তী উদ্ধার সম্পদ বরাদ্দ অপ্টিমাইজ করা
१. সর্বাধিক সম্ভাবনা অনুমান (MLE): পদ্ধতি সরাসরি এবং সহজ, কিন্তু পরামিতি অনুমানের অনিশ্চয়তা পেতে পারে না २. মার্কভ চেইন মন্টে কার্লো (MCMC): বেইজিয়ান পদ্ধতি হিসাবে পরামিতি পশ্চাৎ বিতরণ পেতে পারে, কিন্তু নিম্নলিখিত সমস্যা রয়েছে:
'inlabru' পদ্ধতি সমন্বিত নেস্টেড ল্যাপ্লেস আনুমানিকতা (INLA) এর উপর ভিত্তি করে তৈরি, যা পারে:
१. পদ্ধতিগত অবদান: ETAS মডেলের ভূমিকম্পন পূর্বাভাসে INLA-ভিত্তিক inlabru পদ্ধতি প্রয়োগ করা, ঐতিহ্যবাহী MCMC পদ্ধতির গণনা দক্ষতা সমস্যা সমাধান করা
२. পরামিতি পারস্পরিক সম্পর্ক আবিষ্কার: সিস্টেমেটিক পরীক্ষার মাধ্যমে ETAS মডেলে পরামিতিগুলির মধ্যে পারস্পরিক প্রভাব সম্পর্ক প্রকাশ করা, বিশেষত α এবং K পরামিতির শক্তিশালী সম্পর্ক
३. গণনা দক্ষতা উন্নতি: প্রমাণ করা যে α বা K পরামিতি সত্য মানে স্থির করা মডেল ফিটিং সময় অর্ধেকেরও বেশি হ্রাস করতে পারে
४. পূর্বাভাস ক্ষমতা যাচাইকরণ: বাস্তব ভূমিকম্পন ডেটার উপর ভিত্তি করে পূর্বাভাস মূল্যায়ন কাঠামো প্রতিষ্ঠা করা, নির্ভরযোগ্য পূর্বাভাসের জন্য প্রয়োজনীয় ন্যূনতম ডেটা প্রয়োজনীয়তা নির্ধারণ করা
५. তাত্ত্বিক অন্তর্দৃষ্টি: ভাল পরামিতি পশ্চাৎ বিতরণ পাওয়ার জন্য প্রয়োজনীয় শর্তাবলী স্পষ্ট করা, ব্যবহারিক প্রয়োগের জন্য নির্দেশনা প্রদান করা
ইনপুট: ঐতিহাসিক ভূমিকম্পন ঘটনা ক্রম, যা ঘটনার সময় t, অবস্থান s এবং মাত্রা m অন্তর্ভুক্ত করে আউটপুট: ভবিষ্যত সময় উইন্ডোতে ভূমিকম্পন ঘটনার সম্ভাব্যতা বিতরণ পূর্বাভাস সীমাবদ্ধতা: স্ব-উত্তেজিত পয়েন্ট প্রক্রিয়ার ETAS মডেল কাঠামোর উপর ভিত্তি করে
ETAS মডেল Hawkes প্রক্রিয়ার উপর ভিত্তি করে, এর তীব্রতা ফাংশন:
λ_ETAS(t,m|H_t) = (μ + Σ K·e^(α(m_h-M_0))·(t-t_h/c + 1)^(-p))·π(m)
যেখানে:
Taylor সম্প্রসারণ এবং বিনিং কৌশল ব্যবহার করে লগ সম্ভাবনা আনুমানিক করা:
L(θ|H) = -Λ_0(T_1,T_2) - Σ Σ Λ_i(t_j^(bi), t_{j+1}^(bi)) + Σ log λ(t_i|H_{t_i})
१. প্রাথমিক পরামিতি θ_0 ব্যবহার করে Taylor সিরিজ আনুমানিক লগ সম্ভাবনা ফাংশন পান २. লগ প্রাথমিকের সাথে যোগ করুন, পশ্চাৎ বিতরণের মোড θ_1* গণনা করুন ३. লাইন সার্চের মাধ্যমে পরামিতি আপডেট করুন: θ* = αθ_0 + (1-α)θ_1* ४. সংমিশ্রণ পর্যন্ত পুনরাবৃত্তি করুন (পার্থক্য ১% এর কম)
१. বিনিং কৌশল অপ্টিমাইজেশন: স্ব-অভিযোজিত বিনিং ব্যবহার করা, বাম দিকে সংকীর্ণ বিনিং নির্ভুলতা বৃদ্ধি করে, ডান দিকে প্রশস্ত বিনিং গণনা পরিমাণ হ্রাস করে २. বিপরীত সম্ভাব্যতা সমন্বিত রূপান্তর: গাউসিয়ান বিতরণ নমুনা লক্ষ্য প্রাথমিক বিতরণ নমুনায় রূপান্তরিত করা ३. রৈখিক আনুমানিকতা: Taylor সম্প্রসারণের মাধ্যমে দক্ষ সম্ভাবনা ফাংশন আনুমানিকতা বাস্তবায়ন করা
१. Aquila ভূমিকম্পন: সংশ্লেষিত ক্যাটালগ তৈরি এবং পরামিতি বিশ্লেষণের জন্য ব্যবহৃত २. Amatrice ভূমিকম্পন (২০१६-२०१७ ইতালি): তিনটি প্রধান কম্পন সহ বাস্তব ভূমিকম্পন ক্রম অন্তর্ভুক্ত করে ३. সংশ্লেষিত ক্যাটালগ: Aquila ভূমিকম্পন পরামিতির উপর ভিত্তি করে তৈরি, নিয়ন্ত্রিত পরীক্ষার জন্য ব্যবহৃত
१. গণনা পরীক্ষা (N-test): পূর্বাভাসিত ঘটনা সংখ্যার নির্ভুলতা মূল্যায়ন করা
δ_2 = |{N_j|N_j ≤ N_obs, j=1,...,m}|/m
२. ক্রমাগত র্যাঙ্কিং সম্ভাব্যতা স্কোর (CRPS): নির্ভুলতা এবং নির্ভুলতা একসাথে মূল্যায়ন করা
S(F,N_true) = Σ(F(k) - I(N_true ≤ k))²
প্রাথমিক বিতরণ:
প্রাথমিক মান: μ₀=0.5, K₀=0.1, α₀=1, c₀=0.1, p₀=1.1
মূল আবিষ্কার:
গণনা দক্ষতা উন্নতি:
| মডেল কনফিগারেশন | ক্যাটালগ १ | ক্যাটালগ २ | ক্যাটালগ ३ |
|---|---|---|---|
| সম্পূর্ণ মডেল | 41.73s | 46.51s | 44.92s |
| K স্থির | 26.00s | 20.95s | 21.47s |
| α স্থির | 18.25s | 7.75s | 15.15s |
Amatrice ভূমিকম্পন পূর্বাভাস ফলাফল:
ঐতিহাসিক প্রধান কম্পন সংখ্যার প্রভাব: १. কোন প্রধান কম্পন নেই: μ ছাড়া সমস্ত পরামিতির পশ্চাৎ বিতরণ অত্যন্ত খারাপ २. একটি প্রধান কম্পন: মৌলিক ব্যবহারযোগ্য পূর্বাভাস পেতে পারে, কিন্তু নির্ভুলতা সীমিত ३. দুটি প্রধান কম্পন: পূর্বাভাস কর্মক্ষমতা উল্লেখযোগ্যভাবে উন্নত করে, বিশেষত α পরামিতির পশ্চাৎ বিতরণ
পরামিতি অনুমানের সহজ-কঠিন ক্রম: १. μ (সবচেয়ে সহজ): অল্প সংখ্যক ঘটনা অনুমান করা যায় २. c, p (মধ্যম): আরও বেশি ভূমিকম্পন ঘটনা প্রয়োজন, শান্ত ক্যাটালগে অনুমান করা যায় ३. α, K (সবচেয়ে কঠিন): কমপক্ষে একটি সম্পূর্ণ প্রধান কম্পন-পরবর্তী কম্পন ক্রম প্রয়োজন; আরও বেশি প্রধান কম্পন α এর পশ্চাৎ বিতরণ উল্লেখযোগ্যভাবে উন্নত করতে পারে
१. পরামিতি সম্পর্ক: α-K এবং c-p শক্তিশালী সম্পর্ক রয়েছে, যৌথভাবে বিবেচনা করা প্রয়োজন २. ডেটা প্রয়োজনীয়তা: নির্ভরযোগ্য পূর্বাভাসের জন্য কমপক্ষে একটি সম্পূর্ণ প্রধান কম্পন-পরবর্তী কম্পন ক্রম প্রয়োজন ३. পূর্বাভাস উন্নতি: ঐতিহাসিক প্রধান কম্পন সংখ্যা বৃদ্ধি পূর্বাভাস কর্মক্ষমতা উল্লেখযোগ্যভাবে উন্নত করতে পারে ४. গণনা অপ্টিমাইজেশন: প্রাথমিক তথ্যের যুক্তিসঙ্গত ব্যবহার গণনা দক্ষতা উল্লেখযোগ্যভাবে উন্নত করতে পারে
१. ETAS মডেল উন্নয়ন: Ogata (१९८८) এর মূল মডেল থেকে বিভিন্ন সম্প্রসারণ সংস্করণ পর্যন্ত २. পরামিতি অনুমান পদ্ধতি: MLE থেকে MCMC থেকে INLA-ভিত্তিক পদ্ধতির বিবর্তন ३. ভূমিকম্পন পূর্বাভাস মূল্যায়ন: গণনা পরীক্ষা, CRPS এবং অন্যান্য মূল্যায়ন সূচকের উন্নয়ন
१. গণনা দক্ষতা ঐতিহ্যবাহী MCMC পদ্ধতির চেয়ে উল্লেখযোগ্যভাবে ভাল २. সিস্টেমেটিকভাবে পরামিতিগুলির মধ্যে পারস্পরিক সম্পর্ক বিশ্লেষণ করা ३. ব্যবহারিক পূর্বাভাস মূল্যায়ন কাঠামো প্রদান করা
१. পদ্ধতির কার্যকারিতা: inlabru পদ্ধতি ETAS মডেলে চমৎকার কর্মক্ষমতা প্রদর্শন করে, গণনা দক্ষতা উল্লেখযোগ্যভাবে উন্নত করে २. পরামিতি সম্পর্ক: α এবং K এর পারস্পরিক নির্ভরতা মডেল কর্মক্ষমতা প্রভাবিত করার মূল কারণ ३. পূর্বাভাস ক্ষমতা: মডেল নির্ভরযোগ্য পরবর্তী কম্পন পূর্বাভাস ক্ষমতা রাখে, পূর্বাভাস গুণমান ঐতিহাসিক ডেটা সমৃদ্ধতার সাথে উন্নত হয় ४. ডেটা প্রয়োজনীয়তা: নির্ভরযোগ্য পূর্বাভাসের জন্য কমপক্ষে একটি প্রধান কম্পন-পরবর্তী কম্পন ক্রম প্রয়োজন
१. ডেটা সীমাবদ্ধতা: সংশ্লেষিত ক্যাটালগ শুধুমাত্র একক প্রধান কম্পন অন্তর্ভুক্ত করে, সিদ্ধান্তের সর্বজনীনতা সীমিত করে २. ঘটনা মধ্যবর্তী সময় সংজ্ঞা: শুধুমাত্র সংলগ্ন ঘটনা মধ্যবর্তী সময় বিবেচনা করে, প্রধান কম্পন সমষ্টির অভ্যন্তরীণ কাঠামো বিবেচনা করে না ३. গণনা চ্যালেঞ্জ: প্রধান কম্পন সংখ্যা বৃদ্ধির সাথে মডেল ফিটিং সময় উল্লেখযোগ্যভাবে বৃদ্ধি পায়, সংমিশ্রণ কঠিন ४. বৈশিষ্ট্য বিশ্লেষণ অপর্যাপ্ত: প্রধান কম্পন সংখ্যা ছাড়া, অন্যান্য ভূমিকম্পন বৈশিষ্ট্যের পশ্চাৎ বিতরণে প্রভাব এখনও সম্পূর্ণভাবে অন্বেষণ করা হয়নি
१. পরামিতি আপডেট কৌশল: বৃদ্ধিশীল শেখার পদ্ধতি বিকাশ করা, পূর্ববর্তী প্রধান কম্পনের পশ্চাৎ বিতরণ পরবর্তীটির প্রাথমিক হিসাবে ব্যবহার করা २. বহু-প্রধান কম্পন পরীক্ষা: আরও বেশি প্রধান কম্পন সহ সংশ্লেষিত ক্যাটালগে পরীক্ষা পুনরাবৃত্তি করা ३. ঘটনা সমষ্টি বিশ্লেষণ: প্রধান কম্পন সমষ্টির অভ্যন্তরীণ এবং মধ্যবর্তী ঘটনা মধ্যবর্তী সময় বিতরণ গবেষণা করা ४. বৈশিষ্ট্য প্রভাব গবেষণা: বিভিন্ন ভূমিকম্পন বৈশিষ্ট্যের মডেল কর্মক্ষমতায় প্রভাব সিস্টেমেটিকভাবে বিশ্লেষণ করা
१. পদ্ধতি উদ্ভাবনী: INLA পদ্ধতি ETAS মডেলে সফলভাবে প্রয়োগ করা, ব্যবহারিক গণনা সমস্যা সমাধান করা २. পরীক্ষা ডিজাইন সম্পূর্ণ: তিনটি পরীক্ষা বিভিন্ন কোণ থেকে পদ্ধতির কর্মক্ষমতা সম্পূর্ণভাবে মূল্যায়ন করা ३. পরামিতি বিশ্লেষণ গভীর: ETAS মডেল পরামিতিগুলির মধ্যে পারস্পরিক সম্পর্ক সিস্টেমেটিকভাবে প্রকাশ করা ४. ব্যবহারিক মূল্য উচ্চ: ব্যবহারিক ভূমিকম্পন পূর্বাভাসের জন্য স্পষ্ট ডেটা প্রয়োজনীয়তা নির্দেশনা প্রদান করা ५. ফলাফল বিশ্বাসযোগ্যতা উচ্চ: বাস্তব ভূমিকম্পন ডেটা ব্যবহার করে যাচাইকরণ, সিদ্ধান্ত ব্যবহারিক তাৎপর্য রয়েছে
१. তাত্ত্বিক বিশ্লেষণ অপর্যাপ্ত: পরামিতি সম্পর্কের তাত্ত্বিক ব্যাখ্যার অভাব २. পরীক্ষা পরিসীমা সীমিত: প্রধানত ইতালি ভূমিকম্পন ডেটার উপর ভিত্তি করে, ভৌগোলিক প্রতিনিধিত্ব অপর্যাপ্ত ३. পদ্ধতি তুলনা অপর্যাপ্ত: অন্যান্য আধুনিক বেইজিয়ান পদ্ধতির সাথে বিস্তারিত তুলনা নেই ४. অনিশ্চয়তা পরিমাণীকরণ: পূর্বাভাস অনিশ্চয়তার বিশ্লেষণ এখনও যথেষ্ট গভীর নয়
१. একাডেমিক অবদান: ভূমিকম্পন বিজ্ঞান এবং পরিসংখ্যান ক্রস-ডিসিপ্লিনারি ক্ষেত্রে নতুন পদ্ধতিবিদ্যা প্রদান করা २. ব্যবহারিক মূল্য: সরাসরি ব্যবহারিক ভূমিকম্পন পূর্বাভাস এবং ঝুঁকি মূল্যায়নে প্রয়োগ করা যায় ३. পুনরুৎপাদনযোগ্যতা: খোলা উৎস R প্যাকেজ inlabru এর উপর ভিত্তি করে, ভাল পুনরুৎপাদনযোগ্যতা রয়েছে ४. সম্প্রসারণ সম্ভাবনা: পদ্ধতি অন্যান্য স্ব-উত্তেজিত পয়েন্ট প্রক্রিয়া প্রয়োগ ক্ষেত্রে সম্প্রসারণ করা যায়
१. পরবর্তী কম্পন পূর্বাভাস: কম্পনের পরে স্বল্পমেয়াদী পরবর্তী কম্পন কার্যকলাপ পূর্বাভাস २. ভূমিকম্পন ঝুঁকি মূল্যায়ন: ঐতিহাসিক ডেটার উপর ভিত্তি করে ভূমিকম্পন ঝুঁকি মূল্যায়ন ३. জরুরি প্রতিক্রিয়া পরিকল্পনা: কম্পন-পরবর্তী উদ্ধারের জন্য বৈজ্ঞানিক সিদ্ধান্ত সহায়তা প্রদান করা ४. ভূমিকম্পন বিজ্ঞান গবেষণা: ভূমিকম্পন কার্যকলাপের সময়-স্থান বিতরণ নিয়ম বোঝার জন্য ব্যবহার করা
প্রধান তথ্যসূত্র অন্তর্ভুক্ত করে:
সামগ্রিক মূল্যায়ন: এটি ভূমিকম্পন পূর্বাভাস ক্ষেত্রে গুরুত্বপূর্ণ ব্যবহারিক মূল্য সহ একটি পেপার, ETAS মডেলের গণনা দক্ষতা সমস্যা সফলভাবে সমাধান করে, এবং সিস্টেমেটিক পরীক্ষার মাধ্যমে মূল্যবান পরামিতি বিশ্লেষণ এবং পূর্বাভাস মূল্যায়ন প্রদান করে। যদিও তাত্ত্বিক গভীরতা এবং পরীক্ষা বিস্তৃতিতে উন্নতির জায়গা রয়েছে, তবে এর পদ্ধতির ব্যবহারিকতা এবং সিদ্ধান্তের কার্যকারিতা এটিকে এই ক্ষেত্রের একটি গুরুত্বপূর্ণ অবদান করে তোলে।