تحليل هيلبرت الترددي واستخلاص الغلاف: كشف الشروخ المجهرية قبل انطلاق إنذارات الاهتزاز التقليدية

ينبع — في مجمع التكرير والبتروكيماويات المطل على ساحل البحر الأحمر، تواصل مضخة تبريد وضخ فحم الكوك عالي الضغط—التي يديرها محرك كهربائي ضخم بقدرة 2,500 كيلوواط—عملها ضمن دورة تشغيلية مستمرة وشاقة. على الغلاف الحاضن لكراسي التحميل من جهة المحرك، يبث مقياس تسارع كهرضغطي صناعي تدفقاً مستمراً من إشارات التسارع اللحظية إلى أجهزة القياس والتشخيص. وعندما يفحص فني الموثوقية طيف سرعة الاهتزاز التقليدي المحول بتحويل فورييه السريع (FFT) على شاشة غرفة التحكم، تبدو الآلة في حالة ميكانيكية مثالية لا تشوبها شائبة؛ فالمعدل الكلي للاهتزاز لا يتجاوز 1.4 مليمتر في الثانية، وتستقر القمة الطيفية الرئيسة بأمان عند تردد الدوران البالغ 29.8 هرتز (1X)، معبرة عن الاتزان الميكانيكي الطبيعي. وبمعايير كود ISO 10816 التقليدية، تعمل المضخة في النطاق التشغيلي الآمن والممتاز (المنطقة A).

ولكن بعد ثمانٍ وأربعين ساعة فقط، تعرض كرسي التحميل في المضخة لانهيار كارثي مفاجئ في مجراه الداخلي؛ حيث ولّد التآكل الميكانيكي احتكاكاً حرارياً هائلاً أدى إلى لحام بكرات التحميل بالمسار الدوار وانشطار فاصل الربط الميكانيكي وتوقف وحدة المعالجة بالكامل بصورة طارئة. وكشف الفحص الجنائي بعد تفكيك الآلة عن تعرض المعدن لإجهاد تلامسي تراكمي عميق؛ إذ نشأت شروخ مجهرية تحت السطح منذ شهور طويلة، وظلت تتمدد حتى اخترقت السطح وتسببت في تشظي رقائق معدنية صلبة داخل حوض زيت التزييت. وهنا تردد السؤال الهندسي الحاسم في أرجاء المصنع: كيف يمكن لخلل ميكانيكي مدمر أن ينمو ويتفاقم دون أن يترك أدنى أثر على أطياف الاهتزاز التقليدية حتى لحظة الانهيار؟

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

ونظراً لضآلة زمن هذه النبضة التصادمية، فإن طاقتها الحركية الفيزيائية تتشتت عبر نطاق ترددي هائل يمتد من عشرة كيلوهرتز إلى ما بعد مائة كيلوهرتز. وعندما يحسب مهندسو الصيانة طيف فورييه التقليدي لسرعة الاهتزاز منخفضة التردد (والذي يقتطع عادة عند 1,000 أو 2,000 هرتز لرصد ترددات الدوران الأساسية)، تتلاشى طاقة هذه النبضات الميكروسكوبية لتصل إلى سعة شبه معدومة، وتغرق كلياً تحت القمم الطيفية العالية للاهتزازات الهيكلية الناتجة عن عدم الاتزان الميكانيكي واختلال المحاذاة. إن محاولة رصد شرخ مجهري تحت السطح في كرسي تحميل عبر طيف الاهتزاز التقليدي تعادل رياضياً محاولة الاستماع لدقات ساعة يد صغيرة داخل عنبر غلايات هادر بالضجيج.

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

تسمى هذه الظاهرة الفيزيائية بـ "تعديل السعة" (Amplitude Modulation)؛ حيث يعمل الرنين الهيكلي عالي التردد كموجة حاملة تناظرية، تتعدل سعتها دورياً في كل مرة تصطدم فيها البكرة بالعيب الميكانيكي. ولعزل التردد الحقيقي للعطل، تنفذ الخوارزميات الرقمية عملية استخلاص الغلاف عبر مسار رياضي ثلاثي المراحل.

أولاً، تُمرر إشارة الاهتزاز الخام واسعة النطاق عبر مرشح تمرير نطاقي رقمي عالي الرتبة (Bandpass Filter) يُضبط بدقة حول منطقة الرنين الهيكلي (مثل نطاق 5 إلى 15 كيلوهرتز)، مما يلغي كلياً اهتزازات عدم الاتزان واختلال المحاذاة والاضطرابات الهيدروليكية منخفضة التردد.

ثانياً، تشتق الخوارزمية الإشارة التحليلية z(t) بتطبيق تحويل هيلبرت الرياضي على الإشارة المصفاة، حيث يُحسب الغلاف الزمني الرياضي e(t)—الذي يمثل السعة اللحظية لتذبذبات الرنين الرنانة—باعتباره القيمة المطلقة للمتجه المركب للإشارة التحليلية.

ثالثاً، يُطبق تحويل فورييه السريع (FFT) عالي الدقة ليس على إشارة الاهتزاز الخام الأصلية، بل مباشرة على هذا الغلاف الزمني المستخلص e(t).

