محاسبات سطح سنج راداری FMCW در مهندسی پیشرفته

محاسبات سطح سنج راداری FMCW بخشی از طراحی مهندسی سیستم‌های اندازه‌گیری سطح است که به تحلیل فیزیک انتشار امواج الکترومغناطیسی، پردازش سیگنال بازتابی و تبدیل Distance به Level می‌پردازد.

در یک Radar Level Transmitter، مقدار اندازه‌گیری‌شده تنها به زمان رفت‌وبرگشت موج وابسته نیست. هندسه مخزن، نوع Antenna، وضعیت سطح محصول، شرایط گاز بالای مخزن و عملکرد زنجیره الکترونیکی نیز می‌توانند بر نتیجه نهایی اثر بگذارند.

در سطح Fundamentals، اصول عملکرد فناوری توضیح داده می‌شود و در سطح Professional Engineering، انتخاب و نصب تجهیز مورد بررسی قرار می‌گیرد.

اما در Advanced Engineering، هدف تحلیل کمّی منابع خطا، بررسی فرضیات مدل، تشکیل Uncertainty Budget و تدوین معیارهای Design Verification است.

این کتابچه از مجموعه NICO Technical Learning Series، روابط مهندسی Radar Level Measurement را از فیزیک FMCW تا FMEA و راستی‌آزمایی طراحی بررسی می‌کند.

تمام نتایج عددی حاصل از مدل‌های آموزشی هستند. مقادیری مانند فرکانس، Beam Angle، Signal Margin، Range Resolution و Expanded Uncertainty به هیچ مدل تجاری NICO نسبت داده نشده‌اند.

زنجیره اندازه‌گیری در سطح سنج راداری FMCW

زنجیره اندازه‌گیری FMCW شامل چند مرحله متوالی است.

ابتدا Chirp Generator یک سیگنال با فرکانس متغیر تولید می‌کند.

سیگنال از طریق Antenna به سمت سطح محصول منتشر می‌شود و بخشی از انرژی آن پس از برخورد با سطح بازتاب می‌شود.

موج بازگشتی در Receiver با سیگنال مرجع ترکیب شده و Beat Signal ایجاد می‌شود.

پس از نمونه‌برداری توسط ADC، پردازش FFT برای استخراج مؤلفه‌های فرکانسی انجام می‌شود.

از فرکانس‌های استخراج‌شده، موقعیت Echo و Distance محاسبه می‌شود.

سپس با استفاده از Reference Point، Empty Distance و Full Distance مقدار Level و خروجی مربوطه تعیین می‌شود.

زنجیره عملکرد به‌صورت زیر است:

Chirp Generator → Antenna → Product Echo → Mixer → ADC → FFT → Distance → Level → 4–20 mA / HART

این ترتیب نشان می‌دهد که خطای نهایی می‌تواند از انتشار موج، بازتاب سطح، انتخاب Echo، پارامترهای هندسی یا تبدیل خروجی ایجاد شود.

معادلات FMCW و محاسبه Beat Frequency

در فناوری Frequency-Modulated Continuous Wave، فرکانس موج ارسالی در مدت زمان مشخصی به‌صورت کنترل‌شده تغییر می‌کند.

اگر پهنای باند Sweep برابر B و مدت آن T باشد، Sweep Slope از رابطه زیر به دست می‌آید:

S = B / T

موج بازتابی با تأخیر زمانی نسبت به موج ارسالی دریافت می‌شود.

برای فاصله R و با فرض سرعت انتشار c، تأخیر رفت‌وبرگشت برابر است با:

τ = 2R / c

در یک Sweep خطی، Beat Frequency از ضرب Sweep Slope در تأخیر به دست می‌آید:

fb = S × τ

در نتیجه:

fb = 2RS / c

و فاصله اندازه‌گیری‌شده برابر خواهد بود با:

R = c × fb / (2S)

در مثال آموزشی کتابچه، پهنای باند 4 GHz و زمان Sweep برابر 1 ms فرض شده‌اند.

