本文共 1095 字,大约阅读时间需要 3 分钟。
简单做个备份
VTI介质ADCIGs中剩余曲率根据角度的变化可以进一步分析得出反演策略,根据相速度公式求解的深度残差关于Thomsen参数的偏导数:
#include#include #include float deltaz(float v, float e, float d, float theta){ float vp=2000; float beta=0.0; float t=1; theta *= 3.1415926/180.0; return (v-vp)*t/(2.0*cos(theta)*cos(beta)) *sqrt(0.5 + e*sin(theta-beta)*sin(theta-beta) + 0.5*sqrt( pow(1+2*e*sin(theta-beta)*sin(theta-beta),2.0 ) - 2*(e-d)*sin(2*theta-2*beta)*sin(2*theta-2*beta)) );}int main(int argc, char **argv){ float diff = 0.1; float theta, dzv, dze, dzd; int i; FILE *fp1 = fopen("dz.txt","w"); for(i=0;i<90;i++){ theta = i*1.0; dzv = deltaz(1800-diff, 0.5, 0.3, theta) - deltaz(1800, 0.5, 0.3, theta); dze = deltaz(1800, 0.5+diff, 0.3, theta) - deltaz(1800, 0.5, 0.3, theta); dzd = deltaz(1800, 0.5, 0.3+diff, theta) - deltaz(1800, 0.5, 0.3, theta); printf("%f, %f, %f\n",dzv,dze,dzd); fprintf(fp1,"%f %f %f\n",dzv,dze,dzd); }}
得出关系曲线: