2025-11-24T05:19:17.751702

Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry

Bobrovnikov, Jones, Prasanna et al.
We discuss the following inverse problem: given the run-up data of a tsunami wave, can we recover its initial shape? We study this problem within the framework of the non-linear shallow water equations, a model widely used to study tsunami propagation and inundation. Previously, it has been demonstrated that in the case of infinite sloping bathymetry, it is possible to recover the initial water displacement and velocity from shoreline readings \cite{Rybkin23,Rybkin24,Rybkin25}. We consider a finite sloping bathymerty. We show that it is possible to recover boundary conditions (water displacement and velocity) on a virtual buoy from the shoreline data. Further, we discuss stitching together the shallow water equations and the Boussinesq equation in a more complex piece-wise sloping bathymetry in order to recover the initial conditions, while incorporating the dispersion to our model.
academic

अरैखिक तरंग का पुनर्निर्माण: बॉय डेटा से तटीय डेटा और टुकड़ों में ढलान वाली समुद्र तल के लिए सुनामी व्युत्क्रम समस्या में अनुप्रयोग

मूल जानकारी

  • पेपर ID: 2510.14177
  • शीर्षक: Reconstruction of the non-linear wave at a buoy from shoreline data and applications to the tsunami inverse problem for piece-wise sloping bathymetry
  • लेखक: Oleksandr Bobrovnikov¹, Madison Jones², Shriya Prasanna³, Josiah Smith², Alexei Rybkin¹, Efim Pelinovsky⁴'⁵
  • लेखक संस्थान: ¹अलास्का फेयरबैंक्स विश्वविद्यालय, ²कोलोराडो बोल्डर विश्वविद्यालय, ³वाशिंगटन विश्वविद्यालय, ⁴HSE विश्वविद्यालय, ⁵अनुप्रयुक्त भौतिकी संस्थान
  • वर्गीकरण: math.AP (गणितीय अनुप्रयुक्त विश्लेषण)
  • प्रकाशन तिथि: 17 अक्टूबर 2025
  • पेपर लिंक: https://arxiv.org/abs/2510.14177

सारांश

यह पेपर एक महत्वपूर्ण व्युत्क्रम समस्या पर विचार करता है: सुनामी तरंग के तटीय अभिलेख डेटा को देखते हुए, क्या हम इसके प्रारंभिक आकार को पुनः प्राप्त कर सकते हैं? अनुसंधान अरैखिक उथले जल समीकरणों के ढांचे में किया गया है, जो सुनामी प्रसार और बाढ़ के अध्ययन में व्यापक रूप से उपयोग किए जाते हैं। लेखकों ने सीमित ढलान वाली समुद्र तल स्थलाकृति पर विचार किया है और प्रदर्शित किया है कि तटीय डेटा से आभासी बॉय पर सीमा शर्तें (जल स्तर विस्थापन और वेग) पुनः प्राप्त की जा सकती हैं। आगे चलकर अधिक जटिल टुकड़ों में ढलान वाली समुद्र तल स्थलाकृति में उथले जल समीकरणों को Boussinesq समीकरणों के साथ जोड़ने पर विचार किया गया है, ताकि प्रारंभिक शर्तें पुनः प्राप्त की जा सकें और मॉडल में फैलाव प्रभाव शामिल किए जा सकें।

अनुसंधान पृष्ठभूमि और प्रेरणा

समस्या की महत्ता

  1. आपदा खतरा: 1964 का अलास्का भूकंप, 2004 की हिंद महासागर सुनामी, 2011 का पूर्वी जापान भूकंप और अन्य प्रमुख सुनामी आपदाओं ने भारी जनहानि और आर्थिक नुकसान का कारण बना, जिससे सुनामी चेतावनी और आपदा न्यूनीकरण की महत्ता स्पष्ट हुई।
  2. चेतावनी प्रणाली की आवश्यकता: DART (गहरे समुद्र मूल्यांकन और सुनामी रिपोर्टिंग) प्रणाली के विकास के साथ, बॉय प्लेसमेंट को निर्देशित करने के लिए सैद्धांतिक ढांचे की आवश्यकता है, ताकि सुनामी चेतावनी की सटीकता और समयबद्धता में सुधार हो सके।
  3. गणितीय मॉडलिंग चुनौती: सुनामी प्रसार में जटिल अरैखिक प्रभाव शामिल होते हैं, विशेषकर तटीय क्षेत्र में, जिसके लिए तरंग व्यवहार का वर्णन करने के लिए सटीक गणितीय मॉडल की आवश्यकता होती है।