بنابراین Sweep Slope برابر 4 × 10¹² Hz/s است.

برای فاصله فرضی 10 متر، تأخیر رفت‌وبرگشت تقریباً 0.0667 میکروثانیه و Beat Frequency حدود 266.9 kHz محاسبه شده است.

این مقادیر صرفاً برای نمایش نحوه تبدیل فرکانس به فاصله هستند.

پهنای باند، مدت Sweep و الگوریتم پردازش واقعی باید از اطلاعات رسمی مدل موردنظر استخراج شوند.

محاسبه Range Resolution در سطح سنج راداری FMCW

یکی از مفاهیم مهم در پردازش سیگنال FMCW، Range Resolution یا توانایی تفکیک دو Echo نزدیک به یکدیگر است.

در مدل ساده ارائه‌شده در کتابچه، Range Resolution با پهنای باند Sweep رابطه معکوس دارد:

ΔR = c / (2B)

در این رابطه، ΔR وضوح فاصله و B پهنای باند فرکانسی است.

برای پهنای باند فرضی 4 GHz، مقدار وضوح نظری حدود 37.5 میلی‌متر محاسبه می‌شود.

در پردازش واقعی، استفاده از Window Function می‌تواند باعث افزایش پهنای Main Lobe شود.

با ضریب آموزشی 1.5، وضوح عملی مدل به حدود 56.2 میلی‌متر می‌رسد.

نکته مهم، تفاوت میان Resolution و Measurement Accuracy است.

Resolution به توانایی تفکیک دو بازتاب نزدیک اشاره دارد، درحالی‌که Accuracy به نزدیکی نتیجه اندازه‌گیری با مقدار مرجع مربوط می‌شود.

بنابراین، مقدار 37.5 یا 56.2 میلی‌متر در این مثال، دقت تضمین‌شده هیچ سطح‌سنج راداری نیست.

نقش Wavelength و فرکانس در محاسبات سطح سنج راداری

طول موج یا Wavelength از نسبت سرعت انتشار به فرکانس به دست می‌آید.

λ = c / f

در کتابچه، فرکانس‌های نمونه 6، 26 و 80 GHz برای مقایسه رفتار فیزیکی موج بررسی شده‌اند.

طول موج متناظر در خلأ به‌ترتیب حدود 50.0، 11.5 و 3.75 میلی‌متر است.

با کاهش طول موج و ثابت ماندن ابعاد مؤثر Antenna، امکان ایجاد پرتو متمرکزتر فراهم می‌شود.

در مقابل، اندازه الکتریکی لایه‌های Condensation و ذرات نسبت به طول موج می‌تواند افزایش پیدا کند.

ازاین‌رو، انتخاب فرکانس بالاتر لزوماً به معنای عملکرد بهتر در همه کاربردها نیست.

این تصمیم باید با توجه به Beam Geometry، وضعیت Foam، Condensation، Nozzle و شرایط واقعی فرآیند انجام شود.

تحلیل Antenna و محاسبه Beam Angle

Antenna نقش اصلی در انتشار و دریافت موج الکترومغناطیسی را بر عهده دارد.

الگوی تابش Antenna معمولاً دارای Main Lobe و Side Lobes است.

در بسیاری از مدل‌های تحلیلی، Beam Angle بر اساس پهنای لوب اصلی در سطح Half-Power یا نقطه −3 dB تعریف می‌شود.

بااین‌حال، تعریف دقیق زاویه پرتو به مدارک سازنده وابسته است.

در کتابچه، یک تقریب برای Antenna با دهانه دایره‌ای معرفی شده است:

θ ≈ 70 × λ / D

که در آن θ زاویه بر حسب درجه، λ طول موج و D قطر مؤثر دهانه Antenna است.

برای قطر فرضی 75 میلی‌متر، زاویه تقریبی در فرکانس‌های 6، 26 و 80 GHz به‌ترتیب 46.6، 10.8 و 3.5 درجه به دست آمده است.

این رابطه یک مدل تقریبی است و نباید به‌جای Radiation Pattern واقعی سازنده استفاده شود.

