LEVEL 3 — ADVANCED ENGINEERING

مهندسی پیشرفته دبی جرمی و حجمی — جبران چگالی، ضریب تراکم‌پذیری و تحلیل عدم‌قطعیت RSS

Advanced Engineering of Mass and Volumetric Flow — Density Compensation and Uncertainty Analysis

این راهنما ادامه راهنمای سطح ۱ و سطح ۲ همین مجموعه است: وقتی دبی جرمی از ترکیب یک اندازه‌گیری حجمی/سرعتی با چگالی به‌طور غیرمستقیم استخراج می‌شود، چگونه می‌توان مدل چگالی درست را انتخاب کرد، خطای نادیده‌گرفتن ضریب تراکم‌پذیری Z را کمّی کرد، عدم‌قطعیت را با روش RSS ترکیب کرد، و با FMEA و یک ماتریس راستی‌آزمایی طراحی، تصمیم را پیش از اجرا تأیید کرد — با چهار مثال مهندسی کامل.

شماره سندNICO-EDU-FLW-L3-010-FA بازنگریRev00 انتشارشهریور ۱۴۰۵ مدت مطالعه26 min

مقدمه

راهنمای سطح ۱ همین مجموعه ‎(‎NICO-EDU-FLW-L1-008-FA‎)‎ تفاوت مفهومی دبی جرمی و دبی حجمی را معرفی کرد، و راهنمای سطح ۲ ‎(‎NICO-EDU-FLW-L2-009-FA‎)‎ معیارهای انتخاب و نصب هفت خانواده فناوری اندازه‌گیری دبی را پوشش داد. این سند یک سؤال متفاوت را پاسخ می‌دهد: وقتی دبی جرمی به‌طور غیرمستقیم — از ترکیب یک اندازه‌گیری حجمی/سرعتی با یک مقدار چگالی — استخراج می‌شود، چگونه می‌توان این استخراج را با مدل چگالی درست، معادله درست و عدم‌قطعیت کمّی‌شده و قابل‌دفاع انجام داد؛ و چه زمانی اندازه‌گیری مستقیم جرم با Coriolis از نظر مهندسی توجیه دارد؟

NOTE — این سند گسترش سطح ۲ با واژگان سنگین‌تر نیست

هر بخش این راهنما یا یک رابطه کمّی جدید مشتق می‌کند یا یک مکانیزم خطای قابل‌محاسبه معرفی می‌کند که در سطح ۲ عمداً پوشش داده نشد: مدل چگالی مایع و گاز، اثر Compressibility Factor، انضباط کمّی Base Conditions، فیزیک Coriolis، مدل خطای غیرجمع‌پذیر، ترکیب RSS، تحلیل حساسیت، FMEA و ماتریس راستی‌آزمایی طراحی.

پس از مطالعه این راهنما خواهید توانست:
  1. دبی جرمی، دبی حجمی واقعی و دبی حجمی استاندارد/نرمال را به‌طور دقیق و بدون ابهام تعریف کنید.
  2. چگالی مایع را با مدل خطی دما و چگالی گاز را با مدل ایده‌آل و واقعی (Compressibility Factor) محاسبه کنید.
  3. خطای ناشی از نادیده‌گرفتن ضریب تراکم‌پذیری Z را در یک گاز متراکم کمّی کنید.
  4. اصل فیزیکی Coriolis برای اندازه‌گیری مستقیم جرم و چگالی جانبی را توضیح دهید.
  5. عدم‌قطعیت دبی جرمی غیرمستقیم را با روش RSS ترکیب کنید — نه با جمع ساده.
  6. گزینه اندازه‌گیری مستقیم را در برابر گزینه غیرمستقیم با چگالی ثابت مقایسه کنید.
  7. حالات خرابی مرتبط با چگالی و دبی جرمی غیرمستقیم را با مکانیزم و اقدام اصلاحی فهرست کنید.
  8. یک ماتریس راستی‌آزمایی طراحی کامل کنید و مشخص کنید چه مواردی هنوز باز هستند.

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

DESIGN BASIS — داده‌های ورودی الزامی
  1. سیال فرآیند و فاز. مایع تک‌فاز، گاز، یا بخار — هرکدام مدل چگالی متفاوتی دارند.
  2. محدوده فشار و دمای عملیاتی واقعی. نه یک نقطه تکی، بلکه کل بازه سرویس — چگالی در این بازه چقدر تغییر می‌کند؟
  3. ترکیب گاز (در صورت گاز/بخار). جرم مولی مؤثر M و منبع همبستگی ضریب تراکم‌پذیری Z.
  4. روش اندازه‌گیری دبی انتخابی. مستقیم (Coriolis) یا غیرمستقیم (حجمی/سرعتی + چگالی).
  5. منبع مقدار چگالی در روش غیرمستقیم. چگالی‌سنج برخط، مقدار ثابت از برگه داده، یا محاسبه از فشار/دما.
  6. شرایط مرجع (Base Conditions) الزامی گزارش. اگر خروجی، دبی حجمی استاندارد/نرمال باشد.
  7. عدم‌قطعیت هر مؤلفه اندازه‌گیری. از برگه داده تأییدشده — نه فرض این سند.
  8. معیار پذیرش عدم‌قطعیت ترکیبی دبی جرمی. الزام پروژه یا قرارداد اندازه‌گیری (مثلاً برای تحویل مالی).
  9. الزامات نگهداری و دسترسی برای بازرسی/کالیبراسیون دوره‌ای.

جدول ۱ — نمادها

نمادمعناواحد
m_dotدبی جرمیkg/h
Q_actualدبی حجمی واقعی در شرایط خطm³/h
Q_stdدبی حجمی استاندارد/نرمال در شرایط مرجعNm³/h, Sm³/h
rhoچگالی سیالkg/m³
T، T_refدمای مطلق واقعی و مرجعK, °C
P، P_refفشار مطلق واقعی و مرجعPa, bar
betaضریب انبساط حرارتی حجمی مایع1/°C
Mجرم مولی مؤثر گازkg/kmol
Rثابت جهانی گازهاJ/(kmol·K)
Zضریب تراکم‌پذیری (Compressibility Factor)بدون بعد
u_Q، u_rho، u_mdotعدم‌قطعیت نسبی دبی حجمی، چگالی و دبی جرمی ترکیبی%