मौजूदा विधियों की सीमाएं

मौजूदा सुनामी व्युत्क्रम विधियों में मुख्य रूप से चार दृष्टिकोण शामिल हैं:

  1. भूकंप और तरंग रूप डेटा को जोड़ने वाली विधि
  2. केवल तरंग रूप डेटा का उपयोग करके पुनर्निर्माण विधि
  3. तलछट वितरण का उपयोग करके व्युत्क्रम विधि
  4. बाढ़ डेटा पर आधारित व्युत्क्रम विधि

लेकिन ये विधियां सीमित ढलान वाली समुद्र तल स्थलाकृति को संभालने में सीमाएं रखती हैं, विशेषकर अरैखिक प्रभावों और फैलाव प्रभावों के समन्वित प्रभाव पर विचार करते समय।

अनुसंधान प्रेरणा

इस पेपर की प्रेरणा लेखकों के पिछले अनंत ढलान वाली समुद्र तल स्थलाकृति के मामले में किए गए कार्य को विस्तारित करना है, अधिक यथार्थवादी सीमित ढलान वाली समुद्र तल स्थलाकृति के लिए लागू व्युत्क्रम सिद्धांत विकसित करना है, और बॉय प्लेसमेंट को निर्देशित करने के लिए एक सैद्धांतिक ढांचा प्रदान करना है।

मुख्य योगदान

  1. सैद्धांतिक विस्तार: सुनामी व्युत्क्रम सिद्धांत को अनंत ढलान वाली समुद्र तल से सीमित ढलान वाली समुद्र तल स्थलाकृति तक विस्तारित किया, पूर्ण विश्लेषणात्मक समाधान प्रदान किया।
  2. एल्गोरिदम विकास: तटीय अभिलेख डेटा से आभासी बॉय पर सीमा शर्तें पुनः प्राप्त करने के लिए एक पूर्ण एल्गोरिदम प्रवाह प्रस्तावित किया।
  3. मॉडल जोड़ना: अरैखिक उथले जल समीकरणों (NSWE) को Boussinesq समीकरणों के साथ जोड़ने की विधि विकसित की, जो टुकड़ों में ढलान वाली समुद्र तल स्थलाकृति में अरैखिक और फैलाव प्रभावों को एक साथ मानती है।
  4. संख्यात्मक सत्यापन: विभिन्न भौतिक यथार्थवादी परिदृश्यों के माध्यम से विधि की प्रभावशीलता का सत्यापन किया, जिसमें विभिन्न प्रकार की प्रारंभिक तरंग आकृतियां शामिल हैं।
  5. कम्प्यूटेशनल दक्षता: O(N³) जटिलता के साथ उच्च-दक्ष एल्गोरिदम कार्यान्वयन प्रदान किया, जो तेजी से सुनामी विशेषता गणना का समर्थन करता है।

विधि विवरण

कार्य परिभाषा

इनपुट: तटीय अभिलेख डेटा R(t) आउटपुट: आभासी बॉय पर जल स्तर विस्थापन η(L,t) और वेग u(L,t) बाधाएं: गैर-टूटती तरंग स्थिति, Carrier-Greenspan रूपांतरण की वैधता

मुख्य गणितीय ढांचा

1. अरैखिक उथले जल समीकरण प्रणाली

आयामहीन इकाइयों में, प्रणाली (1) है:

∂tη + u∂x(x + η) + (x + η)∂xu = 0  (द्रव्यमान संरक्षण)
∂tu + u∂xu + ∂xη = 0              (गति संरक्षण)