در مدل ساده کتابچه، نصف شدن Beam Angle می‌تواند باعث بهبود حدود 12 dB در بودجه Echo شود.

اما این مزیت تنها زمانی قابل استفاده است که جهت پرتو، هندسه سطح و شرایط نصب با Beam باریک سازگار باشند.

محاسبه Beam Footprint و فاصله ایمن از دیواره

Beam Footprint ناحیه‌ای است که پرتو Radar در فاصله مشخص روی سطح محصول پوشش می‌دهد.

در تقریب مخروطی ساده، قطر این ناحیه از رابطه زیر محاسبه می‌شود:

Dfp = 2 × d × tan(α / 2)

در این رابطه d فاصله Antenna تا سطح و α زاویه کامل پرتو است.

در مثال آموزشی کتابچه، برای Beam Angle برابر 8 درجه و فاصله 10 متر، قطر Beam Footprint حدود 1.40 متر محاسبه شده است.

با همین فاصله و زاویه 4 درجه، قطر تقریبی 0.70 متر خواهد بود.

اگر شعاع Beam Footprint برابر 0.70 متر و حاشیه مهندسی فرضی 0.30 متر باشد، حداقل فاصله پیشنهادی در این مثال تا دیواره یا مانع برابر 1.00 متر خواهد بود.

این مقدار یک Engineering Recommendation مبتنی بر داده‌های فرضی است.

در طراحی واقعی، Beam Angle، Beam Clearance و فاصله مجاز از دیواره باید بر اساس Datasheet و Installation Manual مدل انتخابی تأیید شوند.

تحلیل Dielectric Reflection در سطح سنج راداری

قدرت بازتاب موج از سطح محصول به Dielectric Constant و ویژگی‌های الکترومغناطیسی ماده وابسته است.

در یک مدل ساده با برخورد عمود موج به سطح تخت و محیط بدون اتلاف، کسر توان بازتابیده از رابطه زیر محاسبه می‌شود:

Γ² = [(1 − √εr) / (1 + √εr)]²

در این رابطه εr ثابت دی‌الکتریک نسبی محصول است.

در جدول آموزشی کتابچه، برای εr برابر 2 مقدار بازتاب حدود 2.94 درصد محاسبه شده است.

در مقابل، برای مقدار مرجع فرضی εr برابر 80، کسر بازتاب حدود 63.82 درصد است.

اختلاف توان بازتابی این دو حالت در مدل آموزشی حدود 13.4 dB خواهد بود.

این تحلیل نشان می‌دهد که کاهش Dielectric Constant می‌تواند باعث تضعیف Product Echo شود.

بااین‌حال، مقدار واقعی Dielectric Constant در فرکانس کاری Radar به دما، ترکیب و ویژگی‌های ماده وابسته است.

مقدار 80 برای آب در کتابچه یک مرجع ساده‌شده است و نباید بدون بررسی در محاسبات واقعی فرکانس‌های راداری استفاده شود.

محاسبه Echo Budget و Signal Margin

Signal Margin فاصله سطح Product Echo تا آستانه تشخیص سیگنال است.

برای بررسی این پارامتر، کتابچه از یک مدل نسبی بر اساس Reference Signal Margin، ضریب بازتاب، فاصله و تلفات استفاده می‌کند.

رابطه ساده‌شده به‌صورت زیر است:

M = M0 + 10log10(Γ²/Γw²) − 20log10(d/dref) − Af

در این رابطه M حاشیه نهایی، M0 حاشیه مرجع، Γ² نسبت توان بازتابی محصول، d فاصله اندازه‌گیری و Af تلفات اضافی است.

در مثال فرضی مخزن هیدروکربن با εr برابر 2، حاشیه مرجع 50 dB در فاصله 5 متر در نظر گرفته شده است.

برای فاصله 14.65 متر و افت Condensation برابر 3 dB، حاشیه نهایی در مدل 8 درجه حدود 24.3 dB محاسبه شده است.