تعریف دقیق Measurand

پیش‌شرط هر تحلیل عدم‌قطعیت، تعریف بدون ابهام کمیت اندازه‌گیری‌شونده (Measurand) است (STANDARD DEFINITION — سازگار با واژگان ISO 4006:1991). دبی جرمی یک کمیت پایسته است — مستقل از فشار و دمای محل اندازه‌گیری — اما دبی حجمی چنین نیست:

Q_actual = m_dot / rho_actual
Q_std = m_dot / rho(P_ref, T_ref)

همان m_dot، بسته به اینکه کدام چگالی — واقعی یا مرجع — در مخرج قرار گیرد، به دو عدد حجمی متفاوت تبدیل می‌شود. گزارش یک عدد حجمی «نرمال» یا «استاندارد» بدون بیان صریح P_ref و T_ref، از نظر علمی ناقص است.

MASS FLOWm_dot [kg/h] - invariantACTUAL VOLUME FLOWQ_actual = m_dot / rho(P,T)_actualSTANDARD VOLUME FLOWQ_std = m_dot / rho(P_ref,T_ref)same m_dot; two different volumetric numbers depending on WHICH density is usedreporting a volumetric number without stating actual vs. standard conditions is incompleteMEASURAND MAP
شکل ۱ — نقشه Measurand؛ دبی جرمی پایسته، دبی حجمی واقعی و استاندارد هر دو وابسته به چگالی

چگالی مایع — مدل خطی با دما

در بازه دمایی متعارف صنعتی و دور از نقطه تغییر فاز، چگالی اغلب مایعات با یک مدل خطی ساده (ENGINEERING CALCULATION) به‌خوبی تقریب زده می‌شود:

rho(T) ≈ rho(T_ref) × [1 − beta × (T − T_ref)]

این مدل برای بازه دمایی وسیع یا نزدیک نقطه تغییر فاز معتبر نیست. ضریب beta باید از برگه داده سیال واقعی پروژه تأیید شود — REQUIRES SOURCE VERIFICATION، نه مقدار عمومی این سند.

rhotemperature TT_refLIQUID DENSITY vs TEMPERATURE — LINEAR MODELrho(T) approx rho(T_ref) * (1 - beta * (T - T_ref)) - valid over a moderate T range only
شکل ۲ — چگالی مایع در برابر دما؛ مدل خطی با ضریب انبساط beta

چگالی گاز — قانون گاز ایده‌آل و ضریب تراکم‌پذیری

برای گاز، مدل ایده‌آل نقطه شروع است (STANDARD DEFINITION — ترمودینامیک کلاسیک):

rho_ideal = P × M / (R × T)

در فشار پایین و دور از شرایط بحرانی، این مدل تقریب خوبی است. اما در فشار بالا — دقیقاً همان محدوده‌ای که بسیاری از خطوط انتقال گاز صنعتی در آن کار می‌کنند — انحراف رفتار واقعی گاز از مدل ایده‌آل دیگر قابل چشم‌پوشی نیست.

rhopressure Pideal gas: rho_ideal = P*M/(R*T)real gas: rho_real = P*M/(Z*R*T)GAS DENSITY — IDEAL vs REAL (Z EFFECT)deviation grows with pressure; Z must come from a verified correlation (e.g. AGA Report No. 8), not assumed = 1
شکل ۳ — چگالی گاز؛ انحراف مدل واقعی از مدل ایده‌آل با افزایش فشار

معادله عمومی چگالی گاز واقعی

انحراف رفتار واقعی با یک ضریب اصلاحی — ضریب تراکم‌پذیری (Compressibility Factor) Z — وارد معادله می‌شود:

rho_real = P × M / (Z × R × T)

Z نزدیک فشار پایین به ۱ میل می‌کند و در فشار بالا یا نزدیک شرایط بحرانی می‌تواند به‌طور محسوس کمتر یا بیشتر از ۱ شود. مقدار دقیق Z باید از یک همبستگی استاندارد تأییدشده — مانند AGA Report No. 8 یا API MPMS Chapter 14.2 برای گاز طبیعی — محاسبه شود، نه فرض گرفته شود (REFERENCE؛ محاسبه Z در این سند صرفاً مفهومی است).

ENGINEERING TIP — چه زمانی Z را نمی‌توان نادیده گرفت

یک قاعده تقریبی: هرچه فشار خط بالاتر و دمای عملیاتی نزدیک‌تر به دمای بحرانی گاز باشد، فاصله Z از ۱ بیشتر است. مثال کارشده B (بخش ۱۳) نشان می‌دهد یک Z ≈ 0.92 نمایشی در فشار متوسط صنعتی می‌تواند به‌تنهایی خطایی در حدود ۸٪ در دبی جرمی محاسبه‌شده ایجاد کند — بزرگ‌تر از بسیاری منابع خطای «رسمی» دیگر.

انضباط Base Conditions

گزارش دبی نرمال/استاندارد بدون بیان صریح دما و فشار مرجع (Base Conditions)، یک عدد ناقص علمی است (NICO GUIDANCE — بسط کمّی همان اصل سطح ۱). صنایع و کشورهای مختلف مراجع متفاوتی به کار می‌برند — نمونه‌های رایج 0°C / 1 atm، 15°C / 1 atm یا 15.56°C / 1 atm بسته به صنعت — و این سند هیچ مرجع خاصی را به‌عنوان «استاندارد جهانی» معرفی نمی‌کند.

INDUSTRY Areference: 0 degC / 1 atmINDUSTRY Breference: 15 degC / 1 atmQ_std,AQ_std,Bsame physical gas stream -> two DIFFERENT standard-volume numbersunless the reference (P_ref, T_ref) is stated, a normal/standard volume figure is not comparableBASE CONDITIONS DISCIPLINE
شکل ۴ — همان جریان فیزیکی گاز، دو مرجع متفاوت، دو عدد دبی حجمی استاندارد متفاوت
NOTE — انتقال مرجع بین دو گزارش

وقتی یک عدد Q_std باید از یک مرجع به مرجع دیگر تبدیل شود، تبدیل باید از طریق m_dot (رابطه ۲) انجام شود، نه با یک ضریب ساده تقریبی — در غیر این صورت خطای دور از انتظار وارد گزارش تحویل/صورت‌حساب می‌شود.

