简单做个备份
VTI介质ADCIGs中剩余曲率根据角度的变化可以进一步分析得出反演策略,根据相速度公式求解的深度残差关于Thomsen参数的偏导数:
#include<stdio.h> #include<stdlib.h> #include<math.h> 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); } }
得出关系曲线: