---

4.1.9 Характеристики измерительного комплекса (измерение температуры)

SIZMTD - номер варианта состава измерительного комплекса:

1) 0 - термометр;

2) 1 - измерительный преобразователь и средство измерений с линейной функцией преобразования;

3) 2 - измерительный преобразователь и средство измерений с квадратичной функцией преобразования;

4) 3 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

5) 4 - измерительный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с квадратичной функцией преобразования;

6) 5 - измерительный преобразователь, средство измерений с квадратичной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования;

7) 6 - измерительный преобразователь, вторичный преобразователь, средство измерений с линейной функцией преобразования и планиметр (или интегратор) с линейной функцией преобразования.

ETA[J] - массив погрешностей измерительного комплекса (таблицы 11, 12, 13).

ЕТА[0] - абсолютная погрешность термометра, °С.

Таблица 11 - Наименование погрешности измерительного преобразователя и ее номер (J) в массиве ETA[J]

Таблица 12 - Наименование погрешности средства измерений или вторичного преобразователя и ее номер (J) в массиве ETA[J]

Таблица 13 - Наименование погрешности средства измерений и ее номер (J) в массиве ETA[J]

ЕТА[32] - систематическая составляющая погрешности планиметра (или интегратора), %.

ЕТА[33] - случайная составляющая погрешности планиметра (или интегратора), %. ЕТА[34] - верхний предел измерения температуры, °С.

ЕТА[35] - нижний предел измерения температуры, °С.

4.1.10 Характеристики измерительного комплекса (измерение плотности природного газа при стандартных условиях)

EROSA[J] - массив погрешностей измерительного комплекса (таблица 14).

Таблица 14 - Наименование погрешности измерительного преобразователя и ее номер (J) в массиве EROSA[J]

EROSA[6] - относительная погрешность средства измерений, %.

EROSA[7] - EROSA[22] - относительные погрешности измерения концентраций компонентов природного газа (см. 3.1.5), %.

4.1.11 ErrTau - погрешность определения интервала времени TauAV, %.

4.1.12 DHD - абсолютная погрешность хода приводного механизма диаграммы, мин.

4.1.13 Параметры эксплуатации (измеряемые за определенный промежуток времени - месяц, год и т.д.)

T1LO, Т1НО - нижнее и верхнее значения температуры окружающей среды, определяемые измерительным преобразователем, °С.

T2LO, Т2НО - нижнее и верхнее значения температуры окружающей среды, определяемые средством измерений, °С.

UL, UH - нижнее и верхнее значения напряжения питания, В.

TLG, THG - нижнее и верхнее значения температуры среды, °С.

PLG, PHG - нижнее и верхнее значения давления среды, бар.

DPL, DPH - нижнее и верхнее значения перепада давления на сужающем устройстве, бар.

ROSL, ROSH - нижнее и верхнее значения плотности природного газа при стандартных условиях, кг/м³.

YIM[J,1] - верхние значения концентраций компонентов природного газа (1 ≤ J ≤ 16, см. 3.1.5).

YIM[J,2] - нижние значения концентраций компонентов природного газа (1 ≤ J ≤ 16, cм. 3.1.5).

RoL, RoH - нижнее и верхнее значения плотности природного газа при рабочих условиях, кг/м³.

4.1.14 Характеристика измеряемых параметров:

IfConstRo = 0, если плотность природного газа при стандартных условиях - непрерывно измеряемая величина; IfConstRo = 1, если плотность природного газа при стандартных условиях - условно-постоянная величина;

IfConstY = 0, если концентрации компонентов природного газа - непрерывно измеряемые величины; IfConstY = 1, если концентрации компонентов природного газа - условно-постоянные величины;

IfConstP = 0, если давление среды - непрерывно измеряемая величина; IfConstP = 1, если давление среды - условно-постоянная величина;

IfConstT = 0, если температура среды - непрерывно измеряемая величина; IfConstT = 1, если температура среды - условно-постоянная величина;

IfConstDp = 0, если перепад давления - непрерывно измеряемая величина; IfConstDp = 1, если перепад давления - условно-постоянная величина.

4.1.15 Характеристика измерительного комплекса (измерение плотности природного газа при рабочих условиях)

ErowA[J] - массив погрешностей измерительного комплекса (таблицы 15, 16).

Таблица 15 - Наименование погрешности измерительного преобразователя и ее номер (J) в массиве ErowA[J]

Таблица 16 - Наименование погрешности корнеизвлекающего средства измерений и ее номер (J) в массиве ERowA[J]

ERowA[21] - верхний предел измерения плотности газа при рабочих условиях, кг/м³. ERowA[22] - нижний предел измерения плотности газа при рабочих условиях, кг/м³.

4.2 Выходные данные

RO, KAPPA, MU, QC, VC, VM, HS[1], HS[2] - см. 3.2.

HSV[1] - высшая теплота сгорания (энергосодержание) среды, МДж.

HSV[2] - низшая теплота сгорания (энергосодержание) среды, МДж.