فیزیک عمیق‌تر Coriolis — رابطه اختلاف فاز و دبی جرمی

اصل فیزیکی Coriolis کلاسیک و شناخته‌شده است (ENGINEERING CALCULATION ساده‌شده، نه معادله کالیبراسیون سازنده): یک یا دو لوله در حال ارتعاش، سیال در حال جریان را حمل می‌کنند. نیروی Coriolis ناشی از ترکیب حرکت ارتعاشی لوله و حرکت خطی سیال داخل آن، باعث پیچش نامتقارن دو بازوی ورودی و خروجی لوله می‌شود؛ این پیچش به‌صورت یک اختلاف فاز/زمان کوچک بین دو حسگر سرعت روی دو بازو قابل اندازه‌گیری است:

dt = k × m_dot

k ثابت هندسی/مکانیکی لوله است — تابعی از جنس، قطر، طول و شکل لوله — و اختصاصی هر سازنده و هر مدل است (MANUFACTURER-SPECIFIC DATA). رابطه کالیبراسیون واقعی هر Coriolis تجاری، فرمول اختصاصی سازنده است؛ رابطه بالا فقط تناسب مستقیم فیزیکی را نشان می‌دهد، نه دقت کالیبراسیون واقعی هیچ تجهیزی.

no-flow: tube twists symmetrically, no phase shiftinlet leg sensoroutlet leg sensorCORIOLIS TUBE TWIST -> PHASE SHIFTflowing mass generates Coriolis force on the vibrating tube -> inlet/outlet legs twist out of phase -> dt
شکل ۵ — پیچش لوله Coriolis تحت جریان و اختلاف فاز حاصل بین دو بازو

چگالی جانبی از فرکانس ارتعاش Coriolis

علاوه بر دبی جرمی، یک Coriolis معمولاً چگالی سیال را نیز به‌عنوان خروجی جانبی گزارش می‌کند (CONFIRMED کیفی، ثابت‌های دقیق MANUFACTURER-SPECIFIC DATA). این خروجی از رابطه بین فرکانس طبیعی ارتعاش لوله و جرم کل سیستم (لوله + سیال داخل آن) استخراج می‌شود: افزایش چگالی سیال، جرم نوسانی را افزایش می‌دهد و فرکانس طبیعی را کاهش می‌دهد. این خروجی چگالی، وقتی در دسترس و معتبر باشد، می‌تواند خود جایگزین یا مکمل مقدار چگالی فرضی در یک روش غیرمستقیم شود.

natural freq. ffluid density rhoDENSITY FROM CORIOLIS VIBRATION FREQUENCYtube mass increases with fluid density -> natural vibration frequency drops (qualitative relation)
شکل ۶ — رابطه کیفی فرکانس طبیعی ارتعاش لوله Coriolis و چگالی سیال

مدل خطای ترکیبی دبی جرمی غیرمستقیم

وقتی دبی جرمی به‌طور غیرمستقیم از ضرب یک دبی حجمی اندازه‌گیری‌شده در یک مقدار چگالی به دست می‌آید:

m_dot = Q_actual × rho

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

VOLUMETRIC METERu_QDENSITY VALUEu_rhom_dot= Q * rhoINDIRECT MASS FLOWu_mdot/mdot = RSS(u_Q, u_rho)COMBINED ERROR MODEL — INDIRECT MASS FLOWtwo independent measurement uncertainties propagate into the derived mass-flow uncertainty
شکل ۷ — مدل خطای ترکیبی دبی جرمی غیرمستقیم؛ دو منبع مستقل به یک عدم‌قطعیت نهایی ترکیب می‌شوند

بودجه عدم‌قطعیت — اصل RSS

روش‌شناسی صحیح ترکیب دو منبع خطای آماری مستقل، جذر مجموع مربعات (Root-Sum-Square، RSS) است — نه جمع خطی (NORMATIVE REQUIREMENT — چارچوب JCGM 100:2008 / ISO/IEC Guide 98-3:2008 (GUM)؛ دقیقاً همان اصلی که در راهنمای سطح ۳ خوشه فشار ‎(‎NICO-EDU-PRS-L3-005-FA‎)‎ برای ترکیب منابع خطای ترانسمیتر به‌کار رفت):

(u_mdot / m_dot)² = (u_Q / Q)² + (u_rho / rho)²

این رابطه فرض می‌کند دو منبع خطا آماری مستقل‌اند. در صورت وجود همبستگی بین آن‌ها (مثلاً هر دو از یک حسگر مشترک تأثیر می‌گیرند)، این فرض باید بازبینی شود — یک مورد ASSUMPTION در محدوده این سند.

u_Qu_rhou_mdot (RSS)u_mdot = sqrt(u_Q^2 + u_rho^2)the two error vectors combine like perpendicular components, not by simple additionRSS UNCERTAINTY COMBINATION (GEOMETRIC VIEW)
شکل ۸ — نمایش هندسی ترکیب RSS؛ دو بردار خطای عمود بر هم، نه جمع خطی

مثال کارشده A — مایع، انتشار عدم‌قطعیت

آب صنعتی، اندازه‌گیری غیرمستقیم دبی جرمی WORKED EXAMPLE

داده (ASSUMPTION — مقادیر نوعی برای این مثال). آب در دمای نزدیک محیط: rho = 998 kg/m³، دبی حجمی واقعی اندازه‌گیری‌شده Q_actual = 50 m³/h. عدم‌قطعیت مرجع فلومتر حجمی (MANUFACTURER-SPECIFIC DATA نمایشی): u_Q = 0.50%؛ عدم‌قطعیت چگالی‌سنج برخط: u_rho = 0.30%.

گام ۱ — دبی جرمی (رابطه ۷).

m_dot = 50 × 998 = 49900 kg/h

گام ۲ — ترکیب عدم‌قطعیت (رابطه ۸).

u_mdot / m_dot = sqrt(0.50² + 0.30²) ≈ 0.583%

نتیجه: u_mdot ≈ 0.583% × 49900 ≈ 291 kg/h. جمع خطی ساده (0.50 + 0.30 = 0.80%) روشی نادرست بود و عدد محافظه‌کارانه‌تری غیرواقعی می‌داد — ترکیب RSS عدد کوچک‌تر و درست‌تر 0.583% را می‌دهد (جدول ۲).

