انقسام الموجة القصية

انقسام الموجة القصية ، أو ما يُعرف أيضًا بالانكسار المزدوج الزلزالي ، هو ظاهرة تحدث عندما تدخل موجة قصية مستقطبة وسطًا غير متجانس . تنقسم الموجة القصية الساقطة إلى موجتين قصيتين مستقطبتين. يُستخدم انقسام الموجة القصية عادةً كأداة لاختبار تباين الخواص في منطقة معينة. تعكس هذه القياسات درجة تباين الخواص، وتؤدي إلى فهم أفضل لكثافة الشقوق واتجاهها أو محاذاة البلورات في تلك المنطقة. ^{[ 1 ]} على سبيل المثال، يمكن تشبيه تباين الخواص في منطقة معينة بصندوق أسود ، وقياسات انقسام الموجة القصية بنظرة على ما بداخل هذا الصندوق.

قد تدخل موجة قص ساقطة وسطًا غير متجانس من وسط متجانس نتيجةً لتغير في اتجاه أو خصائص الوسط. فعندما تدخل موجة قص مستقطبة وسطًا جديدًا غير متجانس، تنقسم إلى موجتين قص. تكون إحدى هاتين الموجتين أسرع من الأخرى وموازية للشقوق أو البلورات في الوسط. أما الموجة الثانية فتكون أبطأ من الأولى، وأحيانًا تكون متعامدة مع كل من الموجة الأولى والشقوق أو البلورات في الوسط. توفر الفروق الزمنية الملحوظة بين موجتي القص البطيئة والسريعة معلومات حول كثافة الشقوق في الوسط. ويشير اتجاه موجة القص السريعة إلى اتجاه الشقوق في الوسط.

عند رسمها باستخدام مخططات الاستقطاب، يمكن تحديد وصول موجات القص المنقسمة من خلال التغيرات المفاجئة في اتجاه حركة الجسيمات.

في مادة متجانسة ذات تباين ضعيف، تنقسم موجة القص الساقطة إلى موجتين شبه قصيتين باستقطابين متعامدين تقريبًا، تصلان إلى جهاز الاستقبال في نفس الوقت تقريبًا. أما في القشرة الأرضية العميقة والوشاح العلوي ، فتنقسم موجات القص عالية التردد تمامًا إلى موجتين قصيتين منفصلتين باستقطابين مختلفين ، مع وجود فارق زمني بينهما قد يصل إلى بضع ثوانٍ. ^{[ 2 ]}

أجرى هيس ^{[ 3 ]} (1964) أولى القياسات لتغيرات السرعة السمتية لموجة P في الأحواض المحيطية . وقد اختيرت هذه المنطقة لهذه الدراسة لأن الأحواض المحيطية تتكون من صخور كبيرة ومتجانسة نسبيًا. لاحظ هيس، من تجارب سابقة لقياس السرعة الزلزالية باستخدام بلورات الأوليفين ، أنه إذا كان للبلورات حتى أدنى توجيه إحصائي، فسيكون ذلك واضحًا للغاية في السرعات الزلزالية المسجلة باستخدام الانكسار الزلزالي. تم اختبار هذا المفهوم باستخدام مقاطع الانكسار الزلزالي من منطقة صدع ميندوسينو . وجد هيس أن موجات الضغط البطيئة تنتشر عموديًا على مستوى الانزلاق، وأن مركبة السرعة الأعلى موازية له. واستنتج أنه يمكن تسجيل بنية الأحواض المحيطية بسرعة وفهمها بشكل أفضل باستخدام هذه التقنيات.