يكشف طيف الغلاف الناتج عن تحول تشخيصي مذهل؛ إذ تتلاشى ضوضاء الرنين عالي التردد كلياً، ليظهر مشهد طيفي نقي ومميز بنسبة إشارة إلى ضوضاء بالغة الارتفاع، تسيطر عليه بدقة الترددات الحركية المميزة لهندسة كرسي التحميل. يستطيع مهندسو الموثوقية فورياً تحديد تردد مرور البكرات على المجرى الخارجي (BPFO)، أو المجرى الداخلي (BPFI)، أو تردد دوران البكرات حول نفسها (BSF)، وتردد عيوب القفص الحاضن (FTF).

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

تحل الأسس الرياضية لتحليل الغلاف الترددي (Envelope Demodulation) أحد أعقد تحديات معالجة الإشارات في تشخيص المعدات الدوارة؛ والمتمثل في استخلاص نبضات الصدمات الميكانيكية المجهرية منخفضة الطاقة والمدفونة تحت ضوضاء الاهتزاز الهيكلي ودوران الآلات الصاخب. فعندما تصطدم إحدى البكرات المتدحرجة بعيب مجهري داخل المجرى المعدني، تولد الصدمة طاقة ميكانيكية لحظية تثير الترددات الرنينية الطبيعية العالية لكرسي التحميل، والتي تقع عادة في النطاق الترددي بين خمسة وعشرين كيلوهرتز. غير أن وتيرة تكرار هذه الصدمات—المحكومة بالسرعة الدورانية والأبعاد الهندسية للمحمل—تحدث عند ترددات منخفضة تقع بين خمسين وأربعمائة هرتز؛ حيث تؤدي اهتزازات الدوران الرئيسة واختلال المحاذاة وتشابك التروس إلى حجب بصمة العطل تماماً في التحليل الطيفي التقليدي السريع (FFT).

تعزل تقنية رنين الترددات العالية (HFRT) هذه النبضات المجهرية عبر خط معالجة رياضي صارم يتألف من ثلاث مراحل متتابعة؛ أولاً: تمر الإشارة الزمنية الخام عبر مرشح تمرير نطاقي رقمي أو تناظري مضبوط بدقة على نطاق الرنين الهيكلي للغلاف المعدني، مما يخمد كافة الظواهر الدورانية المنخفضة التردد. ثانياً: تخضع الإشارة المصفاة $x(t)$ لتحويل هيلبرت الرياضي لإنشاء الإشارة التحليلية المركبة: $z(t) = x(t) + imathcal{H}{x(t)} = A(t) e^{iphi(t)}$؛ حيث يُعرف تحويل هيلبرت بتكامل كوشي للقيمة الأساسية: $mathcal{H}{x(t)} = rac{1}{pi} ext{P.V.} int_{-infty}^{infty} rac{x( au)}{t - au} , d au$. ويستخلص الغلاف السعوي اللحظي $A(t) = sqrt{x^2(t) + [mathcal{H}{x(t)}]^2}$ منحنى التعديل السعوي اللحظي متخلصاً كلياً من الموجة الحاملة عالية التردد. وأخيراً، يؤدي تطبيق تحويل فورييه السريع على الغلاف المستخلص $A(t)$ إلى إظهار ترددات الصدمات الحركية بدقة ونقاء طيفي استثنائيين.

يتطابق طيف الغلاف المستخلص مباشرة مع المعادلات الحركية الحاكمة لترددات عيوب محامل الدحرجة: تردد مرور البكرات على المجرى الخارجي ($BPFO = rac{n}{2} f_r [1 - rac{d}{D}cosalpha]$)، وتردد مرور البكرات على المجرى الداخلي ($BPFI = rac{n}{2} f_r [1 + rac{d}{D}cosalpha]$)، وتردد دوران البكرة حول محورها ($BSF = rac{D}{2d} f_r [1 - (rac{d}{D}cosalpha)^2]$)، وتردد دوران القفص ($FTF = rac{1}{2} f_r [1 - rac{d}{D}cosalpha]$)؛ حيث يمثل $n$ عدد البكرات، و$f_r$ سرعة دوران العمود، و$d$ قطر البكرة، و$D$ قطر دائرة الخطوة، و$alpha$ زاوية التلامس. ومن خلال تتبع قمم التردد وحزم التردد الجانبية المتزامنة مع سرعة الدوران، يحدد مهندسو الموثوقية العيب الميكانيكي بدقة متناهية ومعدل تمدد الشروخ المجهرية.

تأكدت هذه الدقة التحليلية في وحدة تقطير النفط الخام بمصفاة ينبع؛ حيث خضعت مضخة قاع برج التقطير الحيوية التي يديرها محرك كهربائي بقدرة 400 كيلووات لمراقبة مستمرة. وظلت سرعة الاهتزاز الإجمالية مستقرة بأمان عند 2.1 مليمتر/ثانية وفق معيار ISO 10816-3، في حين كشف تحليل هيلبرت للغلاف الترددي عن قمة طيفية صاعدة بدقة عند 142.8 هرتز، متطابقة مع تردد العيب الخارجي BPFO لمحمل التلامس الزاوي SKF 7314، مصحوبة بحزم جانبية تدل على بدء التآكل بمعدل 1.2 ميكرون يومياً. وسمح هذا التشخيص الاستباقي بجدولة تغيير المحمل قبل اثنين وأربعين يوماً من تفتت المعدن، متفادياً انهيار كرسي التحميل وتمزق مانع التسرب واشتعال حريق بترولي قدرت خسائره بـ 1.6 مليون دولار.

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