جدول ۲ — بودجه عدم‌قطعیت مثال کارشده A (% از m_dot)

مؤلفهسهممبنا
فلومتر حجمی (u_Q)±0.50%مرجع سازنده نمایشی
چگالی‌سنج برخط (u_rho)±0.30%مرجع سازنده نمایشی
ترکیب RSS (u_mdot)±0.583%رابطه ۸ — design_calc.py Case A
0.50%volumetric meter0.30%density valueRSS total = 0.58%UNCERTAINTY BUDGET — INDIRECT MASS FLOW (CASE A)component values from design_calc.py Case A; RSS total is smaller than a simple sum
شکل ۹ — بودجه عدم‌قطعیت مثال کارشده A و ترکیب RSS آن

مثال کارشده B — گاز متراکم، تصحیح Z

گاز طبیعی فشار متوسط، خطای نادیده‌گرفتن Z WORKED EXAMPLE

داده (ASSUMPTION — مقادیر نمایشی). P = 51 bara، T = 40°C (313.15 K)، M ≈ 17.4 kg/kmol (گاز طبیعی خط لوله، نوعی)، Z ≈ 0.92 (نمایشی — REQUIRES VERIFICATION در برابر AGA Report No. 8)، دبی حجمی واقعی اندازه‌گیری‌شده Q_actual = 800 m³/h.

گام ۱ — چگالی ایده‌آل و واقعی (روابط ۴ و ۵).

rho_ideal = 51e5 × 17.4 / (8314 × 313.15) ≈ 34.085 kg/m³
rho_real = rho_ideal / Z ≈ 37.048 kg/m³

گام ۲ — دبی جرمی واقعی در برابر دبی جرمی محاسبه‌شده با فرض غلط Z=1.

m_dot,true = 800 × 37.048 ≈ 29639 kg/h m_dot,wrong = 800 × 34.085 ≈ 27268 kg/h

نتیجه: نادیده‌گرفتن Z در این نمونه دبی جرمی را حدود ۸.۰٪ کمتر از مقدار واقعی نشان می‌دهد. این یک بایاس سیستماتیک است — نه یک مؤلفه عدم‌قطعیت تصادفی — و با ترکیب RSS قابل «میانگین‌گیری» نیست؛ باید در مدل چگالی اصلاح شود، نه در بودجه عدم‌قطعیت لحاظ شود.

COMMON MISTAKE — درمان بایاس سیستماتیک به‌عنوان عدم‌قطعیت تصادفی

خطای نادیده‌گرفتن Z، خطای دمای Wet Leg، یا استفاده از یک چگالی ثابت قدیمی برای گاز/بخاری که شرایطش تغییر کرده، همگی بایاس سیستماتیک‌اند — همان جهت خطا در هر بار اندازه‌گیری تکرار می‌شود. اشتباه رایج طراحی، افزودن این نوع خطا به بودجه عدم‌قطعیت RSS (که برای منابع تصادفی/مستقل ساخته شده) به‌جای اصلاح مستقیم مدل است. اقدام درست: مدل چگالی (Z، دما، مرجع) را اصلاح کنید؛ سپس عدم‌قطعیت باقیمانده را با RSS بودجه‌بندی کنید.

مثال کارشده C — Coriolis مستقیم در برابر غیرمستقیم

مقایسه عدم‌قطعیت دو معماری اندازه‌گیری برای همان سرویس WORKED EXAMPLE

داده (MANUFACTURER-SPECIFIC DATA نمایشی / ASSUMPTION). Coriolis (مستقیم): عدم‌قطعیت نمایشی سازنده u_mdot ≈ 0.15% از قرائت. گزینه غیرمستقیم: فلومتر Vortex با u_Q = 1.00%، به‌همراه یک مقدار چگالی ثابت فرضی که خطای احتمالی آن نسبت به چگالی واقعی متغیر فرآیند (نه یک عدم‌قطعیت کالیبراسیون تصادفی، بلکه یک باند ریسک بایاس) حدود ±2.00% برآورد می‌شود.

گام ۱ — ترکیب RSS برای گزینه غیرمستقیم (رابطه ۸).

u_indirect = sqrt(1.00² + 2.00²) ≈ 2.24%

نتیجه: در این مثال نمایشی، عدم‌قطعیت مسیر غیرمستقیم (≈2.24%) حدود ۱۵ برابر عدد Coriolis مستقیم است — و این فاصله تقریباً به‌طور کامل از فرض چگالی ثابت ناشی می‌شود، نه از خود فلومتر حجمی. این محاسبه توجیه مهندسی برای انتخاب اندازه‌گیری مستقیم جرم را در سرویس‌هایی با چگالی متغیر نشان می‌دهد — بدون ادعای برتری مطلق Coriolis در همه شرایط (بخش ۱۸).

مثال کارشده D — بخار، اثر ترکیبی دما و فشار

بخار فوق‌گرم، رانش دما بدون به‌روزرسانی چگالی WORKED EXAMPLE

داده (ASSUMPTION — REQUIRES VERIFICATION در برابر جداول بخار IAPWS-IF97). چگالی نمایشی بخار فوق‌گرم در 200°C: rho ≈ 4.40 kg/m³. دمای فرآیند به 220°C افزایش می‌یابد (فشار عملاً ثابت)، اما مقدار چگالی استفاده‌شده در محاسبه دبی جرمی به‌روزرسانی نمی‌شود. دبی حجمی واقعی ثابت: Q_actual = 300 m³/h.

گام ۱ — چگالی تقریبی در دمای جدید (تقریب آموزشی، نه معادله بخار واقعی).

rho(220°C) ≈ 4.40 × (473.15/493.15) ≈ 4.222 kg/m³

گام ۲ — دبی جرمی واقعی در برابر دبی جرمی با چگالی رانش‌نکرده.

m_dot,true = 300 × 4.222 ≈ 1266 kg/h m_dot,stale = 300 × 4.40 ≈ 1320 kg/h