در مدل فرضی با زاویه 4 درجه، نتیجه به حدود 36.3 dB می‌رسد.

این مقادیر نتایج یک Echo Budget آموزشی هستند، نه Signal Margin اندازه‌گیری‌شده یا تضمین‌شده یک دستگاه واقعی.

در طراحی نهایی، Gain، Noise Floor، Detection Threshold و الزامات عملکرد باید بر اساس اطلاعات سازنده و آزمون واقعی بررسی شوند.

تحلیل Tank Bottom Echo در مایعات کم‌بازتاب

در سطح پایین مخزن، Product Echo ممکن است با Tank Bottom Echo تداخل پیدا کند.

هنگامی که موج بخشی از مسیر خود را در مایع طی می‌کند، سرعت انتشار به ویژگی دی‌الکتریک مایع وابسته خواهد بود.

در مدل ساده کتابچه، فاصله ظاهری Echo کف از رابطه زیر محاسبه می‌شود:

db = (E − h) + √εr × h

که در آن E فاصله مرجع تا کف، h ارتفاع مایع و db فاصله ظاهری بازتاب کف است.

در مثال فرضی با E برابر 15.15 متر، ارتفاع سطح 0.5 متر و εr برابر 2، فاصله ظاهری Echo کف حدود 15.36 متر محاسبه شده است.

در این شرایط، فاصله Product Echo تا Echo کف حدود 0.71 متر است.

در مدل ساده و بدون اتلاف کتابچه، بازتاب کف فلزی می‌تواند حدود 15.1 dB از بازتاب سطح قوی‌تر باشد.

این وضعیت نشان می‌دهد که قوی‌ترین Echo الزاماً مربوط به سطح واقعی محصول نیست.

بنابراین، استفاده از Echo کف به‌عنوان مبنای تعیین Empty Distance بدون درنظرگرفتن شرایط انتشار و هندسه مخزن صحیح نیست.

Multiple Reflection و False Echo در اندازه‌گیری راداری

Multiple Reflection زمانی ایجاد می‌شود که موج الکترومغناطیسی بیش از یک بار میان سطح محصول و اجسام بازتابنده حرکت کند.

در یک مدل ساده از بازتاب‌های تکراری میان سطح و سقف، Echoهای اضافی ممکن است در مضارب فاصله اصلی ظاهر شوند.

dk = k × R

در این رابطه k مرتبه بازتاب و R فاصله واقعی محصول است.

برای مثال، اگر Product Echo در فاصله 7 متر باشد، ممکن است بازتاب مرتبه دوم در فاصله ظاهری 14 متر مشاهده شود.

چنین بازتابی در صورت انتخاب اشتباه توسط الگوریتم دستگاه می‌تواند یک Level ثابت ولی نادرست تولید کند.

این نوع خطا اهمیت زیادی دارد؛ زیرا ممکن است سیستم کنترل مقدار پایداری نمایش دهد، درحالی‌که Echo انتخاب‌شده واقعی نیست.

در این شرایط، بررسی Echo Curve در چند سطح و مقایسه با مرجع مستقل اهمیت دارد.

False Echo Mapping می‌تواند برای مدیریت بازتاب‌های ثابت استفاده شود، اما نمی‌تواند تمام مشکلات هندسی یا هم‌پوشانی Echoها را برطرف کند.

تحلیل الکترومغناطیسی Nozzle و Ring Echo

Nozzle تنها یک اتصال مکانیکی نیست و می‌تواند بر رفتار انتشار موج نیز اثر بگذارد.

در تحلیل پیشرفته، یک Nozzle فلزی ممکن است مانند Waveguide عمل کند.

در کتابچه، فرکانس قطع مود پایه TE11 برای یک لوله دایره‌ای از رابطه زیر بررسی شده است:

fc = 1.8412 × c / (π × Dn)

در این رابطه Dn قطر داخلی Nozzle است.

برای قطر فرضی 100 میلی‌متر، مقدار فرکانس قطع حدود 1.757 GHz محاسبه شده است.

این محاسبه به‌تنهایی به معنای تأیید سازگاری Nozzle با یک مدل راداری نیست.

در Nozzle بلند یا در شرایطی که Antenna داخل نازل قرار گرفته است، بازتاب‌های تکراری می‌توانند Ring Echo ایجاد کنند.

در مدل هندسی ساده کتابچه:

dring,k = k × 2 × Ln

برای طول مؤثر فرضی 0.40 متر، Echoهای حلقه‌ای در فاصله‌های 0.8، 1.6 و 2.4 متر ظاهر می‌شوند.

اگر Full Distance برابر 1.15 متر باشد، این بازتاب‌ها باید از نظر تداخل با سطح پر بررسی شوند.

بنابراین، تحلیل Nozzle نباید فقط به شرایط مخزن تقریباً خالی محدود شود.

تأثیر کجی Nozzle و محور Antenna بر Echo

اگر محور Antenna نسبت به عمود بر سطح کج باشد، بخش مهمی از انرژی بازتابی ممکن است به سمت گیرنده بازنگردد.

در مدل هندسی کتابچه، انحراف جانبی تقریبی بازتاب از رابطه زیر محاسبه می‌شود:

s ≈ 2 × d × sin(β)

در این رابطه β زاویه انحراف محور است.

برای فاصله 10 متر، کجی یک درجه در مدل ساده حدود 0.35 متر انحراف جانبی ایجاد می‌کند.

کتابچه همچنین از یک مدل گاوسی برای بررسی Pointing Loss استفاده می‌کند:

Ap ≈ 24 × (β/α)²

برای کجی 3 درجه و Beam Angle برابر 8 درجه، افت تقریبی 3.4 dB محاسبه شده است.

با همان کجی و Beam Angle برابر 4 درجه، افت حدود 13.5 dB خواهد بود.

این مقایسه نشان می‌دهد که Beam باریک‌تر می‌تواند نسبت به خطای جهت‌گیری حساس‌تر باشد.

مدل گاوسی برای تحلیل آموزشی به کار رفته و در طراحی واقعی باید با Radiation Pattern و داده‌های سازنده راستی‌آزمایی شود.

هندسه سطح جامد و محاسبه خطای سطح در سیلو

در اندازه‌گیری سطح مواد جامد، شکل سطح محصول می‌تواند باعث اختلاف میان ارتفاع نقطه‌ای و ارتفاع نماینده حجم شود.

در مثال کتابچه، سطح ماده دانه‌ای به شکل مخروط با زاویه ریزش مشخص در نظر گرفته شده است.

اگر شعاع سیلو R و زاویه ریزش φ باشد، ارتفاع مخروط از رابطه زیر به دست می‌آید:

hc = R × tan(φ)

ارتفاع معادل حجمی مخروط:

heq = hc / 3

در فاصله شعاعی r از محور، ارتفاع محلی در مدل مخروطی برابر است با:

hr = hc × (1 − r/R)

برای شعاع فرضی 2 متر و زاویه ریزش 35 درجه، ارتفاع مخروط 1.40 متر و ارتفاع معادل حجمی حدود 0.467 متر محاسبه شده است.

اگر Antenna در فاصله یک متر از محور قرار گیرد، ارتفاع محلی 0.70 متر است.

در این مثال، اختلاف ارتفاع نقطه‌ای و ارتفاع معادل حجمی حدود 0.233 متر خواهد بود.

این اختلاف نشان می‌دهد که Measurement Accuracy بالا در یک نقطه لزوماً به معنای Accuracy مناسب برای برآورد Inventory Volume نیست.

در سیلوها، شکل سطح، زاویه ریزش، محل پرشدن و تعریف کمیت موردنیاز باید در طراحی لحاظ شوند.

تحلیل Echo Curve در سه سطح برای Commissioning

بررسی Echo Curve فقط در یک سطح ممکن است برای شناسایی تمام تداخل‌های احتمالی کافی نباشد.

در مطالعه آموزشی کتابچه، Echo Curve برای سه وضعیت Nearly Empty،