2. Carrier-Greenspan रूपांतरण

महत्वपूर्ण रैखिकीकरण रूपांतरण (5):

φ(σ,τ) = u(x,t)
σ = x + η(x,t)
ψ(σ,τ) = η(x,t) + u²(x,t)/2
τ = t - u(x,t)

यह रूपांतरण अरैखिक प्रणाली (1) को रैखिक प्रणाली (6) में परिवर्तित करता है:

∂τψ + σ∂σφ + φ = 0
∂τφ + ∂σψ = 0

3. अग्रगामी समस्या समाधान

नए चर ρ² = σ/σL और फलन θ(ρ,τ) = ψ(σ,τ) - ψb(τ) का परिचय देकर, Fourier अपघटन का उपयोग करते हुए:

θ(τ,ρ) = Σ cn(τ)J₀(jnρ)

जहां J₀ पहली तरह का बेसेल फलन है, jn J₀ का n-वां शून्य है।

4. व्युत्क्रम एल्गोरिदम मुख्य सूत्र

तटीय समीकरण (20):

ψsh(τ) = Σ (bn/√an) ∫₀τ sin(√an(τ-ξ))ψ''b(ξ)dξ + ψb(τ)

व्युत्क्रम सूत्र (24):

ψb(τ) = L⁻¹[Lψsh/(s²Σ bn/(an+s²) + 1)⁻¹]

पूर्ण व्युत्क्रम एल्गोरिदम प्रवाह

  1. डेटा रूपांतरण: CGT का उपयोग करके R(t) को ψsh(τ) में रूपांतरित करें
  2. सीमा पुनर्प्राप्ति: सूत्र (24) के माध्यम से ψsh(τ) से ψb(τ) की गणना करें
  3. क्षेत्र पुनर्निर्माण: सूत्र (15,16) का उपयोग करके पूरे डोमेन में ψ(σ,τ) और φ(σ,τ) की गणना करें
  4. व्युत्क्रम रूपांतरण: भौतिक सीमा शर्तें u(L,t) और η(L,t) पुनः प्राप्त करने के लिए व्युत्क्रम CGT निष्पादित करें

तकनीकी नवाचार बिंदु

  1. विश्लेषणात्मक समाधान विधि: संख्यात्मक विधियों के विपरीत, पूर्ण विश्लेषणात्मक व्युत्क्रम समाधान प्रदान करता है, उच्च कम्प्यूटेशनल दक्षता और नियंत्रणीय सटीकता के साथ।
  2. लाप्लास रूपांतरण तकनीक: कम्प्यूटेशनल दक्षता बढ़ाने के लिए लाप्लास रूपांतरण और व्युत्क्रम तीव्र फूरियर रूपांतरण (iFFT) का कुशलतापूर्वक उपयोग।
  3. मॉडल जोड़ने की रणनीति: NSWE को Boussinesq समीकरणों के साथ जोड़ने का नवाचारी तरीका, अरैखिक और फैलाव प्रभावों को एक साथ संभालता है।

प्रायोगिक सेटअप

संख्यात्मक सत्यापन परिदृश्य

लेखकों ने सत्यापन के लिए चार अलग-अलग प्रारंभिक विस्थापन परिदृश्य तैयार किए:

  1. द्विशिखर तरंग आकृति: η(x,0) = 0.005 sech(x-16) - 0.003e^(-(x-13)²)
  2. एकाकी तरंग आकृति: η(x,0) = 0.005 cosh⁻²(2x-6)
  3. गॉसियन तरंग आकृति: η(x,0) = 0.005e^(-(x-7)²)
  4. बहु-शिखर तरंग आकृति: η(x,0) = 0.005e^(-2(x-6)²) + 0.003e^(-(x-10)²)

मूल्यांकन संकेतक

  • सटीकता तुलना: पुनः प्राप्त ψb(τ) की तुलना सटीक मान ψᵉb(τ) से
  • सीमा शर्त पुनर्प्राप्ति: u(L,t) और η(L,t) की पुनर्प्राप्ति सटीकता का सत्यापन
  • कम्प्यूटेशनल दक्षता: एल्गोरिदम रनटाइम और जटिलता को मापना

कार्यान्वयन विवरण

  • सीमा सेटिंग: L = 1 (आयामहीन इकाई)
  • डेटा बिंदु संख्या: N = 1500 डेटा बिंदु
  • काटने का पैरामीटर: बेसेल फलन शून्य को n = 500 तक काटा गया
  • कम्प्यूटिंग प्लेटफॉर्म: Google Colab वातावरण

प्रायोगिक परिणाम

मुख्य परिणाम

1. व्युत्क्रम सटीकता सत्यापन

चित्र 4 और 5 से देखा जा सकता है कि सभी चार परीक्षण परिदृश्यों के लिए, पुनः प्राप्त सीमा शर्तें ψb(τ) और φb(τ) सटीक समाधान के साथ अत्यंत सुसंगत हैं, जो एल्गोरिदम की प्रभावशीलता को सत्यापित करता है।

2. कम्प्यूटेशनल प्रदर्शन

  • रनटाइम: 1500 डेटा बिंदुओं की गणना Google Colab पर 76 सेकंड का समय लेती है
  • जटिलता सत्यापन: समय जटिलता सैद्धांतिक विश्लेषण के O(N³) के अनुरूप है
  • फिटिंग परिणाम: समय और डेटा बिंदु संख्या के न्यूनतम वर्ग फिटिंग की ढलान 2.82 है, जो सैद्धांतिक मान 3 के करीब है

3. भौतिक सीमा शर्त पुनर्प्राप्ति

चित्र 7 चार परिदृश्यों में पुनः प्राप्त भौतिक सीमा शर्तें u(L,t) और η(L,t) दिखाता है, परिणाम दर्शाते हैं कि एल्गोरिदम आभासी बॉय पर तरंग क्षेत्र की जानकारी को सटीकता से पुनर्निर्माण कर सकता है।

Boussinesq समीकरण जोड़ने का प्रयोग

1. द्विगुण एकाकी तरंग समाधान सत्यापन

Hirota द्विगुण एकाकी तरंग समाधान को सीमा शर्त के रूप में उपयोग करते हुए, NSWE और Boussinesq समीकरणों की जोड़ने की विधि को सफलतापूर्वक सत्यापित किया। प्रयोग पैरामीटर:

  • q₁ = 0.1, q₂ = √0.1
  • विभिन्न समय पैरामीटर t₁, t₂ संयोजन

2. प्रारंभिक शर्त पुनर्निर्माण

चित्र 13 सुनामी घटना समय (tevent = -100s) पर प्रारंभिक शर्त पुनर्निर्माण परिणाम दिखाता है:

  • सटीक शिखर स्थिति: x₁ = 119.55, x₂ = 123.63
  • पुनः प्राप्त शिखर स्थिति: x₁ = 119.50, x₂ = 123.64
  • त्रुटि न्यूनतम है, जो विधि की सटीकता को सत्यापित करता है

रैखिक उथले जल समीकरण तुलना

चित्र 16 समान सीमा शर्तों के तहत Boussinesq मॉडल और रैखिक उथले जल समीकरण मॉडल की प्रारंभिक स्थिति की तुलना दिखाता है, फैलाव प्रभाव को उजागर करता है।

संबंधित कार्य

सुनामी व्युत्क्रम अनुसंधान दिशाएं

  1. भूकंप-तरंग रूप संयुक्त व्युत्क्रम: Yokota et al. (2011), Yue et al. (2014)
  2. शुद्ध तरंग रूप डेटा व्युत्क्रम: Ho et al. (2019), Yamanaka and Tanioka (2024)
  3. तलछट व्युत्क्रम: Tang and Weiss (2015), Ioki and Tanioka (2016)
  4. बाढ़ डेटा व्युत्क्रम: Lee et al. (2021), Piatanesi et al. (1996)