نتیجه: عدم به‌روزرسانی چگالی پس از رانش دما، خطایی حدود ۴.۲٪+ در دبی جرمی گزارش‌شده ایجاد می‌کند. برای سیالات تراکم‌پذیر/فوق‌گرم، یک مقدار چگالی ثابت تنها برای یک نقطه (P,T) معتبر است؛ بدون جبران برخط، رانش دما یا فشار مستقیماً و بدون اصلاح وارد دبی جرمی گزارش‌شده می‌شود.

تحلیل حساسیت و نمودار Tornado

مقایسه مرتبه‌بزرگی منابع خطای کمّی‌شده در بخش‌های ۱۲ تا ۱۵ نشان می‌دهد کدام‌یک باید در بازبینی طراحی اولویت بگیرد:

جدول ۳ — تحلیل حساسیت (بزرگی خطا/عدم‌قطعیت برای ورودی داده‌شده)

منبعنوعبزرگیرتبه
چگالی ثابت فرضی به‌جای برخط (Case C)بایاس ریسک±2.00%بالا
نادیده‌گرفتن Z در گاز متراکم (Case B)بایاس سیستماتیک≈8.0%بسیار بالا
فلومتر حجمی u_Q (Case A/C)عدم‌قطعیت تصادفی0.50–1.00%متوسط
چگالی‌سنج برخط u_rho (Case A)عدم‌قطعیت تصادفی±0.30%پایین
fixed-density assumption (Case C)Z-factor ignored (Case B)volumetric meter u_Qonline density u_rhoSENSITIVITY — RELATIVE MAGNITUDE OF ERROR SOURCEStornado-style ranking; magnitudes illustrative, from design_calc.py Cases A-C
شکل ۱۰ — نمودار Tornado؛ بایاس سیستماتیک ناشی از فرض چگالی غالب‌ترین منبع خطاست

سناریوی میدانی

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

«یک ایستگاه اندازه‌گیری گاز فروش، دبی حجمی خط را با یک Vortex اندازه می‌گرفت و برای تبدیل به دبی جرمی/انرژی از یک چگالی ثابت — تعیین‌شده در زمان راه‌اندازی — استفاده می‌کرد. پس از چند ماه، مقایسه با یک اندازه‌گیری مرجع مستقل نشان داد صورت‌حساب محاسبه‌شده به‌طور پیوسته کمتر از مقدار واقعی است.»

گام ۱ — چرا چگالی ثابت زمان راه‌اندازی کافی نبود

فشار خط تغذیه در طول سال، متناسب با تقاضای شبکه، در بازه‌ای گسترده نوسان می‌کرد. رابطه ۵ نشان می‌دهد چگالی واقعی گاز به‌طور مستقیم به فشار وابسته است؛ یک مقدار ثابت فقط برای فشار لحظه راه‌اندازی معتبر بود.

گام ۲ — چه چیزی نادیده گرفته شده بود

علاوه بر فشار متغیر، ضریب تراکم‌پذیری Z نیز در فشار بالاتر خط از ۱ فاصله می‌گرفت — دقیقاً مکانیزم Worked Case B (بخش ۱۳). ترکیب این دو اثر — فشار متغیر نادیده‌گرفته‌شده و Z نادیده‌گرفته‌شده — بایاس منفی پایدار در صورت‌حساب ایجاد کرده بود.

گام ۳ — اقدام اصلاحی و چرا کافی بود

جایگزینی چگالی ثابت با یک محاسبه برخط از فشار و دمای لحظه‌ای خط، به‌همراه یک همبستگی Z تأییدشده (بخش ۶)، بایاس سیستماتیک را عملاً حذف کرد؛ عدم‌قطعیت باقیمانده سپس با روش RSS (بخش ۱۱) در گواهی اندازه‌گیری مستند شد.

مطالعه سه‌گزینه‌ای طراحی

برای یک سرویس با چگالی متغیر، سه معماری قابل مقایسه‌اند — هیچ‌کدام برتری مطلق ندارد (NICO GUIDANCE؛ تصمیم باید مبتنی بر داده واقعی سرویس باشد، نه رتبه‌بندی عمومی این سند):

DIRECT MASSCoriolisINDIRECT, ONLINE RHOVolumetric + densitometerINDIRECT, FIXED RHOVolumetric + assumed densityeach option trades installed cost, complexity and uncertainty differentlyno single option is universally correct - see Section 17 for the trade-off discussionTHREE-OPTION DESIGN COMPARISON
شکل ۱۱ — سه گزینه طراحی برای اندازه‌گیری دبی جرمی یک سرویس با چگالی متغیر

جدول ۴ — مقایسه سه گزینه طراحی

گزینهعدم‌قطعیت نمونهپیچیدگی/هزینه نصبریسک اصلی
مستقیم — Coriolisکم (Case C)بالاتر (تجهیز اختصاصی)افت فشار، محدودیت اندازه خط، هزینه اولیه
غیرمستقیم — حجمی + چگالی‌سنج برخطمتوسط (Case A)متوسط (دو تجهیز)نگهداری/کالیبراسیون دو حلقه مستقل
غیرمستقیم — حجمی + چگالی ثابتبالا در چگالی متغیر (Case C)پایین‌ترینبایاس سیستماتیک نامرئی (بخش ۱۳)

تحلیل حالات خرابی (FMEA)

Failure modePhysical mechanismReading behaviourDiagnostic evidenceEngineering correctionFMEA LOGIC CHAINeach failure mode in the FMEA table follows this same five-column chain
شکل ۱۲ — زنجیره منطقی FMEA استفاده‌شده در جدول ۵

جدول ۵ — حالات خرابی، مکانیزم فیزیکی، رفتار قرائت، شواهد تشخیصی و اصلاح مهندسی