ركزت دراسة أندو ^{[ 4 ]} (1980) على تحديد تباين الموجات القصية في الوشاح العلوي للأرض . وتناولت هذه الدراسة انقسام الموجات القصية المسجلة بالقرب من منطقة تشوبو البركانية في اليابان . وباستخدام محطات رصد الزلازل عن بُعد حديثة الإنشاء ، تمكن الباحثون من تسجيل وصول كل من الموجات الأولية (P) والثانوية (S) من الزلازل التي وقعت على عمق يصل إلى 260 كيلومترًا تحت المنطقة البركانية. وتجعل أعماق هذه الزلازل هذه المنطقة مثالية لدراسة بنية الوشاح العلوي. وقد لاحظ الباحثون وصول موجتين قصيتين متميزتين باستقطابين مختلفين (شمال-جنوب، سريع، وشرق-غرب، بطيء) بفارق زمني يبلغ حوالي 0.7 ثانية. وخلصت الدراسة إلى أن الانقسام لم يكن ناتجًا عن مصدر الزلزال، بل عن مسار الموجات في طريقها إلى أجهزة قياس الزلازل . واستُخدمت بيانات من محطات أخرى قريبة لتحديد مصدر التباين الزلزالي. وجد أن التباين يتوافق مع المنطقة الواقعة أسفل المنطقة البركانية مباشرةً، وافترض أن ذلك ناتج عن بلورات موجهة في حجرة صهارة عميقة الجذور . إذا احتوت حجرة الصهارة على شوائب بيضاوية الشكل موجهة تقريبًا باتجاه الشمال والجنوب، فسيكون اتجاه السرعة القصوى أيضًا باتجاه الشمال والجنوب، مما يفسر وجود الانكسار المزدوج الزلزالي . 

اقترح كرامبين ^{[ 5 ]} (1980) نظرية التنبؤ بالزلازل باستخدام قياسات انقسام الموجات القصية. تستند هذه النظرية إلى حقيقة أن الشقوق الدقيقة بين حبيبات أو بلورات الصخور تتسع أكثر من المعتاد عند مستويات الإجهاد العالية. وبعد انحسار الإجهاد، تعود هذه الشقوق إلى مواقعها الأصلية. تُعرف هذه الظاهرة، المتمثلة في اتساع الشقوق وانغلاقها استجابةً لتغيرات الإجهاد، بالتمدد . ولأن بصمات انقسام الموجات القصية تعتمد على كلٍ من اتجاه الشقوق الدقيقة (عموديًا على اتجاه الإجهاد السائد) ووفرة الشقوق، فإن هذه البصمة تتغير بمرور الوقت لتعكس تغيرات الإجهاد في المنطقة. وبمجرد تحديد بصمات منطقة ما، يمكن استخدامها للتنبؤ بالزلازل القريبة التي تحمل نفس البصمات.

أقرّ كرامبين ^{[ 6 ]} (1981) لأول مرة بظاهرة انقسام موجات القصّ المتراصفة في القشرة الأرضية . وقد استعرض النظرية السائدة، وحدّث المعادلات لفهم انقسام موجات القصّ بشكل أفضل، وقدّم بعض المفاهيم الجديدة. وأثبت كرامبين إمكانية تطوير حلول لمعظم المسائل غير المتناحية. فإذا أمكن صياغة حلٍّ مماثل لحالة متناحية، يُمكن التوصل إلى حلٍّ للحالة غير المتناحية بإجراء المزيد من الحسابات. ويُعدّ التحديد الصحيح لاستقطاب الموجات الجسمية والسطحية مفتاحًا لتحديد درجة عدم التناحية. ويمكن تبسيط نمذجة العديد من المواد ثنائية الطور باستخدام ثوابت المرونة غير المتناحية . ويمكن إيجاد هذه الثوابت من خلال مراجعة البيانات المسجلة. وقد لوحظ ذلك في العديد من المناطق حول العالم. ^{[ 7 ]}

يمكن تفسير اختلاف سرعات انتشار موجتي القص بمقارنة استقطابهما مع الاتجاه السائد للتباين في المنطقة. ويمكن استخدام التفاعلات بين الجسيمات الدقيقة المكونة للمواد الصلبة والسائلة كمثال لكيفية انتقال الموجة عبر وسط ما. تتميز المواد الصلبة بجسيمات مترابطة بإحكام شديد تنقل الطاقة بسرعة وكفاءة عاليتين. أما في السوائل، فتكون الجسيمات أقل ترابطًا، وعادةً ما يستغرق انتقال الطاقة وقتًا أطول. ويعود ذلك إلى أن الجسيمات تقطع مسافة أطول لنقل الطاقة من جسيم إلى آخر. إذا استقطبت موجة القص بالتوازي مع الشقوق في هذا الوسط المتباين الخواص، فقد تبدو كالموجة الزرقاء الداكنة في الشكل 1. تؤثر هذه الموجة على الجسيمات كما لو كانت الطاقة تنتقل عبر مادة صلبة، وستكون سرعتها عالية نظرًا لتقارب الحبيبات. إذا وُجدت موجة قصّ مستقطبة عموديًا على الشقوق المملوءة بالسوائل أو بلورات الأوليفين المتطاولة الموجودة في الوسط، فإنها ستؤثر على هذه الجسيمات كما تؤثر على الجسيمات المكونة للسائل أو الغاز . تنتقل الطاقة ببطء أكبر عبر الوسط، وتكون سرعتها أبطأ من سرعة موجة القصّ الأولى. يعتمد التأخير الزمني بين وصول موجات القصّ على عدة عوامل، منها درجة التباين والمسافة التي تقطعها الموجات إلى محطة التسجيل. تتميز الأوساط ذات الشقوق الأوسع والأكبر بتأخير زمني أطول من الأوساط ذات الشقوق الصغيرة أو حتى المغلقة. يستمر انقسام موجة القصّ حتى يصل تباين سرعة موجة القصّ إلى حوالي 5.5%. ^{[ 7 ]}