EQR - случайная составляющая относительной погрешности расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

EQS - систематическая составляющая относительной погрешности расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

EQ1 - полная относительная погрешность расчета расхода среды для измерительного трубопровода с порядковым номером NNit.

EVC - полная относительная погрешность расчета количества среды для измерительного трубопровода с порядковым номером NNit за время TauAV.

ЕН[1] - полная относительная погрешность расчета высшей теплоты сгорания среды для измерительного трубопровода с порядковым номером NNit.

ЕН[2] - полная относительная погрешность расчета низшей теплоты сгорания среды для измерительного трубопровода с порядковым номером Nnit.

4.3 Листинг модуля расчета погрешности определения расхода и количества жидкостей и газов

4.3.1 Типы используемых переменных:

FL: text; NNIT, IERR, SIZMDPD, SIZMPD, SIZMTD: byte; T1HO, T1LO, UH, UL, PHG, PLG, THG, TLG, DPH, DPL, T2HO, T2LO, ROSH, ROSL, EDPS, EDPR, EDP, EPS, EPR, EP, ETS, ETR, ET, EROSS, EROSR, EROS, EROWS, EROWR, EROW, EQS, EQR, EQ1, EVC, DD20, DT20, DD, DT, ALFAT, ALFASU, RSH, RN, TAUP, TAUAV, ERRTAU, DHD, RO, KAPPA, MU, VM, ROSC, QC, VCV, VC, RoH, RoL, RoM: real; RocStr: string[6]; LRA: array [1..10] of real; DAC, TAC: array [1..3] of real; EPA: err1; YIM: err2; HS, HSV, EH: err3; NMETHKA, NSUBA, NSUZA, NVARYA, VARPA, VarRoA, AMOUNTR, SODSUA: err7; NRA: err8; EDPA: err9; ETA: err10; EROSA, ERowA; err11;

type err1 = array [1..26] of real; err2 = array [1..16,1..2] of real; err3 = array [1..2] of real; err7 = array [1..30] of byte; err8 = array [1..10] of byte; err9 = array [1..35] of rsal; err10 = array [0..35] of real; err11 = array [1..22] of real;

4.3.2 Листинг модуля ERRQSNX

Unit ERRQSNX;

Interface

{----------------------------------------------------------------------}

Procedure QS;

Implementation uses DOS, CRT;

{----------------------------------------------------------------------}

var

Zc: real;

const

Bi:array( 1.. 16] of real=(0.0436.0.0894.0.1288,0.1783,0.1703,0.2345,

0.2168,0.2846,0.3521,0.4278,0.0173,0.0728,

0.1,0.0,0.02,-0.0051);

Procedure DP(E:err9; DPQreal; var EDPSY, EDPRA, EDP : real);

forward;

Procedure P(E:err1; PM:real; var EPSY, EPRA, EP : real);

forward;

Procedure T(E : err10; TA:real; var ETSY, ETRA, ET : real);

forward;

Procedure ROS(E:err11;RIM : err1; var EROSS, EROSR, EROS : real);

forward;

Procedure C(BET,RE,DP,KSH,KK:real; var ERQreal; var IERR:byte);

forward;

Procedure EPSI(BET,DP,P:real; var EEPS:real);

forward;

Procedure EPH(YR:err1; T,P,ROS:real; var ERO,EKAP,EVIS:real);

forward;

Procedure EPHP(N:byte; T,P:real; var ERO,EKAP,EVIS:real);

forward;

Procedure ROW(E:err11; var EROWS, EROWR, EROW : real); forward;

Procedure QCalc(NRQ:byte; T.P,DP,YA,YY,DDO,DTO:real; YR:err1;

var BET,RE.KSH,KK,ROC,QC,VCV:real; var VCC:real; var VMC.ROX,

KAPPAX,MUX:real; var HS:err3); forward;

Procedure HSP(RIM:err1; YA,YY:real; var ERRH1,ERRH2:real);

forward;

function r_(A, R: real): real; forward;

{--------------------------------------------------------------------------}

{Основная процедура: выполняет диспетчерские функции, а также расчет коэффициентов влияния измеряемых параметров и погрешностей определения расхода и количества среды}

Procedure QS;

const

EDD = 0.07; EDT = 0.40;

label

1,3;

var

I,J,NRQ : byte;

RIM,RIMX,DVY,TETYI : err1; DVR,YA,YY,DVDP,DVPG,DVTG,T1,T2 :

real;

DVAY,HSX : err3;

BET,RE,KSH,KK,QC1,VCV1,VM1,QC2,VCV2,VM2: real;

VC1,VC2 : real;

ERC,EEPS,ERO,EKAP,EVIS,ERRH1,ERRH2:real;

DPMD,DPM1,DPM2,TETDP,DDD,DD01,DD02,TETDD,DDT,DT01,DT02,TETDT,

DKAP,KAP1,KAP2,TETKAP,DMU,MU1,MU2,TETMU,TETC,TETEPS,TETRO,

