baimananopti.cpp

石油公司网管系统优化系统软件源码,很有价值的.

//---------------------------------------------------------------------------

#include <vcl.h>
#pragma hdrstop

#include "BaiMaNanOpti.h"
#include <math.h>
//---------------------------------------------------------------------------
#pragma package(smart_init)
#pragma resource "*.dfm"
TTBaiMaNanOpti *TBaiMaNanOpti;
//---------------------------------------------------------------------------
__fastcall TTBaiMaNanOpti::TTBaiMaNanOpti(TComponent* Owner)
    : TForm(Owner)
{
    Terminated=false;
}
//---------------------------------------------------------------------------
void __fastcall TTBaiMaNanOpti::TimerTimer(TObject *Sender)
{
    //Timer->Enabled=false;
    Execute();
    ModalResult=mrOk;
}
//---------------------------------------------------------------------------
double TTBaiMaNanOpti::GetYaJiang(double G,double T,double Pz,
        double QiYouBi,double HanShuiLv,double d,double L,double GaoCha)
{
    double Q0,ql,Ql,Qg,Rs,B,QiMiDu,YouMiDu,Pr,_P,HL0,HL,TiJiHanYe,TiJiHanQi,HunHeliuSu,A;
    double QiNianDu,YouNianDu,YeNianDu,TiJiHanShui,HunHeNianDu,HunHeMiDu,ShuiMiDu,YeMiDu;
    double Z,C,S0,MoZuXiShu,ShuiPingYaJiang,GaoChengYaJiang,QiZheSuanLiuSu,Gl,YaJiang;
    double Tr,Prr,HunHeLiuSu,Re,Fe,Trr;
    L*=1e-3;
    Pr=Pz+1000;;
    int i;
    for(i=0;i<100;i++)
    {
        _P=(Pr+Pz)/2;
        //气液密度
        QiMiDu=WuXing.GetQiMiDu(T,_P);
        if(QiMiDu==0)
        {
            Application->MessageBox("天然气密度计算错误！计算中止！", "错误", MB_OK);
            Terminate();
            return 0;
        }
        YouMiDu=WuXing.GetYouMiDu(T);
        //气液流量
        Q0=G/YouMiDu;
        ql=Q0/(1-HanShuiLv);
        Rs=WuXing.GetRs(YouMiDu,_P,T);
        if(Rs>QiYouBi)
        {
            Application->MessageBox("天然气溶解系数大于气油比！计算中止！", "错误", MB_OK);
            Terminate();
            return 0;
        }
        B=WuXing.GetB(Rs,YouMiDu,T);
        Ql=ql*B;
        Trr=(T+273.15)/264.87;
        Prr=(_P)/(43.960*1e5);
        Z=WuXing.GetZ(Prr,Trr);
        if(Z==0)
        {
            Application->MessageBox("天然气压缩系数计算错误！计算中止！", "错误", MB_OK);
            Terminate();
            return 0;
        }
//        Qg=ql*(QiYouBi*G/Q0-Rs)*0.101*T*Z/(_P*20);
        Qg=ql*(QiYouBi-Rs)*101325*(T+273.15)*Z/(_P*273.15);
        //体积含液
        TiJiHanYe=Ql/(Ql+Qg);
        TiJiHanQi=1-TiJiHanYe;
        ShuiMiDu=WuXing.GetShuiMiDu(T);
        //液相密度
        YeMiDu=ShuiMiDu*HanShuiLv+YouMiDu*(1-HanShuiLv);
        YeMiDu=(YeMiDu+QiXiangDuiMiDu*1.205*Rs)/B;
        //混合流速
        HunHeLiuSu=4*(Ql+Qg)/(M_PI*d*d);
        //气油粘度
        QiNianDu=WuXing.GetQiNianDu(QiMiDu,T);
        YouNianDu=WuXing.GetYouNianDu(T);
        if(YouNianDu==0)
        {
            Application->MessageBox("粘度计算错误，计算终止！", "错误", MB_OK);
            Terminate();
            return 0;
        }
        //液相粘度
        TiJiHanShui=(HanShuiLv/ShuiMiDu)/(HanShuiLv/ShuiMiDu+(1-HanShuiLv)/YouMiDu);
        if(TiJiHanShui<0.7)
        {
            YeNianDu=YouNianDu*exp(4*TiJiHanShui);
        }
        else
        {
            YeNianDu=YouNianDu*exp(6.5*TiJiHanShui);
        }
        //混合粘度
        HunHeNianDu=QiNianDu*TiJiHanQi+YeNianDu*TiJiHanYe;
        HL=TiJiHanYe;
        //循环确定截面含液率
        do
        {
            HL0=HL;
            //混合物密度
            HunHeMiDu=YeMiDu*TiJiHanYe*TiJiHanYe/HL0+QiMiDu*TiJiHanQi*TiJiHanQi/(1-HL0);
            //雷诺数
            Re=d*HunHeLiuSu*HunHeMiDu/HunHeNianDu;
            //查图求截面含液率
            HL=WuXing.GetHL(Re,TiJiHanYe);
            if(HL==-1)
            {
                Application->MessageBox("持液率计算错误，计算终止！", "错误", MB_OK);
                Terminate();
                return 0;
            }
        }while(fabs(HL-HL0)/HL>0.05);
        //阻力系数
        C=log(TiJiHanYe);
        S0=1.281+0.478*C+0.444*C*C+0.094*C*C*C+0.00843*C*C*C*C;
        C=1-C/S0;
        MoZuXiShu=C*(0.0056+0.5/pow(Re,0.32));
        //计算水平压降
        ShuiPingYaJiang=MoZuXiShu*HunHeLiuSu*HunHeLiuSu*L*HunHeMiDu*1e-3/(2*d);
        //计算因高程产生的压降
        A=M_PI*d*d/4;
        QiZheSuanLiuSu=Qg/A;
        if(QiZheSuanLiuSu<=15)
        {
            Fe=1/(1+1.0785*pow(QiZheSuanLiuSu,1.006));
        }
        else
        {
            Gl=YeMiDu*Ql;
            Fe=3.175*1e-5*sqrt(Gl/A)/pow(QiZheSuanLiuSu,0.7);
        }
        GaoChengYaJiang=Fe*YeMiDu*g*GaoCha/1e6;
        //总压将
        YaJiang=ShuiPingYaJiang+GaoChengYaJiang;
        YaJiang*=1e6;
        //判断精度
        if(fabs(Pr-Pz-YaJiang)/YaJiang<0.05)
        {
            Pr=Pz+YaJiang;
            break;
        }
        Pr=Pz+YaJiang;
    }
    return Pr;
}

