博客
关于我
强烈建议你试试无所不能的chatGPT,快点击我
VTI介质Thomsen参数关于入射角的敏感性曲线
阅读量:2032 次
发布时间:2019-04-28

本文共 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); }}

得出关系曲线:

你可能感兴趣的文章
Linux - Shell
查看>>
MySQL 如何执行关联查询
查看>>
算法运行时间n3,n2, n, nlogn对比
查看>>
java 中的suppressedException
查看>>
Rope --高效字符串处理数据结构
查看>>
Gap buffer -- 一个数据结构为可编辑的文本
查看>>
X/Open DTP模型,两阶段提交,JTA接口定义
查看>>
MySql binlog 日志
查看>>
spring 事务管理
查看>>
Redis集群
查看>>
线程中实现不可中断的任务
查看>>
世界城市时间计算
查看>>
MD5算法
查看>>
Redis基本命令-- 集合(Set)
查看>>
Hessian原理分析
查看>>
elastic-job --介绍
查看>>
顺序表——数组实现(C++代码)
查看>>
Java常用的输入输出语句
查看>>
码农常用的十个网站
查看>>
怎么理解用代码来实现数据结构
查看>>