حالت خرابیمکانیزم فیزیکیرفتار قرائت مورد انتظارشواهد تشخیصیاصلاح مهندسی
نادیده‌گرفتن Z در محاسبه چگالی گازاستفاده از مدل ایده‌آل در فشار بالابایاس منفی پایدار در دبی جرمیمقایسه با اندازه‌گیری مرجع مستقل یا Coriolisبه‌کارگیری همبستگی Z تأییدشده (بخش ۶)
چگالی ثابت رانش‌نکرده با تغییر P/T فرآیندعدم به‌روزرسانی مقدار چگالی پس از تغییر شرایطبایاس تدریجی هم‌جهت با رانش P/Tهمبستگی خطا با سابقه فشار/دمای خطمحاسبه چگالی برخط از P و T لحظه‌ای
خطای کالیبراسیون یا رانش چگالی‌سنج برخطکثیف‌شدگی، رسوب یا رانش سنسوربایاس یا نویز افزایش‌یابنده در چگالی گزارش‌شدهکالیبراسیون دوره‌ای در برابر نمونه مرجعسرویس/کالیبراسیون دوره‌ای طبق رویه سازنده
تغییر ترکیب گاز نسبت به مقدار M فرض‌شدهتغییر منبع تغذیه یا اختلاط جریان‌هابایاس متغیر، همبسته با منبع/زمانآنالیز کروماتوگرافی دوره‌ای گازبه‌روزرسانی M طبق آنالیز واقعی، یا آنالایزر آنلاین
جریان دوفازی ناخواسته در خط مایع/گازحباب گاز در مایع یا قطرات مایع در گازنوسان یا خطای بزرگ غیرقابل‌پیش‌بینینوسان شدید و ناپایدار قرائتجداساز دوفازی، بازبینی نقطه نمونه‌برداری
فرسایش یا رسوب داخل لوله Coriolisسایش مکانیکی یا رسوب‌گذاری روی سطح داخلی لولهZero shift یا خطای دبی جرمی تدریجیمقایسه Zero فعلی با Zero کارخانه‌ایبازرسی/تعویض طبق برنامه نگهداری پیشگیرانه
تنش نصب روی لوله Coriolisعدم تراز یا تنش مکانیکی از لوله‌کشی مجاورZero shift ثابت پس از نصب/تعمیرZero Verification بلافاصله پس از نصبپشتیبانی مکانیکی صحیح طبق دستورالعمل نصب سازنده
کثیف‌شدگی یا خوردگی حسگر Vortexرسوب یا خوردگی روی جسم مولد گردابافت حساسیت یا از دست‌رفتن سیگنال در دبی پایینعدم تشخیص گرداب زیر یک آستانه دبی مشخصبازرسی/تمیزکاری دوره‌ای، انتخاب متریال مقاوم
خطای مرجع Base Conditions نامنطبق بین دو سامانهدو سیستم گزارش با P_ref/T_ref متفاوت مقایسه می‌شونداختلاف پایدار بین دو گزارش بدون خطای واقعی سنسوربازبینی مستندات پیکربندی هر سیستمیکسان‌سازی یا تبدیل صریح مرجع (بخش ۷)
کیفیت/رطوبت بخار متفاوت از فرض خشک اشباعوجود قطرات مایع در جریان بخارخطای مثبت یا منفی غیرقابل‌پیش‌بینی در دبی جرمیاندازه‌گیری کیفیت بخار در صورت شکجداساز رطوبت، یا انتخاب فناوری کمتر حساس به کیفیت
مقایسه نادرست دبی حجمی استاندارد بین دو صنعت/قراردادعدم تبدیل صریح بین دو Base Conditions متفاوتاختلاف عددی بدون خطای واقعی اندازه‌گیریبازبینی توافق‌نامه مرجع اندازه‌گیری در قراردادتبدیل از طریق m_dot طبق رابطه ۱/۲ (بخش ۷)

ماتریس راستی‌آزمایی طراحی

RequirementDesign inputCalculated valueAcceptance basisStatus + evidenceDESIGN VERIFICATION MATRIX — RECORD STRUCTUREone row of this chain is recorded per verified design requirement
شکل ۱۳ — ساختار ثبت ماتریس راستی‌آزمایی طراحی (جدول ۶)

جدول ۶ — الزام، ورودی طراحی، مقدار محاسبه‌شده، مبنای پذیرش، وضعیت و شواهد

الزامورودی طراحیمقدار محاسبه‌شدهمبنای پذیرشوضعیتمورد باز
عدم‌قطعیت دبی جرمی غیرمستقیم (مایع)u_Q=0.50%، u_rho=0.30%≈0.583%< 1.0%قبول (این مثال)—
خطای نادیده‌گرفتن Z در گاز متراکمP=51 bara، Z=0.92 (نمایشی)≈8.0%< 1.0% یا اصلاح مدلرد — نیازمند اصلاح مدلهمبستگی Z تأییدشده (AGA 8)
گزینه مستقیم در برابر غیرمستقیم با چگالی ثابتu_mdot,Coriolis=0.15%، u_indirect≈2.24%نسبت ≈۱۵×تصمیم مهندسی مستندمشروطتأیید داده واقعی سرویس هدف
خطای چگالی رانش‌نکرده در سیال تراکم‌پذیرdT=+20°C در بخار≈4.2%+< 2% یا جبران برخطرد — نیازمند جبران برخطافزودن اندازه‌گیری P/T برخط
مبنای ضریب beta مایع مورد استفاده——برگه داده سیال واقعیدر انتظار راستی‌آزماییتأیید beta از برگه داده پروژه
مبنای همبستگی Z گاز مورد استفاده——AGA Report No. 8 / ترکیب واقعی گازدر انتظار راستی‌آزماییانتخاب همبستگی و تأیید ترکیب گاز پروژه
ENGINEERING NOTE — چرا دو ردیف «رد» طبیعی است

دو ردیف «رد» در این جدول عمداً نشان داده شده‌اند تا ماهیت واقعی ماتریس راستی‌آزمایی را نشان دهند: این ابزار برای کشف خطاهای اصلاح‌نشده ساخته شده، نه برای تأیید خودکار هر طراحی. یک طراحی واقعی که به این نتایج برسد باید مدل چگالی (Z یا جبران برخط) را پیش از اجرا اصلاح کند.

اشتباهات مفهومی رایج پیشرفته

COMMON MISTAKE — استفاده از یک چگالی «کارخانه‌ای» ثابت برای همیشه

یک چگالی اندازه‌گیری‌شده یا محاسبه‌شده در زمان راه‌اندازی، فقط برای شرایط فشار/دمای همان لحظه معتبر است. اشتباه رایج، فرض تلویحی این‌که این عدد برای کل عمر بهره‌برداری ثابت می‌ماند — دقیقاً مکانیزمی که Worked Case B، Worked Case D و سناریوی میدانی (بخش ۱۷) هرکدام با اعداد نشان دادند.