संख्यात्मक विधि विकास

  • परिमित अंतर/परिमित तत्व: Kounadis and Dougalis (2020), Xing (2017)
  • भौतिकी-सूचित तंत्रिका नेटवर्क: Brecht et al. (2025)
  • टुकड़ों में स्थलाकृति मॉडलिंग: Synolakis (1987), Varsoliwala and Singh (2021)

इस पेपर के सापेक्ष लाभ

मौजूदा कार्य की तुलना में, यह पेपर पूर्ण विश्लेषणात्मक व्युत्क्रम विधि प्रदान करता है, संख्यात्मक विधियों की संचयी त्रुटि से बचता है, और पहली बार सीमित ढलान वाली स्थलाकृति की व्युत्क्रम समस्या को व्यवस्थित रूप से संभालता है।

निष्कर्ष और चर्चा

मुख्य निष्कर्ष

  1. सैद्धांतिक पूर्णता: सुनामी व्युत्क्रम सिद्धांत को सीमित ढलान वाली समुद्र तल स्थलाकृति तक सफलतापूर्वक विस्तारित किया, पूर्ण विश्लेषणात्मक समाधान प्रदान किया।
  2. व्यावहारिक मूल्य: एल्गोरिदम DART प्रणाली जैसी सुनामी चेतावनी बॉय के अनुकूल प्लेसमेंट को निर्देशित करने के लिए उपयोग किया जा सकता है, चेतावनी दक्षता में सुधार करता है।
  3. मॉडल एकीकरण: NSWE और Boussinesq समीकरणों की जोड़ने को सफलतापूर्वक कार्यान्वित किया, अधिक जटिल स्थलाकृति में सुनामी मॉडलिंग के लिए नया मार्ग प्रदान किया।

सीमाएं

  1. गैर-टूटती तरंग धारणा: एल्गोरिदम केवल गैर-टूटती तरंग स्थिति में लागू होता है, जब Jacobian निर्धारक शून्य हो तो विधि विफल हो जाती है।
  2. एक-आयामी मॉडल सीमा: वर्तमान में केवल एक-आयामी स्थिति पर विचार किया गया है, वास्तविक सुनामी प्रसार त्रि-आयामी घटना है।
  3. स्थलाकृति सरलीकरण: हालांकि टुकड़ों में स्थलाकृति पर विचार किया गया है, लेकिन यह वास्तविक समुद्र तल स्थलाकृति से बहुत सरल है।
  4. फैलाव प्रभाव सीमा: Boussinesq समीकरण की प्रयोज्यता कुछ मामलों में अभी भी विवादास्पद है।

भविष्य की दिशाएं

  1. बहु-आयामी विस्तार: द्वि-आयामी और त्रि-आयामी व्युत्क्रम सिद्धांत विकसित करना
  2. जटिल स्थलाकृति: बहु-खंड ढलान वाली समुद्र तल स्थलाकृति की व्युत्क्रम समस्या को संभालना
  3. टूटती तरंग मानदंड: स्वीकार्य तटीय डेटा के लिए स्पष्ट मानदंड स्थापित करना
  4. व्यावहारिक अनुप्रयोग: वास्तविक सुनामी अवलोकन डेटा के साथ सत्यापन

गहन मूल्यांकन

शक्तियां

  1. गणितीय कठोरता: कठोर गणितीय सिद्धांत पर आधारित, पूर्ण विश्लेषणात्मक समाधान प्रदान करता है, संख्यात्मक विधियों की त्रुटि संचय से बचता है।
  2. कम्प्यूटेशनल दक्षता: O(N³) जटिलता सुनामी चेतावनी की वास्तविक समय आवश्यकता के तहत स्वीकार्य है, विशेषकर एक-आयामी मॉडल की कम्प्यूटेशनल गति लाभ स्पष्ट है।
  3. मजबूत नवाचार: पहली बार सीमित ढलान वाली स्थलाकृति की सुनामी व्युत्क्रम समस्या को व्यवस्थित रूप से हल किया, और विभिन्न समीकरणों की जोड़ने को नवाचारी रूप से कार्यान्वित किया।
  4. पर्याप्त सत्यापन: विभिन्न तरंग आकृति परिदृश्यों के संख्यात्मक सत्यापन के माध्यम से, विधि की स्थिरता और सटीकता को प्रमाणित किया।

कमियां

  1. प्रयोज्यता सीमा: गैर-टूटती तरंग धारणा मजबूत अरैखिक स्थितियों में विधि के अनुप्रयोग को सीमित करती है।
  2. स्थलाकृति मॉडल सरलीकरण: वास्तविक समुद्र तल स्थलाकृति टुकड़ों में ढलान से कहीं अधिक जटिल है, विधि के व्यावहारिक अनुप्रयोग प्रभाव को सत्यापन की आवश्यकता है।
  3. वास्तविक डेटा सत्यापन की कमी: केवल संश्लेषित डेटा का उपयोग करके सत्यापन किया गया है, वास्तविक सुनामी अवलोकन डेटा के साथ तुलना की कमी है।
  4. सीमा शर्त संगतता: हालांकि संगतता शर्तों का उल्लेख किया गया है, लेकिन व्यावहारिक अनुप्रयोग में इसके प्रभाव पर विस्तार से चर्चा नहीं की गई है।

प्रभाव

  1. सैद्धांतिक योगदान: सुनामी व्युत्क्रम सिद्धांत के लिए महत्वपूर्ण गणितीय आधार प्रदान करता है, विशेषकर विश्लेषणात्मक विधियों के क्षेत्र में।
  2. व्यावहारिक मूल्य: सुनामी चेतावनी प्रणाली के बॉय प्लेसमेंट के लिए सैद्धांतिक मार्गदर्शन प्रदान करता है, आपदा न्यूनीकरण में महत्वपूर्ण महत्व रखता है।
  3. पद्धतिगत महत्व: Carrier-Greenspan रूपांतरण का नवाचारी अनुप्रयोग अन्य अरैखिक तरंग समस्याओं के लिए संदर्भ प्रदान करता है।

प्रयोज्य परिदृश्य

  1. सुनामी चेतावनी प्रणाली डिजाइन: DART बॉय जैसे निगरानी उपकरणों के अनुकूल प्लेसमेंट को निर्देशित करना
  2. ऐतिहासिक सुनामी पुनर्निर्माण: प्राचीन सुनामी तलछट डेटा के आधार पर ऐतिहासिक घटनाओं का पुनर्निर्माण
  3. सैद्धांतिक अनुसंधान: सुनामी प्रसार तंत्र अनुसंधान के लिए गणितीय उपकरण प्रदान करना
  4. इंजीनियरिंग अनुप्रयोग: तटीय सुरक्षा इंजीनियरिंग डिजाइन के लिए संदर्भ

संदर्भ

पेपर में संबंधित साहित्य के समृद्ध संदर्भ हैं, मुख्य रूप से:

  • Carrier and Greenspan (1958): CGT रूपांतरण का मूल कार्य
  • Rybkin et al. (2021, 2023, 2024, 2025): लेखक दल का श्रृंखला कार्य
  • Synolakis (1987): टुकड़ों में स्थलाकृति मॉडलिंग का शास्त्रीय कार्य
  • Hirota (1973): Boussinesq समीकरण एकाकी तरंग समाधान का महत्वपूर्ण योगदान

समग्र मूल्यांकन: यह गणितीय भौतिकी का एक उच्च गुणवत्ता वाला पेपर है, जो सुनामी व्युत्क्रम सिद्धांत में महत्वपूर्ण योगदान देता है। विधि में कठोर गणितीय आधार और अच्छी कम्प्यूटेशनल दक्षता है, सुनामी चेतावनी और आपदा न्यूनीकरण के लिए महत्वपूर्ण सैद्धांतिक और व्यावहारिक मूल्य रखता है। हालांकि कुछ सीमाएं हैं, लेकिन यह क्षेत्र के आगे विकास के लिए एक मजबूत आधार प्रदान करता है।