التفسير الرياضي (نظرية الأشعة) ^{[ 8 ]}

يمكن كتابة معادلة الحركة في الإحداثيات الديكارتية المستطيلة على النحو التالي:

${\displaystyle {\frac {\partial }{\partial x_{i}}}\left[c_{ijkl}{\frac {\partial U_{k}}{\partial x_{l}}}\right]=\rho {\frac {\partial ^{2}U_{j}}{\partial t^{2}}}}$

1

حيث يمثل t الزمن،${\displaystyle \rho }$الكثافة ،​${\displaystyle U_{j}}$ يمثل هذا المكون من متجه الإزاحة U ، و${\displaystyle c_{ijkl}}$يمثل موتر المرونة . ويمكن وصف جبهة الموجة بالمعادلة

${\displaystyle t=\tau \left(x_{i}\right)}$

2

يمكن التعبير عن حل المعادلة ( 1 ) على شكل متسلسلة أشعة

${\displaystyle U_{k}\left(x_{i},t\right)=\sum _{n=0}^{\infty }U_{k}^{\left(n\right)}\left(x_{i}\right)f_{n}\left(t-\tau \left(x_{i}\right)\right)}$

3

حيث الدالة${\displaystyle f_{n}\left(\vartheta \right)}$ يفي بالعلاقة

${\displaystyle df_{n+1}\left(\vartheta \right)/d\vartheta =f_{n}\left(\vartheta \right)}$

4

استبدل ( 3 ) في ( 1 )،

${\displaystyle N\left(U^{\left(n\right)}\right)-M\left(U^{\left({n-1}\right)}\right)+L\left(U^{\left({n-2}\right)}\right)=0}$

5

حيث تُعطى معاملات المتجهات N وM وL بالصيغة التالية:

${\displaystyle {\begin{cases}N_{j}\left(U^{\left(n\right)}\right)=\Gamma _{jk}U_{k}^{\left(n\right)}-U_{j}^{\left(n\right)}\\M_{j}\left(U^{\left(n\right)}\right)=p_{i}~a_{ijkl}{\frac {\partial U_{k}^{\left(n\right)}}{\partial x_{l}}}+\rho ^{-1}{\frac {\partial }{\partial x_{i}}}\left(\rho ~a_{ijkl}~p_{l}U_{k}^{\left(n\right)}\right)\\L_{j}\left(U^{\left(n\right)}\right)=\rho ^{-1}{\frac {\partial }{\partial x_{i}}}\left(\rho ~a_{ijkl}{\frac {\partial U_{k}^{\left(n\right)}}{\partial x_{l}}}\right)\end{cases}}}$

6

أين

${\displaystyle \Gamma _{jk}=p_{i}~p_{l}~a_{ijkl},\quad a_{ijkl}=c_{ijkl}/\rho ,\quad p_{i}={\frac {\partial \tau }{\partial x_{i}}}}$

7

للطلب الأول${\displaystyle n=0}$، لذا${\displaystyle U^{\left(-1\right)}=U^{\left(-2\right)}=0}$ وبالتالي ، لم يتبق سوى المكون الأول من المعادلة ( 5 ).

${\displaystyle N_{j}\left(U^{\left(0\right)}\right)=\Gamma _{jk}U_{k}^{\left(0\right)}-U_{j}^{\left(0\right)}==\left(\Gamma _{jk}-\delta _{jk}\right)U_{k}^{\left(0\right)}==0}$