void TTBaiMaNanOpti::BuJuJieGuo()
{
    WellUnit *Well,*CurrentWell,*TempWell,*TempPreWell,*BaseWell,*PreWell;
    WellHead *CurrentHead=Head,*TempHead,*CurrentTemp;
    int TempInt,i;
    AnsiString TempStr;
    TFileStream *File;
    //总体结果
    ADOResult->TableName="over";
    ADOResult->Open();
    ADOResult->First();
    //油井位置
//    ADOResult->MoveBy(1);
    //总耗钢量
    ADOResult->Edit();
    ADOResult->FieldByName("耗钢量")->AsFloat=BestUnit->Fitness;
    //总输量
    ADOResult->FieldByName("流量")->AsFloat=Head->YouLiuLiang;
    //增压点个数
    TempInt=0;
    CurrentHead=Head->Next;
    while(CurrentHead!=0)
    {
        TempHead=CurrentHead->Sub;
        while(TempHead!=0)
        {
            TempInt++;
            TempHead=TempHead->Sub;
        }
        CurrentHead=CurrentHead->Next;
    }
    CurrentHead=Head->Sub;
    if(CurrentHead!=0)
    {
        CurrentHead=CurrentHead->Next;
    }
    while(CurrentHead!=0)
    {
        TempHead=CurrentHead->Sub;
        while(TempHead!=0)
        {
            TempInt++;
            TempHead=TempHead->Sub;
        }
        CurrentHead=CurrentHead->Next;
    }
    ADOResult->FieldByName("增压点")->AsInteger=TempInt;
    //转油站个数
    TempInt=0;
    CurrentHead=Head->Next;
    while(CurrentHead!=0)
    {
        if(!CurrentHead->Fixed)
        {
            TempInt++;
        }
        CurrentHead=CurrentHead->Next;
    }
    CurrentHead=Head->Sub;
    if(CurrentHead!=0)
    {
        CurrentHead=CurrentHead->Next;
    }
    while(CurrentHead!=0)
    {
        if(!CurrentHead->Fixed)
        {
            TempInt++;
        }
        CurrentHead=CurrentHead->Next;
    }
    ADOResult->FieldByName("转油站")->AsInteger=TempInt;
    ADOResult->Post();
    ADOResult->Close();
    File=new TFileStream("Nan.pos",fmCreate|fmOpenWrite|fmShareExclusive);
    //编号
    TempStr=Head->No;
    TempInt=TempStr.Length();
    File->Write(&TempInt,sizeof(int));
    File->Write(TempStr.c_str(),TempStr.Length());
    //新建？
    if(Head->Fixed)
    {
        TempInt=1;
    }
    else
    {
        TempInt=0;
    }
    File->Write(&TempInt,sizeof(int));
    //联合站坐标
    File->Write(&(BestUnit->Gens[2*GenNum-2]),sizeof(double));
    File->Write(&(BestUnit->Gens[2*GenNum-1]),sizeof(double));
    CurrentHead=Head->Next;
    i=0;
    while(CurrentHead!=0)
    {
        //标识位，0-油井，1-增压点，2-转油站，3-联合站
        TempInt=2;
        File->Write(&TempInt,sizeof(int));
        //转油站编号
        TempStr=CurrentHead->No;
        //编号所占字节数
        TempInt=TempStr.Length();
        File->Write(&TempInt,sizeof(int));
        File->Write(TempStr.c_str(),TempStr.Length());
        //新建？
        if(CurrentHead->Fixed)
        {
            TempInt=1;
        }
        else
        {
            TempInt=0;
        }
        File->Write(&TempInt,sizeof(int));
        Well=CurrentHead->Well;
        while(Well!=0)
        {
            //标识位
            TempInt=0;
            File->Write(&TempInt,sizeof(int));
            //油井编号
            TempInt=Well->No.Length();
            File->Write(&TempInt,sizeof(int));
            File->Write(Well->No.c_str(),Well->No.Length());
            Well=Well->Next;
        }
        TempHead=CurrentHead->Sub;
        while(TempHead!=0)
        {
            //标识位，0-油井，1-增压点，2-转油站，3-联合站
            TempInt=1;
            File->Write(&TempInt,sizeof(int));
            //增压点编号
            TempStr=TempHead->No;
            //编号所占字节数
            TempInt=TempStr.Length();
            File->Write(&TempInt,sizeof(int));
            File->Write(TempStr.c_str(),TempStr.Length());
            //坐标
            TempHead->X=BestUnit->Gens[2*i];
            TempHead->Y=BestUnit->Gens[2*i+1];
            i++;
            Well=TempHead->Well;
            while(Well!=0)
            {
                //标识位
                TempInt=0;
                File->Write(&TempInt,sizeof(int));
                //油井编号
                TempInt=Well->No.Length();
                File->Write(&TempInt,sizeof(int));
                File->Write(Well->No.c_str(),Well->No.Length());
                Well=Well->Next;
            }
            TempHead=TempHead->Sub;
        }
        CurrentHead->X=BestUnit->Gens[2*i];
        CurrentHead->Y=BestUnit->Gens[2*i+1];
        i++;
        CurrentHead=CurrentHead->Next;
    }
    Head->X=BestUnit->Gens[2*i];
    Head->Y=BestUnit->Gens[2*i+1];
    i++;
    CurrentHead=Head->Sub;
    if(CurrentHead!=0)
    {
        //标识位
        TempInt=3;
        File->Write(&TempInt,sizeof(int));
        //编号
        TempStr=CurrentHead->No;
        TempInt=TempStr.Length();
        File->Write(&TempInt,sizeof(int));
        File->Write(TempStr.c_str(),TempStr.Length());
        //新建？
        if(CurrentHead->Fixed)
        {
            TempInt=1;
        }
        else
        {
            TempInt=0;
        }
        File->Write(&TempInt,sizeof(int));
        //联合站坐标
//        CurrentHead->X=BestUnit->Gens[2*i];
//        CurrentHead->Y=BestUnit->Gens[2*i+1];
//        i++;
        CurrentHead=CurrentHead->Next;
        while(CurrentHead!=0)
        {
            //标识位，0-油井，1-增压点，2-转油站，3-联合站
            TempInt=2;
            File->Write(&TempInt,sizeof(int));
            //转油站编号
            TempStr=CurrentHead->No;
            //编号所占字节数
            TempInt=TempStr.Length();
            File->Write(&TempInt,sizeof(int));
            File->Write(TempStr.c_str(),TempStr.Length());
            //新建？
            if(CurrentHead->Fixed)
            {
                TempInt=1;
            }
            else
            {
                TempInt=0;
            }
            File->Write(&TempInt,sizeof(int));
            Well=CurrentHead->Well;
            while(Well!=0)
            {
                //标识位
                TempInt=0;
                File->Write(&TempInt,sizeof(int));
                //油井编号
                TempInt=Well->No.Length();
                File->Write(&TempInt,sizeof(int));
                File->Write(Well->No.c_str(),Well->No.Length());
                Well=Well->Next;
            }
            TempHead=CurrentHead->Sub;
            while(TempHead!=0)
            {
                //标识位，0-油井，1-增压点，2-转油站，3-联合站
                TempInt=1;
                File->Write(&TempInt,sizeof(int));
                //增压点编号
                TempStr=TempHead->No;
                //编号所占字节数
                TempInt=TempStr.Length();
                File->Write(&TempInt,sizeof(int));
                File->Write(TempStr.c_str(),TempStr.Length());
                //坐标
                TempHead->X=BestUnit->Gens[2*i];
                TempHead->Y=BestUnit->Gens[2*i+1];
                i++;
                Well=TempHead->Well;
                while(Well!=0)