自然电位测电阻率.cpp

利用sp求地层水电阻率原码

#include<iostream.h>                  
#include<math.h>
class well{
public:
	double ssp,Rmf,Rw,T1;
private:
	double  K, Rfw, RmfN, RmfeN, RweN, RwN;
public:
     well(double ssp=0.0,double Rmf=0.0,double T1=0.0);
	 double Jisuan(double ssp,double Rmf,double T1)
	{   K=70.7*(273+T1)/298;
		Rfw=pow(10,-ssp/K);
		RmfN=(1.8*T1+39)*Rmf/(1.8*24+39);
		if(RmfN>0.1) RmfeN=0.85*RmfN;
		else RmfeN=(146*RmfN-5)/(337*RmfN+77);
		RweN=RmfeN/Rfw;
		if(RweN>0.12) RwN=-0.58+pow(10,0.69*RweN-0.24);
		else RwN=(77*RweN+5)/(146-337*RweN);
		Rw=82*RwN/(1.8*T1+39);
		return Rw;
};
     ~well();
	 void print()
	 {cout<<"地层水电阻率：Rw="<<Rw<<"Ω·m"<<endl;};
};	
well::well(double ssp,double Rmf,double T1)  { ssp=0.0;Rmf=0.0;T1=0.0;};	
well::~well() {cout<<"Thanks for your trip!"<<endl;};
main()                                                              
{ 	well well1;                                                       
	double ssp,Rmf,T1;                                                 
	cout<<"请输入ssp："<<endl;                                      
	cin>>ssp;
cout<<"请输入地层温度（华氏）T1："<<endl;
    cin>>T1;
	cout<<"请输入地层温度下泥浆滤液电阻率Rmf："<<endl;
	cin>>Rmf;
	T1=(T1-32)/1.8;
    well1.Jisuan(ssp,Rmf,T1);
    well1.print();
    return 0;                                
}