8

للحصول على حل المعادلة ( 8 ) ، القيم الذاتية والمتجهات الذاتية للمصفوفة${\displaystyle \Gamma _{jk}}$هناك حاجة إلى ذلك،

${\displaystyle Det\left(\Gamma _{jk}-G\delta _{jk}\right)=0}$

9

والتي يمكن إعادة كتابتها على النحو التالي

${\displaystyle G^{3}-PG^{2}+QG-R=0}$

9

حيث القيم${\displaystyle P,Q}$و${\displaystyle R}$هي ثوابت المصفوفة المتناظرة${\displaystyle \Gamma _{jk}}$المصفوفة ​${\displaystyle \Gamma _{jk}}$له ثلاثة متجهات ذاتية:${\displaystyle g_{1},~g_{2},~g_{3}}$، والتي تتوافق مع ثلاث قيم ذاتية لـ ${\displaystyle G_{1},~G_{2},}$و${\displaystyle ~G_{3}}$.

• بالنسبة للأوساط المتجانسة،${\displaystyle G_{1}=\alpha ^{2}p_{i}p_{i}}$يتوافق مع الموجة الانضغاطية و ${\displaystyle G_{2}=G_{3}=\beta ^{2}p_{i}p_{i}}$يتوافق ذلك مع موجتي القص اللتين تتحركان معًا.
• بالنسبة للأوساط غير المتجانسة،${\displaystyle G_{2}\neq G_{3}}$يشير هذا إلى أن موجتي القص قد انقسمتا.

قياس معلمات انقسام الموجة القصية

النمذجة

^{[ 9 ]} في وسط متجانس متساوي الخواص، يمكن كتابة دالة موجة القص على النحو التالي

${\displaystyle u\left(\omega \right)=A~w\left(\omega \right)\exp \left[-i\omega T_{0}\right]\cdot {\hat {p}}}$

10

حيث A هي السعة المركبة ،${\displaystyle w\left(\omega \right)}$هي دالة الموجة (نتيجة تحويل فورييه لدالة زمن المصدر)، و${\displaystyle {\hat {p}}}$هو متجه وحدة حقيقي يشير إلى اتجاه الإزاحة ويقع في المستوى العمودي على اتجاه الانتشار . ويمكن تمثيل عملية انقسام موجة القص بتطبيق عامل الانقسام على دالة موجة القص.

${\displaystyle \Gamma =\exp \left[i\omega \delta t/2\right]{\hat {f}}{\hat {f}}+\exp \left[-i\omega \delta t/2\right]{\hat {s}}{\hat {s}}}$

11

أين${\displaystyle {\hat {f}}}$و${\displaystyle {\hat {s}}}$هي متجهات ذاتية لمصفوفة الاستقطاب ذات قيم ذاتية تتوافق مع سرعتي موجتي القص. شكل الموجة المنقسمة الناتج هو

${\displaystyle u_{s}\left(\omega \right)=A~w\left(\omega \right)\exp \left[-i\omega T_{0}\right]\Gamma \left(\phi ,\delta t\right)\cdot {\hat {p}}}$

12

أين ${\displaystyle \delta t}$هو التأخير الزمني بين موجات القص البطيئة والسريعة و${\displaystyle \phi }$هي الزاوية بين استقطاب موجة القص الساقطة${\displaystyle {\hat {p}}}$واستقطاب موجة القص السريعة${\displaystyle {\hat {f}}}$يمكن تقدير هذين المعيارين بشكل فردي من تسجيلات الزلازل متعددة المكونات.

تم استخدام قياسات انقسام الموجة القصية لاستكشاف التنبؤ بالزلازل ، ولرسم خرائط شبكات الكسور الناتجة عن التكسير عالي الضغط للخزانات .

بحسب كرامبين ^{[ 5 ]،} يمكن استخدام قياسات انقسام الموجات القصية لرصد مستويات الإجهاد في باطن الأرض. ومن المعروف أن الصخور القريبة من المناطق المعرضة للزلازل تُظهر تمددًا . وينتج انقسام الموجات القصية عن الموجات الزلزالية التي تنتقل عبر وسط يحتوي على شقوق أو بلورات موجهة. ويمكن دراسة التغيرات في قياسات انقسام الموجات القصية خلال الفترة التي تسبق وقوع الزلزال الوشيك، وذلك لفهم توقيت الزلزال وموقعه. وقد تُلاحظ هذه الظواهر على بُعد مئات الكيلومترات من مركز الزلزال.

