☆─────────────────────────────────────☆
tongxl (但使龙城飞将在，不叫胡马度阴山) 于 (Thu Jan 20 20:58:51 2022) 提到:
其实机缘巧合，我自己因为是油田的子弟，对油田有一种发自内心的热爱和责任感。也是这份责任感，让我选择了今天的道路。既然有人质疑我科研水平。我也就说最近两年疫情期间成果吧。国内最近这30年吧，一直吹波动方程成像，大学就是同济那边马在田王华忠，院所就是中科院地球物理所李幼铭那边。做的东西，真的是沐猴而冠，关键一个个胆子都很大，敢说敢骗。
一个波动方程的数值解，就是抄抄开源的seismic unix。可是问题来了，有限差分数值解，第一个问题，频散（色散)怎么处理？空间的，时间的？第二个问题，边界条件怎么处理？这么多人搞了这么多边界，效果如何？pml，海绵，abs？等等等。第三个问题，近场不连续，怎么处理？这三个问题，我可以负责任的说，本来是高校院所该去推公式、想办法解决的问题没有一个人去做。除了抄抄海外文章，灌灌水。
这两年吧，色散问题，我这处理的，空间完全可以做到nyquist，时间方向误差不超过万分之一。边界，新算法简单有效，计算量低，高低频率都可以处理的干干净净。至于说，近场，也有了好的解析解，至少结果和实际数据可以几乎完美匹配。说了这些，什么目的呢？
科研，本来就是为了解决问题。但是很多看似愚蠢的问题，正是启发我们思维的过程。人们也不愿意把自己的曾经的愚蠢暴露在大庭广众之下。克服了这些，才能真的去解决问题，做到真的把论文写在祖国大地上。
另外一个问题，就是长春地院的风气问题。实在不愿意提起来。忽然想起来，这个单位第一个校长叫李四光，居然能够，也敢说大庆油田是他发现的。算是开了新中国不要脸的先河。大庆油田的发现，是黄汲青、谢家荣以及捷克地震队的功劳，和李四光半毛钱关系都没有。有这样的校长，后面的弟子们欺世盗名也就可以理解了
有 1 位用户评价了这篇文章：
[+5] freevictor: 就爱看你的帖子
☆─────────────────────────────────────☆
macas (ma) 于 (Fri Jan 21 08:52:36 2022) 提到:
想起《荀子·性恶》中的这句话：“故圣人化性而起伪，伪起而生礼义，礼义生而制法度。”
真相如何也许并没有那么重要，社会和谐才有利于民生
【 在 tongxl 的大作中提到: 】
: 其实机缘巧合，我自己因为是油田的子弟，对油田有一种发自内心的热爱和责任感。也是这份责任感，让我选择了今天的道路。既然有人质疑我科研水平。我也就说最近两年疫情期间成果吧。国内最近这30年吧，一直吹波动方程成像，大学就是同济那边马在田王华忠，院所就是中科院地球物理所李幼铭那边。做的东西，真的是沐猴而冠，关键一个个胆子都很大，敢说敢骗。
: 一个波动方程的数值解，就是抄抄开源的seismic unix。可是问题来了，有限差分数值解，第一个问题，频散（色散)怎么处理？空间的，时间的？第二个问题，边界条件怎么处理？这么多人搞了这么多边界，效果如何？pml，海绵，abs？等等等。第三个问题，近场不连续，怎么处理？这三个问题，我可以负责任的说，本来是高校院所该去推公式、想办法解决的问题没有一个人去做。除了抄抄海外文章，灌灌水。
: 这两年吧，色散问题，我这处理的，空间完全可以做到nyquist，时间方向误差不超过万分之一。边界，新算法简单有效，计算量低，高低频率都可以处理的干干净净。至于说，近场，也有了好的解析解，至少结果和实际数据可以几乎完美匹配。说了这些，什么目的呢？
: ...................
☆─────────────────────────────────────☆
Myarmor (余哀) 于 (Fri Jan 21 11:55:27 2022) 提到:
纯粹好奇问一下，一个PDE的数值解，有那么难，需要那么多人去做吗？
【 在 tongxl 的大作中提到: 】
: 其实机缘巧合，我自己因为是油田的子弟，对油田有一种发自内心的热爱和责任感。也是这份责任感，让我选择了今天的道路。既然有人质疑我科研水平。我也就说最近两年疫情期间成果吧。国内最近这30年吧，一直吹波动方程成像，大学就是同济那边马在田王华忠，院所就是中科院地球物理所李幼铭那边。做的东西，真的是沐猴而冠，关键一个个胆子都很大，敢说敢骗。
: 一个波动方程的数值解，就是抄抄开源的seismic unix。可是问题来了，有限差分数值解，第一个问题，频散（色散)怎么处理？空间的，时间的？第二个问题，边界条件怎么处理？这么多人搞了这么多边界，效果如何？pml，海绵，abs？等等等。第三个问题，近场不连续，怎么处理？这三个问题，我可以负责任的说，本来是高校院所该去推公式、想办法解决的问题没有一个人去做。除了抄抄海外文章，灌灌水。
: 这两年吧，色散问题，我这处理的，空间完全可以做到nyquist，时间方向误差不超过万分之一。边界，新算法简单有效，计算量低，高低频率都可以处理的干干净净。至于说，近场，也有了好的解析解，至少结果和实际数据可以几乎完美匹配。说了这些，什么目的呢？
: ...................
☆─────────────────────────────────────☆
tongxl (但使龙城飞将在，不叫胡马度阴山) 于 (Fri Jan 21 12:08:49 2022) 提到:
真有，现在都没到尽头，还是有很多问题没解决。而且各自有各自的问题。
疫情刚开始，我甚至翻出来warlfram的自动元胞机算法测试了一下，发现这玩意最大问题是
空间阶数只有2阶，太低，要做对，网格太小了，吃不消。
每种算法，做透了，和你看到的论文差异巨大。有的是真的愚蠢，做不对。有的是故意把人往坑里带
【 在 Myarmor 的大作中提到: 】
: 纯粹好奇问一下，一个PDE的数值解，有那么难，需要那么多人去做吗？
:
☆─────────────────────────────────────☆
Myarmor (余哀) 于 (Fri Jan 21 12:22:07 2022) 提到:
“空间阶数只有2阶”， 作何解？
几乎所有物理方程都是2阶吧。
3阶以上系统不稳定。
【 在 tongxl 的大作中提到: 】
: 真有，现在都没到尽头，还是有很多问题没解决。而且各自有各自的问题。
: 疫情刚开始，我甚至翻出来warlfram的自动元胞机算法测试了一下，发现这玩意最大问题是
: 空间阶数只有2阶，太低，要做对，网格太小了，吃不消。
: ....................
☆─────────────────────────────────────☆
glandetian (glandetian) 于 (Fri Jan 21 15:20:19 2022) 提到:
对您解决色散问题的方法特别感兴趣。在我的认识里，要么减小步长增加网格，要么用高精度的方法。不知道您是怎么做的。
还有吸收边界条件的问题。不用PML等传统方法，您是怎么处理的呢？
☆─────────────────────────────────────☆
haipai312 (haipai312) 于 (Fri Jan 21 15:22:46 2022) 提到:
☆─────────────────────────────────────☆
glandetian (glandetian) 于 (Fri Jan 21 15:41:29 2022) 提到:
楼主说的“2阶”是指数值算法的空间收敛阶，不是原微分方程的阶数。
2阶，大概可以理解为：如果步长缩小到原来的一半，那么误差变为原来的四分之一 ----- 二分之一的2次方。
【 在 Myarmor 的大作中提到: 】
: “空间阶数只有2阶”， 作何解？
: 几乎所有物理方程都是2阶吧。
: 3阶以上系统不稳定。
☆─────────────────────────────────────☆
tongxl (但使龙城飞将在，不叫胡马度阴山) 于 (Fri Jan 21 16:13:27 2022) 提到:
减小步长真的不行，我早年也是被国内的专家误导了，步长再小，传的时间足够长，还是会出来。
这个方向上，我知道的第一个生产相关的研究的人叫周铮铮，他推导出了二阶的时间频散的相位变化
做了第一个工作，叫做tvps（时变相位校正），是周铮铮在ion期间的工作。戴南荀把源代码偷到了物探局之后
发的文章（不厚道）。2013年，鬼子有开始校正振幅的，后来statoil的徐升、王美霞发了一篇cos函数校正的
这个研究算是成熟了。这些朋友们发文章的可以说。徐升王美霞那篇文章的可以把时间二阶控制在千分之几的范畴了
说通俗点就是数据加载前做振幅相位矫正，传播后做矫正。不影响计算效率
空间频散，spectrum的算法可以做到100%nyquist，里面技巧搞计算数学很多人已经搞定了。
PML的没问题，这个是传统公式实现有问题。
【 在 glandetian 的大作中提到: 】
: 对您解决色散问题的方法特别感兴趣。在我的认识里，要么减小步长增加网格，要么用高精度的方法。不知道您是怎么做的。
: 还有吸收边界条件的问题。不用PML等传统方法，您是怎么处理的呢？
☆─────────────────────────────────────☆
glandetian (glandetian) 于 (Fri Jan 21 16:16:57 2022) 提到:
受教了，非常感谢！
【 在 tongxl 的大作中提到: 】
: 减小步长真的不行，我早年也是被国内的专家误导了，步长再小，传的时间足够长，还是会出来。
: 这个方向上，我知道的第一个生产相关的研究的人叫周铮铮，他推导出了二阶的时间频散的相位变化
: 做了第一个工作，叫做tvps（时变相位校正），是周铮铮在ion期间的工作。戴南荀把源代码偷到了物探局之后
: ...................
☆─────────────────────────────────────☆
koleasu (koleasu) 于 (Fri Jan 21 20:03:48 2022) 提到:
哎 偏微分方程忘光了 一直心心念念想自己动手编数值解