TETYA.TETYY,

DROS,ROS1,ROS2,TETROS,DYA,DYY,YA1,YY1,YA2,YY2,R01,R02,DTMG,

TETT.DPMG,

TETP,EYR:real;

DRIM : real;

begin { QS }

ROSC:=0.0; DVR:=0.0; for I := 1 to 2 do DVAY[I]:=0.0;

for I := 1 to 16 do DVY[I]:=0.0;

if NSUBA[NNIT] = 0 then begin

ifNVARYA[NNIT] = 1 then begin

ROSC := 2.0*ROSH*ROSL/(ROSH+ROSL);

if IfConstRo <> 0 then DVR:=100.0"(ROSH-ROSL)/(ROSH+ROSL);

YA:=0.5*(YIM[11,1]+YIM[11,2]);

YY:=0.5*(YIM[12,1]+YIM[12,2]);

for J := 11 to 12 do begin

I:=J-10;

if((YIM[J,1] <> 0) and (IfConstY <> 0)) or

((YIM[J,2] <> 0) and (IfConstY <> 0)) then

DVAY[I]:=100.0*(YIM[J,1]-YIM(J,2])/(YIM[J,1]+YIM[J,2])end;

end else

for I := 1 to 16 do begin RIM[I] := 0.5*(YIM[I,1]+YIM[I,2]);

if((YIM[I,1] <> 0) and (IfConstY <> 0)) or

((YIM[I,2] <> 0) and (IfConstY <> 0)) then

DVY[I]:=100.0*(YIM[I,1]-YIM[I,2])/(YIM[I,1]+YIM[I,2]);

RIMX[I]:=RIM[I]

end;

end; {endif}

PMG sqr(0.5*(sqrt(PHG)+sqrt(PLG)));

DVPG := 0.0; ifIfConstP <> 0 then DVPG:=100.0*(PHG-PLG)/(PHG+PLG);

T1:=THG+273.15; T2:=TLG+273.15;

TMG .= 4.0*Tl*T2/sqr(sqrt(Tl)+sqrt(T2))-273.15;

DVTG := 0.0; ifIfConstT <> 0 then DVTG:=100.0*(T1-T2)/(T1+T2);

DPM:=sqr(0.5*(sqrt(DPH)+sqrt(DPL)));

DVDP := 0.0; ifIfConstDp <> 0 then DVDP:=100.0*(DPH-DPL)/(DPH+DPL);

if(NSubA[NNit]=0) and (VarRoA[NNit]=l) then

RoM:=sqr(0.5*(sqrt(RoH)+sqrt(RoL)));

NRQ:=1;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QC,VCV,VC,VM,RO,KAPPA,MU,HS);

C(BET,RE,DPH,KSH,KK, ERCJERR);

if IERR = 0 then begin DP(EDPA,DPM,EDPS,EDPR,EDP);

P(EPA,PMG,EPS, EPR, EP);

T(ETA, TMG, ETS, ETR, ET);

ROS(EROSA, RIM, EROSS, EROSR, EROS);

EPSI(BET.DPM,PMG,EEPS);

EPH(RIM,TMG,PMG,ROSC, ERO,EKAP,EVIS);

if (NSubA[NNIT]=0 and (VarRoA[NNit]=l) then begin

ROW(EROWA, EROWS, EROWR, EROW); ERO := EROWS

end

else EROWR := 0.0;

NRQ:=0;

TETDP:=0.0;

ifEDP <> 0.0 then begin

DPMD:=0.001*EDP*DPM; DPM1:=DPM+DPMD; DPM2:=DPM-DPMD;

QCalc(NRQ,TMG,PMG,DPMl,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QC1,VCV1,VC1,VM1,RO,KAPPA,MU,HS);

QCalc(NRQ,TMG,PMG,DPM2,YA,YY,DD20,DT20,RIM,BET,RE,KSH,KK,ROSC,

QC2,VCV2,VC2,VM2,RO.KAPPA,MU,HS);

TETDP:=DPM*(QC1-QC2)/(2.0*DPMD)/QC;

end:

DDD:=0.001*EDD*DD20: DD01:=DD20+DDD; DD02:=DD20-DDD;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD01,DT20,RIM,BET,RE,KSH,KK,ROSC,

QC1,VCV1,VC1,VM1,RO,KAPPA.MU,HS):

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD02,DT20,RIM,BET,R£,KSH,KK,ROSC,

QC2,VCV2,VC2,VM2,RO,KAPPA,MU,HS);

TETDD:=DD20*(QC1-QC2)/(2.0*DDD)/QC;

DDT:=0.001*EDT*DT20; DT01:=DT20+DDT; DT02:=DT20-DDT;

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT01,RIM,BET,RE,KSH,KK,ROSC,

QC 1,VCV1,VC1,VM1,RO, KAPPA, MU,HS);

QCalc(NRQ,TMG,PMG,DPM,YA,YY,DD20,DT02,RIM,BET,RE,KSH,KK,ROSC,