تستخدم صناعة النفط قياسات انقسام الموجات القصية لرسم خرائط الشقوق في جميع أنحاء مكامن الهيدروكربونات . وحتى الآن، تُعد هذه الطريقة الأفضل للحصول على معلومات مباشرة عن شبكة الشقوق الموجودة في مكامن الهيدروكربونات . ^{[ 10 ]} يرتبط أفضل إنتاج في الحقل بمنطقة تحتوي على العديد من الشقوق الصغيرة المفتوحة، مما يسمح بتدفق مستمر للهيدروكربونات . تُسجل قياسات انقسام الموجات القصية وتُحلل للحصول على درجة التباين في جميع أنحاء المكمن. وعمومًا، تُعد المنطقة ذات أعلى درجة من التباين أفضل مكان للحفر لأنها تحتوي على أكبر عدد من الشقوق المفتوحة. ^{[ 11 ]}

تُوفّر قياسات انقسام الموجات القصية معلومات دقيقة وشاملة حول منطقة معينة. مع ذلك، توجد قيود يجب مراعاتها عند تسجيل أو تحليل هذه القياسات. تشمل هذه القيود حساسية الموجات القصية، وتغير انقسامها بتغير زاوية السقوط والسمت، واحتمالية انقسامها عدة مرات في وسط غير متجانس، وربما في كل مرة يتغير فيها اتجاهها. ^{[ 15 ]}

يُعدّ انقسام الموجات القصية شديد الحساسية للتغيرات الطفيفة في ضغط المسام في قشرة الأرض. وللكشف بنجاح عن درجة التباين في منطقة ما، يجب أن يكون هناك عدد كافٍ من الموجات الواصلة موزعة زمنيًا بشكل جيد. فقلة عدد الموجات لا تكفي لرصد التغير حتى لو كانت من أشكال موجية متشابهة. ^{[ 7 ]} يتغير انقسام الموجات القصية بتغير كل من زاوية السقوط واتجاه الانتشار. وما لم تُعرض هذه البيانات في إسقاط قطبي، فإن طبيعتها ثلاثية الأبعاد لا تنعكس وقد تكون مضللة. ^{[ 7 ]} قد ينتج انقسام الموجات القصية عن أكثر من طبقة واحدة متباينة الخواص تقع في أي مكان بين المصدر ومحطة الاستقبال. تتميز قياسات انقسام الموجات القصية بدقة أفقية عالية، ولكن بدقة رأسية ضعيفة جدًا. ^{[ 16 ]} تختلف استقطابات الموجات القصية في جميع أنحاء الكتلة الصخرية. لذلك، قد تكون الاستقطابات المرصودة خاصة بالبنية القريبة من السطح، وليست بالضرورة ممثلة للبنية محل الاهتمام. ^{[ 17 ]}

نظراً لطبيعة موجات القص المنقسمة، فإنها عند تسجيلها في مخططات الزلازل ثلاثية المكونات النموذجية ، تُنتج إشارات معقدة للغاية. إذ تتشتت الاستقطابات والتأخيرات الزمنية بشكل كبير وتختلف اختلافاً واسعاً زمانياً ومكانياً. وبسبب هذا التباين في الإشارة، يسهل إساءة تفسير وصول موجات القص واستقطابها. ^{[ 18 ]} فيما يلي شرح لبعض المفاهيم الخاطئة الشائعة المرتبطة بموجات القص، ويمكن الاطلاع على مزيد من المعلومات في كرامبين وبيكوك (2008). ^{[ 7 ]}

• تكون استقطابات موجات القص المنقسمة متعامدة. ^{[ 7 ] [ 18 ]}

تتميز موجات القص التي تنتشر على طول مسار الشعاع بسرعة جماعية باستقطابات متعامدة فقط في اتجاهات محددة. أما استقطابات الموجات الجسمية فهي متعامدة في جميع اتجاهات سرعة الطور ، إلا أن هذا النوع من الانتشار يصعب رصده أو تسجيله عمومًا.

