◇◇新语丝(www.xys.org)(xys.dxiong.com)(xys.dropin.org)(xys-reader.org)◇◇ 再析《基于可公度方法的川滇地区地震趋势研究》 zeroyear 　　关于《基于可公度方法的川滇地区地震趋势研究》一文，网友们已经纷纷指 出其假设的灾害的可预测的周期性根本不存在，也举了些反例说明：符合其所谓 可公度特性的年份未必发生地震，不符合的也未必不震。这个显然连历史都“预 测”很不准的“模型”，还敢拿来预测未来？但是，这论文似乎还就是会给人们 留了这样一个印象：也许其假设是很不可靠的，甚至完全荒唐，但若不强求其都 能测准，并允许其适当测漏，那总算也是有其一套粗略理论的。这么一来，举些 反例虽然能够简单地推翻其预测，却难免让人以为反例的数量也可能是很有限的， 论文中的数字游戏多少还有些神秘色彩未被揭开。因此，我就趁周末时间想找找 更多的反例，来彻底弄弄明白这里的所谓可公度方法究竟价值何在，算是再抬举 它一次。好在这个文章实在与地震专业无关，其所谓推理也只需用初等数学知识 就能对付，验证起来其实毫不费力。下面让我们先直奔其“精髓”： 　　作者在2.1节的“三元可公度法预测”中浩浩荡荡地列出75个组合（其中有 不少重复的）之后称道，“通过计算可以看出，得出的结果与实际的年份一致， 且每一个结果都能写出三组以上的三元可公度式子，由此说明，川滇地区的强震 表现出了极好的可公度特性。”我不知道所谓可公度特性的好坏究竟是怎么衡量 的，但是作者既然认为若每一个结果都能写得出三组以上的式子那就能算是“极 好”了，因此给我的感觉就是：每个三组足以说明问题。有道是，三就是多嘛！ 然后，作者列举出分别与2007和2008相关的几个组合以说明预测2007和2008的合 理性，其中与2007相关的2组，与2008相关的3组，与2007相关的这2组后面还诡 异地跟了1组本该是与2006相关的组合。这就更让我以为能举出3组与2008相关的 就已经很说明2008的问题了。相比起来，与2007相关的似乎只能举出2组，那个 与2006相关的大概就是用来近似给2007凑数的，这样似乎也就体现出了2007的 “灾难信号”比2008的弱。再接下来，作者在第2.2节的“四元可公度法预测” 中列了4个组合以说明12是个神奇的数字，最后又在第2.3节的“五元可公度法预 测”中列出了与2008相关的6个组合。 　　初看此文时，这让我感觉能列出6个组合简直就该是很罕见的情况了。也许 作者在这里并没有玩障眼法的意图，而只是有些笔误和有些解释不足，但是我的 确有这种被误导的感觉。而相比其文字，作者给出的数据显得如此洋洋洒洒，一 下子还真容易让人以为这大震年号之间确有玄机。现在就让我们也用颇占篇幅的 数据，用三个表格来彻底揭了这所谓可公度方法预测地震的唬人之皮。 　　一、无关痛痒的“三组以上” 　　不管真正的“三元可公度法”是个什么东西，作者在论文中无非就是用了 Y1=Y3+Y4-Y2这样一个式子（Y即YEAR。不妨将25或者26个年号作为一个数列或者 集合，Y1、Y2、Y3、Y4表示其中任意四个年号），这就等价于Y1+Y2=Y3+Y4，即： 要求两个年号相加之和等于另外两个年号相加之和。如果不动手亲自算一算，我 们直观上恐怕不大知道满足这个要求的组合数到底会是在哪个数量级，由于作者 只列三个，这就难免让人以为能找出十个恐怕也不容易。但果然是不算不知道！ 其实与这25或者26个年号的每一个相关的组合数均是大大地超过三个的。比如， 在25个年号的数列中，与Y1=X1=1937相关的就有如下组合（其中的数字为序号）： [1+09=04+06][1+11=02+09][1+11=03+08][1+13=03+10] [1+13=04+09][1+14=04+12][1+15=03+13][1+15=05+10] [1+15=07+08][1+16=02+14][1+16=06+09][1+17=06+10] [1+18=02+15][1+19=07+09][1+20=04+14][1+20=08+10] [1+21=03+16][1+21=08+11][1+22=04+19][1+22=07+13] [1+23=03+20][1+23=08+13][1+23=10+11][1+24=07+14] [1+24=09+13][1+25=05+19][1+25=06+17][1+25=08+14] 　　由于作者采用的数据从1913始到1996年止，这跨度为83年的25个年号之间， 至多隔7年就会发生大地震，那么我们不妨按照论文发表的时间，朝前后各推10 年，看看1996年（即X25）后的这20年的任何一年若大震，则取这年为X26，形成 的每一种由26个年号构成的数列的所谓可公度特性究竟会有多大差别，由此来看 看作者慧眼识出来的这个2008究竟有多么特别。下表即列出这26个年号的数据， 最后一行的0000年号随列号变动，列号97至16即表示1997至2016这20年，前25行 表示与这25个年号中的每一个相关的组合数，最后一行则表示与所在列对应的年 号相关的组合数。比如，与2008相关的组合有42个，对应的意思是：如果预测 1996之后直到2008年才发生地震，那么取X26=2008后整个数列的所谓可公度特性 就体现在08这一列。从这些都远远地大于三的数字里面，作者究竟是怎么看出平 淡无奇的2007和2008的“灾害信号”比较强的！ 表一： 0000 97 98 99 00 01 02 03 04 05 06 07 08 09 10 11 12 13 14 15 16 1913 30 30 35 31 31 33 29 31 31 31 31 33 33 31 31 31 31 30 31 33 1917 41 43 40 42 41 41 41 43 43 41 41 41 41 40 41 43 41 44 39 40 1923 43 46 46 43 44 44 43 42 44 46 43 46 41 42 42 40 44 44 43 42 1925 44 41 40 41 41 39 41 43 40 43 39 40 39 37 41 41 40 39 39 41 1933 53 55 51 50 51 48 53 54 53 51 50 53 49 52 53 51 51 51 52 52 1936 52 53 52 55 54 54 51 52 53 50 55 55 52 52 53 53 52 53 52 51 1941 60 57 57 59 55 58 59 58 58 57 58 56 57 55 56 56 56 58 55 55 1942 59 61 61 57 61 61 58 60 61 60 59 60 57 58 56 57 61 59 57 56 1948 55 57 57 56 59 60 56 58 55 55 58 56 54 55 57 56 54 57 54 54 1950 57 57 53 54 53 52 54 55 55 55 52 52 53 53 51 54 54 52 52 53 1952 61 62 62 62 57 58 61 58 60 61 59 61 58 60 59 58 58 59 60 59 1955 55 57 55 59 58 58 58 55 56 55 56 58 58 56 58 56 56 57 56 54 1960 64 64 63 65 62 63 63 61 63 64 64 65 62 62 61 60 62 64 61 63 1967 54 57 56 57 56 60 55 58 57 54 57 55 53 57 52 54 55 56 56 53 1970 57 57 61 55 57 55 55 55 58 54 58 56 54 56 54 53 56 55 54 55 1971 62 62 62 63 57 60 59 60 60 62 57 61 59 58 60 59 59 60 57 58 1973 57 56 55 53 55 58 52 55 51 53 53 55 51 55 56 53 53 52 51 53 1974 51 51 52 50 50 51 55 47 50 47 48 49 52 49 52 49 48 50 47 47 1976 54 56 55 55 58 53 53 53 57 51 55 54 54 52 54 51 55 54 53 53 1979 59 64 55 58 58 57 57 60 59 56 57 59 54 56 55 57 56 58 54 57 1981 54 52 53 58 50 53 55 53 53 54 51 51 52 55 49 51 50 50 51 53 1988 45 47 44 45 44 45 46 47 44 45 51 43 45 46 45 46 47 45 44 44 1989 58 55 56 54 56 54 56 55 56 55 55 60 53 56 55 55 56 57 55 54 1995 50 51 48 50 50 51 53 49 51 48 50 48 49 50 50 49 50 54 47 49 1996 41 43 44 41 43 44 44 45 42 44 41 43 41 43 42 44 41 43 48 39 0000 44 50 43 43 39 42 41 41 42 36 38 42 29 34 33 31 34 39 28 28 二、莫名其妙的“周期” 　　论文中的“四元可公度法”不过就是将其“三元可公度法”的Y1+Y2=Y3+Y4 替换成Y1+Y2=Y3+Y4+12，我们不妨再扩展一步，即用Y1+Y2=Y3+Y4+T来衡量。下 表第一行，即1996那行，表示没有预测时，即只有25个年号，取各个T值所对应 的所有不重复的组合数（所谓重复的组合即如，Y1=Y3+Y4+T-Y2是与Y1相关的， 但同时也是与Y2相关的，即满足Y2=Y3+Y4+T-Y1。若T=0，则这还同时是与Y3和Y4 相关的），其它行则表示预测了对应的年发生地震后的“可公度特性”。该表的 第1列为年号，第2列至第13列依次表示T取0，2，5，8，10，12，15，19，25， 30，40，50（这里填充的数字0是为了文字排版的对齐。注意，若取T=0，实质就 退化为所谓“三元可公度法”，所以T=0对应的列的每一项，即是表一的某一列 的总和再减去其中重复组合数的结果）。比如，若预测1996年之后直到2008年发 生地震，那么取X26=2008来计算，得到所有满足Y1+Y2=Y3+Y4+12的不重复的组合 数为625个（论文中列出了其中的区区4个组合）。从这些数字里面，作者又是为 何偏偏就选择了T=12？ 表二： 0000 000 002 005 008 010 012 015 019 025 030 040 050 1996 296 584 574 532 577 562 554 492 506 524 429 392 1997 340 671 665 624 667 642 638 576 583 600 509 458 1998 346 664 663 619 660 653 643 572 580 602 508 455 1999 339 667 658 616 656 647 646 577 581 604 498 463 2000 339 676 660 623 662 634 633 567 578 596 497 450 2001 335 668 647 609 656 646 639 571 579 595 501 463 2002 338 668 656 609 657 652 633 572 591 597 510 459 2003 337 665 666 609 659 636 624 566 580 594 505 454 2004 337 662 646 600 659 639 629 565 577 588 504 463 2005 338 663 651 610 649 637 631 576 581 593 501 450 2006 332 667 646 621 642 634 625 563 572 591 501 455 2007 334 655 647 603 648 631 624 555 576 598 489 463 2008 338 654 649 603 656 625 624 565 579 591 495 464 2009 325 655 654 598 645 631 624 561 575 588 500 454 2010 330 657 644 603 646 633 617 556 573 587 489 457 2011 329 647 638 598 641 627 613 560 579 584 491 457 2012 327 657 639 599 640 632 621 556 570 591 494 461 2013 330 645 645 599 640 620 622 552 558 592 486 456 2014 335 643 628 589 645 626 615 554 567 587 482 459 2015 324 645 634 592 638 626 616 545 564 583 486 456 2016 324 652 632 601 635 621 609 552 562 589 483 458 三、“元数”高了反而暴露出自相矛盾 　　论文中的“五元可公度法”不过就是将“三元”的扩展成用 Y1+Y2+Y3=Y4+Y5+Y6来衡量。下面的表三就是与“三元”对应的那个表一的“五 元”版，考虑到帖子的宽度有限，拆分成A和B两部分列出如下： 表三A： 0000 1997 1998 1999 2000 2001 2002 2003 2004 2005 2006 1913 2719 2712 2677 2694 2711 2685 2702 2694 2680 2678 1917 2789 2758 2793 2753 2764 2766 2746 2731 2727 2742 1923 3020 2984 2990 3002 2994 2982 2975 2989 2974 2940 1925 3073 3059 3071 3064 3052 3073 3045 3000 3037 3026 1933 3220 3187 3216 3210 3194 3215 3162 3148 3157 3158 1936 3304 3281 3288 3256 3261 3246 3258 3241 3231 3235 1941 3327 3333 3336 3297 3333 3292 3275 3268 3274 3268 1942 3352 3322 3322 3327 3294 3285 3303 3277 3262 3265 1948 3414 3404 3406 3377 3384 3350 3362 3347 3356 3345 1950 3478 3440 3466 3472 3451 3443 3447 3407 3388 3403 1952 3453 3417 3418 3413 3444 3425 3382 3401 3385 3361 1955 3501 3482 3500 3453 3451 3443 3431 3429 3432 3410 1960 3446 3438 3442 3417 3434 3406 3394 3403 3387 3360 1967 3517 3480 3502 3450 3473 3430 3437 3417 3427 3431 1970 3480 3461 3436 3460 3454 3470 3437 3437 3407 3425 1971 3416 3393 3390 3371 3419 3385 3401 3359 3366 3339 1973 3435 3421 3452 3446 3439 3402 3428 3393 3430 3400 1974 3452 3448 3440 3451 3455 3424 3387 3437 3412 3418 1976 3407 3385 3393 3373 3354 3388 3371 3344 3329 3357 1979 3343 3270 3358 3315 3314 3308 3294 3265 3283 3290 1981 3335 3305 3302 3259 3300 3260 3273 3253 3242 3256 1988 3226 3203 3215 3190 3223 3176 3167 3171 3165 3171 1989 3106 3143 3135 3133 3107 3132 3094 3100 3101 3091 1995 2998 2979 3012 2975 2981 2959 2927 2972 2964 2975 1996 3018 3009 2994 3006 2998 2979 2956 2951 2988 2980 0000 2993 2890 2938 2860 2884 2812 2758 2714 2708 2692 表三B： 0000 2007 2008 2009 2010 2011 2012 2013 2014 2015 2016 1913 2681 2649 2657 2669 2653 2659 2643 2660 2657 2615 1917 2733 2734 2728 2732 2728 2689 2702 2680 2708 2704 1923 2958 2941 2969 2949 2946 2953 2915 2915 2912 2915 1925 3033 3006 3025 3032 2984 2989 2977 2966 2986 2972 1933 3151 3121 3153 3119 3100 3104 3091 3079 3080 3064 1936 3176 3184 3206 3190 3182 3164 3157 3137 3147 3146 1941 3248 3250 3248 3247 3230 3218 3213 3169 3204 3191 1942 3246 3243 3258 3238 3259 3224 3173 3188 3203 3198 1948 3300 3307 3340 3316 3283 3301 3271 3240 3267 3241 1950 3406 3380 3386 3379 3362 3344 3329 3318 3330 3323 1952 3359 3343 3357 3323 3319 3320 3311 3287 3275 3279 1955 3399 3373 3374 3376 3361 3346 3343 3297 3318 3315 1960 3350 3342 3352 3336 3335 3336 3307 3282 3301 3273 1967 3404 3407 3416 3374 3415 3394 3343 3326 3338 3328 1970 3378 3395 3409 3365 3389 3370 3334 3343 3353 3316 1971 3367 3312 3334 3329 3293 3311 3288 3261 3298 3277 1973 3386 3361 3398 3337 3345 3344 3320 3325 3320 3307 1974 3389 3375 3369 3376 3351 3353 3345 3302 3341 3316 1976 3323 3330 3321 3310 3301 3307 3261 3253 3282 3252 1979 3266 3251 3286 3253 3269 3231 3225 3198 3240 3188 1981 3254 3226 3246 3204 3227 3232 3231 3172 3196 3179 1988 3104 3139 3148 3119 3113 3108 3093 3086 3116 3086 1989 3075 3025 3077 3048 3056 3044 3024 3008 3028 3016 1995 2944 2972 2944 2925 2934 2930 2909 2866 2927 2897 1996 2989 2943 2988 2943 2944 2948 2939 2916 2877 2936 0000 2611 2549 2625 2525 2503 2471 2376 2282 2368 2294 　　乍一看可能感觉这些数据都实在大得出乎意料，对此作者怎么会只列出6个 组合来说明问题呢！是不是我搞错了？为了直观起见，不妨再占些篇幅列出一个 例子的小部分清单。比如，如果固定Y2，取Y2=X2=1917，那么满足要求的与 Y1=X1=1913相关的组合就已有如下的如此之多（第一列[]中的数字为年号数列的 序号，第二列为Y1+Y2+Y3，第三列为Y4+Y5+Y6，每一列都满足Y1=1913的所谓可 公度，即X1=Y4+Y5+Y6-X2-Y3）： [1+2+13=03+04+08] 1913+1917+1960=5790 1923+1925+1942=5790 [1+2+14=03+05+07] 1913+1917+1967=5797 1923+1933+1941=5797 [1+2+15=03+04+11] 1913+1917+1970=5800 1923+1925+1952=5800 [1+2+15=03+06+07] 1913+1917+1970=5800 1923+1936+1941=5800 [1+2+15=04+05+08] 1913+1917+1970=5800 1925+1933+1942=5800 [1+2+16=03+06+08] 1913+1917+1971=5801 1923+1936+1942=5801 [1+2+17=03+04+12] 1913+1917+1973=5803 1923+1925+1955=5803 [1+2+17=04+06+08] 1913+1917+1973=5803 1925+1936+1942=5803 [1+2+18=03+05+09] 1913+1917+1974=5804 1923+1933+1948=5804 [1+2+19=03+05+10] 1913+1917+1976=5806 1923+1933+1950=5806 [1+2+19=03+07+08] 1913+1917+1976=5806 1923+1941+1942=5806 [1+2+19=04+05+09] 1913+1917+1976=5806 1925+1933+1948=5806 [1+2+20=03+06+10] 1913+1917+1979=5809 1923+1936+1950=5809 [1+2+20=04+06+09] 1913+1917+1979=5809 1925+1936+1948=5809 [1+2+21=03+05+12] 1913+1917+1981=5811 1923+1933+1955=5811 [1+2+21=03+06+11] 1913+1917+1981=5811 1923+1936+1952=5811 [1+2+21=04+06+10] 1913+1917+1981=5811 1925+1936+1950=5811 [1+2+21=05+06+08] 1913+1917+1981=5811 1933+1936+1942=5811 [1+2+22=03+04+15] 1913+1917+1988=5818 1923+1925+1970=5818 [1+2+22=04+05+13] 1913+1917+1988=5818 1925+1933+1960=5818 [1+2+22=04+07+11] 1913+1917+1988=5818 1925+1941+1952=5818 [1+2+23=03+04+16] 1913+1917+1989=5819 1923+1925+1971=5819 [1+2+23=03+06+13] 1913+1917+1989=5819 1923+1936+1960=5819 [1+2+23=03+07+12] 1913+1917+1989=5819 1923+1941+1955=5819 [1+2+23=04+08+11] 1913+1917+1989=5819 1925+1942+1952=5819 [1+2+23=05+06+10] 1913+1917+1989=5819 1933+1936+1950=5819 [1+2+23=06+07+08] 1913+1917+1989=5819 1936+1941+1942=5819 [1+2+24=03+08+13] 1913+1917+1995=5825 1923+1942+1960=5825 [1+2+24=03+10+11] 1913+1917+1995=5825 1923+1950+1952=5825 [1+2+24=04+05+14] 1913+1917+1995=5825 1925+1933+1967=5825 [1+2+24=04+09+11] 1913+1917+1995=5825 1925+1948+1952=5825 [1+2+24=05+08+10] 1913+1917+1995=5825 1933+1942+1950=5825 [1+2+24=06+07+09] 1913+1917+1995=5825 1936+1941+1948=5825 [1+2+25=03+05+15] 1913+1917+1996=5826 1923+1933+1970=5826 [1+2+25=03+06+14] 1913+1917+1996=5826 1923+1936+1967=5826 [1+2+25=03+09+12] 1913+1917+1996=5826 1923+1948+1955=5826 [1+2+25=04+07+13] 1913+1917+1996=5826 1925+1941+1960=5826 [1+2+25=05+07+11] 1913+1917+1996=5826 1933+1941+1952=5826 [1+2+25=06+08+09] 1913+1917+1996=5826 1936+1942+1948=5826 　　如果我给出的数据无误，那么由表三B可见，论文中列举的仅仅是与2008相 关的2549个组合中的区区6个，而被作者认为“灾害信号”较之为弱从而根本没 有在文中列出示例组合的2007所相关的组合数却是较之更大的2611个。总体上 2007那列的数据也并不比2008那列的小。茫茫然全部都是这么大的组合数，不管 预测哪一年会地震，所谓可公度特性难道不都是极其极其极其的“极好”吗？ 　　通过如上三个表的数据，我们容易发现其中的微弱规律也不过是来自年号的 奇偶间隔，以及距边界年号（即1913、1996和2016）的距离，这是数列本身的特 点，与地震发生与否何干？就算灾害的可预测的周期性是存在的，恐怕用这么个 “可公度方法”也是根本无法将之找出来的。玩弄数字游戏的作者（们），如网 友们所言，估计就是事先认定了2008，再手工尝试着找出相关的一些组合列出来 唬唬人，很可能竟然连最基本的计算机程序都根本没有去用，就如此无知无畏地 “计算”出2008年的“灾害信号”尤其强！ (XYS20080518) ◇◇新语丝(www.xys.org)(xys.dxiong.com)(xys.dropin.org)(xys-reader.org)◇◇