سه اشتباه مفهومی دیگر که در بازبینی طراحی سطح پیشرفته باید صریحاً بررسی شوند:

حدود اعتبار و رجیستر فرض‌ها

LIMITS OF APPLICABILITY — حدود اعتبار روش این سند
  1. مدل خطی چگالی مایع (رابطه ۳) فقط در بازه دمایی متعارف و دور از نقطه تغییر فاز معتبر است. برای بازه دمایی وسیع یا نزدیک تبخیر/انجماد، یک مدل غیرخطی لازم است.
  2. مقدار Z در Worked Case B نمایشی و صرفاً برای نشان‌دادن مکانیزم خطاست. REQUIRES SOURCE VERIFICATION در برابر یک همبستگی استاندارد مانند AGA Report No. 8 برای ترکیب گاز واقعی پروژه.
  3. رابطه ساده‌شده Coriolis (رابطه ۶) یک تقریب آموزشی است، نه معادله کالیبراسیون واقعی هیچ سازنده‌ای. هیچ عدد این بخش نباید به‌عنوان مشخصه NICO یا سازنده دیگر نقل شود.
  4. ترکیب RSS (رابطه ۸) فرض استقلال آماری دو منبع خطا را دارد. در صورت همبستگی واقعی بین دو منبع، این فرض باید بازبینی و روش ترکیب اصلاح شود.
  5. تقریب چگالی بخار در Worked Case D صرفاً یک مدل آموزشی ساده‌شده (نسبت مستقیم با دما) است. REQUIRES VERIFICATION در برابر جداول بخار IAPWS-IF97 برای هر کاربرد واقعی.
  6. اعداد جدول ۲، ۳ و ۴ نمایشی‌اند. مقادیر واقعی باید از برگه داده تجهیز و داده میدانی پروژه بازسازی شوند.

چک‌لیست مهندسی طراحی

DEFINE MEASURANDm_dot / Q_actual / Q_stdMODEL DENSITYliquid beta or gas P,T,ZDIRECT ORINDIRECT MASS FLOW?INDIRECTCOMBINE u_Q, u_rhoRSS (Sec. 10)BUILD UNCERTAINTYBUDGETCOMPLETE VERIFICATIONMATRIX (Sec. 19)
شکل ۱۴ — مسیر کامل از تعریف Measurand تا تکمیل ماتریس راستی‌آزمایی

چک‌لیست بازبینی طراحی

توالی راستی‌آزمایی راه‌اندازی

جمع‌بندی

این سند تعریف دقیق Measurand را نقطه شروع قرار داد و از آن، مدل چگالی مایع و گاز واقعی (با ضریب تراکم‌پذیری Z)، انضباط کمّی Base Conditions، فیزیک عمیق‌تر Coriolis، مدل خطای ترکیبی دبی جرمی غیرمستقیم و ترکیب عدم‌قطعیت به روش RSS را مشتق کرد. چهار مثال کارشده نشان دادند که بزرگ‌ترین منبع خطا در عمل اغلب یک بایاس سیستماتیک قابل‌اصلاح — نادیده‌گرفتن Z یا استفاده از چگالی ثابت رانش‌نکرده — است، نه عدم‌قطعیت تصادفی تجهیز. تحلیل حساسیت، مطالعه سه‌گزینه‌ای طراحی، FMEA و ماتریس راستی‌آزمایی، مسیر کامل از «طراحی مفهومی» تا «تأیید عددی قبل از اجرا» را می‌بندند. مقادیر Z، beta و چگالی بخار این سند نمایشی‌اند و باید پیش از تأیید نهایی طراحی، با داده واقعی پروژه جایگزین شوند.

جمع‌بندی نکات کلیدی
  • دبی جرمی پایسته است؛ دبی حجمی واقعی و استاندارد هر دو وابسته به چگالی‌اند و باید با مرجع صریح گزارش شوند.
  • نادیده‌گرفتن ضریب تراکم‌پذیری Z در گاز متراکم می‌تواند خطایی بزرگ‌تر از بسیاری منابع «رسمی» دیگر ایجاد کند.
  • بایاس سیستماتیک (Z، چگالی ثابت رانش‌نکرده) باید در مدل اصلاح شود؛ در بودجه RSS جذب نمی‌شود.
  • ترکیب صحیح منابع خطای مستقل RSS است، نه جمع خطی.
  • انتخاب مستقیم (Coriolis) در برابر غیرمستقیم باید با مقایسه عدم‌قطعیت مستند، نه صرفاً هزینه اولیه، تصمیم‌گیری شود.
  • ماتریس راستی‌آزمایی باید واقعی باشد؛ یک ردیف «رد» به معنای بازنگری مدل است، نه شکست فرآیند.

پرسش‌های متداول

