No processo de perfilagem de poços de petróleo e gás, as propriedades entre duas camadas axiais podem apresentar diferenças. Estas descontinuidades das propriedades algumas vezes trazem dificuldades na representação do comportamento dos campos electromagnéticos. Nas análises eletromagnéticas de perfilagem de poço de petróleo, o Método de Casamento de Modos (NMM) vem sendo utilizado a partir da combinação de técnicas numéricas com analíticas e da resolução das equações de Maxwell em meios heterogêneos, cujos campos eletromagnéticos, na direção axial, são representados via decomposição espectral, e, na direção radial, por meio de solução das equações diferenciais ordinárias. Assim, com o objetivo de explorar a representação dos campos eletromagnéticos sobre as interfaces axiais via NMM, este trabalho propõe a utilização de funções B-Splines cúbicas na expansão dos campos, na direção axial, bem como a incorporação da multiplicidade de nós, que permitem uma melhor representação do comportamento dos campos nas interfaces axiais. O algoritmo implementado foi validado nas suas diversas etapas ao comparar os autovalores da representação modal com os obtidos analiticamente em problemas canônicos, e os resultados finais comparados com exemplos apresentados na literatura e calculados com o Método de Diferenças Finitas no Tempo (FDTD). O método NMM, combinado às funções B-Splines cúbicas, foi aplicado a vários perfis de poços de petroléo que serviram de exemplo neste estudo. Um estudo comparativo mostrou que o uso das funções B-Splines cúbicas com multiplicidade dos nós sobre as descontinuidades permitiu uma redução no número total de nós nas expansões dos campos eletromagnéticos.