• تكون استقطابات موجات القص المنقسمة ثابتة، موازية للشقوق، أو عمودية على مراكز الانتشار. ^{[ 7 ] [ 18 ]}

حتى عند انتشارها عبر الشقوق المتوازية أو العمودية على مراكز الانتشار أو الموازية للشقوق، فإن استقطاب موجات القص سيتغير دائمًا في ثلاثة أبعاد مع زاوية السقوط والسمت داخل نافذة موجة القص.

• يقل تباين الخواص في الشقوق دائمًا مع العمق حيث يتم إغلاق الشقوق المملوءة بالسوائل بفعل الضغط الليثوستاتيكي . ^{[ 7 ] [ 18 ]}

لا يصح هذا القول إلا إذا أُزيل السائل الموجود في الشقوق بطريقة ما. قد يتم ذلك عن طريق الامتصاص الكيميائي، أو التصريف، أو التدفق إلى السطح. مع ذلك، فإن هذه الحالات نادرة نسبيًا، وهناك أدلة تدعم وجود سوائل في الأعماق. تشمل هذه الأدلة بيانات من بئر كولا العميق، ووجود موصلية عالية في القشرة الأرضية السفلى.

• يمكن تحسين نسب الإشارة إلى الضوضاء لانقسام الموجات القصية فوق الزلازل الصغيرة عن طريق التجميع. ^{[ 7 ] [ 18 ]}

يُعدّ تجميع البيانات الزلزالية من مسح انعكاسي مفيدًا لأنها جُمعت باستخدام مصدر يمكن التنبؤ به والتحكم فيه. أما عندما يكون المصدر غير مُتحكم فيه وغير قابل للتنبؤ، فإن تجميع البيانات يؤدي فقط إلى تدهور الإشارة. ولأن تأخيرات زمن موجات القص المسجلة واستقطابها تختلف في زاوية سقوطها وسمت انتشارها الراديوي ، فإن تجميع هذه الموجات الواصلة سيؤدي إلى تدهور الإشارة وانخفاض نسبة الإشارة إلى الضوضاء، مما ينتج عنه رسم بياني مشوش ويصعب تفسيره في أحسن الأحوال. ^{[ 7 ]}

يتطور فهمنا لانقسام الموجات القصية وكيفية الاستفادة المثلى من القياسات باستمرار. ومع تحسن معرفتنا في هذا المجال، ستظهر حتماً طرق أفضل لتسجيل هذه القياسات وتفسيرها، وستتاح فرص أكبر لاستخدام البيانات. ويجري حالياً تطوير هذه التقنية لاستخدامها في صناعة النفط والتنبؤ بالزلازل والانفجارات البركانية .

استُخدمت قياسات انقسام الموجات القصية بنجاح للتنبؤ بالعديد من الزلازل. وبفضل المعدات المتطورة ومحطات التسجيل المتباعدة بكثافة أكبر، تمكّنا من دراسة التغيرات المميزة لانقسام الموجات القصية أثناء الزلازل في مناطق مختلفة. تتغير هذه السمات بمرور الوقت لتعكس مقدار الإجهاد الموجود في المنطقة. بعد تسجيل ودراسة العديد من الزلازل، تصبح سمات انقسام الموجات القصية قبيل وقوع الزلزال معروفة جيدًا، ويمكن استخدامها للتنبؤ بالزلازل المستقبلية. يمكن ملاحظة هذه الظاهرة نفسها قبل ثوران بركاني، ويُستنتج أنه يمكن التنبؤ بها بالطريقة نفسها.

تستخدم صناعة النفط منذ سنوات قياسات انقسام الموجات القصية المسجلة فوق مكامن الهيدروكربونات للحصول على معلومات قيّمة حول هذه المكامن. ويجري تحديث المعدات باستمرار للكشف عن صور جديدة ومعلومات إضافية. ^{[ 7 ]}

• معهد ألفريد فيجنر للأبحاث القطبية والبحرية (AWI) (ألمانيا)
• انقسام الموجات القصية في برنامج ماتلاب (فرنسا)
• الكثير من الصور الزلزالية المثيرة للاهتمام (جامعة ولاية أريزونا)
• معلومات عن المواد الصلبة والسائلة والغازية
• مجموعة انقسام الموجات القصية في الوسائط غير المتجانسة (جامعة غوته فرانكفورت، ألمانيا)