山野莽夫

  • 归档
    • 随笔
    • 建站资源
    • 分享
    • 代码
  • 地球物理学
    • 专业课
    • 概念解释
  • 计算机
  • 互联网
  • 教程
  • 规划
  • 实验室
    • 珍藏的软件
    • 贴吧云签到
    • A1账号自助申请
山野莽夫
小学生的挣扎的点点滴滴
  1. 首页
  2. 地球物理学
  3. 专业课
  4. 正文

地下电偶极子电场空间分布

2016年7月9日 4876点热度 1人点赞 2条评论

实习项目一 地下电偶极子电场空间分布

一、实习目的:

掌握地下不同倾角电偶极子地面电场分布特征和规律。

二、要求:

1、推导并计算地下电偶极子的空间分布;

2、熟悉使用相关软件绘出地下不同倾角电偶极子地面电场强度、电位场分;

3、熟悉不同参数偶极子地面场分布规律及特征。

三、实验过程

1、实验计算原理



图1 地下电偶极子结构图

根据示意图确立地表电位、电场强度计算公式:




2、试验参数

1)计算区域坐标范围__(-800,800)________________

2)电偶极子几何参数:倾角0、30、60、90、120、150、180(取值范围0-180度,至少取7个);偶极距长__200______;埋深___300_____。

3) 电偶极子物理参数:电流___1.2_____;介质电阻率___200_____。

3、程序设计及代码(可附打印源代码)

Matlab代码

p=200;I=1.2;pai=3.141592653;L=200;h=300;

E=zeros(161,161);

U=zeros(161,161);

PU=zeros(161);

PE=zeros(161);


for i=0:30:180


for x=-800:10:800


for y=-800:10:800

m=(x+810)/10;

n=(y+810)/10;

s=-1*x*I*p/(2*pai)*power(power((x+L/2*cos(i/180.0*pai)),2.0)+y*y+power((h-L/2*sin(i/180.0*pai)),2),(-1.5))...

+x*I*p/(2*pai)*power(power((x-L/2*cos(i/180.0*pai)),2.0)+y*y+power((h+L/2*sin(i/180.0*pai)),2),(-1.5)) ;

E(m,n)=s;


end


end

x=-800:10:800;

y=-800:10:800;

subplot(4,2,(i/30+1));

contour(x,y,E,100);


end

figure


for i=0:30:180


for x=-800:10:800


for y=-800:10:800

m=(x+810)/10;

n=(y+810)/10;

ss=I*p/(2*pai)/power(power((x+L/2*cos(i/180.0*pai)),2)+y*y+power((h-L/2*sin(i/180.0*pai)),2),(0.5))...

-I*p/(2*pai)/power(power((x-L/2*cos(i/180.0*pai)),2)+y*y+power((h+L/2*sin(i/180.0*pai)),2),(0.5)) ;

U(m,n)=ss;


end


end

x=-800:10:800;

y=-800:10:800;

subplot(4,2,(i/30+1));

contour(x,y,U,100);


end

figure

for i=0:30:180

for x=-800:10:800

m=(x+810)/10;

PU(m)=I*p/(2*pai)*L*(x*cos(i/180.0*pai)-h*sin(i/180.0*pai))/(power((h*h+x*x),1.5));

end

subplot(4,2,(i/30+1));

x=-800:10:800;

plot(x,PU);

end

figure


for i=0:30:180

for x=-800:10:800

m=(x+810)/10;

PE(m)=-I*p/(2*pai)*L*((h*h-2*x*x)*cos(i/180.0*pai)+3*h*x*sin(i/180.0*pai))/(power((h*h+x*x),2.5));

end

subplot(4,2,(i/30+1));

x=-800:10:800;

plot(x,PE);

end

C语言代码

#include<stdio.h>

#include<stdlib.h>

#include<math.h>

int main()