چرا دبی حجمی استاندارد بدون ذکر Base Conditions یک عدد ناقص است؟
چون همان دبی جرمی، بسته به مرجع فشار/دمای انتخابی، به اعداد حجمی متفاوتی تبدیل می‌شود (رابطه ۲). بدون بیان صریح P_ref و T_ref، مقایسه یا تحویل بر مبنای این عدد ممکن نیست (بخش ۷).
چرا نادیده‌گرفتن ضریب تراکم‌پذیری Z مهم است؟
چون در فشار بالا، انحراف رفتار واقعی گاز از مدل ایده‌آل می‌تواند به‌تنهایی خطایی چند درصدی — در مثال کارشده B حدود ۸٪ — در دبی جرمی محاسبه‌شده ایجاد کند؛ عددی بزرگ‌تر از بسیاری منابع خطای «رسمی» دیگر (بخش ۱۳).
تفاوت بایاس سیستماتیک و عدم‌قطعیت تصادفی در این سند چیست؟
بایاس سیستماتیک (مثل خطای Z یا چگالی رانش‌نکرده) همیشه در یک جهت تکرار می‌شود و باید در مدل اصلاح شود. عدم‌قطعیت تصادفی (مثل u_Q یا u_rho) حول مقدار واقعی پراکنده است و با روش RSS ترکیب می‌شود. آمیختن این دو، تصمیم طراحی نادرست ایجاد می‌کند (بخش ۱۳، Common Mistake).
چرا ترکیب RSS به‌جای جمع ساده استفاده می‌شود؟
چون دو منبع خطای مستقل مانند دو بردار عمود بر هم عمل می‌کنند (شکل ۸)؛ جمع ساده همیشه محافظه‌کارانه‌تر از واقعیت آماری است و می‌تواند به رد بیش‌ازحد یک طراحی مناسب بینجامد (بخش ۱۱).
آیا Coriolis همیشه بهترین انتخاب برای دبی جرمی است؟
نه لزوماً. در سرویس با چگالی نسبتاً پایدار و چگالی‌سنج برخط قابل‌اعتماد، یک روش غیرمستقیم می‌تواند عدم‌قطعیت قابل‌قبولی داشته باشد. Coriolis زمانی مزیت روشنی دارد که چگالی به‌طور محسوس نوسان دارد یا نمی‌توان آن را به‌طور مستقل و مطمئن اندازه گرفت (بخش ۱۴، جدول ۴).
آیا خروجی چگالی جانبی یک Coriolis همیشه قابل‌اتکاست؟
اصل کیفی آن (فرکانس طبیعی وابسته به چگالی) شناخته‌شده است، اما دقت عددی آن خروجی، مشخصه‌ای اختصاصی هر سازنده و مدل است. استفاده مهندسی از آن باید بر مبنای برگه داده تأییدشده همان مدل باشد، نه فرض عمومی (بخش ۹).
اگر ماتریس راستی‌آزمایی یک ردیف «رد» نشان دهد، مرحله بعدی چیست؟
یعنی طراحی فعلی یکی از الزامات را برآورده نمی‌کند — نه اینکه سند اشتباه است. جدول ۶ نشان داد این می‌تواند به معنی نیاز به همبستگی Z تأییدشده یا افزودن جبران برخط چگالی باشد. تصمیم باید مستند و در همان جدول ثبت شود (بخش ۲۰).
آیا اعداد نمایشی این سند برای طراحی یا خرید واقعی قابل استفاده‌اند؟
نه به‌تنهایی. جدول ۶ به‌صراحت وضعیت هر الزام را «قبول»، «مشروط»، «رد» یا «در انتظار راستی‌آزمایی» اعلام می‌کند. این سند روش و ساختار محاسبه را نشان می‌دهد؛ تأیید نهایی نیازمند جایگزینی مقادیر نمایشی با داده واقعی پروژه و برگه داده تجهیز است.

واژه‌نامه

جدول ۷ — واژه‌نامه

اصطلاحتوضیح
Mass Flow (دبی جرمی)جرم سیال عبوری از یک مقطع در واحد زمان؛ کمیتی پایسته و مستقل از شرایط.
Actual Volume Flow (دبی حجمی واقعی)حجم سیال عبوری در شرایط فشار و دمای واقعی محل اندازه‌گیری.
Standard/Normal Volume Flow (دبی حجمی استاندارد/نرمال)معادل حجمی دبی جرمی، محاسبه‌شده با چگالی در یک شرایط مرجع تعریف‌شده.
Base Conditionsفشار و دمای مرجع مورد استفاده برای محاسبه دبی حجمی استاندارد/نرمال.
Compressibility Factor (ضریب تراکم‌پذیری)ضریب اصلاحی Z که انحراف رفتار گاز واقعی از مدل ایده‌آل را نشان می‌دهد.
Coriolisاصل اندازه‌گیری مستقیم دبی جرمی بر پایه نیروی Coriolis روی لوله در حال ارتعاش.
RSS (Root-Sum-Square)روش صحیح ترکیب عدم‌قطعیت منابع خطای آماری مستقل؛ جذر مجموع مربعات.
FMEAتحلیل حالات و اثرات خرابی؛ ابزار سیستماتیک شناسایی مکانیزم، رفتار و اصلاح خرابی‌های محتمل.
بایاس سیستماتیکخطایی که همواره در یک جهت تکرار می‌شود و باید در مدل اصلاح شود، نه در بودجه عدم‌قطعیت تصادفی.
ماتریس راستی‌آزمایی طراحیسند ثبت الزام، ورودی طراحی، مقدار محاسبه‌شده، مبنای پذیرش و وضعیت هر بند طراحی.

استانداردها و منابع

  1. ISO 4006:1991 — Measurement of fluid flow in closed conduits — Vocabulary and symbols مبنای واژگان Measurand، دبی جرمی و دبی حجمی در بخش‌های ۳ و ۷ (STANDARD DEFINITION).
  2. AGA Report No. 8 — Compressibility Factors of Natural Gas and Other Related Hydrocarbon Gases مرجع عمومی محاسبه Z برای گاز طبیعی؛ مقدار Z این سند نمایشی است و از این همبستگی محاسبه نشده — REFERENCE، نه CONFIRMED.
  3. API MPMS Chapter 14.2 — Compressibility Factors of Natural Gas and Other Related Hydrocarbon Gases مرجع API معادل AGA Report No. 8؛ عنوان و کاربرد عمومی، بدون نقل عدد اختصاصی در این سند.
  4. JCGM 100:2008 / ISO/IEC Guide 98-3:2008 (GUM) — Evaluation of measurement data — Guide to the expression of uncertainty in measurement روش‌شناسی ترکیب RSS به‌کاررفته در بخش ۱۱ و همه مثال‌های کارشده (NORMATIVE REQUIREMENT برای روش‌شناسی).
  5. Steam-table order-of-magnitude reference for superheated steam density near 200°C (general engineering reference) مبنای rho ≈ 4.40 kg/m³ در Worked Case D، بخش ۱۵ — ASSUMPTION، نیازمند تأیید در برابر IAPWS-IF97.
  6. General classical-mechanics reference for the Coriolis effect (general engineering reference) مبنای شکل کیفی رابطه ۶ در بخش ۸؛ ثابت k اختصاصی هر سازنده و هر مدل است.
نسخه کامل PDF فایل A4 با چهارده شکل فنی، دو جدول عدم‌قطعیت و حساسیت، جدول مقایسه سه گزینه طراحی، جدول FMEA ده‌ردیفه، ماتریس راستی‌آزمایی طراحی، یک سناریوی میدانی، چهار مثال مهندسی کامل با اعداد کامل و یک واژه‌نامه — مناسب بازبینی مهندسی و پیوست به مدارک پروژه.
دانلود NICO-EDU-FLW-L3-010-FA
مطالعه بعدی