{

    float p=200,I=1.2,pai=3.141592653,L=200,h=300,E,U,PE,PU,Ey;

    int i,j,x,y;

    i=30;

    FILE*fp;FILE *fp1,*fp2,*fp3;

    fp=fopen("xxx.txt","w");

fp1=fopen("www.txt","w");

fp2=fopen("ooo.txt","w");

fp3=fopen("kkk.txt","w");

    for(x=-800;x<=800;x=x+10)

    {

for(y=-800;y<=800;y=y+10)

     {

E=-1*x*I*p/(2*pai)*pow(pow((x+L/2*cos(i/180.0*pai)),2.0)+y*y+pow((h-L/2*sin(i/180.0*pai)),2),(-1.5))

             +x*I*p/(2*pai)*pow(pow((x-L/2*cos(i/180.0*pai)),2.0)+y*y+pow((h+L/2*sin(i/180.0*pai)),2),(-1.5)) ;

    //    printf("%d\t%d\t%.12f\n",x,y,E);

fprintf(fp,"%d\t%d\t%.12f\n",x,y,E);

     }

    }

for(x=-800;x<=800;x=x+10)

    {

for(y=-800;y<=800;y=y+10)

     {

Ey=-1*y*I*p/(2*pai)*pow(pow((x+L/2*cos(i/180.0*pai)),2.0)+y*y+pow((h-L/2*sin(i/180.0*pai)),2),(-1.5))

             +y*I*p/(2*pai)*pow(pow((x-L/2*cos(i/180.0*pai)),2.0)+y*y+pow((h+L/2*sin(i/180.0*pai)),2),(-1.5)) ;

    //    printf("%d\t%d\t%.12f\n",x,y,Ey);

fprintf(fp3,"%d\t%d\t%.12f\n",x,y,Ey);

     }

    }

for(x=-800;x<=800;x=x+10)

    {

for(y=-800;y<=800;y=y+10)

     {

U=I*p/(2*pai)/pow(pow((x+L/2*cos(i/180.0*pai)),2)+y*y+pow((h-L/2*sin(i/180.0*pai)),2),(0.5))

             -I*p/(2*pai)/pow(pow((x-L/2*cos(i/180.0*pai)),2)+y*y+pow((h+L/2*sin(i/180.0*pai)),2),(0.5)) ;

         //printf("%d\t%d\t%f\n",x,y,U);

fprintf(fp1,"%d\t%d\t%f\n",x,y,U);

     }

    }

    for(x=-800;x<=800;x=x+10)

    {

for(i=0;i<=180;i=i+15)

        {

        PU=I*p/(2*pai)*L*(x*cos(i/180.0*pai)-h*sin(i/180.0*pai))/(pow((h*h+x*x),1.5));

        PE=-I*p/(2*pai)*L*((h*h-2*x*x)*cos(i/180.0*pai)+3*h*x*sin(i/180.0*pai))/(pow((h*h+x*x),2.5));

        if(i==0)

        {

        printf("%d\t%.10f\t%.10f",x,PU,PE);

fprintf(fp2,"%d\t%.10f\t%.10f",x,PU,PE);

        }

        else

        {

     // printf("\t%.10f\t%.10f",PU,PE);

fprintf(fp2,"\t%.10f\t%.10f",PU,PE);

        }

        }

// printf("\n");

        fprintf(fp2,"\n");

    }

    fclose(fp);

    fclose(fp1);

fclose(fp2);

    return 0;

}

4、电场图(可附录打印图,要求给出图例及相关计算参数)

1)电位图主剖面图

070916_0713_4.jpg

2)电场强度剖面图


3)电位平面图


4)电场强度图平面图


四、实验结论及认识

本次试验成图使用matlab,在质量方面不及surfer的效果。电场强度随测点远离偶极源时衰减很快(与x的三次方成反比),而电位则衰减较慢(与x的平方成反比),但均较一个点电源的电场衰减快。电场强度曲线也不对称,在偶极子的右上方,由于电流的集中,电场强度出现主极值。在电位曲线的两个极值点处,电场强度等于零,在零值点两侧有两个次极值点,它们与向下流的电流回流到地表的两个集流区对应。如果在地表绘出电场强度的等值线,便会看到在偶极子的右上方有明显之负心,负心两侧有两个正心,三者成直线。

电位曲线的特点是到处均为负值,在x=0处有负极值,在x=±h/2处为电位曲线的拐点,分别对应E曲线负的和正的极值。

标签: 地球电磁学 电偶极子 电场分布
最后更新:2018年4月6日

小菜菜

菜鸟

打赏 点赞
下一篇 >

文章评论

  • admire

    cdut的??

    2020年11月23日
    回复
    • 小菜菜

      @admire 当年是的嘞

      2020年11月24日
      回复
  • razz evil exclaim smile redface biggrin eek confused idea lol mad twisted rolleyes wink cool arrow neutral cry mrgreen drooling persevering
    取消回复

    此站点使用Akismet来减少垃圾评论。了解我们如何处理您的评论数据。

    联系方式

    QQ群 | TG群 | 邮箱

    最新 热点 随机
    最新 热点 随机
    Azure Student 微软云 学生订阅 免费12个月用量避坑注意点集合 MP3音频文件格式详细解析 python按固定采样点个数分割wav格式音频 愉快使用谷歌免费人工智能平台colab,训练你的神经网络模型,为你的学术生活添砖加瓦 华为云版轻量应用服务器-云耀云服务器简单体验评测 Cloudflare 免费CDN自定义节点ip之自选cloudflare 高速节点ip工具分享
    wordpress用户角色介绍以及如何最简单的修改wordpress用户不同角色的详细权限 大明王朝1566云盘资源 宝塔docker管理器的使用(以docker搭建临时邮件forsaken-mail为例) 十八届五中全会课件 中梯法的电阻率正演模拟 二层电测深电阻率正演模拟
    标签聚合
    虚拟机 onedrive c语言 ppt 宝塔面板 wordpress 地震学程序 模板
    最近评论
    小菜菜 发布于 6 个月前(11月24日) 这玩意已经废了,成收割工具了,不能再用了。
    eamon 发布于 7 个月前(11月07日) 我一年不用了才发现这个休眠管理费每月15,一共扣了我135元,然后我消费还消费不了,我宁愿消费掉也不...
    magic 发布于 11 个月前(07月03日) 请问账号不注销会有什么影响吗?
    magic 发布于 11 个月前(07月01日) 我想问一下 如果不注销账号就留着会怎么样
    qwp6601 发布于 12 个月前(06月04日) 有没有方法改为bing

    COPYRIGHT © 2021 shanyemangfu.com. ALL RIGHTS RESERVED.

    Theme Kratos Made By Seaton Jiang

    蜀ICP备15031791号-2